Single nucleotide polymorphism (SNP) selection for genotype-phenotype association studies

With hundreds of single nucleotide polymorphisms (SNPs) in a candidate gene and millions of SNPs across the genome, selecting an informative subset of SNPs to maximize the ability to detect genotype-phenotype association is of great interest and importance. In addition, with a large number of SNPs, analytic methods are needed that allow investigators to control the false positive rate resulting from large numbers of SNP genotype-phenotype analyses. This dissertation uses simulated data to explore methods for selecting SNPs for genotype-phenotype association studies. I examined the pattern of linkage disequilibrium (LD) across a candidate gene region and used this pattern to aid in localizing a disease-influencing mutation. The results indicate that the r2 measure of linkage disequilibrium is preferred over the common D′ measure for use in genotype-phenotype association studies. Using step-wise linear regression, the best predictor of the quantitative trait was not usually the single functional mutation. Rather it was a SNP that was in high linkage disequilibrium with the functional mutation. Next, I compared three strategies for selecting SNPs for application to phenotype association studies: based on measures of linkage disequilibrium, based on a measure of haplotype diversity, and random selection. The results demonstrate that SNPs selected based on maximum haplotype diversity are more informative and yield higher power than randomly selected SNPs or SNPs selected based on low pair-wise LD. The data also indicate that for genes with small contribution to the phenotype, it is more prudent for investigators to increase their sample size than to continuously increase the number of SNPs in order to improve statistical power. When typing large numbers of SNPs, researchers are faced with the challenge of utilizing an appropriate statistical method that controls the type I error rate while maintaining adequate power. We show that an empirical genotype based multi-locus global test that uses permutation testing to investigate the null distribution of the maximum test statistic maintains a desired overall type I error rate while not overly sacrificing statistical power. The results also show that when the penetrance model is simple the multi-locus global test does as well or better than the haplotype analysis. However, for more complex models, haplotype analyses offer advantages. The results of this dissertation will be of utility to human geneticists designing large-scale multi-locus genotype-phenotype association studies. ^

Recognizing a Single-Issue Spatial Election

A single-issue spatial election is a voter preference profile derived from an arrangement of candidates and voters on a line, with each voter preferring the nearer of each pair of candidates. We provide a polynomial-time algorithm that determines whether a given preference profile is a single-issue spatial election and, if so, constructs such an election. This result also has preference representation and mechanism design applications.

Multiscale modeling of the plastic behaviour in single crystal tungsten: from atomistic to crystal plasticity simulations

