Breakthroughs in genetic studies of ankylosing spondylitis

Ankylosing spondylitis (AS), the prototypic seronegative arthropathy, is known to be highly heritable, with >90% of the risk of developing the disease determined genetically. As with most common heritable diseases, progress in identifying the genes involved using family-based or candidate gene approaches has been slow. The recent development of the genome-wide association study approach has revolutionized genetic studies of such diseases. Early studies in ankylosing spondylitis have produced two major breakthroughs in the identification of genes contributing roughly one third of the population attributable risk of the disease, and pointing directly to a potential therapy. These exciting findings highlight the potential of future more comprehensive genetic studies of determinants of disease risk and clinical manifestations, and are the biggest advance in our understanding of the causation of the disease since the discovery of the association with HLA-B27.

Influence of LRP5 polymorphisms on normal variation in BMD

Genetic studies based on cohorts with rare and extreme bone phenotypes have shown that the LRP5 gene is an important genetic modulator of BMD. Using family-based and case-control approaches, this study examines the role of the LRP5 gene in determining normal population variation of BMD and describes significant association and suggestive linkage between LRP5 gene polymorphisms and BMD in >900 individuals with a broad range of BMD. Introduction: Osteoporosis is a common, highly heritable condition determined by complex interactions of genetic and environmental etiologies. Genetic factors alone can account for 50-80% of the interindividual variation in BMD. Mutations in the LRP5 gene on chromosome 11q12-13 have been associated with rare syndromes characterized by extremely low or high BMD, but little is known about the contribution of this gene to the development of osteoporosis and determination of BMD in a normal population. Materials and Methods: To examine the entire spectrum of low to high BMD, 152 osteoporotic probands, their families (597 individuals), and 160 women with elevated BMD (T score > 2.5) were recruited. BMD at the lumbar spine, femoral neck, and hip were measured in each subject using DXA. Results: PAGE sequencing of the LRP5 gene revealed 10 single nucleotide polymorphisms (SNPs), 8 of which had allele frequencies of >5%, in exons 8, 9, 10, 15, and 18 and in introns 6, 7, and 21. Within families, a strong association was observed between an SNP at nucleotide C171346A in intron 21 and total hip BMD (p < 1 × 10-5 in men only, p = 0.0019 in both men and women). This association was also observed in comparisons of osteoporotic probands and unrelated elevated BMD in women (p = 0.03), along with associations with markers in exons 8 (C135242T, p = 0.007) and 9 (C141759T, p = 0.02). Haplotypes composed of two to three of the SNPs G121513A, C135242T, G138351A, and C141759T were strongly associated with BMD when comparing osteoporotic probands and high BMD cases (p < 0.003). An SNP at nucleotide C165215T in exon 18 was linked to BMD at the lumbar spine, femoral neck, and total hip (parametric LOD scores = 2.8, 2.5, and 2.2 and nonparametric LOD scores = 0.3, 1.1, and 2.2, respectively) but was not genetically associated with BMD variation. Conclusion: These results show that common LRP5 polymorphisms contribute to the determination of BMD in the general population.

The ATHEROMA (Atorvastatin Therapy: Effects on Reduction of Macrophage Activity) Study. Evaluation using ultrasmall superparamagnetic iron oxide-enhanced magnetic resonance imaging in carotid disease

Objectives: The aim of this study was to evaluate the effects of low-dose (10 mg) and high-dose (80 mg) atorvastatin on carotid plaque inflammation as determined by ultrasmall superparamagnetic iron oxide (USPIO)-enhanced carotid magnetic resonance imaging (MRI). The hypothesis was that treatment with 80 mg atorvastatin would demonstrate quantifiable changes in USPIO-enhanced MRI-defined inflammation within the first 3 months of therapy. Background: Preliminary studies indicate that USPIO-enhanced MRI can identify macrophage infiltration in human carotid atheroma in vivo and hence may be a surrogate marker of plaque inflammation. Methods: Forty-seven patients with carotid stenosis >40% on duplex ultrasonography and who demonstrated intraplaque accumulation of USPIO on MRI at baseline were randomly assigned in a balanced, double-blind manner to either 10 or 80 mg atorvastatin daily for 12 weeks. Baseline statin therapy was equivalent to 10 mg of atorvastatin or less. The primary end point was change from baseline in signal intensity (ΔSI) on USPIO-enhanced MRI in carotid plaque at 6 and 12 weeks. Results: Twenty patients completed 12 weeks of treatment in each group. A significant reduction from baseline in USPIO-defined inflammation was observed in the 80-mg group at both 6 weeks (ΔSI 0.13; p = 0.0003) and at 12 weeks (ΔSI 0.20; p < 0.0001). No difference was observed with the low-dose regimen. The 80-mg atorvastatin dose significantly reduced total cholesterol by 15% (p = 0.0003) and low-density lipoprotein cholesterol by 29% (p = 0.0001) at 12 weeks. Conclusions: Aggressive lipid-lowering therapy over a 3-month period is associated with significant reduction in USPIO-defined inflammation. USPIO-enhanced MRI methodology may be a useful imaging biomarker for the screening and assessment of therapeutic response to "anti-inflammatory" interventions in patients with atherosclerotic lesions. (Effects of Atorvastatin on Macrophage Activity and Plaque Inflammation Using Magnetic Resonance Imaging [ATHEROMA]; NCT00368589).

