INTERACTING WITH LOCAL AND REMOTE DATA RESPOSITORIES USING THE stashR PACKAGE

The stashR package (a Set of Tools for Administering SHared Repositories) for R implements a simple key-value style database where character string keys are associated with data values. The key-value databases can be either stored locally on the user's computer or accessed remotely via the Internet. Methods specific to the stashR package allow users to share data repositories or access previously created remote data repositories. In particular, methods are available for the S4 classes localDB and remoteDB to insert, retrieve, or delete data from the database as well as to synchronize local copies of the data to the remote version of the database. Users efficiently access information from a remote database by retrieving only the data files indexed by user-specified keys and caching this data in a local copy of the remote database. The local and remote counterparts of the stashR package offer the potential to enhance reproducible research by allowing users of Sweave to cache their R computations for a research paper in a localDB database. This database can then be stored on the Internet as a remoteDB database. When readers of the research paper wish to reproduce the computations involved in creating a specific figure or calculating a specific numeric value, they can access the remoteDB database and obtain the R objects involved in the computation.

Analysis of gene order data supports vertical inheritance of the leukotoxin operon and genome rearrangements in the 5' flanking region in genus Mannheimia

BACKGROUND: The Mannheimia subclades belong to the same bacterial genus, but have taken divergent paths toward their distinct lifestyles. For example, M. haemolytica + M. glucosida are potential pathogens of the respiratory tract in the mammalian suborder Ruminantia, whereas M. ruminalis, the supposed sister group, lives as a commensal in the ovine rumen. We have tested the hypothesis that vertical inheritance of the leukotoxin (lktCABD) operon has occurred from the last common ancestor of genus Mannheimia to any ancestor of the diverging subclades by exploring gene order data. RESULTS: We examined the gene order in the 5' flanking region of the leukotoxin operon and found that the 5' flanking gene strings, hslVU-lapB-artJ-lktC and xylAB-lktC, are peculiar to M. haemolytica + M. glucosida and M. granulomatis, respectively, whereas the gene string hslVU-lapB-lktC is present in M. ruminalis, the supposed sister group of M. haemolytica + M. glucosida, and in the most ancient subclade M. varigena. In M. granulomatis, we found remnants of the gene string hslVU-lapB-lktC in the xylB-lktC intergenic region. CONCLUSION: These observations indicate that the gene string hslVU-lapB-lktC is more ancient than the hslVU-lapB-artJ-lktC and xylAB-lktC gene strings. The presence of (remnants of) the ancient gene string hslVU-lapB-lktC among any subclades within genus Mannheimia supports that it has been vertically inherited from the last common ancestor of genus Mannheimia to any ancestor of the diverging subclades, thus reaffirming the hypothesis of vertical inheritance of the leukotoxin operon. The presence of individual 5' flanking regions in M. haemolytica + M. glucosida and M. granulomatis reflects later genome rearrangements within each subclade. The evolution of the novel 5' flanking region in M. haemolytica + M. glucosida resulted in transcriptional coupling between the divergently arranged artJ and lkt promoters. We propose that the chimeric promoter have led to high level expression of the leukotoxin operon which could explain the increased potential of certain M. haemolytica + M. glucosida strains to cause a particular type of infection.

Electroelasticity of charged black branes

We present the first-order corrected dynamics of fluid branes carrying higher-form charge by obtaining the general form of their equations of motion to pole-dipole order in the absence of external forces. Assuming linear response theory, we characterize the corresponding effective theory of stationary bent charged (an)isotropic fluid branes in terms of two sets of response coefficients, the Young modulus and the piezoelectric moduli. We subsequently find large classes of examples in gravity of this effective theory, by constructing stationary strained charged black brane solutions to first order in a derivative expansion. Using solution generating techniques and bent neutral black branes as a seed solution, we obtain a class of charged black brane geometries carrying smeared Maxwell charge in Einstein-Maxwell-dilaton gravity. In the specific case of ten-dimensional space-time we furthermore use T-duality to generate bent black branes with higher-form charge, including smeared D-branes of type II string theory. By subsequently measuring the bending moment and the electric dipole moment which these geometries acquire due to the strain, we uncover that their form is captured by classical electroelasticity theory. In particular, we find that the Young modulus and the piezoelectric moduli of our strained charged black brane solutions are parameterized by a total of 4 response coefficients, both for the isotropic as well as anisotropic cases.

