Microbial catabolic activities are naturally selected by metabolic energy harvest rate

The fundamental trade-off between yield and rate of energy harvest per unit of substrate has been largely discussed as a main characteristic for microbial established cooperation or competition. In this study, this point is addressed by developing a generalized model that simulates competition between existing and not experimentally reported microbial catabolic activities defined only based on well-known biochemical pathways. No specific microbial physiological adaptations are considered, growth yield is calculated coupled to catabolism energetics and a common maximum biomass-specific catabolism rate (expressed as electron transfer rate) is assumed for all microbial groups. Under this approach, successful microbial metabolisms are predicted in line with experimental observations under the hypothesis of maximum energy harvest rate. Two microbial ecosystems, typically found in wastewater treatment plants, are simulated, namely: (i) the anaerobic fermentation of glucose and (ii) the oxidation and reduction of nitrogen under aerobic autotrophic (nitrification) and anoxic heterotrophic and autotrophic (denitrification) conditions. The experimentally observed cross feeding in glucose fermentation, through multiple intermediate fermentation pathways, towards ultimately methane and carbon dioxide is predicted. Analogously, two-stage nitrification (by ammonium and nitrite oxidizers) is predicted as prevailing over nitrification in one stage. Conversely, denitrification is predicted in one stage (by denitrifiers) as well as anammox (anaerobic ammonium oxidation). The model results suggest that these observations are a direct consequence of the different energy yields per electron transferred at the different steps of the pathways. Overall, our results theoretically support the hypothesis that successful microbial catabolic activities are selected by an overall maximum energy harvest rate

How Chaotic is the Stadium Billiard? A Semiclassical Analysis

The impression gained from the literature published to date is that the spectrum of the stadium billiard can be adequately described, semiclassically, by the Gutzwiller periodic orbit trace formula together with a modified treatment of the marginally stable family of bouncing ball orbits. I show that this belief is erroneous. The Gutzwiller trace formula is not applicable for the phase space dynamics near the bouncing ball orbits. Unstable periodic orbits close to the marginally stable family in phase space cannot be treated as isolated stationary phase points when approximating the trace of the Green function. Semiclassical contributions to the trace show an $\hbar$ - dependent transition from hard chaos to integrable behavior for trajectories approaching the bouncing ball orbits. A whole region in phase space surrounding the marginal stable family acts, semiclassically, like a stable island with boundaries being explicitly $\hbar$-dependent. The localized bouncing ball states found in the billiard derive from this semiclassically stable island. The bouncing ball orbits themselves, however, do not contribute to individual eigenvalues in the spectrum. An EBK-like quantization of the regular bouncing ball eigenstates in the stadium can be derived. The stadium billiard is thus an ideal model for studying the influence of almost regular dynamics near marginally stable boundaries on quantum mechanics.Comment: 27 pages, 6 figures, submitted to J. Phys.

Approaching the ideal elastic limit of metallic glasses

The ideal elastic limit is the upper bound to the stress and elastic strain a material can withstand. This intrinsic property has been widely studied for crystalline metals, both theoretically and experimentally. For metallic glasses, however, the ideal elastic limit remains poorly characterized and understood. Here we show that the elastic strain limit and the corresponding strength of submicron-sized metallic glass specimens are about twice as high as the already impressive elastic limit observed in bulk metallic glass samples, in line with model predictions of the ideal elastic limit of metallic glasses. We achieve this by employing an in situ transmission electron microscope tensile deformation technique. Furthermore, we propose an alternative mechanism for the apparent 'work hardening' behaviour observed in the tensile stress–strain curves

Zircon ages in granulite facies rocks: decoupling from geochemistry above 850 °C?

Granulite facies rocks frequently show a large spread in their zircon ages, the interpretation of which raises questions: Has the isotopic system been disturbed? By what process(es) and conditions did the alteration occur? Can the dates be regarded as real ages, reflecting several growth episodes? Furthermore, under some circumstances of (ultra-)high-temperature metamorphism, decoupling of zircon U–Pb dates from their trace element geochemistry has been reported. Understanding these processes is crucial to help interpret such dates in the context of the P–T history. Our study presents evidence for decoupling in zircon from the highest grade metapelites (> 850 °C) taken along a continuous high-temperature metamorphic field gradient in the Ivrea Zone (NW Italy). These rocks represent a well-characterised segment of Permian lower continental crust with a protracted high-temperature history. Cathodoluminescence images reveal that zircons in the mid-amphibolite facies preserve mainly detrital cores with narrow overgrowths. In the upper amphibolite and granulite facies, preserved detrital cores decrease and metamorphic zircon increases in quantity. Across all samples we document a sequence of four rim generations based on textures. U–Pb dates, Th/U ratios and Ti-in-zircon concentrations show an essentially continuous evolution with increasing metamorphic grade, except in the samples from the granulite facies, which display significant scatter in age and chemistry. We associate the observed decoupling of zircon systematics in high-grade non-metamict zircon with disturbance processes related to differences in behaviour of non-formula elements (i.e. Pb, Th, U, Ti) at high-temperature conditions, notably differences in compatibility within the crystal structure

Inverse bifurcation analysis: application to simple gene systems