En una planta de fusión, los materiales en contacto con el plasma así como los materiales de primera pared experimentan condiciones particularmente hostiles al estar expuestos a altos flujos de partículas, neutrones y grandes cargas térmicas. Como consecuencia de estas diferentes y complejas condiciones de trabajo, el estudio, desarrollo y diseño de estos materiales es uno de los más importantes retos que ha surgido en los últimos años para la comunidad científica en el campo de los materiales y la energía. Debido a su baja tasa de erosión, alta resistencia al sputtering, alta conductividad térmica, muy alto punto de fusión y baja retención de tritio, el tungsteno (wolframio) es un importante candidato como material de primera pared y como posible material estructural avanzado en fusión por confinamiento magnético e inercial. Sin embargo, el tiempo de vida del tungsteno viene controlado por diversos factores como son su respuesta termo-mecánica en la superficie, la posibilidad de fusión y el fallo por acumulación de helio. Es por ello que el tiempo de vida limitado por la respuesta mecánica del tungsteno (W), y en particular su fragilidad, sean dos importantes aspectos que tienes que ser investigados. El comportamiento plástico en materiales refractarios con estructura cristalina cúbica centrada en las caras (bcc) como el tungsteno está gobernado por las dislocaciones de tipo tornillo a escala atómica y por conjuntos e interacciones de dislocaciones a escalas más grandes. El modelado de este complejo comportamiento requiere la aplicación de métodos capaces de resolver de forma rigurosa cada una de las escalas. El trabajo que se presenta en esta tesis propone un modelado multiescala que es capaz de dar respuestas ingenieriles a las solicitudes técnicas del tungsteno, y que a su vez está apoyado por la rigurosa física subyacente a extensas simulaciones atomísticas. En primer lugar, las propiedades estáticas y dinámicas de las dislocaciones de tipo tornillo en cinco potenciales interatómicos de tungsteno son comparadas, determinando cuáles de ellos garantizan una mayor fidelidad física y eficiencia computacional. Las grandes tasas de deformación asociadas a las técnicas de dinámica molecular hacen que las funciones de movilidad de las dislocaciones obtenidas no puedan ser utilizadas en los siguientes pasos del modelado multiescala. En este trabajo, proponemos dos métodos alternativos para obtener las funciones de movilidad de las dislocaciones: un modelo Monte Cario cinético y expresiones analíticas. El conjunto de parámetros necesarios para formular el modelo de Monte Cario cinético y la ley de movilidad analítica son calculados atomísticamente. Estos parámetros incluyen, pero no se limitan a: la determinación de las entalpias y energías de formación de las parejas de escalones que forman las dislocaciones, la parametrización de los efectos de no Schmid característicos en materiales bcc,etc. Conociendo la ley de movilidad de las dislocaciones en función del esfuerzo aplicado y la temperatura, se introduce esta relación como ecuación de flujo dentro de un modelo de plasticidad cristalina. La predicción del modelo sobre la dependencia del límite de fluencia con la temperatura es validada experimentalmente con ensayos uniaxiales en tungsteno monocristalino. A continuación, se calcula el límite de fluencia al aplicar ensayos uniaxiales de tensión para un conjunto de orientaciones cristalográticas dentro del triángulo estándar variando la tasa de deformación y la temperatura de los ensayos. Finalmente, y con el objetivo de ser capaces de predecir una respuesta más dúctil del tungsteno para una variedad de estados de carga, se realizan ensayos biaxiales de tensión sobre algunas de las orientaciones cristalográficas ya estudiadas en función de la temperatura.-------------------------------------------------------------------------ABSTRACT ----------------------------------------------------------Tungsten and tungsten alloys are being considered as leading candidates for structural and functional materials in future fusion energy devices. The most attractive properties of tungsten for the design of magnetic and inertial fusion energy reactors are its high melting point, high thermal conductivity, low sputtering yield and low longterm disposal radioactive footprint. However, tungsten also presents a very low fracture toughness, mostly associated with inter-granular failure and bulk plasticity, that limits its applications. As a result of these various and complex conditions of work, the study, development and design of these materials is one of the most important challenges that have emerged in recent years to the scientific community in the field of materials for energy applications. The plastic behavior of body-centered cubic (bcc) refractory metals like tungsten is governed by the kink-pair mediated thermally activated motion of h¿ (\1 11)i screw dislocations on the atomistic scale and by ensembles and interactions of dislocations at larger scales. Modeling this complex behavior requires the application of methods capable of resolving rigorously each relevant scale. The work presented in this thesis proposes a multiscale model approach that gives engineering-level responses to the technical specifications required for the use of tungsten in fusion energy reactors, and it is also supported by the rigorous underlying physics of extensive atomistic simulations. First, the static and dynamic properties of screw dislocations in five interatomic potentials for tungsten are compared, determining which of these ensure greater physical fidelity and computational efficiency. The large strain rates associated with molecular dynamics techniques make the dislocation mobility functions obtained not suitable to be used in the next steps of the multiscale model. Therefore, it is necessary to employ mobility laws obtained from a different method. In this work, we suggest two alternative methods to get the dislocation mobility functions: a kinetic Monte Carlo model and analytical expressions. The set of parameters needed to formulate the kinetic Monte Carlo model and the analytical mobility law are calculated atomistically. These parameters include, but are not limited to: enthalpy and energy barriers of kink-pairs as a function of the stress, width of the kink-pairs, non-Schmid effects ( both twinning-antitwinning asymmetry and non-glide stresses), etc. The function relating dislocation velocity with applied stress and temperature is used as the main source of constitutive information into a dislocation-based crystal plasticity framework. We validate the dependence of the yield strength with the temperature predicted by the model against existing experimental data of tensile tests in singlecrystal tungsten, with excellent agreement between the simulations and the measured data. We then extend the model to a number of crystallographic orientations uniformly distributed in the standard triangle and study the effects of temperature and strain rate. Finally, we perform biaxial tensile tests and provide the yield surface as a function of the temperature for some of the crystallographic orientations explored in the uniaxial tensile tests.

A pair of adjacent glucocorticoid response elements regulate expression of two mouse metallothionein genes

