Disentangling the role of phenotypic plasticity and genetic divergence in contemporary ecotype formation during a biological invasion

The occurrence of contemporary ecotype formation through adaptive divergence of populations within the range of an invasive species typically requires standing genetic variation but can be facilitated by phenotypic plasticity. The relative contributions of both of these to adaptive trait differentiation have rarely been simultaneously quantified in recently diverging vertebrate populations. Here we study a case of intraspecific divergence into distinct lake and stream ecotypes of threespine stickleback that evolved in the past 140 years within the invasive range in Switzerland. Using a controlled laboratory experiment with full-sib crosses and treatments mimicking a key feature of ecotypic niche divergence, we test if the phenotypic divergence that we observe in the wild results from phenotypic plasticity or divergent genetic predisposition. Our experimental groups show qualitatively similar phenotypic divergence as those observed among wild adults. The relative contribution of plasticity and divergent genetic predisposition differs among the traits studied, with traits related to the biomechanics of feeding showing a stronger genetic predisposition, whereas traits related to locomotion are mainly plastic. These results implicate that phenotypic plasticity and standing genetic variation interacted during contemporary ecotype formation in this case.

Analysis of Phenotypic Variation in Childhood Wheezing Disorders

Recurrent wheezing or asthma is a common problem in children that has increased considerably in prevalence in the past few decades. The causes and underlying mechanisms are poorly understood and it is thought that a numb er of distinct diseases causing similar symptoms are involved. Due to the lack of a biologically founded classification system, children are classified according to their observed disease related features (symptoms, signs, measurements) into phenotypes. The objectives of this PhD project were a) to develop tools for analysing phenotypic variation of a disease, and b) to examine phenotypic variability of wheezing among children by applying these tools to existing epidemiological data. A combination of graphical methods (multivariate co rrespondence analysis) and statistical models (latent variables models) was used. In a first phase, a model for discrete variability (latent class model) was applied to data on symptoms and measurements from an epidemiological study to identify distinct phenotypes of wheezing. In a second phase, the modelling framework was expanded to include continuous variability (e.g. along a severity gradient) and combinations of discrete and continuo us variability (factor models and factor mixture models). The third phase focused on validating the methods using simulation studies. The main body of this thesis consists of 5 articles (3 published, 1 submitted and 1 to be submitted) including applications, methodological contributions and a review. The main findings and contributions were: 1) The application of a latent class model to epidemiological data (symptoms and physiological measurements) yielded plausible pheno types of wheezing with distinguishing characteristics that have previously been used as phenotype defining characteristics. 2) A method was proposed for including responses to conditional questions (e.g. questions on severity or triggers of wheezing are asked only to children with wheeze) in multivariate modelling.ii 3) A panel of clinicians was set up to agree on a plausible model for wheezing diseases. The model can be used to generate datasets for testing the modelling approach. 4) A critical review of methods for defining and validating phenotypes of wheeze in children was conducted. 5) The simulation studies showed that a parsimonious parameterisation of the models is required to identify the true underlying structure of the data. The developed approach can deal with some challenges of real-life cohort data such as variables of mixed mode (continuous and categorical), missing data and conditional questions. If carefully applied, the approach can be used to identify whether the underlying phenotypic variation is discrete (classes), continuous (factors) or a combination of these. These methods could help improve precision of research into causes and mechanisms and contribute to the development of a new classification of wheezing disorders in children and other diseases which are difficult to classify.

Variable phenotypic expressivity in inbred retinal degeneration mouse lines: A comparative study of C3H/HeOu and FVB/N rd1 mice

