Computational Techniques for Haplotype Inference and for Local Alignment Significance

This thesis which consists of an introduction and four peer-reviewed original publications studies the problems of haplotype inference (haplotyping) and local alignment significance. The problems studied here belong to the broad area of bioinformatics and computational biology. The presented solutions are computationally fast and accurate, which makes them practical in high-throughput sequence data analysis. Haplotype inference is a computational problem where the goal is to estimate haplotypes from a sample of genotypes as accurately as possible. This problem is important as the direct measurement of haplotypes is difficult, whereas the genotypes are easier to quantify. Haplotypes are the key-players when studying for example the genetic causes of diseases. In this thesis, three methods are presented for the haplotype inference problem referred to as HaploParser, HIT, and BACH. HaploParser is based on a combinatorial mosaic model and hierarchical parsing that together mimic recombinations and point-mutations in a biologically plausible way. In this mosaic model, the current population is assumed to be evolved from a small founder population. Thus, the haplotypes of the current population are recombinations of the (implicit) founder haplotypes with some point--mutations. HIT (Haplotype Inference Technique) uses a hidden Markov model for haplotypes and efficient algorithms are presented to learn this model from genotype data. The model structure of HIT is analogous to the mosaic model of HaploParser with founder haplotypes. Therefore, it can be seen as a probabilistic model of recombinations and point-mutations. BACH (Bayesian Context-based Haplotyping) utilizes a context tree weighting algorithm to efficiently sum over all variable-length Markov chains to evaluate the posterior probability of a haplotype configuration. Algorithms are presented that find haplotype configurations with high posterior probability. BACH is the most accurate method presented in this thesis and has comparable performance to the best available software for haplotype inference. Local alignment significance is a computational problem where one is interested in whether the local similarities in two sequences are due to the fact that the sequences are related or just by chance. Similarity of sequences is measured by their best local alignment score and from that, a p-value is computed. This p-value is the probability of picking two sequences from the null model that have as good or better best local alignment score. Local alignment significance is used routinely for example in homology searches. In this thesis, a general framework is sketched that allows one to compute a tight upper bound for the p-value of a local pairwise alignment score. Unlike the previous methods, the presented framework is not affeced by so-called edge-effects and can handle gaps (deletions and insertions) without troublesome sampling and curve fitting.

Inventories and significance of the genetic resources of an African mahogany species (Khaya senegalensis (Desr.) A. Juss.) assembled and further developed in Australia.

The forest tree species Khaya senegalensis (Desr.) A. Juss. occurs in a belt across 20 African countries from Senegal-Guinea to Sudan-Uganda where it is a highly important resource. However, it is listed as Vulnerable (IUCN 2015-3). Since introduction in northern Australia around 1959, the species has been planted widely, yielding high-value products. The total area of plantations of the species in Australia exceeds 15,000 ha, mostly planted in the Northern Territory since 2006, and includes substantial areas across 60-70 woodlots and industrial plantations established in north-eastern Queensland since the early-1990s and during 2005-2007 respectively. Collaborative conservation and tree improvement by governments began in the Northern Territory and Queensland in 2001 based on provenance and other trials of the 1960s-1970s. This work has developed a broad base of germplasm in clonal seed orchards, hedge gardens and trials (clone and progeny). Several of the trials were established collaboratively on private land. Since the mid-2000s, commercial growers have introduced large numbers of provenance-bulk and individual-tree seedlots to establish industrial plantations and trials, several of the latter in collaboration with the Queensland Government. Provenance bulks (>140) and families (>400) from 17 African countries are established in Australia, considered the largest genetic base of the species in a single country outside Africa. Recently the annual rate of industrial planting of the species in Australia has declined, and R&D has been suspended by governments and reduced by the private sector. However, new commercial plantings in the Northern Territory and Queensland are proposed. In domesticating a species, the strategic importance of a broad genetic base is well known. The wide range of first- and advanced-generation germplasm of the species established in northern Australia and documented in this paper provides a sound basis for further domestication and industrial plantation and woodlot expansion, when investment conditions are favourable

