Advances in random matrix theory, zeta functions, and sphere packing

Over four hundred years ago, Sir Walter Raleigh asked his mathematical assistant to find formulas for the number of cannonballs in regularly stacked piles. These investigations aroused the curiosity of the astronomer Johannes Kepler and led to a problem that has gone centuries without a solution: why is the familiar cannonball stack the most efficient arrangement possible? Here we discuss the solution that Hales found in 1998. Almost every part of the 282-page proof relies on long computer verifications. Random matrix theory was developed by physicists to describe the spectra of complex nuclei. In particular, the statistical fluctuations of the eigenvalues (“the energy levels”) follow certain universal laws based on symmetry types. We describe these and then discuss the remarkable appearance of these laws for zeros of the Riemann zeta function (which is the generating function for prime numbers and is the last special function from the last century that is not understood today.) Explaining this phenomenon is a central problem. These topics are distinct, so we present them separately with their own introductory remarks.

Generation of a high-titer retroviral vector capable of expressing high levels of the human beta-globin gene.

Retrovirus-mediated gene transfer into hematopoietic cells may provide a means of treating both inherited and acquired diseases involving hematopoietic cells. Implementation of this approach for disorders resulting from mutations affecting the beta-globin gene (e.g., beta-thalassemia and sickle cell anemia), however, has been hampered by the inability to generate recombinant viruses able to efficiently and faithfully transmit the necessary sequences for appropriate gene expression. We have addressed this problem by carefully examining the interactions between retroviral and beta-globin gene sequences which affect vector transmission, stability, and expression. First, we examined the transmission properties of a large number of different recombinant proviral genomes which vary both in the precise nature of vector, beta-globin structural gene, and locus control region (LCR) core sequences incorporated and in the placement and orientation of those sequences. Through this analysis, we identified one specific vector, termed M beta 6L, which carries both the human beta-globin gene and core elements HS2, HS3, and HS4 from the LCR and faithfully transmits recombinant proviral sequences to cells with titers greater than 10(6) per ml. Populations of murine erythroleukemia (MEL) cells transduced by this virus expressed levels of human beta-globin transcript which, on a per gene copy basis, were 78% of the levels detected in an MEL-derived cell line, Hu11, which carries human chromosome 11, the site of the beta-globin locus. Analysis of individual transduced MEL cell clones, however, indicated that, while expression was detected in every clone tested (n = 17), the levels of human beta-globin treatment varied between 4% and 146% of the levels in Hu11. This clonal variation in expression levels suggests that small beta-globin LCR sequences may not provide for as strict chromosomal position-independent expression of beta-globin as previously suspected, at least in the context of retrovirus-mediated gene transfer.

Critical properties of the four-state commutative random permutation glassy Potts model in three and four dimensions

We investigate the critical properties of the four-state commutative random permutation glassy Potts model in three and four dimensions by means of Monte Carlo simulations and a finite-size scaling analysis. By using a field programmable gate array, we have been able to thermalize a large number of samples of systems with large volume. This has allowed us to observe a spin-glass ordered phase in d=4 and to study the critical properties of the transition. In d=3, our results are consistent with the presence of a Kosterlitz-Thouless transition, but also with different scenarios: transient effects due to a value of the lower critical dimension slightly below 3 could be very important.

Análise da estrutura populacional e do desequilíbrio de ligação de um painel de acessos de sorgo: uma abordagem usando teoria da coalescência