Heterotic T-folds with a small number of neutral moduli

We discuss non-geometric supersymmetric heterotic string models in D=4, in the framework of the free fermionic construction. We perform a systematic scan of models with four a priori left-right asymmetric Z2 projections and shifts. We analyze some 220 models, identifying 18 inequivalent classes and addressing variants generated by discrete torsions. They do not contain geometrical or trivial neutral moduli, apart from the dilaton. However, we show the existence of flat directions in the form of exactly marginal deformations and identify patterns of symmetry breaking where product gauge groups, realized at level one, are broken to their diagonal at higher level. We also describe an “inverse Gepner map” from Heterotic to Type II models that could be used, in certain non geometric settings, to define “effective” topological invariants.

Dirichlet-Neumann inverse spectral problem for a star graph of Stieltjes strings

We solve two inverse spectral problems for star graphs of Stieltjes strings with Dirichlet and Neumann boundary conditions, respectively, at a selected vertex called root. The root is either the central vertex or, in the more challenging problem, a pendant vertex of the star graph. At all other pendant vertices Dirichlet conditions are imposed; at the central vertex, at which a mass may be placed, continuity and Kirchhoff conditions are assumed. We derive conditions on two sets of real numbers to be the spectra of the above Dirichlet and Neumann problems. Our solution for the inverse problems is constructive: we establish algorithms to recover the mass distribution on the star graph (i.e. the point masses and lengths of subintervals between them) from these two spectra and from the lengths of the separate strings. If the root is a pendant vertex, the two spectra uniquely determine the parameters on the main string (i.e. the string incident to the root) if the length of the main string is known. The mass distribution on the other edges need not be unique; the reason for this is the non-uniqueness caused by the non-strict interlacing of the given data in the case when the root is the central vertex. Finally, we relate of our results to tree-patterned matrix inverse problems.

Altered transmission of HOX and apoptotic SNPs identify a potential common pathway for clubfoot.

Clubfoot is a common birth defect that affects 135,000 newborns each year worldwide. It is characterized by equinus deformity of one or both feet and hypoplastic calf muscles. Despite numerous study approaches, the cause(s) remains poorly understood although a multifactorial etiology is generally accepted. We considered the HOXA and HOXD gene clusters and insulin-like growth factor binding protein 3 (IGFBP3) as candidate genes because of their important roles in limb and muscle morphogenesis. Twenty SNPs from the HOXA and HOXD gene clusters and 12 SNPs in IGFBP3 were genotyped in a sample composed of non-Hispanic white and Hispanic multiplex and simplex families (discovery samples) and a second sample of non-Hispanic white simplex trios (validation sample). Four SNPs (rs6668, rs2428431, rs3801776, and rs3779456) in the HOXA cluster demonstrated altered transmission in the discovery sample, but only rs3801776, located in the HOXA basal promoter region, showed altered transmission in both the discovery and validation samples (P = 0.004 and 0.028). Interestingly, HOXA9 is expressed in muscle during development. An SNP in IGFBP3, rs13223993, also showed altered transmission (P = 0.003) in the discovery sample. Gene-gene interactions were identified between variants in HOXA, HOXD, and IGFBP3 and with previously associated SNPs in mitochondrial-mediated apoptotic genes. The most significant interactions were found between CASP3 SNPS and variants in HOXA, HOXD, and IGFBP3. These results suggest a biologic model for clubfoot in which perturbation of HOX and apoptotic genes together affect muscle and limb development, which may cause the downstream failure of limb rotation into a plantar grade position.

Characteristics of a spina bifida population including North American Caucasian and Hispanic individuals.