Effects of added shelter and stocking density on growth of sleepy cod Oxyeleotris lineolatus in ponds

Sleepy cod Oxyeleotris lineolatus is a species of freshwater goby in demand in Australian markets by consumers of Asian origin. It is related to marble goby Oxyeleotris marmoratus, the most expensive freshwater food fish in Asia, which is cultured throughout southeast Asia in ponds and cages. The performance of sleepy cod in culture conditions was investigated to assess the viability of farming them in northern Australia. Sleepy cod fingerlings (62.8 +/- 0.8 mm total length and 2.56 +/- 0.095 g) were stocked into experimental ponds at 32,857 fish/ha, and grown out for 8 mo. Shelter was provided in each of three replicate ponds and was absent in three control ponds. The provision of shelter in juvenile growout was found to be of no benefit, although fish in ponds provided with shelter weighed slightly more per unit length than fish in ponds without shelter. Cannibalism was not a problem in growout, and survival was close to 100%. After the shelter trial was completed, fish were graded into large and small classes (three replicates of each), and grown out without shelter at the same density for 158 d. Following that, fish were again graded, and the largest 30% retained from growout at a density of 8,857 fish/ha (large, 198 +/-6.44 g) or 10,000 fish/ha (small, 48.9 +/-1.27 g). These were grown out for 188 d. Growth of selected stock at low densities was slower than earlier growth rates, although smaller fish gained weight more rapidly than larger fish. Growth rates were better than the only published data for marble goby. Further investigation into high density culture and different genotypes of sleepy cod needs to be undertaken to determine the viability of pond culture.

Row spacing and planting density effects on the growth and yield of sugarcane. 1. Responses in fumigated and non-fumigated soil.

It has been reported that high-density planting of sugarcane can improve cane and sugar yield through promoting rapid canopy closure and increasing radiation interception earlier in crop growth. It is widely known that the control of adverse soil biota through fumigation (removes soil biological constraints and improves soil health) can improve cane and sugar yield. Whether the responses to high-density planting and improved soil health are additive or interactive has important implications for the sugarcane production system. Field experiments established at Bundaberg and Mackay, Queensland, Australia, involved all combinations of 2-row spacings (0.5 and 1.5 m), two planting densities (27 000 and 81 000 two-eyed setts/ha), and two soil fumigation treatments (fumigated and non-fumigated). The Bundaberg experiment had two cultivars (Q124, Q155), was fully irrigated, and harvested 15 months after planting. The Mackay experiment had one cultivar (Q117), was grown under rainfed conditions, and harvested 10 months after planting. High-density planting (81 000 setts/ha in 0.5-m rows) did not produce any more cane or sugar yield at harvest than low-density planting (27 000 setts/ha in 1.5-m rows) regardless of location, crop duration (15 v. 10 months), water supply (irrigated v. rainfed), or soil health (fumigated v. non-fumigated). Conversely, soil fumigation generally increased cane and sugar yields regardless of site, row spacing, and planting density. In the Bundaberg experiment there was a large fumigation x cultivar x density interaction (P<0.01). Cultivar Q155 responded positively to higher planting density in non-fumigated soil but not in fumigated soil, while Q124 showed a negative response to higher planting density in non-fumigated soil but no response in fumigated soil. In the Mackay experiment, Q117 showed a non-significant trend of increasing yield in response to increasing planting density in non-fumigated soil, similar to the Q155 response in non-fumigated soil at Bundaberg. The similarity in yield across the range of row spacings and planting densities within experiments was largely due to compensation between stalk number and stalk weight, particularly when fumigation was used to address soil health. Further, the different cultivars (Q124 and Q155 at Bundaberg and Q117 at Mackay) exhibited differing physiological responses to the fumigation, row spacing, and planting density treatments. These included the rate of tiller initiation and subsequent loss, changes in stalk weight, and propensity to lodging. These responses suggest that there may be potential for selecting cultivars suited to different planting configurations.

