Sample size considerations and augmentation of computer experiments

Computer Experiments, consisting of a number of runs of a computer model with different inputs, are now common-place in scientific research. Using a simple fire model for illustration some guidelines are given for the size of a computer experiment. A graph is provided relating the error of prediction to the sample size which should be of use when designing computer experiments. Methods for augmenting computer experiments with extra runs are also described and illustrated. The simplest method involves adding one point at a time choosing that point with the maximum prediction variance. Another method that appears to work well is to choose points from a candidate set with maximum determinant of the variance covariance matrix of predictions.

Some practical issues in the design and analysis of computer experiments

Deterministic computer simulations of physical experiments are now common techniques in science and engineering. Often, physical experiments are too time consuming, expensive or impossible to conduct. Complex computer models or codes, rather than physical experiments lead to the study of computer experiments, which are used to investigate many scientific phenomena of this nature. A computer experiment consists of a number of runs of the computer code with different input choices. The Design and Analysis of Computer Experiments is a rapidly growing technique in statistical experimental design. This thesis investigates some practical issues in the design and analysis of computer experiments and attempts to answer some of the questions faced by experimenters using computer experiments. In particular, the question of the number of computer experiments and how they should be augmented is studied and attention is given to when the response is a function over time.

The effect of different levels of inbreeding on self-incompatibility and inbreeding depression in pineapple

A high level of heterozygosity in domesticated pineapple is one of the main obstacles hindering the efficient development of new varieties. Selfing has been proposed as a strategy to minimise this heterozygosity through the development of parentals with a greater level of homozygosity. Selfing and a range of lesser levels of inbreeding were evaluated for their effects on seed production and inbreeding depression in the early growth of seedlings. Selfing produced few seeds, and very few viable seedlings. The paternal backcross and several half-sib combinations exhibited minimal effects on seed development or early growth and in some cases were similar to the outcross. Sibcrosses were generally unsuccessful. The highest inbreeding coefficient that was not associated with severe inbreeding depression was approximately 0.25. The effect of inbreeding depression and the level of homozygosity on several quantitative traits including those related to fruit quality within the inbred populations, is now being assessed.

Catalyzed air oxidation for treatment of thiosalt effulents

Continuous slurry reactor runs of two to four weeks duration were carried out for catalyzed air oxidation of thiosalts under a variety of conditions using poly (4-vinylpyridine) - Cu (II) and quaternized poly (4-vinylpyridine) - Cu (II) catalysts. Results obtained indicate that these catalysts have high activity and relatively long-term catalyst stability for thiosalt waste streams of < 1000 ppm thiosalt level. Using 2% (w/w) slurries of the poly (4-vinylpyridine) Cu (II) catalyst, effective oxidation of 700 ppm S2O32− influent to an effluent of < 100 ppm total thio-salts can be carried out continuously for at least one month when operating at 20 to 30°C with solution flow rates of$˜1l/h and aeration of 1300 XXX/h using a two-stage reactor system comprised of 12 l reactors. At higher thiosalt influent levels (i.e. > 1600 ppm) increased reaction temperatures enable depletion to < 100 ppm thiosalt effluent levels for up to one week of continuous operation. The catalysts deactivate much more readily at these higher influent levels as a result of greater copper losses and appreciable adsorption of S2O32− and S4O62−. The behaviour of continuous slurry reactors employed in the experimental studies, by use of batch reaction data for the poly (4-vinylpyridine) Cu (II) catalyzed oxidation of thiosalts, can be modelled successfully. Quaternized poly (4-vinylpyridine) Cu (II) catalyst has good long-term stability and copper losses are very low. The poly (4-vinylpyridine) Cu (II) catalyst, however, is susceptible to appreciable oxidation of the polymer matrix on long-term usage. This oxidation of the polymer matrix results in a substantial loss in the activity of the regenerated catalyst.

NUVIEW:software for display and interactive manipulation of nucleic acid models

The NUVIEW software package allows skeletal models of any double helical nucleic acid molecule to be displayed out a graphics monitor and to apply various rotations, translations and scaling transformations interactively, through the keyboard. The skeletal model is generated by connecting any pair of representative points, one from each of the bases in the basepair. In addition to the above mentioned manipulations, the base residues can be identified by using a locator and the distance between any pair of residues can be obtained. A sequence based color coded display allows easy identification of sequence repeats, such as runs of Adenines. The real time interactive manipulation of such skeletal models for large DNA/RNA double helices, can be used to trace the path of the nucleic acid chain in three dimensions and hence get a better idea of its topology, location of linear or curved regions, distances between far off regions in the sequence etc. A physical picture of these features will assist in understanding the relationship between base sequence, structure and biological function in nucleic acids.

Consequences of Supersymmetry in the Early Universe