PURPOSE Recent advances in optogenetics and gene therapy have led to promising new treatment strategies for blindness caused by retinal photoreceptor loss. Preclinical studies often rely on the retinal degeneration 1 (rd1 or Pde6b(rd1)) retinitis pigmentosa (RP) mouse model. The rd1 founder mutation is present in more than 100 actively used mouse lines. Since secondary genetic traits are well-known to modify the phenotypic progression of photoreceptor degeneration in animal models and human patients with RP, negligence of the genetic background in the rd1 mouse model is unwarranted. Moreover, the success of various potential therapies, including optogenetic gene therapy and prosthetic implants, depends on the progress of retinal degeneration, which might differ between rd1 mice. To examine the prospect of phenotypic expressivity in the rd1 mouse model, we compared the progress of retinal degeneration in two common rd1 lines, C3H/HeOu and FVB/N. METHODS We followed retinal degeneration over 24 weeks in FVB/N, C3H/HeOu, and congenic Pde6b(+) seeing mouse lines, using a range of experimental techniques including extracellular recordings from retinal ganglion cells, PCR quantification of cone opsin and Pde6b transcripts, in vivo flash electroretinogram (ERG), and behavioral optokinetic reflex (OKR) recordings. RESULTS We demonstrated a substantial difference in the speed of retinal degeneration and accompanying loss of visual function between the two rd1 lines. Photoreceptor degeneration and loss of vision were faster with an earlier onset in the FVB/N mice compared to C3H/HeOu mice, whereas the performance of the Pde6b(+) mice did not differ significantly in any of the tests. By postnatal week 4, the FVB/N mice expressed significantly less cone opsin and Pde6b mRNA and had neither ERG nor OKR responses. At 12 weeks of age, the retinal ganglion cells of the FVB/N mice had lost all light responses. In contrast, 4-week-old C3H/HeOu mice still had ERG and OKR responses, and we still recorded light responses from C3H/HeOu retinal ganglion cells until the age of 24 weeks. These results show that genetic background plays an important role in the rd1 mouse pathology. CONCLUSIONS Analogous to human RP, the mouse genetic background strongly influences the rd1 phenotype. Thus, different rd1 mouse lines may follow different timelines of retinal degeneration, making exact knowledge of genetic background imperative in all studies that use rd1 models.

Cryptic invasion drives phenotypic changes in central European threespine stickleback

Cryptic invasions are commonly associated with genetic changes of the native species or genetic lineage that the invaders replace. Phenotypic shifts resulting from cryptic invasions are less commonly reported given the relative paucity of historical specimens that document such phenotypic changes. Here, I study such a case in two populations of threespine stickleback from central Europe, comparing contemporary patterns of gene flow with phenotypic changes between historical and contemporary population samples. I find gene flow from an invasive lineage to be associated with significant phenotypic changes, where the degree of phenotypic change corresponds with the level of gene flow that a population receives. These findings underline the utility of combining genetic approaches with phenotypic data to estimate the impact of gene flow in systems where anthropogenic alterations have removed former geographic barriers promoting cryptic invasions.

Phenotypic plasticity is a negative, though weak, predictor of the commonness of 105 grassland species

Aim The usual hypothesis about the relationship between niche breadth and range size posits that species with the capacity to use a wider range of resources or to tolerate a greater range of environmental conditions should be more widespread. In plants, broader niches are often hypothesized to be due to pronounced phenotypic plasticity, and more plastic species are therefore predicted to be more common. We examined the relationship between the magnitude of phenotypic plasticity in five functional traits, mainly related to leaves, and several measures of abundance in 105 Central European grassland species. We further tested whether mean values of traits, rather than their plasticity, better explain the commonness of species, possibly because they are pre-adapted to exploiting the most common resources. Location Central Europe. Methods In a multispecies experiment with 105 species we measured leaf thickness, leaf greenness, specific leaf area, leaf dry matter content and plant height, and the plasticity of these traits in response to fertilization, waterlogging and shading. For the same species we also obtained five measures of commonness, ranging from plot-level abundance to range size in Europe. We then examined whether these measures of commonness were associated with the magnitude of phenotypic plasticity, expressed as composite plasticity of all traits across the experimental treatments. We further estimated the relative importance of trait plasticity and trait means for abundance and geographical range size. Results More abundant species were less plastic. This negative relationship was fairly consistent across several spatial scales of commonness, but it was weak. Indeed, compared with trait means, plasticity was relatively unimportant for explaining differences in species commonness. Main conclusions Our results do not indicate that larger phenotypic plasticity of leaf morphological traits enhances species abundance. Furthermore, possession of a particular trait value, rather than of trait plasticity, is a more important determinant of species commonness.