LRP5 regulates human body fat distribution by modulating adipose progenitor biology in a dose- and depot-specific fashion

Summary Common variants in WNT pathway genes have been associated with bone mass and fat distribution, the latter predicting diabetes and cardiovascular disease risk. Rare mutations in the WNT co-receptors LRP5 and LRP6 are similarly associated with bone and cardiometabolic disorders. We investigated the role of LRP5 in human adipose tissue. Subjects with gain-of-function LRP5 mutations and high bone mass had enhanced lower-body fat accumulation. Reciprocally, a low bone mineral density-associated common LRP5 allele correlated with increased abdominal adiposity. Ex vivo LRP5 expression was higher in abdominal versus gluteal adipocyte progenitors. Equivalent knockdown of LRP5 in both progenitor types dose-dependently impaired β-catenin signaling and led to distinct biological outcomes: diminished gluteal and enhanced abdominal adipogenesis. These data highlight how depot differences in WNT/β-catenin pathway activity modulate human fat distribution via effects on adipocyte progenitor biology. They also identify LRP5 as a potential pharmacologic target for the treatment of cardiometabolic disorders. © 2015 The Authors.

Multistage genome-wide association meta-analyses identified two new loci for bone mineral density

Aiming to identify novel genetic variants and to confirm previously identified genetic variants associated with bone mineral density (BMD), we conducted a three-stage genome-wide association (GWA) meta-analysis in 27 061 study subjects. Stage 1 meta-analyzed seven GWA samples and 11 140 subjects for BMDs at the lumbar spine, hip and femoral neck, followed by a Stage 2 in silico replication of 33 SNPs in 9258 subjects, and by a Stage 3 de novo validation of three SNPs in 6663 subjects. Combining evidence from all the stages, we have identified two novel loci that have not been reported previously at the genome-wide significance (GWS; 5.0 × 10-8) level: 14q24.2 (rs227425, P-value 3.98 × 10-13, SMOC1) in the combined sample of males and females and 21q22.13 (rs170183, P-value 4.15 × 10-9, CLDN14) in the female-specific sample. The two newly identified SNPs were also significant in the GEnetic Factors for OSteoporosis consortium (GEFOS, n 5 32 960) summary results. We have also independently confirmed 13 previously reported loci at the GWS level: 1p36.12 (ZBTB40), 1p31.3 (GPR177), 4p16.3 (FGFRL1), 4q22.1 (MEPE), 5q14.3 (MEF2C), 6q25.1 (C6orf97, ESR1), 7q21.3 (FLJ42280, SHFM1), 7q31.31 (FAM3C, WNT16), 8q24.12 (TNFRSF11B), 11p15.3 (SOX6), 11q13.4 (LRP5), 13q14.11 (AKAP11) and 16q24 (FOXL1). Gene expression analysis in osteogenic cells implied potential functional association of the two candidate genes (SMOC1 and CLDN14) in bone metabolism. Our findings independently confirm previously identified biological pathways underlying bone metabolism and contribute to the discovery of novel pathways, thus providing valuable insights into the intervention and treatment of osteoporosis. © The Author 2013. Published by Oxford University Press.

Chemically Programmed Ultrahigh Density Two-Dimensional Semiconductor Superlattice Array

Designing an ultrahigh density linear superlattice array consisting of periodic blocks of different semiconductors in the strong confinement regime via a direct synthetic route remains an unachieved challenge in nanotechnology. We report a general synthesis route for the formulation of a large-area ultrahigh density superlattice array that involves adjoining multiple units of ZnS rods by prolate US particles at the tips. A single one-dimensional wire is 300-500 nm long and consists of periodic quantum wells with a barrier width of 5 nm provided by ZnS and a well width of 1-2 nm provided by CdS, defining a superlattice structure. The synthesis route allows for tailoring of ultranarrow laserlike emissions (fwhm approximate to 125 meV) originating from strong interwell energy dispersion along with control of the width, pitch, and registry of the superlattice assembly. Such an exceptional high-density superlattice array could form the basis of ultrahigh density memories in addition to offering opportunities for technological advancement in conventional heterojunction-based device applications.