BACKGROUND: Meningomyelocele (MM) is a common human birth defect. MM is a disorder of neural development caused by contributions from genes and environmental factors that result in the NTD and lead to a spectrum of physical and neurocognitive phenotypes. METHODS: A multidisciplinary approach has been taken to develop a comprehensive understanding of MM through collaborative efforts from investigators specializing in genetics, development, brain imaging, and neurocognitive outcome. Patients have been recruited from five different sites: Houston and the Texas-Mexico border area; Toronto, Canada; Los Angeles, California; and Lexington, Kentucky. Genetic risk factors for MM have been assessed by genotyping and association testing using the transmission disequilibrium test. RESULTS: A total of 509 affected child/parent trios and 309 affected child/parent duos have been enrolled to date for genetic association studies. Subsets of the patients have also been enrolled for studies assessing development, brain imaging, and neurocognitive outcomes. The study recruited two major ethnic groups, with 45.9% Hispanics of Mexican descent and 36.2% North American Caucasians of European descent. The remaining patients are African-American, South and Central American, Native American, and Asian. Studies of this group of patients have already discovered distinct corpus callosum morphology and neurocognitive deficits that associate with MM. We have identified maternal MTHFR 667T allele as a risk factor for MM. In addition, we also found that several genes for glucose transport and metabolism are potential risk factors for MM. CONCLUSIONS: The enrolled patient population provides a valuable resource for elucidating the disease characteristics and mechanisms for MM development.

Gene expression in sea urchin development: Study of {\it LpS1\/} gene regulation and cloning and characterization of the {\it SpKrox}-{\it 1\/} gene

Cell differentiation are associated with activation of cell lineage-specific genes. The $LpS{\it 1}\beta$ gene of Lytechinus pictus is activated at the late cleavage stage. $LpS{\it 1}\beta$ transcripts accumulate exclusively in aboral ectoderm lineages. Previous studies demonstrated two G-string DNA-elements, proximal and distal G-strings, which bind to an ectoderm-enriched nuclear factor. In order to define the cis-elements which control positive expression of the $LpS{\it 1}\beta$ gene, the regulatory region from $-$108 to +17 bp of the $LpS{\it 1}\beta$ gene promoter was characterized. The ectoderm G-string factor binds to a G/C-rich region larger than the G-string itself and the binding of the G-string factor requires sequences immediately downstream from the G-string. These downstream sequences are essential for full promoter activity. In addition, only 108 bp of $LpS{\it 1}\beta\ 5\sp\prime$ flanking DNA drives $LpS{\it 1}\beta$ gene expression in aboral ectoderm/mesenchyme cells. Therefore, for positive control of $LpS{\it 1}\beta$ gene expression, two regions of 5$\sp\prime$ flanking DNA are required: region I from base pairs $-$762 to $-$511, and region II, which includes the G/C-rich element, from base pairs $-$108 to $-$61. A mesenchyme cell repressor element is located within region I.^ DNA-binding proteins play key roles in determination of cell differentiation. The zinc finger domain is a DNA-binding domain present in many transcription factors. Based on homologies in zinc fingers, a zinc finger-encoding gene, SpKrox-1, was cloned from S. purpuratus. The putative SpKrox-1 protein has all structural characteristics of a transcription factor: four zinc fingers for DNA binding; acidic domain for transactivation; basic domain for nuclear targeting; and leucine zipper for dimerization. SpKrox-1 RNA transcripts showed a transient expression pattern which correlates largely with early embryonic development. The spatial expression of SpKrox-1 mRNA was distributed throughout the gastrula and larva ectodermal wall. However, SpKrox-1 was not expressed in pigment cells. The SpKrox-1 gene is thus a marker of a subset of SMCs or ectoderm cells. The structural features, and the transient temporal and restricted spatial expression patterns suggest that SpKrox-1 plays a role in a specific developmental event. ^

Cost and power analysis of linkage disequilibrium methods for locating genes influencing quantitative traits in humans

