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.

Experimental Evidence of an Eco-evolutionary Feedback during Adaptive Divergence

Differences in how organisms modify their environment can evolve rapidly and might influence adaptive population divergence [1, 2]. In a common garden experiment in aquatic mesocosms, we found that adult stickleback from a recently diverged pair of lake and stream populations had contrasting effects on ecosystem metrics. These modifications were caused by both genetic and plastic differences between populations and were sometimes comparable in magnitude to those caused by the presence/ absence of stickleback. Lake and streamfish differentially affected the biomass of zooplankton and phytoplankton, the concentration of phosphorus, and the abundance of several prey (e.g., copepods) and non-prey (e.g., cyanobacteria) species. The adult mediated effects on mesocosm ecosystems influenced the survival and growth of a subsequent generation of juvenile stickleback reared in the same mesocosms. The prior presence of adults decreased the overall growth rate of juveniles, and the prior presence of stream adults lowered overall juvenile survival. Among the survivors, lake juveniles grew faster than co-occurring stream juveniles, except in mesocosm ecosystems previously modified by adult lake fish that were reared on plankton. Overall, our results provide evidence for reciprocal interactions between ecosystem dynamics and evolutionary change (i.e., eco-evolutionary feedbacks) in the early stages of adaptive population divergence.

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.

Low but contrasting neutral genetic differentiation shaped by winter temperature in European great tits

Gene flow is usually thought to reduce genetic divergence and impede local adaptation by homogenising gene pools between populations. However, evidence for local adaptation and phenotypic differentiation in highly mobile species, experiencing high levels of gene flow, is emerging. Assessing population genetic structure at different spatial scales is thus a crucial step towards understanding mechanisms underlying intraspecific differentiation and diversification. Here, we studied the population genetic structure of a highly mobile species – the great tit Parus major – at different spatial scales. We analysed 884 individuals from 30 sites across Europe including 10 close-by sites (< 50 km), using 22 microsatellite markers. Overall we found a low but significant genetic differentiation among sites (FST = 0.008). Genetic differentiation was higher, and genetic diversity lower, in south-western Europe. These regional differences were statistically best explained by winter temperature. Overall, our results suggest that great tits form a single patchy metapopulation across Europe, in which genetic differentiation is independent of geographical distance and gene flow may be regulated by environmental factors via movements related to winter severity. This might have important implications for the evolutionary trajectories of sub-populations, especially in the context of climate change, and calls for future investigations of local differences in costs and benefits of philopatry at large scales.

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.

Phylogenetic analysis of receptor FgfrL1 shows divergence of the C-terminal end in rodents

FGFRL1 is a member of the fibroblast growth factor receptor (FGFR) family. Similar to the classical receptors FGFR1-FGFR4, it contains three extracellular Ig-like domains and a single transmembrane domain. However, it lacks the intracellular tyrosine kinase domain that would be required for signal transduction, but instead contains a short intracellular tail with a peculiar histidine-rich motif. This motif has been conserved during evolution from mollusks to echinoderms and vertebrates. Only the sequences of FgfrL1 from a few rodents diverge at the C-terminal region from the canonical sequence, as they appear to have suffered a frameshift mutation within the histidine-rich motif. This mutation is observed in mouse, rat and hamster, but not in the closely related rodents mole rat (Nannospalax) and jerboa (Jaculus), suggesting that it has occurred after branching of the Muridae and Cricetidae from the Dipodidae and Spalacidae. The consequence of the frameshift is a deletion of a few histidine residues and an extension of the C-terminus by about 40 unrelated amino acids. A similar frameshift mutation has also been observed in a human patient with a craniosynostosis syndrome as well as in several patients with colorectal cancer and bladder tumors, suggesting that the histidine-rich motif is prone to mutation. The reason why this motif was conserved during evolution in most species, but not in mice, is not clear.

Rapid divergence in physiological and life-history traits between northern and southern populations of the British introduced neo-species, Senecio squalidus

Alien plants provide a unique opportunity to study evolution in novel environments, but relatively little is known about the extent to which they become locally adapted to different environments across their new range. Here, we compare northern and southern populations of the introduced species Senecio squalidus in Britain; S. squalidus has been in southern Britain for approximately 200 years and reached Scotland only about 50 years ago. We conducted common garden experiments at sites in the north and south of the species’ range in Britain. We also conducted glasshouse and growth chamber experiments to test the hypothesis that southern genotypes flower later, are more drought-tolerant, germinate and establish better at warmer temperatures, and are less sensitive to cold stress than their more northern counterparts. Results from the common garden experiments are largely consistent with the hypothesis of rapid adaptive divergence of populations of the species within the introduced range, with genotypes typically showing a home-site advantage. Results from the glasshouse and growth chamber experiments demonstrate adaptive divergence in ability to tolerate drought stress and high temperatures, as well as in phenology. In particular, southern genotypes were more tolerant of dry conditions and high temperatures and they flowered later than northern genotypes. Our results show that rapid local adaptation can occur in alien species, and they have implications for our understanding of the ecological genetics of range expansion of introduced weeds.

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.

On the theory of high and low pressure areas: The significance of divergence in pressure areas

This is the edited and translated version of the article by Richard Scherhag “Zur Theorie der Hoch- und Tiefdruckgebiete. Die Bedeutung der Divergenz in Druckfeldern” (On the theory of high and low pressure areas: The significance of divergence in pressure areas), which was published in Meteorologische Zeitschrift 51, 129–138.

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 ^