Genetic Mechanisms Underlying Autism Spectrum Disorders

Autism is a childhood-onset developmental disorder characterized by deficits in reciprocal social interaction, verbal and non-verbal communication, and dependence on routines and rituals. It belongs to a spectrum of disorders (autism spectrum disorders, ASDs) which share core symptoms but show considerable variation in severity. The whole spectrum affects 0.6-0.7% of children worldwide, inducing a substantial public health burden and causing suffering to the affected families. Despite having a very high heritability, ASDs have shown exceptional genetic heterogeneity, which has complicated the identification of risk variants and left the etiology largely unknown. However, recent studies suggest that rare, family-specific factors contribute significantly to the genetic basis of ASDs. In this study, we investigated the role of DISC1 (Disrupted-in-schizophrenia-1) in ASDs, and identified association with markers and haplotypes previously associated with psychiatric phenotypes. We identified four polymorphic micro-RNA target sites in the 3 UTR of DISC1, and showed that hsa-miR-559 regulates DISC1 expression in vitro in an allele-specific manner. We also analyzed an extended autism pedigree with genealogical roots in Central Finland reaching back to the 17th century. To take advantage of the beneficial characteristics of population isolates to gene mapping and reduced genetic heterogeneity observed in distantly related individuals, we performed a microsatellite-based genome-wide screen for linkage and linkage disequilibrium in this pedigree. We identified a putative autism susceptibility locus on chromosome 19p13.3 and obtained further support for previously reported loci at 1q23 and 15q11-q13. To follow-up these findings, we extended our study sample from the same sub-isolate and initiated a genome-wide analysis of homozygosity and allelic sharing using high-density SNP markers. We identified a small number of haplotypes shared by different subsets of the genealogically connected cases, along with convergent biological pathways from SNP and gene expression data, which highlighted axon guidance molecules in the pathogenesis of ASDs. In conclusion, the results obtained in this thesis show that multiple distinct genetic variants are responsible for the ASD phenotype even within single pedigrees from an isolated population. We suggest that targeted resequencing of the shared haplotypes, linkage regions, and other susceptibility loci is essential to identify the causal variants. We also report a possible micro-RNA mediated regulatory mechanism, which might partially explain the wide-range neurobiological effects of the DISC1 gene.

Morphogenesis of Mammalian Endoplasmic Reticulum and Golgi Apparatus throughout the Cell Cycle

The endoplasmic reticulum (ER) and the Golgi apparatus are organelles that produce, modify and transport proteins and lipids and regulate Ca2+ environment within cells. Structurally they are composed of sheets and tubules. Sheets may take various forms: intact, fenestrated, single or stacked. The ER, including the nuclear envelope, is a single continuous network, while the Golgi shows only some level of connectivity. It is often unclear, how different morphologies correspond to particular functions. Previous studies indicate that the structures of the ER and Golgi are dynamic and regulated by fusion and fission events, cytoskeleton, rate of protein synthesis and secretion, and specific structural proteins. For example, many structural proteins shaping tubular ER have been identified, but sheet formation is much more unclear. In this study, we used light and electron microscopy to study morphological changes of the ER and Golgi in mammalian cells. The proportion, type, location and dynamics of ER sheets and tubules were found to vary in a cell type or cell cycle stage dependent manner. During interphase, ER and Golgi structures were demonstrated to be regulated by p37, a cofactor of the fusion factor p97, and microtubules, which also affected the localization of the organelles. Like previously shown for the Golgi, the ER displayed a tendency for fenestration and tubulation during mitosis. However, this shape change did not result in ER fragmentation as happens to Golgi, but a continuous network was retained. The activity of p97/p37 was found to be important for the reassembly of both organelles after mitosis. In EM images, ER sheet membranes appear rough, since they contain attached ribosomes, whereas tubular membranes appear smooth. Our studies revealed that structural changes of the ER towards fenestrated and tubular direction correlate with loss of ER-bound ribosomes and vice versa. High and low curvature ER membranes have a low and high density of ribosomes, respectively. To conclude, both ER and Golgi architecture depend on fusion activity of p97/p37. ER morphogenesis, particularly of the sheet shape, is intimately linked to the density of membrane bound ribosomes.