Molecular Breeding for Complex Adaptive Traits: How Integrating Crop Ecophysiology and Modelling Can Enhance Efficiency

Progress in crop improvement is limited by the ability to identify favourable combinations of genotypes (G) and management practices (M) in relevant target environments (E) given the resources available to search among the myriad of possible combinations. To underpin yield advance we require prediction of phenotype based on genotype. In plant breeding, traditional phenotypic selection methods have involved measuring phenotypic performance of large segregating populations in multi-environment trials and applying rigorous statistical procedures based on quantitative genetic theory to identify superior individuals. Recent developments in the ability to inexpensively and densely map/sequence genomes have facilitated a shift from the level of the individual (genotype) to the level of the genomic region. Molecular breeding strategies using genome wide prediction and genomic selection approaches have developed rapidly. However, their applicability to complex traits remains constrained by gene-gene and gene-environment interactions, which restrict the predictive power of associations of genomic regions with phenotypic responses. Here it is argued that crop ecophysiology and functional whole plant modelling can provide an effective link between molecular and organism scales and enhance molecular breeding by adding value to genetic prediction approaches. A physiological framework that facilitates dissection and modelling of complex traits can inform phenotyping methods for marker/gene detection and underpin prediction of likely phenotypic consequences of trait and genetic variation in target environments. This approach holds considerable promise for more effectively linking genotype to phenotype for complex adaptive traits. Specific examples focused on drought adaptation are presented to highlight the concepts.

Optimal adaptive sequential spatial sampling of soil using pair-copulas

A spatial sampling design that uses pair-copulas is presented that aims to reduce prediction uncertainty by selecting additional sampling locations based on both the spatial configuration of existing locations and the values of the observations at those locations. The novelty of the approach arises in the use of pair-copulas to estimate uncertainty at unsampled locations. Spatial pair-copulas are able to more accurately capture spatial dependence compared to other types of spatial copula models. Additionally, unlike traditional kriging variance, uncertainty estimates from the pair-copula account for influence from measurement values and not just the configuration of observations. This feature is beneficial, for example, for more accurate identification of soil contamination zones where high contamination measurements are located near measurements of varying contamination. The proposed design methodology is applied to a soil contamination example from the Swiss Jura region. A partial redesign of the original sampling configuration demonstrates the potential of the proposed methodology.

Fourier transform of the collagen triple-helical structure and its significance

The Fourier transforms of the collagen molecular structure have been calculated taking into consideration various side chain atoms, as well as the presence of bound water molecules. There is no significant change in the calculated intensity distribution on including the side chain atoms of non-imino-acid residues. Taking into account the presence of about two bound water molecules per tripeptide unit, the agreement with the observed x-ray pattern is slightly improved. Fourier transforms have also been calculated for the detailed molecular geometries proposed from other laboratories. It is found that there are no major differences between them, as compared to our structure, either in the positions of peak intensity or in the intensity distribution. Hence it is not possible to judge the relative merits of the various molecular geometries for the collagen triple helix from a comparison of the calculated transforms with the meagre data available from its x-ray fibre pattern. It is also concluded that the collagen molecular structure should be regarded as a somewhat flexible chain structure, capable of adapting itself to the requirements of the different side groups which occur in each local region.

Assessing the significance of climate and community factors on urban water demand

Ensuring adequate water supply to urban areas is a challenging task due to factors such as rapid urban growth, increasing water demand and climate change. In developing a sustainable water supply system, it is important to identify the dominant water demand factors for any given water supply scheme. This paper applies principal components analysis to identify the factors that dominate residential water demand using the Blue Mountains Water Supply System in Australia as a case study. The results show that the influence of community intervention factors (e.g. use of water efficient appliances and rainwater tanks) on water demand are among the most significant. The result also confirmed that the community intervention programmes and water pricing policy together can play a noticeable role in reducing the overall water demand. On the other hand, the influence of rainfall on water demand is found to be very limited, while temperature shows some degree of correlation with water demand. The results of this study would help water authorities to plan for effective water demand management strategies and to develop a water demand forecasting model with appropriate climatic factors to achieve sustainable water resources management. The methodology developed in this paper can be adapted to other water supply systems to identify the influential factors in water demand modelling and to devise an effective demand management strategy.