A estrutura populacional e o desequilíbrio de ligação são dois processos fundamentais para estudos evolutivos e de mapeamento associativo. Tradicionalmente, ambos têm sido investigados por meio de métodos clássicos comumente utilizados. Tais métodos certamente forneceram grandes avanços no entendimento dos processos evolutivos das espécies. No entanto, em geral, nenhum deles utiliza uma visão genealógica de forma a considerar eventos genéticos ocorridos no passado, dificultando a compreensão dos padrões de variação observados no presente. Uma abordagem que possibilita a investigação retrospectiva com base no atual polimorfismo observado é a teoria da coalescência. Assim, o objetivo deste trabalho foi analisar, com base na teoria da coalescência, a estrutura populacional e o desequilíbrio de ligação de um painel mundial de acessos de sorgo (Sorghum bicolor). Para tanto, análises de mutação, migração com fluxo gênico e recombinação foram realizadas para cinco regiões genômicas relacionadas à altura de plantas e maturidade (Dw1, Dw2, Dw4, Ma1 e Ma3) e sete populações previamente selecionadas. Em geral, elevado fluxo gênico médio (Μ = m/μ = 41,78 − 52,07) foi observado entre as populações considerando cada região genômica e todas elas simultaneamente. Os padrões sugeriram intenso intercâmbio de acessos e história evolutiva específica para cada região genômica, mostrando a importância da análise individual dos locos. A quantidade média de migrantes por geração (Μ) não foi simétrica entre pares recíprocos de populações, de acordo com a análise individual e simultânea das regiões. Isso sugere que a forma pela qual as populações se relacionaram e continuam interagindo evolutivamente não é igual, mostrando que os métodos clássicos utilizados para investigar estrutura populacional podem ser insatisfatórios. Baixas taxas médias de recombinação (ρL = 2Ner = 0,030 − 0,246) foram observadas utilizando o modelo de recombinação constante ao longo da região. Baixas e altas taxas médias de recombinação (ρr = 2Ner = 0,060 − 3,395) foram estimadas utilizando o modelo de recombinação variável ao longo da região. Os métodos tradicional (r2) e via coalescência (E[r2 rhomap]) utilizados para a estimação do desequilíbrio de ligação mostraram resultados próximos para algumas regiões genômicas e populações. No entanto, o r2 sugeriu padrões descontínuos de desequilíbrio em várias ocasiões, dificultando o entendimento e a caracterização de possíveis blocos de associação. O método via coalescência (E[r2 rhomap]) forneceu resultados que pareceram ter sido mais consistentes, podendo ser uma estratégia eventualmente importante para um refinamento dos padrões não-aleatórios de associação. Os resultados aqui encontrados sugerem que o mapeamento genético a partir de um único pool gênico pode ser insuficiente para detectar associações causais importantes para características quantitativas em sorgo.

On the time to reach a critical number of infections in epidemic models with infective and susceptible immigrants.

In this paper we examine the time T to reach a critical number K0 of infections during an outbreak in an epidemic model with infective and susceptible immigrants. The underlying process X, which was first introduced by Ridler-Rowe (1967), is related to recurrent diseases and it appears to be analytically intractable. We present an approximating model inspired from the use of extreme values, and we derive formulae for the Laplace-Stieltjes transform of T and its moments, which are evaluated by using an iterative procedure. Numerical examples are presented to illustrate the effects of the contact and removal rates on the expected values of T and the threshold K0, when the initial time instant corresponds to an invasion time. We also study the exact reproduction number Rexact,0 and the population transmission number Rp, which are random versions of the basic reproduction number R0.

Electrochemical generation of oxygen-containing groups in an ordered microporous zeolite-templated carbon

Zeolite templated carbon (ZTC) was electrochemically oxidized under various conditions, and its chemistry and structural evolution were compared to those produced by conventional chemical oxidation. In both oxidation methods, a general loss of the original structure regularity and high surface area was observed with increasing amount of oxidation. However, the electrochemical method showed much better controllability and enabled the generation of a large number of oxygen functional groups while retaining the original structure of the ZTC. Unlike chemical treatments, highly microporous carbons with an ordered 3-D structure, high surface area (ranging between 1900 and 3500 m2/g) and a large number of oxygen groups (O = 11,000–3300 μmol/g), have been prepared by the electrochemical method. Some insights into the electrooxidation mechanism of carbon materials are proposed from the obtained polarization curves, using ZTC as a model carbon material.

Boston, Massachusetts Region, Digital Elevation Model with Bathymetry

This raster layer represents surface elevation and bathymetry data for the Boston Region, Massachusetts. It was created by merging portions of MassGIS Digital Elevation Model 1:5,000 (2005) data with NOAA Estuarine Bathymetric Digital Elevation Models (30 m.) (1998). DEM data was derived from the digital terrain models that were produced as part of the MassGIS 1:5,000 Black and White Digital Orthophoto imagery project. Cellsize is 5 meters by 5 meters. Each cell has a floating point value, in meters, which represents its elevation above or below sea level.

Unity stronger than divisions. Ukraine's internal diversity. OSW Point of View Number 40, March 2014