Lahopuumäärän ennustaminen ja kartoitus lentokonelaserkeilauksen avulla

Lahopuun määrästä ja sijoittumisesta ollaan kiinnostuneita paitsi elinympäristöjen monimuotoisuuden, myös ilmakehän hiilen varastoinnin kannalta. Tutkimuksen tavoitteena oli kehittää aluepohjainen laserkeilausdataa hyödyntävä malli lahopuukohteiden paikantamiseksi ja lahopuun määrän estimoimiseksi. Samalla tutkittiin mallin selityskyvyn muuttumista mallinnettavan ruudun kokoa suurennettaessa. Tutkimusalue sijaitsi Itä-Suomessa Sonkajärvellä ja koostui pääasiassa nuorista hoidetuista talousmetsistä. Tutkimuksessa käytettiin harvapulssista laserkeilausdataa sekä kaistoittain mitattua maastodataa kuolleesta puuaineksesta. Aineisto jaettiin siten, että neljäsosa datasta oli käytössä mallinnusta varten ja loput varattiin valmiiden mallien testaamiseen. Lahopuun mallintamisessa käytettiin sekä parametrista että ei-parametrista mallinnusmenetelmää. Logistisen regression avulla erikokoisille (0,04, 0,20, 0,32, 0,52 ja 1,00 ha) ruuduille ennustettiin todennäköisyys lahopuun esiintymiselle. Muodostettujen mallien selittävät muuttujat valittiin 80 laserpiirteen ja näiden muunnoksien joukosta. Mallien selittävät muuttujat valittiin kolmessa vaiheessa. Aluksi muuttujia tarkasteltiin visuaalisesti kuvaamalla ne lahopuumäärän suhteen. Ensimmäisessä vaiheessa sopivimmiksi arvioitujen muuttujien selityskykyä testattiin mallinnuksen toisessa vaiheessa yhden muuttujan mallien avulla. Lopullisessa usean muuttujan mallissa selittävien muuttujien kriteerinä oli tilastollinen merkitsevyys 5 % riskitasolla. 0,20 hehtaarin ruutukoolle luotu malli parametrisoitiin muun kokoisille ruuduille. Logistisella regressiolla toteutetun parametrisen mallintamisen lisäksi, 0,04 ja 1,0 hehtaarin ruutukokojen aineistot luokiteltiin ei-parametrisen CART-mallinnuksen (Classification and Regression Trees) avulla. CARTmenetelmällä etsittiin aineistosta vaikeasti havaittavia epälineaarisia riippuvuuksia laserpiirteiden ja lahopuumäärän välillä. CART-luokittelu tehtiin sekä lahopuustoisuuden että lahopuutilavuuden suhteen. CART-luokituksella päästiin logistista regressiota parempiin tuloksiin ruutujen luokituksessa lahopuustoisuuden suhteen. Logistisella mallilla tehty luokitus parani ruutukoon suurentuessa 0,04 ha:sta(kappa 0,19) 0,32 ha:iin asti (kappa 0,38). 0,52 ha:n ruutukoolla luokituksen kappa-arvo kääntyi laskuun (kappa 0,32) ja laski edelleen hehtaarin ruutukokoon saakka (kappa 0,26). CART-luokitus parani ruutukoon kasvaessa. Luokitustulokset olivat logistista mallinnusta parempia sekä 0,04 ha:n (kappa 0,24) että 1,0 ha:n (kappa 0,52) ruutukoolla. CART-malleilla määritettyjen ruutukohtaisten lahopuutilavuuksien suhteellinen RMSE pieneni ruutukoon kasvaessa. 0,04 hehtaarin ruutukoolla koko aineiston lahopuumäärän suhteellinen RMSE oli 197,1 %, kun hehtaarin ruutukoolla vastaava luku oli 120,3 %. Tämän tutkimuksen tulosten perusteella voidaan todeta, että maastossa mitatun lahopuumäärän ja tutkimuksessa käytettyjen laserpiirteiden yhteys on pienellä ruutukoolla hyvin heikko, mutta vahvistuu hieman ruutukoon kasvaessa. Kun mallinnuksessa käytetty ruutukoko kasvaa, pienialaisten lahopuukeskittymien havaitseminen kuitenkin vaikeutuu. Tutkimuksessa kohteen lahopuustoisuus pystyttiin kartoittamaan kohtuullisesti suurella ruutukoolla, mutta pienialaisten kohteiden kartoittaminen ei onnistunut käytetyillä menetelmillä. Pienialaisten kohteiden paikantaminen laserkeilauksen avulla edellyttää jatkotutkimusta erityisesti tiheäpulssisen laserdatan käytöstä lahopuuinventoinneissa.