Significance of Branching Angles in Streamer Discharges

The phenomenon of branching at specific angles in streamer breakdown studies is found to be more universal than was ever thought. The angles measured in the breakdown of gases show that the coefficient of field distortion, K, lies in the range from 1 to less than 0.1. The values of K, so obtained agree well with those envisaged in the criterion of the streamer mechanism. It is hoped that branching angles observed in various types of breakdown may, possibly, be explained similarly.

Pathogenicity, functional significance and clinical phenotype of mismatch repair gene MSH2 variants found in cancer patients

DNA ja siinä sijaitsevat geenit ohjaavat kaikkea solujen toimintaa. DNA-molekyyleihin kuitenkin kertyy mutaatioita sekä ympäristön vaikutuksen, että solujen oman toiminnan tuloksena. Mikäli virheitä ei korjata, saattaa tuloksena olla solun muuttuminen syöpäsoluksi. Soluilla onkin käytössä useita DNA-virheiden korjausmekanismeja, joista yksi on ns. mismatch repair (MMR). MMR vastaa DNA:n kahdentumisessa syntyvien virheiden korjauksesta. Periytyvät mutaatiot geeneissä, jotka vastaavat MMR-proteiinien rakentamisesta, aiheuttavat ongelmia DNA:n korjauksessa ja altistavat kantajansa periytyvälle ei-polypoottiselle paksusuolisyöpäoireyhtymälle (hereditary nonpolyposis colorectal cancer, HNPCC). Yleisimmin mutatoituneet MMR-geenit ovat MLH1 ja MSH2. HNPCC periytyy vallitsevasti, eli jo toiselta vanhemmalta peritty geenivirhe altistaa syövälle. MMR-geenivirheen kantaja sairastuu syöpään elämänsä aikana suurella todennäköisyydellä, ja sairastumisikä on vain noin 40 vuotta. Syövälle altistavan geenivirheen löytäminen mutaation kantajilta on hyvin tärkeää, sillä säännöllinen seuranta mahdollistaa kehittymässä olevan kasvaimen havaitsemisen ja poistamisen jo aikaisessa vaiheessa. Tämän on osoitettu alentavan syöpäkuolleisuutta merkittävästi. Varma tieto altistuksen alkuperästä on tärkeä myös niille syöpäsuvun jäsenille, jotka eivät kanna kyseistä mutaatiota. Syövälle altistavien mutaatioiden ohella MMR-geeneistä löydetään säännöllisesti muutoksia, jotka ovat normaalia henkilöiden välistä geneettistä vaihtelua, eikä niiden oleteta lisäävän syöpäaltistusta. Altistavien mutaatioiden erottaminen näistä neutraaleista variaatioista on vaikeaa, mutta välttämätöntä altistuneiden tehokkaan seurannan varmistamiseksi. Tässä väitöskirjassa tutkittiin 18:a MSH2 -geenin mutaatiota. Mutaatiot oli löydetty perheistä, joissa esiintyi paljon syöpiä, mutta niiden vaikutus DNA:n korjaustehoon ja syöpäaltistukseen oli epäselvä. Työssä tutkittiin kunkin mutaation vaikutusta MSH2-proteiinin normaaliin toimintaan, ja tuloksia verrattiin potilaiden ja sukujen kliinisiin tietoihin. Tutkituista mutaatiosta 12 aiheutti puutteita MMR-korjauksessa. Nämä mutaatiot tulkittiin syövälle altistaviksi. Analyyseissä normaalisti toimineet 4 mutaatiota eivät todennäköisesti ole syynä syövän syntyyn kyseisillä perheillä. Tulkinta jätettiin avoimeksi 2 mutaation kohdalla. Tutkimuksesta hyötyivät suoraan kuvattujen mutaatioiden kantajaperheet, joiden geenivirheen syöpäaltistuksesta saatiin tietoa, mahdollistaen perinnöllisyysneuvonnan ja seurannan kohdentamisen sitä tarvitseville. Työ selvensi myös mekanismeja, joilla mutatoitunut MSH2-proteiini voi menettää toimintakykynsä.