Ukraine is deeply divided internally, although as a result of the changes that have taken place since its independence, the country’s internal divisions now have less and less to do with territorial divides, and the split into historical ‘sub-Ukraines’ has become less pronounced, especially for the younger generation. Ukraine is not a country of two competing regional identities, one in the west, the other in the east. The western identity, in which the unity of Ukraine is a key value, coexists with the multiple and diverse local patriotisms of the different regions in the east and the south of the country. The present protest movement has consolidated the country’s sense of unity. Its opponents have also been championing the indivisibility of Ukraine, even while they demanded a thorough decentralisation of the country, which was often mistaken for separatism. Russia has been stirring up separatist tendencies in Ukraine, but with little success. Crimea is an exception here, because in most respects it has remained unaffected by the dynamics of the social processes transforming mainland Ukraine – separatist tendencies are indeed deeply rooted in the peninsula.

Ukraine’s wartime nationalism. OSW COMMENTARY NUMBER 179, 28.08.2015

From the Introduction. The Euromaidan revolution and the war in Donbas that followed it have brought about deep changes in the way Ukrainian citizens think and have led to a rapid transformation of how people view the concept of national community. It began to be increasingly common for proponents of nationalism to refer to solidarity with the state instead of referring to the nation understood as an ethnic community, as they had previously done. This is due to the fact that an entire young generation has matured which takes the Ukrainian state for granted, but also because in the Russian-Ukrainian conflict most Russian-speaking citizens of Ukraine have opted for Ukraine.

The shadow of Chornobyl. Ukraine thirty years after the disaster. OSW COMMENTARY NUMBER 210/22.04.2016

Thirty years after the Chornobyl nuclear power plant disaster, its aftermath and consequences are still a permanent element of the economic, environmental and social situation of Ukraine, Belarus and some regions of Russia. Ukraine, to which the scope of this text is limited, experienced the most severe shock because, among other factors, the plant where the accident took place was located just 100 km away from Kyiv. Its consequences have affected the course of political developments in the country, and have become part of the newly-shaped national identity of independent Ukraine. The country bore the huge cost of the clean-up effort but did not give up on nuclear energy, and today nuclear power plants generate more than half of its electricity. The system of social benefits for people recognised as disaster survivors, which was put in place by the Soviet government, has become a huge burden on the country’s budget; if implemented fully, it would account for more than 10% of total public spending, and is therefore being implemented to only a partial extent. This system has reinforced the Ukrainian people’s sense of helplessness and dependence on the state. The disaster has also become part of the ‘victim nation’ blueprint of the Ukrainian national myth, which it has further solidified. The technological and environmental consequences of the disaster, and hence also its economic costs, will persist for centuries, while the social consequences will dissipate as the affected generation passes away. In any case, Chornobyl will remain an important part of the life of the Ukrainian state and society.

Next-generation sequencing of iron-metabolism related genes in Portuguese patients with iron overload: novel pathogenic genetic variants

Objective: In Southern European countries up to one-third of the patients with hereditary hemochromatosis (HH) do not present the common HFE risk genotype. In order to investigate the molecular basis of these cases we have designed a gene panel for rapid and simultaneous analysis of 6 HH-related genes (HFE, TFR2, HJV, HAMP, SLC40A1 and FTL) by next-generation sequencing (NGS). Materials and Methods: Eighty-eight iron overload Portuguese patients, negative for the common HFE mutations, were analysed. A TruSeq Custom Amplicon kit (TSCA, by Illumina) was designed in order to generate 97 amplicons covering exons, intron/exon junctions and UTRs of the mentioned genes with a cumulative target sequence of 12115bp. Amplicons were sequenced in the MiSeq instrument (IIlumina) using 250bp paired-end reads. Sequences were aligned against human genome reference hg19 using alignment and variant caller algorithms in the MiSeq reporter software. Novel variants were validated by Sanger sequencing and their pathogenic significance were assessed by in silico studies. Results: We found a total of 55 different genetic variants. These include novel pathogenic missense and splicing variants (in HFE and TFR2), a very rare variant in IRE of FTL, a variant that originates a novel translation initiation codon in the HAMP gene, among others. Conclusion: The merging of TSCA methodology and NGS technology appears to be an appropriate tool for simultaneous and fast analysis of HH-related genes in a large number of samples. However, establishing the clinical relevance of NGS-detected variants for HH development remains a hard-working task, requiring further functional studies.