Modelling vacuum arcs : from plasma initiation to surface interactions

A better understanding of vacuum arcs is desirable in many of today's 'big science' projects including linear colliders, fusion devices, and satellite systems. For the Compact Linear Collider (CLIC) design, radio-frequency (RF) breakdowns occurring in accelerating cavities influence efficiency optimisation and cost reduction issues. Studying vacuum arcs both theoretically as well as experimentally under well-defined and reproducible direct-current (DC) conditions is the first step towards exploring RF breakdowns. In this thesis, we have studied Cu DC vacuum arcs with a combination of experiments, a particle-in-cell (PIC) model of the arc plasma, and molecular dynamics (MD) simulations of the subsequent surface damaging mechanism. We have also developed the 2D Arc-PIC code and the physics model incorporated in it, especially for the purpose of modelling the plasma initiation in vacuum arcs. Assuming the presence of a field emitter at the cathode initially, we have identified the conditions for plasma formation and have studied the transitions from field emission stage to a fully developed arc. The 'footing' of the plasma is the cathode spot that supplies the arc continuously with particles; the high-density core of the plasma is located above this cathode spot. Our results have shown that once an arc plasma is initiated, and as long as energy is available, the arc is self-maintaining due to the plasma sheath that ensures enhanced field emission and sputtering. The plasma model can already give an estimate on how the time-to-breakdown changes with the neutral evaporation rate, which is yet to be determined by atomistic simulations. Due to the non-linearity of the problem, we have also performed a code-to-code comparison. The reproducibility of plasma behaviour and time-to-breakdown with independent codes increased confidence in the results presented here. Our MD simulations identified high-flux, high-energy ion bombardment as a possible mechanism forming the early-stage surface damage in vacuum arcs. In this mechanism, sputtering occurs mostly in clusters, as a consequence of overlapping heat spikes. Different-sized experimental and simulated craters were found to be self-similar with a crater depth-to-width ratio of about 0.23 (sim) - 0.26 (exp). Experiments, which we carried out to investigate the energy dependence of DC breakdown properties, point at an intrinsic connection between DC and RF scaling laws and suggest the possibility of accumulative effects influencing the field enhancement factor.

Revised density of magnetized nuclear matter at the neutron drip line

We study the onset of the neutron drip in high-density matter in the presence of a magnetic field. It has been found that, for systems having only protons and electrons, in the presence of a magnetic field greater than or similar to 10(15) G, neutronization occurs at a density that is at least an order of magnitude higher compared to that in a nonmagnetic system. In a system with heavier ions, the effect of the magnetic field, however, starts arising at a much higher field, greater than or similar to 10(17) G. These results may have important implications for high-magnetic-field neutron stars and white dwarfs and, in general, in nuclear astrophysics when the system is embedded within a strong magnetic field.

A molecular dynamics study of ambient and high pressure phases of silica: Structure and enthalpy variation with molar volume

Extensive molecular dynamics studies of 13 different silica polymorphs are reported in the isothermal-isobaric ensemble with the Parrinello-Rahman variable shape simulation cell. The van Beest-Kramer-van Santen (BKS) potential is shown to predict lattice parameters for most phases within 2%-3% accuracy, as well as the relative stabilities of different polymorphs in agreement with experiment. Enthalpies of high-density polymorphs - CaCl2-type, alpha-PbO2-type, and pyrite-type for which no experimental data are available as yet, are predicted here. Further, the calculated enthalpies exhibit two distinct regimes as a function of molar volume-for low and medium-density polymorphs, it is almost independent of volume, while for high-pressure phases a steep dependence is seen. A detailed analysis indicates that the increased short-range contributions to enthalpy in the high-density phases arise not only from an increased coordination number of silicon but also shorter Si-O bond lengths. Our results indicate that amorphous phases of silica exhibit better optimization of short-range interactions than crystalline phases at the same density while the magnitude of Coulombic contributions is lower in the amorphous phase. (C) 2014 AIP Publishing LLC.