PHENOTYPIC CHARACTERIZATION OF SPONTANEOUS AND IN VITRO-MAINTAINED AKR LYMPHOMAS

T-cell lymphomas from AKR mice were studied to determine their potential as a model of T-cell differentiation. Homogeneous tumor cell lines have been used as model to study normal lymphocyte subpopulations, including differentiation lineages, functional properties, and the inducibility to maturation. The underlying concept is that each lymphoid tumor represents a monoclonal neoplastic proliferation of a discrete lymphoid subpopulation arrested at a particular differentiation stage.^ Individual tumors were analyzed to determine the extent of intertumor heterogeneity, and to determine whether lymphomas represented different thymocyte subsets, by determining the cell-surface antigenic phenotype, PNA-binding capacity, and terminal deoxynucleotidyl transferase (TdT) activity. Splenic and thymic tumor cells were compared to determine if the particular lymphoid microenvironment influenced T-cell marker expression. Several of the lymphomas were passaged in syngeneic hosts to verify the original tumor phenotype and to assess the stability of the cell surface and TdT phenotype after transplantation.^ Lymphomas were adapted to in vitro culture to determine whether the T-cell phenotype was maintained in the absence of the host microenvironment. Clonal progeny were analyzed and compared with each other and with parent cell lines to determine the extent of intratumor heterogeneity in this lymphoma system. Parent and cloned cell lines were passaged in vivo to determine whether alterations in surface phenotype occurred after transplantation.^ Our investigation has verified that most spontaneous AKR lymphomas phenotypically resemble known T-cell subsets, including both immature and mature thymic subpopulations. The in vitro lines, however, expressed a highly unstable phenotype in culture that included loss of Ly-1 and Ly-2 antigen expression. After transplantation in vivo, the in vitro lines exhibited alterations in phenotype, including re-expression of Ly antigen on some lymphomas. The inducibility of T-cell antigen markers on tumor cell lines passaged in vivo suggests that the in vitro lines may serve as a possible model system to study the molecular events involved in gene expression in the T-cell system. ^

TUMOR PROGRESSION: ANALYSIS OF THE INSTABILITY OF THE METASTATIC PHENOTYPE, SENSITIVITY TO RADIATION AND CHEMOTHERAPY (PHENOTYPIC DRIFT, BREAST CANCER)

The major complications for tumor therapy are (i) tumor spread (metastasis); (ii) the mixed nature of tumors (heterogeneity); and (iii) the capacity of tumors to evolve (progress). To study these tumor characteristics, the rat 13762NF mammary adenocarcinoma was cloned and studied for metastatic properties and sensitivities to therapy (chemotherapy, radiation and hyperthermia). The cell clones were heterogeneous and no correlation between metastatic potential and therapeutic sensitivities was observed. Further, these phenotypes were unstable during passage in vitro; yet, the changes were clone dependent and reproducible using different cryoprotected cell stocks. To understand the phenotypic instability, subclones were isolated from low and high passage cell clones. Each subclone possessed a unique composite phenotype. Again, no apparent correlation was seen between metastatic potential and sensitivity to therapy. The results demonstrated that (1) tumor cells are heterogeneous for multiple phenotypes; (2) tumor cells are unstable for multiple phenotypes; (3) the magnitude, direction and time of occurrence of phenotypic drift is clone dependent; (4) the sensitivity of cell clones to ionizing radiation (gamma or heat) and chemotherapy agents is independent of their metastatic potential; (5) shifts in metastatic potential and sensitivity to therapy may occur simultaneously but are not linked; and (6) tumor cells independently diverge to form several subpopulations with unique phenotypic profiles. ^