Linkage disequilibrium methods can be used to find genes influencing quantitative trait variation in humans. Linkage disequilibrium methods can require smaller sample sizes than linkage equilibrium methods, such as the variance component approach to find loci with a specific effect size. The increase in power is at the expense of requiring more markers to be typed to scan the entire genome. This thesis compares different linkage disequilibrium methods to determine which factors influence the power to detect disequilibrium. The costs of disequilibrium and equilibrium tests were compared to determine whether the savings in phenotyping costs when using disequilibrium methods outweigh the additional genotyping costs.^ Nine linkage disequilibrium tests were examined by simulation. Five tests involve selecting isolated unrelated individuals while four involved the selection of parent child trios (TDT). All nine tests were found to be able to identify disequilibrium with the correct significance level in Hardy-Weinberg populations. Increasing linked genetic variance and trait allele frequency were found to increase the power to detect disequilibrium, while increasing the number of generations and distance between marker and trait loci decreased the power to detect disequilibrium. Discordant sampling was used for several of the tests. It was found that the more stringent the sampling, the greater the power to detect disequilibrium in a sample of given size. The power to detect disequilibrium was not affected by the presence of polygenic effects.^ When the trait locus had more than two trait alleles, the power of the tests maximized to less than one. For the simulation methods used here, when there were more than two-trait alleles there was a probability equal to 1-heterozygosity of the marker locus that both trait alleles were in disequilibrium with the same marker allele, resulting in the marker being uninformative for disequilibrium.^ The five tests using isolated unrelated individuals were found to have excess error rates when there was disequilibrium due to population admixture. Increased error rates also resulted from increased unlinked major gene effects, discordant trait allele frequency, and increased disequilibrium. Polygenic effects did not affect the error rates. The TDT, Transmission Disequilibrium Test, based tests were not liable to any increase in error rates.^ For all sample ascertainment costs, for recent mutations ($<$100 generations) linkage disequilibrium tests were less expensive than the variance component test to carry out. Candidate gene scans saved even more money. The use of recently admixed populations also decreased the cost of performing a linkage disequilibrium test. ^

Possible Demographic Consequences of Intraspecific Shelter Competition Among American Lobsters

Pioneering work by J. Stan Cobb described how habitat architecture and body size scaling affect shelter-related behavior of American lobsters. Subsequent research suggested that shelter availability and competition could set local carrying capacity and demographics for this species. To determine how shelter spacing affects population density, the intensity of intraspecific competition and the distribution of body size for this species, I deployed sets of 10 identically sized artificial shelters spaced at distances of 2.5, 0.5, 1.0, 1.5 and 2.0 meters on otherwise featureless substrate at 10 m depth in mid-coast Maine, U.S.A. Five sets had two parallel strings of five opposing shelters and an additional linear string set 2 to apart without opposing shelters was the most widely separated treatment. Shelters spaced I m apart and closer had higher lobster population densities, more intraspecific competition and higher proportions of empty shelters. Surprisingly, lobsters there were also significantly smaller, declining from 62.7 mm to 50.9 on the carapace (CL) for 2 to linear to 0.25 m spaced shelters, respectively. Nearly all 932 lobsters measured in this study were juvenile (< 90 mm CL) and preharvestable (< 83 mm CL) sized, so mate selection and fishing effects were unlikely. At the scale of the experiment, larger lobsters leave or avoid areas of high lobster population density and intense competition for areas of low population density and relaxed competition (called "demographic diffusion"). Scuba surveys in coastal zones found lobster population densities scale with shelter densities and were highest in boulder habitat where, like the experiment, more than half the shelters were vacant. Fisheries independent scuba and trawl surveys in Maine's shallow coastal zone repeatedly recorded declines of preharvestable, lobsters larger than 60 turn CL in size and increases of those sizes offshore and in deep water. It is possible that this demographic diffusion is driven by behaviors associated with intraspecific shelter competition.

Search for contact interactions and large extra dimensions in dilepton events from pp collisions at sqrts=7 TeV with the ATLAS detector