Synthesis of mouse metallothionein (MT)-I and MT-II is transcriptionally induced by the synthetic glucocorticoid, dexamethasone (DEX) or both in vivo as well as in numerous cell lines. However, the location(s) of a glucocorticoid response element (GRE) has not been described. The observation that a marked MT-I gene, as well as heterologous genes, when placed in the context of 17 kb of flanking sequence from the MT locus, are inducible by DEX and lipopolysaccharide in transgenic mice renewed the search for the GRE. Analysis of a series of deletion constructs from this 17-kb region in cultured cells identified a single 455-bp region that conferred DEX induction on a reporter gene. This 455-bp region contains two GREs that bind to the glucocorticoid receptor as assessed by gel mobility shift. Deletion of this fragment from the 17-kb flanking region eliminates the DEX responsiveness of reporter genes. The two GREs, which are located ≈1 kb upstream of the MT-II gene and ≈7 kb upstream of the MT-I gene, are necessary for induction of both genes and can function independently of elements within the proximal promoter region of either gene.

A novel pair of immunoglobulin-like receptors expressed by B cells and myeloid cells

An Fcα receptor probe of human origin was used to identify novel members of the Ig gene superfamily in mice. Paired Ig-like receptors, named PIR-A and PIR-B, are predicted from sequence analysis of the cDNAs isolated from a mouse splenic library. Both type I transmembrane proteins possess similar ectodomains with six Ig-like loops, but have different transmembrane and cytoplasmic regions. The predicted PIR-A protein has a short cytoplasmic tail and a charged Arg residue in the transmembrane region that, by analogy with the FcαR relative, suggests the potential for association with an additional transmembrane protein to form a signal transducing unit. In contrast, the PIR-B protein has an uncharged transmembrane region and a long cytoplasmic tail containing four potential immunoreceptor tyrosine-based inhibitory motifs. These features are shared by the related killer inhibitory receptors. PIR-A proteins appear to be highly variable, in that predicted peptide sequences differ for seven randomly selected PIR-A clones, whereas PIR-B cDNA clones are invariant. Southern blot analysis with PIR-B and PIR-A-specific probes suggests only one PIR-B gene and multiple PIR-A genes. The PIR-A and PIR-B genes are expressed in B lymphocytes and myeloid lineage cells, wherein both are expressed simultaneously. The characteristics of the highly-conserved PIR-A and PIR-B genes and their coordinate cellular expression suggest a potential regulatory role in humoral, inflammatory, and allergic responses.

Single-molecule protein folding: Diffusion fluorescence resonance energy transfer studies of the denaturation of chymotrypsin inhibitor 2

We report single-molecule folding studies of a small, single-domain protein, chymotrypsin inhibitor 2 (CI2). CI2 is an excellent model system for protein folding studies and has been extensively studied, both experimentally (at the ensemble level) and theoretically. Conformationally assisted ligation methodology was used to synthesize the proteins and site-specifically label them with donor and acceptor dyes. Folded and denatured subpopulations were observed by fluorescence resonance energy transfer (FRET) measurements on freely diffusing single protein molecules. Properties of these subpopulations were directly monitored as a function of guanidinium chloride concentration. It is shown that new information about different aspects of the protein folding reaction can be extracted from such subpopulation properties. Shifts in the mean transfer efficiencies are discussed, FRET efficiency distributions are translated into potentials, and denaturation curves are directly plotted from the areas of the FRET peaks. Changes in stability caused by mutation also are measured by comparing pseudo wild-type CI2 with a destabilized mutant (K17G). Current limitations and future possibilities and prospects for single-pair FRET protein folding investigations are discussed.

Crystal structure of the Holliday junction DNA in complex with a single RuvA tetramer

In the major pathway of homologous DNA recombination in prokaryotic cells, the Holliday junction intermediate is processed through its association with RuvA, RuvB, and RuvC proteins. Specific binding of the RuvA tetramer to the Holliday junction is required for the RuvB motor protein to be loaded onto the junction DNA, and the RuvAB complex drives the ATP-dependent branch migration. We solved the crystal structure of the Holliday junction bound to a single Escherichia coli RuvA tetramer at 3.1-Å resolution. In this complex, one side of DNA is accessible for cleavage by RuvC resolvase at the junction center. The refined junction DNA structure revealed an open concave architecture with a four-fold symmetry. Each arm, with B-form DNA, in the Holliday junction is predominantly recognized in the minor groove through hydrogen bonds with two repeated helix-hairpin-helix motifs of each RuvA subunit. The local conformation near the crossover point, where two base pairs are disrupted, suggests a possible scheme for successive base pair rearrangements, which may account for smooth Holliday junction movement without segmental unwinding.

A single G-to-C change causes human centromere TGGAA repeats to fold back into hairpins.