STRUCTURAL ANALYSIS OF THE TRANSCRIPTION PRODUCTS OF A TEMPERATURE-SENSITIVE PHENOTYPIC MUTANT OF MOLONEY SARCOMA VIRUS--MUSV TS110

Cells infected with a temperature sensitive phenotypic mutant of Moloney sarcoma virus (MuSVts110) exhibit a transformed phenotype at 33('(DEGREES)) and synthesize two virus specific proteins, p85('gag-mos), a gag-mos fusion protein and p58('gag), a truncated gag precursor protein (the gag gene codes for viral structural proteins and mos is the MuSV transforming gene). At 39('(DEGREES)) only p58('gag) is synthesized and the morphology of the cells is similar to uninfected NRK parental cells. Two MuSVts110 specific RNAs are made in MuSVts110-infected cells, one of 4.0 kb in length, the other of 3.5 kb. Previous work indicated that each of these RNAs arose by a single central deletion of parental MuSV genetic material, and that p58('gag) was made by the 4.0 kb RNA and p85('gag-mos) from the 3.5 kb RNA. The objective of my dissertation research was to map precisely the deletion boundaries of both of the MuSVts110 RNAs, and to determine the proper reading frame across both deletion borders. This work succeeded in arriving at the following conclusions: (a) Using S-1 nuclease analysis and primer extension sequencing, it was found that the 4.0 kb MuSVts110 RNA arose by a 1488 base deletion of 5.2 kb parental MuSV genomic RNA. This deletion resulted in an out of frame fusion of the gag and mos genes that resulted in the formation of a "stop" codon which causes termination of translation just beyond the c-terminus of the gag region. Thus, this RNA can only be translated into the truncated gag protein p58('gag). (b) S-1 analysis of RNA from cells cultivated at different temperatures demonstrated that the 4.0 kb RNA was synthesized at all temperatures but that synthesis of the 3.5 kb RNA was temperature sensitive. These observations supported the data derived from blot hybridization experiments the interpretation of which argued for the existence of a single provirus in MuSVts110 infected cells, and hence only a single primary transcript (the 4.0 kb RNA). (c) Analyses similar to those described in (a) above showed that the 3.5 kb RNA was derived from the 4.0 kb MuSVts110 RNA by a further deletion of 431 bases, fusing the gag and mos genes into a continuous reading frame capable of directing synthesis of the p85('gag-mos) protein. These sequence data and the presence of only one MuSVts110-specific provirus, indicate that a splice mechanism is employed to generate the 3.5 kb RNA since the gag and mos genes are observed to be fused in frame in this RNA. . . . (Author's abstract exceeds stipulated maximum length. Discontinued here with permission of author.) UMI ^

Molecular differentiation by PCR-RFLP technique shows that Antarctic fishes of the genus Notothenia, claimed to be two species, are phenotypic varieties with the same genetic pattern (Table 1)

The taxonomy of Antarctic fishes has been predominantly based on morphological characteristics rather than on genetic criteria. A typical example is the Notothenia group, which includes N. coriiceps Richardson, 1844, N. neglecta Nybelin, 1951 and N. rossii Richardson, 1844. The Polymerase Chain Reaction and Restriction Fragment Length Polymorphism (PCR-RFLP) technique was used to determine whether N. coriiceps Richardson, 1844 and N. neglecta Nybelin, 1951 are different or whether they are the same species with morphological, physiological and behavioural variability. N. rossii was used as control. Mitochondrial DNA (mtDNA) was isolated from muscle specimens of N. coriiceps Richardson, 1844, N. neglecta Nybelin, 1951 and N. rossii, which were collected in Admiralty Bay, King George Island. The DNA was used to amplify a fragment (690 base pairs) of the mitochondrial gene coding region of NADH dehydrogenase subunit 2. Further, the amplicon was digested with the following restriction enzymes: DdeI, HindIII and RsaI. The results showed a variation of the digestion pattern of the fragment amplified between N. rossii, and N. coriiceps Richardson, 1844 or N. neglecta Nybelin, 1951. However, no differences were found between N. coriiceps Richardson, 1844 and N. neglecta Nybelin, 1951, on the grounds of the same genetic pattern shown by the two fish.