A search for nonresonant new phenomena, originating from either contact interactions or large extra spatial dimensions, has been carried out using events with two isolated electrons or muons. These events, produced at the LHC in proton-proton collisions at root s = 7 TeV, were recorded by the ATLAS detector. The data sample, collected throughout 2011, corresponds to an integrated luminosity of 4.9 and 5.0 fb(-1) in the e(+)e(-) and mu(+)mu(-) channels, respectively. No significant deviations from the Standard Model expectation are observed. Using a Bayesian approach, 95% confidence level lower limits ranging from 9.0 to 13.9 TeV are placed on the energy scale of llqq contact interactions in the left-left isoscalar model. Lower limits ranging from 2.4 to 3.9 TeV are also set on the string scale in large extra dimension models. After combining these limits with results from a similar search in the diphoton channel, slightly more stringent limits are obtained.

ATLAS search for new phenomena in dijet mass and angular distributions using pp collisions at sqrts=7 TeV

Mass and angular distributions of dijets produced in LHC proton-proton collisions at a centre-of-mass energy root s = 7TeV have been studied with the ATLAS detector using the full 2011 data set with an integrated luminosity of 4.8 fb(-1). Dijet masses up to similar to 4.0TeV have been probed. No resonance-like features have been observed in the dijet mass spectrum, and all angular distributions are consistent with the predictions of QCD. Exclusion limits on six hypotheses of new phenomena have been set at 95% CL in terms of mass or energy scale, as appropriate. These hypotheses include excited quarks below 2.83 TeV, colour octet scalars below 1.86TeV, heavy W bosons below 1.68 TeV, string resonances below 3.61 TeV, quantum black holes with six extra space-time dimensions for quantum gravity scales below 4.11 TeV, and quark contact interactions below a compositeness scale of 7.6 TeV in a destructive interference scenario.

Black probes of Schrödinger spacetimes

We consider black probes of Anti-de Sitter and Schrödinger spacetimes embedded in string theory and M-theory and construct perturbatively new black hole geometries. We begin by reviewing black string configurations in Anti-de Sitter dual to finite temperature Wilson loops in the deconfined phase of the gauge theory and generalise the construction to the confined phase. We then consider black strings in thermal Schrödinger, obtained via a null Melvin twist of the extremal D3-brane, and construct three distinct types of black string configurations with spacelike as well as lightlike separated boundary endpoints. One of these configurations interpolates between the Wilson loop operators, with bulk duals defined in Anti-de Sitter and another class of Wilson loop operators, with bulk duals defined in Schrödinger. The case of black membranes with boundary endpoints on the M5-brane dual to Wilson surfaces in the gauge theory is analysed in detail. Four types of black membranes, ending on the null Melvin twist of the extremal M5-brane exhibiting the Schrödinger symmetry group, are then constructed. We highlight the differences between Anti-de Sitter and Schrödinger backgrounds and make some comments on the properties of the corresponding dual gauge theories.

Crystalline Confinement

We show that exotic phases arise in generalized lattice gauge theories known as quantum link models in which classical gauge fields are replaced by quantum operators. While these quantum models with discrete variables have a finite-dimensional Hilbert space per link, the continuous gauge symmetry is still exact. An efficient cluster algorithm is used to study these exotic phases. The (2+1)-d system is confining at zero temperature with a spontaneously broken translation symmetry. A crystalline phase exhibits confinement via multi stranded strings between chargeanti-charge pairs. A phase transition between two distinct confined phases is weakly first order and has an emergent spontaneously broken approximate SO(2) global symmetry. The low-energy physics is described by a (2 + 1)-d RP(1) effective field theory, perturbed by a dangerously irrelevant SO(2) breaking operator, which prevents the interpretation of the emergent pseudo-Goldstone boson as a dual photon. This model is an ideal candidate to be implemented in quantum simulators to study phenomena that are not accessible using Monte Carlo simulations such as the real-time evolution of the confining string and the real-time dynamics of the pseudo-Goldstone boson.