Phylogeography and Adaptive Divergence of Three-spined Stickleback Populations

Genetic studies on phylogeography and adaptive divergence in Northern Hemisphere fish species such as three-spined stickleback (Gasterosteus aculeatus) provide an excellent opportunity to investigate genetic mechanisms underlying population differentiation. According to the theory, the process of population differentiation results from a complex interplay between random and deterministic processes as well historical factors. The main scope in this thesis was to study how historical factors like the Pleistocene ice ages have shaped the patterns molecular diversity in three-spined stickleback populations in Europe and how this information could be utilized in the conservation genetic context. Furthermore, identifying footprints of natural selection at the DNA level might be used in identifying genes involved in evolutionary change. Overall, the results from phylogeographic studies indicate that the three-spined stickleback has colonized the Atlantic basin relatively recently but constitutes three major evolutionary lineages in Europe. In addition, the colonization of freshwater appears to result from multiple and independent invasions by the marine conspecifics. Molecular data together with morphology suggest that the most divergent freshwater populations are located in the Balkan Peninsula and these populations deserve a special conservation genetic status without warranting further taxonomical classification. In order to investigate the adaptive divergence in Fennoscandian three-spined stickleback populations several approaches were used. First, sequence variability in the Eda-gene, coding for the number of lateral plates, was concordant with the previously observed global pattern. Full plated allele is in high frequencies among marine populations whereas low plated allele dominates in the freshwater populations. Second, a microsatellite based genome scan identified both indications of balancing and directional selection in the three-spined stickleback genome, i.e. loci with unusually similar or unusually different allele frequencies over populations. The directionally selected loci were mainly associated with the adaptation to freshwater. A follow up study conducting a more detailed analysis in a chromosome region containing a putatively selected gene locus identified a fairly large genomic region affected by natural selection. However, this region contained several gene predictions, all of which might be the actual target of natural selection. All in all, the phylogeographic and adaptive divergence studies indicate that most of the genetic divergence has occurred in the freshwater populations whereas the marine populations have remained relatively uniform.

trans gene regulation in adaptive evolution: A genetic algorithm model

This is a continuation of earlier studies on the evolution of infinite populations of haploid genotypes within a genetic algorithm framework. We had previously explored the evolutionary consequences of the existence of indeterminate—“plastic”—loci, where a plastic locus had a finite probability in each generation of functioning (being switched “on”) or not functioning (being switched “off”). The relative probabilities of the two outcomes were assigned on a stochastic basis. The present paper examines what happens when the transition probabilities are biased by the presence of regulatory genes. We find that under certain conditions regulatory genes can improve the adaptation of the population and speed up the rate of evolution (on occasion at the cost of lowering the degree of adaptation). Also, the existence of regulatory loci potentiates selection in favour of plasticity. There is a synergistic effect of regulatory genes on plastic alleles: the frequency of such alleles increases when regulatory loci are present. Thus, phenotypic selection alone can be a potentiating factor in a favour of better adaptation.

Differentiating between pollutants build-up on roads and roofs: Significance of roofs as a stormwater pollutant source

Impervious surfaces in an urban catchment are the primary stormwater pollutant contributing areas. Appropriate treatment of stormwater runoff from these impervious surfaces is essential to safeguard the urban water environment. While urban roads have received significant research attention in this regard, roofs have not been well investigated. Key pollutant processes such as build-up on roads and roofs can be different due to the different surface characteristics. This entails different treatment strategies being needed for road and roofs. The research study characterized roof pollutants build-up by differentiating with road surfaces. It was noted that pollutants are more highly concentrated on particles and particularly finer particles in the case of roof surfaces, compared to road surfaces. Additionally, pollutants built-up on roof surfaces tend to be relatively more variable from one day to another in terms of pollutant loads. These results highlight the significance of roofs as a stormwater pollutant source and the important need for a specific stormwater treatment strategy rather than the application of a combined approach for treating stormwater runoff from both, roads and roofs.

Behavioural and physiological responses to predators of captive-bred Arctic charr: significance of genetics, learning and ontogeny