Currently, we live in an era characterized by the completion and first runs of the LHC accelerator at CERN, which is hoped to provide the first experimental hints of what lies beyond the Standard Model of particle physics. In addition, the last decade has witnessed a new dawn of cosmology, where it has truly emerged as a precision science. Largely due to the WMAP measurements of the cosmic microwave background, we now believe to have quantitative control of much of the history of our universe. These two experimental windows offer us not only an unprecedented view of the smallest and largest structures of the universe, but also a glimpse at the very first moments in its history. At the same time, they require the theorists to focus on the fundamental challenges awaiting at the boundary of high energy particle physics and cosmology. What were the contents and properties of matter in the early universe? How is one to describe its interactions? What kind of implications do the various models of physics beyond the Standard Model have on the subsequent evolution of the universe? In this thesis, we explore the connection between in particular supersymmetric theories and the evolution of the early universe. First, we provide the reader with a general introduction to modern day particle cosmology from two angles: on one hand by reviewing our current knowledge of the history of the early universe, and on the other hand by introducing the basics of supersymmetry and its derivatives. Subsequently, with the help of the developed tools, we direct the attention to the specific questions addressed in the three original articles that form the main scientific contents of the thesis. Each of these papers concerns a distinct cosmological problem, ranging from the generation of the matter-antimatter asymmetry to inflation, and finally to the origin or very early stage of the universe. They nevertheless share a common factor in their use of the machinery of supersymmetric theories to address open questions in the corresponding cosmological models.

Future collider signatures of the possible 750 GeV state

If the recent indications of a possible state I broken vertical bar with mass similar to 750 GeV decaying into two photons reported by ATLAS and CMS in LHC collisions at 13 TeV were to become confirmed, the prospects for future collider physics at the LHC and beyond would be affected radically, as we explore in this paper. Even minimal scenarios for the I broken vertical bar resonance and its gamma gamma decays require additional particles with masses . We consider here two benchmark scenarios that exemplify the range of possibilities: one in which I broken vertical bar is a singlet scalar or pseudoscalar boson whose production and gamma gamma decays are due to loops of coloured and charged fermions, and another benchmark scenario in which I broken vertical bar is a superposition of (nearly) degenerate CP-even and CP-odd Higgs bosons in a (possibly supersymmetric) two-Higgs doublet model also with additional fermions to account for the gamma gamma decay rate. We explore the implications of these benchmark scenarios for the production of I broken vertical bar and its new partners at colliders in future runs of the LHC and beyond, at higher-energy pp colliders and at e (+) e (-) and gamma gamma colliders, with emphasis on the bosonic partners expected in the doublet scenario and the fermionic partners expected in both scenarios.

TUNPOP, a computer simulation model of the yellowfin tuna population and the surface tuna fishery of the Eastern Pacific Ocean

ENGLISH: Mathematical documentation of TUNP0P, an age-structured computer simulation model of the yellowfin tuna population and surface tuna fishery of the eastern Pacific Ocean, is described. Example runs of the model are presented and discussed, and the sensitivity of the model output to changes in various parameters is examined. SPANISH: Se describe la documentación matemática de TUNP0P, un modelo computador de simulación basado en la edad de la población del atún aleta amarilla y de la pesca atunera epipelágíca del Océano Pacífico oriental. Se presentan y se discuten ejemplos de las pasadas del modelo, y se examina la sensibilidad de los resultados de salida con relación a los cambios de varios parámetros. (PDF contains 47 pages.)

Modeling the effect of striped bass (Morone saxatilis) on the population viability of Sacramento River winter-run chinook salmon (Onchorhynchus tshawytscha)

We estimated the impact of striped bass (Morone saxatilis) predation on winter-run chinook salmon (Oncorhynchus tshawytscha) with a Bayesian population dynamics model using striped bass and winter-run chinook salmon population abundance data. Winter-run chinook salmon extinction and recovery probabilities under different future striped bass abundance levels were estimated by simulating from the posterior distribution of model parameters. The model predicts that if the striped bass population declines to 512,000 adults as expected in the absence of stocking, winter-run chinook salmon will have about a 28% chance of quasi-extinction (defined as three consecutive spawning runs of fewer than 200 adults) within 50 years. If stocking stabilizes the striped bass population at 700,000 adults, the predicted quasi-extinction probability is 30%. A more ambitious stocking program that maintains a population of 3 million adult striped bass would increase the predicted quasi-extinction probability to 55%. Extinction probability, but not recovery probability, was fairly insensitive to assumptions about density dependence. We conclude that winter-run chinook salmon face a serious extinction risk without augmentation of the striped bass population and that substantial increases in striped bass abundance could significantly increase the threat to winter-run chi-nook salmon if not mitigated by increasing winter chinook salmon survival in some other way.

Acute toxicity of unionized ammonia to milkfish (Chanos chanos Forsskal) fry