Stickleback F2 generation transgenerational plasticity experiment: egg and body size data

Nongenetic inheritance mechanisms such as transgenerational plasticity (TGP) can buffer populations against rapid environmental change such as ocean warming. Yet, little is known about how long these effects persist and whether they are cumulative over generations. Here, we tested for adaptive TGP in response to simulated ocean warming across parental and grandparental generations of marine sticklebacks. Grandparents were acclimated for two months during reproductive conditioning, whereas parents experienced developmental acclimation, allowing us to compare the fitness consequences of short-term vs. prolonged exposure to elevated temperature across multiple generations. We found that reproductive output of F1 adults was primarily determined by maternal developmental temperature, but carry-over effects from grandparental acclimation environments resulted in cumulative negative effects of elevated temperature on hatching success. In very early stages of growth, F2 offspring reached larger sizes in their respective paternal and grandparental environment down the paternal line, suggesting that other factors than just the paternal genome may be transferred between generations. In later growth stages, maternal and maternal granddam environments strongly influenced offspring body size, but in opposing directions, indicating that the mechanism(s) underlying the transfer of environmental information may have differed between acute and developmental acclimation experienced by the two generations. Taken together, our results suggest that the fitness consequences of parental and grandparental TGP are highly context dependent, but will play an important role in mediating some of the impacts of rapid climate change in this system.

Calculation of GFSR pseudorandom number binary starting matrix,

Mode of access: Internet.

Method for generation of homogeneous multicellular tumor spheroids applicable to a wide variety of cell types

Multicellular tumor spheroids (MCTS) are used as organotypic models of normal and solid tumor tissue. Traditional techniques for generating MCTS, such as growth on nonadherent surfaces, in suspension, or on scaffolds, have a number of drawbacks, including the need for manual selection to achieve a homogeneous population and the use of nonphysiological matrix compounds. In this study we describe a mild method for the generation of MCTS, in which individual spheroids form in hanging drops suspended from a microtiter plate. The method has been successfully applied to a broad range of cell lines and shows nearly 100% efficiency (i.e., one spheroid per drop). Using the hepatoma cell line, HepG2, the hanging drop method generated well-rounded MCTS with a narrow size distribution (coefficient of variation [CV] 10% to 15%, compared with 40% to 60% for growth on nonadherent surfaces). Structural analysis of HepG2 and a mammary gland adenocarcinoma cell line, MCF-7, composed spheroids, revealed highly organized, three-dimensional, tissue-like structures with an extensive extracellular matrix. The hanging drop method represents an attractive alternative for MCTS production, because it is mild, can be applied to a wide variety of cell lines, and can produce spheroids of a homogeneous size without the need for sieving or manual selection. The method has applications for basic studies of physiology and metabolism, tumor biology, toxicology, cellular organization, and the development of bioartificial tissue. (C) 2003 Wiley Periodicals, Inc.

Binomial leap methods for simulating stochastic chemical kinetics

This paper discusses efficient simulation methods for stochastic chemical kinetics. Based on the tau-leap and midpoint tau-leap methods of Gillespie [D. T. Gillespie, J. Chem. Phys. 115, 1716 (2001)], binomial random variables are used in these leap methods rather than Poisson random variables. The motivation for this approach is to improve the efficiency of the Poisson leap methods by using larger stepsizes. Unlike Poisson random variables whose range of sample values is from zero to infinity, binomial random variables have a finite range of sample values. This probabilistic property has been used to restrict possible reaction numbers and to avoid negative molecular numbers in stochastic simulations when larger stepsize is used. In this approach a binomial random variable is defined for a single reaction channel in order to keep the reaction number of this channel below the numbers of molecules that undergo this reaction channel. A sampling technique is also designed for the total reaction number of a reactant species that undergoes two or more reaction channels. Samples for the total reaction number are not greater than the molecular number of this species. In addition, probability properties of the binomial random variables provide stepsize conditions for restricting reaction numbers in a chosen time interval. These stepsize conditions are important properties of robust leap control strategies. Numerical results indicate that the proposed binomial leap methods can be applied to a wide range of chemical reaction systems with very good accuracy and significant improvement on efficiency over existing approaches. (C) 2004 American Institute of Physics.