Phenotypic evolution in a fossil gastropod species lineage: evidence for adaptive radiation?

Detecting speciation in the fossil record is a particularly challenging matter. Palaeontologists are usually confronted with poor preservation and limited knowledge on the palaeoenvironment. Even in the contrary case of adequate preservation and information, the linkage of pattern to process is often obscured by insufficient temporal resolution. Consequently, reliable documentations of speciation in fossils with discussions on underlying mechanisms are rare. Here we present a well-resolved pattern of morphological evolution in a fossil species lineage of the gastropod Melanopsis in the long-lived Lake Pannon. These developments are related to environmental changes, documented by isotope and stratigraphical data. After a long period of stasis, the ancestral species experiences a phenotypic change expressed as shift and expansion of the morphospace. The appearance of several new phenotypes along with changes in the environment is interpreted as adaptive radiation. Lake-level high stands affect distribution and availability of habitats and, as a result of varied functional demands on shell geometry, the distribution of phenotypes. The on-going divergence of the morphospace into two branches argues for increasing reproductive isolation, consistent with the model of ecological speciation. In the latest phase, however, progressively unstable conditions restrict availability of niches, allowing survival of one branch only.

Phenotypic evolution in a venerid bivalve species lineage from the late Middle Miocene Central Paratethys Sea: a multi-approach morphometric analysis

A morphometric analysis was performed for the late Middle Miocene bivalve species lineage of Polititapes tricuspis (Eichwald, 1829) (Veneridae: Tapetini). Specimens from various localities grouped into two stratigraphically successive biozones, i.e. the upper Ervilia Zone and the Sarmatimactra Zone, were investigated using a multi-method approach. A Generalized Procrustes Analysis was computed for fifteen landmarks, covering characteristics of the hinge, muscle scars, and pallial line. The shell outline was separately quantified by applying the Fast Fourier Transform, which redraws the outline by fitting in a combination of trigonometric curves. Shell size was calculated as centroid size from the landmark configuration. Shell thickness, as not covered by either analysis, was additionally measured at the centroid. The analyses showed significant phenotypic differentiation between specimens from the two biozones. The bivalves become distinctly larger and thicker over geological time and develop circular shells with stronger cardinal teeth and a deeper pallial sinus. Data on the paleoenvironmental changes in the late Middle Miocene Central Paratethys Sea suggest the phenotypic shifts to be functional adaptations. The typical habitats for Polititapes changed to extensive, very shallow shores exposed to high wave action and tidal activity. Caused by the growing need for higher mechanical stability, the bivalves produced larger and thicker shells with stronger cardinal teeth. The latter are additionally shifted towards the hinge center to compensate for the lacking lateral teeth and improve stability. The deepening pallial sinus is related to a deeper burrowing habit, which is considered to impede being washed out in the new high-energy settings.

Phenotypic plasticity of coralline algae in a High CO2 world

It is important to understand how marine calcifying organisms may acclimatize to ocean acidification to assess their survival over the coming century. We cultured the cold water coralline algae, Lithothamnion glaciale, under elevated pCO2 (408, 566, 770, and 1024 µatm) for 10 months. The results show that the cell (inter and intra) wall thickness is maintained, but there is a reduction in growth rate (linear extension) at all elevated pCO2. Furthermore a decrease in Mg content at the two highest CO2 treatments was observed. Comparison between our data and that at 3 months from the same long-term experiment shows that the acclimation differs over time since at 3 months, the samples cultured under high pCO2 showed a reduction in the cell (inter and intra) wall thickness but a maintained growth rate. This suggests a reallocation of the energy budget between 3 and 10 months and highlights the high degree plasticity that is present. This might provide a selective advantage in future high CO2 world.

Phenotypic and genetic sources of variability of cavitation resistance in Pinus canariensis

Phenotypic and genetic sources of variability of cavitation resistance in Pinus canariensis