Predation forms one of the main selective forces in nature and in a vast number of prey species the behavioural responses form the main way to avoid predation. World wide numerous captive breeding programs are used to produce fish and other animal species for conservational reintroductions. However, rearing animals in the absence of predators in captivity has been shown to weaken their predator avoidance skills and lead to behavioural divergence between wild and captive-bred populations. In my thesis I studied the effects of predator odour exposures on antipredator behavioural and physiological responses of captive reared Saimaa Arctic charr. This charr population is the most endangered fish population in Finland and a sample of the remaining population has been taken to captive breeding and used for an extensive reintroduction program. Lowered responsiveness to predators is probably one of the major reasons for the poor survival probability of the charr after release into the wild. The main aims of my thesis were to explore the reasons for behavioural phenotypic variation in this charr population and whether naïve charr young could be trained to recognise their natural predators. The predator species in my thesis were burbot (Lota lota) and pikeperch (Sander lucioperca). In my thesis I showed that the captive-bred charr responded to chemical cues from burbot and pikeperch, but the magnitude of responses was linked to the predator species. The burbot odour increased the spatial odour avoidance of the charr young. On the other hand, in the pikeperch treatment charr reduced their relative swimming activity and tended to show more freezing behaviour relative to the burbot treatment. It seems evident that these different responses are related to the different hunting tactics of predator species. Furthermore, I detected wide between-family differences in antipredator responsiveness (i.e. inherited variation in antipredator behaviours) in this captive stock. Detected differences were greater in the response towards burbot than towards pikeperch. These results, in addition to predator-specific antipredator responses, suggest that there is a clear inherited component in antipredator responsiveness in Saimaa charr population and that the detected inherited differences could explain a part of the behavioural phenotypic variation in this population. In my thesis I also found out that both social learning and direct exposure to live predators enhance the antipredator responsiveness of charr young. In addition, I obtained indications that predator odour exposures (i.e. life-skills training) in alevin and fry stages can fine-tune the innate antipredator responsiveness of charr. Thus, all these methods have the potential to enhance the innate antipredator responsiveness of naïve charr young, possibly also improving the post-release survival of these trained individuals in the wild. However, the next logical phase would be to carry out large scale survival studies in the wild to test this hypothesis. Finally, the results of my thesis emphasize that possible long-term life-skills training methods should take into account not only the behavioural but also the physiological effects of training.

A Control Strategy for Reference Wave Adaptive Current Generation

Speed control of ac motors requires variable frequency, variable current, or variable voltage supply. Variable frequency supply can be obtained directly from a fixed frequency supply by using a frequency converter or from a dc source using inverters. In this paper a control technique for reference wave adaptive-current generation by modulating the inverter voltage is explained. Extension of this technique for three-phase induction-motor speed control is briefly explained. The oscillograms of the current waveforms obtained from the experimental setup are also shown.

The tectonic significance of lower Permain successions in the Texas Orocline (Eastern Australia)

The Texas Orocline is a prominent orogenic curvature that developed during the early Permian in the southern New England Orogen. Outliers preserving lower Permian sedimentary successions (Bondonga, Silver Spur, Pikedale, Terrica, Alum Rock and Ashford beds) approximately outline the oroclinal structure, but the tectonic processes responsible for the development of these basinal successions, and their relationships to the Texas Orocline, are unclear. Here we address this shortcoming by providing new U–Pb detrital and primary zircon ages from these successions, as well as detailed stratigraphic and structural data from the largest exposed succession (Bondonga beds). Field observations and U–Pb geochronological data suggest that the lower Permian successions in the Texas Orocline are remnants of a single, formerly larger basin that was deposited after ca 302 Ma. Time constraints for formation of this basin are correlative with constraints from the lower Permian Nambucca Block, which was likely deposited in response to regional back-arc extension during and/or after the development of the Texas Orocline. The conclusion that the lower Permian sedimentary basins in the Texas Orocline belong to this back-arc extensional system supports the suggestion that oroclinal bending in the New England Orogen was primarily controlled by trench retreat and associated overriding-plate extension.