The study was conducted to determine the effects of varying concentrations of ammonia to milkfish fry. Two runs of static 96h bioassays were conducted to determine the median lethal concentration (LC 50) of unionized ammonia (NH₃) to milkfish fry. Test concentrations were based on exploratory 24h and 48h bioassays and were made in three replicates. Reagent grade ammonium chloride (NH₄Cl) was used to adjust the level of unionized ammonia. The 96h median lethal concentration, determined by the Reed Muench method was calculated at 28.029 ppm NH₃ 29.69 ppm. Even at high concentrations of unionized ammonia, most of the fry mortality occurred after 48 to 96 hours exposure. Severe gill damage occurs only at concentrations above 20 ppm, especially above the LC 50. The high LC 50 value obtain shows that milkfish fry has great tolerance to ammonia, that even fry with severely-damaged gills can still recover days after it is returned to favorable culture condition. The result suggest that observed mortalities of milkfish fry under culture conditions are not due to ammonia toxicity.

Characterization of the main light-harvesting chlorophyll a/b-protein complex of green alga, Bryopsis corticulans

The main light-harvesting chlorophyll a/b -protein complex (LHC II) has been isolated directly from thylakoid membranes of shiphonous green alga, Bryopsis corticulans Setch. by using two consecutive runs of anion exchange and gel-filtration chromatography. Monomeric and trimeric subcomplexes of LHC 11 were obtained by using sucrose gradient ultracentrifugation. Pigment analysis by reversed-phase high performance liquid chromatography showed that chlorophyll a (Chl a), chlorophyll b (Chl b), neoxanthin, violaxanthin and siphonaxanthin were involved in LHC 11 from B. corticulans. The properties of electronic transition of monomeric LHC II showed similarities to those of trimeric LHC II. Circular dichroism spectroscopy showed that strong intramolecular interaction of excitonic dipoles between Chl a and between Chl b exist in one LHC II apoprotein, while the intermolecular interaction of these dipoles can be intensified in the trimeric structure. The monomer has high efficient energy transfer from Chl b and siphonaxanthin to Chl a similarly to that of the trimer. Our results suggest that in B. corticulans, LHC II monomer has high ordered pigment organization that play effective physiological function as the trimer, and thus it might be also a functional organization existing in thylakoid membrane of B. corticulans.

Isolation of the main light-harvesting chlorophyll a/b-protein complex from thylakoid membranes of marine alga, Bryopsis corticulans by a direct method

The main chlorophyll a/b light-harvesting complex (LHC 11) has been isolated directly from thylakoid membranes of marine green alga (Bryopsis corticulans Setch.) by two consecutive runs of anion exchange and gel-filtration chromatography. LHC 11 proteins in the membrane extracts treated with 3% n-Octyl-b-D-glucopyranoside (OG) obtained specific binding ability on Q Sepharose column, and thus were isolated from the thylakoid membranes in a highly selective fraction. The monomeric, trimeric and oligomeric subcomplexes of LHC 11 have been obtained by fractionation of the LHC 11 mixes with sucrose density gradient ultracentrifugation. The SDS-PAGE analysis of peptide composition and absorption spectrum showed that LHC 11 monomers, trimers and oligomers prepared through this work were intact and in high purity. Our report is the first to show that it is possible to purify LHC If directly from thylakoid membranes without extensively biochemical purification.

Molecular basis of complete C4 deficiency. A study of three patients.

The highly polymorphic fourth component of human complement (C4) is usually encoded by two genes, C4A and C4B, adjacent to the 21-hydroxylase (21-OH) genes and is also remarkable by the high frequency of the null alleles, C4A*Q0 and C4B*Q0. Complete C4 deficiency is exceptional because this condition appears only in homozygotes for the very rare double-null haplotype C4AQ0,BQ0. This condition in most cases gives rise to systemic lupus erythematosus and an increased susceptibility to infections. The molecular basis for complete C4 deficiency has not yet been established. Therefore we studied the DNA of three previously described C4 deficient patients belonging to unrelated families by restriction fragment length polymorphism analysis using C4 and 21-OH probes. These studies revealed a deletion of the C4B and 21-OHA genes in two patients and no deletion at all in the third patient. Therefore, complete C4 deficiency as a result of homozygosity for the C4AQ0, BQ0 haplotype is not a consequence of a deletion of the C4 genes. The molecular basis of this genetic abnormality is certainly very complex and may vary also from one case to another.

Characterisation of a novel panel of polymorphic microsatellite loci for the liver fluke, Fasciola hepatica, using a next generation sequencing approach