Recently, we established that satellite III (TGGAA)n tandem repeats, which occur at the centromeres of human chromosomes, pair with themselves to form an unusual "self-complementary" antiparallel duplex containing (GGA)2 motifs in which two unpaired guanines from opposite strands intercalate between sheared G.A base pairs. In separate studies, we have also established that the GCA triplet does not form bimolecular (GCA)2 motifs but instead promotes the formation of hairpins containing a GCA-turn motif in which the loop contains a single cytidine closed by a sheared G.A pair. Since TGCAA is the most frequent variant of TGGAA found in satellite III repeats, we reasoned that the potential of this variant to form GCA-turn miniloop fold-back structures might be an important factor in modulating the local structure in natural (TGGAA)n repeats. We report here the NMR-derived solution structure of the heptadecadeoxynucleotide (G)TGGAATGCAATGGAA(C) in which a central TGCAA pentamer is flanked by two TGGAA pentamers. This 17-mer forms a rather unusual and very stable hairpin structure containing eight base pairs in the stem, only four of which are Watson-Crick pairs, and a loop consisting of a single cytidine residue. The stem contains a (GGA)2 motif with intercalative 14G/4G stacking between two sheared G.A base pairs; the loop end of the stem consists of a sheared 8G.10A closing pair with the cytosine base of the 9C loop stacked on 8G. The remarkable stability of this unusual hairpin structure (Tm = 63 degrees C) suggests that it probably plays an important role in modulating the folding of satellite III (TGGAA)n repeats at the centromere.

Probing the interaction between two single molecules: fluorescence resonance energy transfer between a single donor and a single acceptor.

We extend the sensitivity of fluorescence resonance energy transfer (FRET) to the single molecule level by measuring energy transfer between a single donor fluorophore and a single acceptor fluorophore. Near-field scanning optical microscopy (NSOM) is used to obtain simultaneous dual color images and emission spectra from donor and acceptor fluorophores linked by a short DNA molecule. Photodestruction dynamics of the donor or acceptor are used to determine the presence and efficiency of energy transfer. The classical equations used to measure energy transfer on ensembles of fluorophores are modified for single-molecule measurements. In contrast to ensemble measurements, dynamic events on a molecular scale are observable in single pair FRET measurements because they are not canceled out by random averaging. Monitoring conformational changes, such as rotations and distance changes on a nanometer scale, within single biological macromolecules, may be possible with single pair FRET.

Activity of primate inferotemporal neurons related to a sought target in pair-association task.

Visual long-term memory in primates has been assessed by using the pair-association (PA) task, in which a subject retrieves and chooses the paired associate of a cue picture. Our previous studies on single neurons in the anterior inferotemporal (AIT) cortex suggested their roles in representing paired associates in the mind. To test the possibility that the delay activity of AIT neurons is related to a particular picture as a sought target, we devised the PA with color switch (PACS) task. In the PACS task, the necessity for memory retrieval and its initiation time were controlled by a color switch in the middle of the delay period. A control task, in which there is no color switch, corresponds to the conventional delayed matching-to-sample (DMS) task where the monkey chooses the same picture as a cue. We found that AIT neurons started to respond just after the color switch in the PACS task, when the cue-optimal picture's associate was presented as a cue. In contrast, they showed no response change in the DMS task. We confirmed that this effect is not due to the visual response to colors. Furthermore, when the cue-optimal picture was presented as a cue, these neurons showed suppression after the color switch in the PACS task. These results suggest that the activity of AIT neurons mediates gating mechanisms that preferentially pass information about a sought target, even when the sought target is retrieved from long-term memory.

Repair of x-ray-induced DNA double-strand breaks in specific Not I restriction fragments in human fibroblasts: joining of correct and incorrect ends.

An assay that allows measurement of absolute induction frequencies for DNA double-strand breaks (dsbs) in defined regions of the genome and that quantitates rejoining of correct DNA ends has been used to study repair of dsbs in normal human fibroblasts after x-irradiation. The approach involves hybridization of single-copy DNA probes to Not I restriction fragments separated according to size by pulsed-field gel electrophoresis. Induction of dsbs is quantitated from the decrease in the intensity of the hybridizing restriction fragment and an accumulation of a smear below the band. Rejoining of dsbs results in reconstitution of the intact restriction fragment only if correct DNA ends are joined. By comparing results from this technique with results from a conventional electrophoresis assay that detects all rejoining events, it is possible to quantitate the misrejoining frequency. Three Not I fragments on the long arm of chromosome 21 were investigated with regard to dsb induction, yielding an identical induction rate of 5.8 X 10(-3) break per megabase pair per Gy. Correct dsb rejoining was measured for two of these Not I fragments after initial doses of 80 and 160 Gy. The misrejoining frequency was about 25% for both fragments and was independent of dose. This result appears to be representative for the whole genome as shown by analysis of the entire Not I fragment distribution. The correct rejoining events primarily occurred within the first 2 h, while the misrejoining kinetics included a much slower component, with about half of the events occurring between 2 and 24 h. These misrejoining kinetics are similar to those previously reported for production of exchange aberrations in interphase chromosomes.