BACKGROUND: Bifurcation analysis has proven to be a powerful method for understanding the qualitative behavior of gene regulatory networks. In addition to the more traditional forward problem of determining the mapping from parameter space to the space of model behavior, the inverse problem of determining model parameters to result in certain desired properties of the bifurcation diagram provides an attractive methodology for addressing important biological problems. These include understanding how the robustness of qualitative behavior arises from system design as well as providing a way to engineer biological networks with qualitative properties. RESULTS: We demonstrate that certain inverse bifurcation problems of biological interest may be cast as optimization problems involving minimal distances of reference parameter sets to bifurcation manifolds. This formulation allows for an iterative solution procedure based on performing a sequence of eigen-system computations and one-parameter continuations of solutions, the latter being a standard capability in existing numerical bifurcation software. As applications of the proposed method, we show that the problem of maximizing regions of a given qualitative behavior as well as the reverse engineering of bistable gene switches can be modelled and efficiently solved

Electron spin coherence exceeding seconds in high purity silicon

Silicon is undoubtedly one of the most promising semiconductor materials for spin-based information processing devices. Its highly advanced fabrication technology facilitates the transition from individual devices to large-scale processors, and the availability of an isotopically-purified $^{28}$Si form with no magnetic nuclei overcomes what is a main source of spin decoherence in many other materials. Nevertheless, the coherence lifetimes of electron spins in the solid state have typically remained several orders of magnitude lower than what can be achieved in isolated high-vacuum systems such as trapped ions. Here we examine electron spin coherence of donors in very pure $^{28}$Si material, with a residual $^{29}$Si concentration of less than 50 ppm and donor densities of $10^{14-15}$ per cm$^3$. We elucidate three separate mechanisms for spin decoherence, active at different temperatures, and extract a coherence lifetime $T_2$ up to 2 seconds. In this regime, we find the electron spin is sensitive to interactions with other donor electron spins separated by ~200 nm. We apply a magnetic field gradient in order to suppress such interactions and obtain an extrapolated electron spin $T_2$ of 10 seconds at 1.8 K. These coherence lifetimes are without peer in the solid state by several orders of magnitude and comparable with high-vacuum qubits, making electron spins of donors in silicon ideal components of a quantum computer, or quantum memories for systems such as superconducting qubits.Comment: 18 pages, 4 figures, supplementary informatio

A comparison of genomic copy number calls by Partek Genomics Suite, Genotyping Console and Birdsuite algorithms to quantitative PCR

<p>Abstract</p> <p>Background</p> <p>Copy number variants are >1 kb genomic amplifications or deletions that can be identified using array platforms. However, arrays produce substantial background noise that contributes to high false discovery rates of variants. We hypothesized that quantitative PCR could finitely determine copy number and assess the validity of calling algorithms.</p> <p>Results</p> <p>Using data from 29 Affymetrix SNP 6.0 arrays, we determined copy numbers using three programs: Partek Genomics Suite, Affymetrix Genotyping Console 2.0 and Birdsuite. We compared array calls at 25 chromosomal regions to those determined by qPCR and found nearly identical calls in regions of copy number 2. Conversely, agreement differed in regions called variant by at least one method. The highest overall agreement in calls, 91%, was between Birdsuite and quantitative PCR. Partek Genomics Suite calls agreed with quantitative PCR 76% of the time while the agreement of Affymetrix Genotyping Console 2.0 with quantitative PCR was 79%.</p> <p>Conclusions</p> <p>In 38 independent samples, 96% of Birdsuite calls agreed with quantitative PCR. Analysis of three copy number calling programs and quantitative PCR showed Birdsuite to have the greatest agreement with quantitative PCR.</p

The Characterisation of Three Types of Genes that Overlie Copy Number Variable Regions

Background: Due to the increased accuracy of Copy Number Variable region (CNV) break point mapping, it is now possible to say with a reasonable degree of confidence whether a gene (i) falls entirely within a CNV; (ii) overlaps the CNV or (iii) actually contains the CNV. We classify these as type I, II and III CNV genes respectively. Principal Findings: Here we show that although type I genes vary in copy number along with the CNV, most of these type I genes have the same expression levels as wild type copy numbers of the gene. These genes must, therefore, be under homeostatic dosage compensation control. Looking into possible mechanisms for the regulation of gene expression we found that type I genes have a significant paucity of genes regulated by miRNAs and are not significantly enriched for monoallelically expressed genes. Type III genes, on the other hand, have a significant excess of genes regulated by miRNAs and are enriched for genes that are monoallelically expressed. Significance: Many diseases and genomic disorders are associated with CNVs so a better understanding of the different ways genes are associated with normal CNVs will help focus on candidate genes in genome wide association studies

Genome-Wide Analysis of Copy Number Variation in Type 1 Diabetes

Type 1 diabetes (T1D) tends to cluster in families, suggesting there may be a genetic component predisposing to disease. However, a recent large-scale genome-wide association study concluded that identified genetic factors, single nucleotide polymorphisms, do not account for overall familiality. Another class of genetic variation is the amplification or deletion of >1 kilobase segments of the genome, also termed copy number variations (CNVs). We performed genome-wide CNV analysis on a cohort of 20 unrelated adults with T1D and a control (Ctrl) cohort of 20 subjects using the Affymetrix SNP Array 6.0 in combination with the Birdsuite copy number calling software. We identified 39 CNVs as enriched or depleted in T1D versus Ctrl. Additionally, we performed CNV analysis in a group of 10 monozygotic twin pairs discordant for T1D. Eleven of these 39 CNVs were also respectively enriched or depleted in the Twin cohort, suggesting that these variants may be involved in the development of islet autoimmunity, as the presently unaffected twin is at high risk for developing islet autoimmunity and T1D in his or her lifetime. These CNVs include a deletion on chromosome 6p21, near an HLA-DQ allele. CNVs were found that were both enriched or depleted in patients with or at high risk for developing T1D. These regions may represent genetic variants contributing to development of islet autoimmunity in T1D