The liver fluke, Fasciola hepatica is an economically important pathogen of sheep and cattle and has been described by the WHO as a re-emerging zoonosis. Control is heavily reliant on the use of drugs, particularly triclabendazole and as a result resistance has now emerged. The population structure of F. hepatica is not well known, yet it can impact on host-parasite interactions and parasite control with drugs, particularly regarding the spread of triclabendazole resistance. We have identified 2448 potential microsatellites from 83Mb of F. hepatica genome sequence using msatfinder. Thirty-five loci were developed and optimised for microsatellite PCR, resulting in a panel of 15 polymorphic loci, with a range of three to 15 alleles. This panel was validated on genomic DNA from 46 adult F. hepatica; 38 liver flukes sourced from a Northwest abattoir, UK and 8 liver flukes from an established isolate (Shrewsbury; Ridgeway Research). Evidence for null alleles was found at four loci (Fh_1, Fh_8, Fh_13 and Fh_14), which showed markedly higher levels of homozygosity than the remaining 11 loci. Of the 38 liver flukes isolated from cattle livers (n=10) at the abattoir, 37 genotypes were identified. Using a multiplex approach all 15 loci could be amplified from several life cycle stages that typically yield low amounts of DNA, including metacercariae, the infective life cycle stage present on pasture, highlighting the utility of this multiplex microsatellite panel. This study reports the largest panel of microsatellite markers available to date for population studies of F. hepatica and the first multiplex panel of microsatellite markers that can be used for several life cycle stages.

Genetic Landscape of Joubert syndrome in French Canadians

Le syndrome de Joubert est une maladie récessive caractérisée par une malformation congénitale distincte du tronc cérébral et du cervelet, associée à une anomalie des mouvements oculaires (apraxie oculomotrice), une respiration irrégulière, un retard de développement, et une ataxie à la démarche. Au cours de la dernière décennie, plus de 20 gènes responsables ont été identifiés, tous ayant un rôle important dans la structure et la fonction des cils primaires. Ainsi, le syndrome de Joubert est considéré une ciliopathie. Bien que le Syndrome de Joubert ait été décrit pour la première fois dans une famille canadienne-française en 1969, le(s) gène(s) causal demeurait inconnu dans presque tous les cas de syndrome de Joubert recensés en 2010 dans la population canadienne-française, soit début de mon projet doctoral. Nous avons identifié un total de 43 individus canadiens-français (35 familles) atteints du syndrome de Joubert. Il y avait un regroupement de familles dans la région du Bas-Saint-Laurent de la province de Québec, suggérant la présence d'un effet fondateur. L’objectif de ce projet était de caractériser la génétique du syndrome de Joubert dans la population canadienne-française. Notre hypothèse était qu’il existait un effet fondateur impliquant au moins un nouveau gène JBTS. Ainsi, dans un premier temps, nous avons utilisé une approche de cartographie par homozygotie. Cependant, nous n’avons pas identifié de région d’homozygotie partagée parmi les individus atteints, suggérant la présence d’une hétérogénéité génétique ou allélique. Nous avons donc utilisé le séquençage exomique chez nos patients, ce qui représente une approche plus puissante pour l’étude de conditions génétiquement hétérogènes. Nos travaux ont permis l’identification de deux nouveaux gènes responsables du syndrome de Joubert: C5orf42 et TMEM231. Bien que la localisation cellulaire et la fonction de C5orf42 soient inconnus au moment de cette découverte, nos résultats génétiques combinés avec des études ultérieures ont établi un rôle important de C5orf42 dans la structure et la fonction ciliaire, en particulier dans la zone de transition, qui est une zone de transition entre le cil et le reste de la cellule. TMEM231 avait déjà un rôle établi dans la zone de transition ciliaire et son interaction avec d’autres protéines impliquées dans le syndrome de Joubert était connu. Nos études ont également identifié des variants rares délétères chez un patient JBTS dans le gène ciliaire CEP104. Nous proposons donc CEP104 comme un gène candidat JBTS. Nous avons identifié des mutations causales dans 10 gènes, y compris des mutations dans CC2D2A dans 9 familles et NPHP1 dans 3 familles. Au total, nous avons identifié les mutations causales définitives chez 32 des 35 familles étudiées (91% des cas). Nous avons documenté un effet fondateur complexe dans la population canadienne-française avec de multiples mutations récurrentes dans quatre gènes différents (C5orf42, CC2D2A, TMEM231, NPHP1). Au début de ce projet de recherche, l’étiologie génétique était inconnue chez les 35 familles touchées du syndrome de Joubert. Maintenant, un diagnostique moléculaire définitif est identifié chez 32 familles, et probable chez les 3 autres. Nos travaux ont abouti à la caractérisation génétique du syndrome de Joubert dans la population canadienne-française grâce au séquençage exomique, et révèlent la présence d'un effet fondateur complexe avec une l'hétérogénéité allélique et intralocus importante. Ces découvertes ont éclairé la physiologie de cette maladie. Finalement, l’identification des gènes responsables ouvre de nouvelles perspectives diagnostiques ante-natales, et de conseils génétique, très précieuses pour les familles.