Single-exciton spectroscopy of single Mn doped InAs quantum dots

The optical spectroscopy of a single InAs quantum dot doped with a single Mn atom is studied using a model Hamiltonian that includes the exchange interactions between the spins of the quantum dot electron-hole pair, the Mn atom, and the acceptor hole. Our model permits linking the photoluminescence spectra to the Mn spin states after photon emission. We focus on the relation between the charge state of the Mn, A0 or A−, and the different spectra which result through either band-to-band or band-to-acceptor transitions. We consider both neutral and negatively charged dots. Our model is able to account for recent experimental results on single Mn doped InAs photoluminescence spectra and can be used to account for future experiments in GaAs quantum dots. Similarities and differences with the case of single Mn doped CdTe quantum dots are discussed.

Advantage of single-trial models for response to selection in wheat breeding multi-environment trials

An investigation was conducted to evaluate the impact of experimental designs and spatial analyses (single-trial models) of the response to selection for grain yield in the northern grains region of Australia (Queensland and northern New South Wales). Two sets of multi-environment experiments were considered. One set, based on 33 trials conducted from 1994 to 1996, was used to represent the testing system of the wheat breeding program and is referred to as the multi-environment trial (MET). The second set, based on 47 trials conducted from 1986 to 1993, sampled a more diverse set of years and management regimes and was used to represent the target population of environments (TPE). There were 18 genotypes in common between the MET and TPE sets of trials. From indirect selection theory, the phenotypic correlation coefficient between the MET and TPE single-trial adjusted genotype means [r(p(MT))] was used to determine the effect of the single-trial model on the expected indirect response to selection for grain yield in the TPE based on selection in the MET. Five single-trial models were considered: randomised complete block (RCB), incomplete block (IB), spatial analysis (SS), spatial analysis with a measurement error (SSM) and a combination of spatial analysis and experimental design information to identify the preferred (PF) model. Bootstrap-resampling methodology was used to construct multiple MET data sets, ranging in size from 2 to 20 environments per MET sample. The size and environmental composition of the MET and the single-trial model influenced the r(p(MT)). On average, the PF model resulted in a higher r(p(MT)) than the IB, SS and SSM models, which were in turn superior to the RCB model for MET sizes based on fewer than ten environments. For METs based on ten or more environments, the r(p(MT)) was similar for all single-trial models.

Pair breaking and bound states in disordered superconductors

We study the effects of inhomogeneous pairing interactions and impurities in short-coherence-length superconductors. Within the Born approximation, the effects of pairing disorder and magnetic impurities are identical. The T-matrices for pairing disorder sites with and without an impurity give rise to bound states within the BCS (Bardeen-Cooper-Schrieffer) gap, consistent with scanning tunnelling microscopy results on Bi2Sr2CaCu2O8+δ with Zn or Ni impurities.

Search for Pair-Produced Supersymmetric Top Quark Partners with the ATLAS Experiment

Searches for the supersymmetric partner of the top quark (stop) are motivated by natural supersymmetry, where the stop has to be light to cancel the large radiative corrections to the Higgs boson mass. This thesis presents three different searches for the stop at √s = 8 TeV and √s = 13 TeV using data from the ATLAS experiment at CERN’s Large Hadron Collider. The thesis also includes a study of the primary vertex reconstruction performance in data and simulation at √s = 7 TeV using tt and Z events. All stop searches presented are carried out in final states with a single lepton, four or more jets and large missing transverse energy. A search for direct stop pair production is conducted with 20.3 fb−1 of data at a center-of-mass energy of √s = 8 TeV. Several stop decay scenarios are considered, including those to a top quark and the lightest neutralino and to a bottom quark and the lightest chargino. The sensitivity of the analysis is also studied in the context of various phenomenological MSSM models in which more complex decay scenarios can be present. Two different analyses are carried out at √s = 13 TeV. The first one is a search for both gluino-mediated and direct stop pair production with 3.2 fb−1 of data while the second one is a search for direct stop pair production with 13.2 fb−1 of data in the decay scenario to a bottom quark and the lightest chargino. The results of the analyses show no significant excess over the Standard Model predictions in the observed data. Consequently, exclusion limits are set at 95% CL on the masses of the stop and the lightest neutralino.