A finite-element model of mechanosensation by a Pacinian corpuscle cluster in human skin

The Pacinian corpuscle (PC) is the cutaneous mechanoreceptor responsible for sensation of high-frequency (20–1000 Hz) vibrations. PCs lie deep within the skin, often in multicorpuscle clusters with overlapping receptive fields. We developed a finite-element mechanical model of one or two PCs embedded within human skin, coupled to a multiphysics PC model to simulate action potentials elicited by each PC. A vibration was applied to the skin surface, and the resulting mechanical signal was analyzed using two metrics: the deformation amplitude ratio (\({\rho }_{\mathrm{1S}} \), \({\rho }_{\mathrm{2S}} )\) and the phase shift of the vibration (\({\delta }_{\mathrm{S}1}^{\mathrm{mech}} \), \({\delta }_{\mathrm{S}2}^{\mathrm{mech}} )\) between the stimulus and the PC. Our results showed that the amplitude attenuation and phase shift at a PC increased with distance from the stimulus to the PC. Differences in amplitude (\(\rho _{12} )\) and phase shift (\({\delta }_{12}^{\mathrm{mech}} )\) between the two PCs in simulated clusters directly affected the interspike interval between the action potentials elicited by each PC (\({\delta }_{12}^{\mathrm{spike}} )\). While \({\delta }_{12}^{\mathrm{mech}} \) had a linear relationship with \({\delta }_{12}^{\mathrm{spike}} \), \(\rho _{12} \)’s effect on \({\delta }_{12}^{\mathrm{spike}} \) was greater for lower values of \(\rho _{12} \). In our simulations, the separation between PCs and the distance of each PC from the stimulus location resulted in differences in amplitude and phase shift at each PC that caused \({\delta }_{12}^{\mathrm{spike}} \) to vary with PC location. Our results suggest that PCs within a cluster receive different mechanical stimuli which may enhance source localization of vibrotactile stimuli, drawing parallels to sound localization in binaural hearing.     

Potential use of sulfite as a supplemental electron donor for wastewater denitrification

Biological denitrification typically requires the addition of a supplemental electron donor, which can add a significant operating expense to wastewater treatment facilities. Most common electron donors are organic, but reduced inorganic sulfur compounds (RISCs), such as sulfide (HS^?) and elemental sulfur (S⁰), may be more cost-effective. S⁰ is an inexpensive and well characterized electron donor, but it provides slow denitrification rates due to its low solubility. A lesser-known RISC is sulfite (\({\text{SO}}_{3}^{2 - }\)), which can be easily produced from S⁰ by a simple combustion process. Unlike S⁰, \({\text{SO}}_{3}^{2 - }\) is highly soluble, and therefore may provide higher denitrification rates. However, very little is known about microbial denitrification with \({\text{SO}}_{3}^{2 - }\). Also, \({\text{SO}}_{3}^{2 - }\) is a strong reductant that reacts abiotically with oxygen and has toxic effects on microorganisms. This paper reviews \({\text{SO}}_{3}^{2 - }\) in the environment, \({\text{SO}}_{3}^{2 - }\) chemistry, microbiology, toxicity, and its potential use for denitrification. Since \({\text{SO}}_{3}^{2 - }\) is an intermediate in the sulfur oxidation pathway of most sulfur-oxidizing microorganisms, it is an energetic electron donor and it should select for a \({\text{SO}}_{3}^{2 - }\)-oxidizing community. Our review of the literature, as well as our own lab experience, suggests that \({\text{SO}}_{3}^{2 - }\) can effectively serve as an electron donor for denitrification. Further research is needed to determine the kinetics of \({\text{SO}}_{3}^{2 - }\)-based denitrification, its toxic threshold for sulfur-oxidizing microorganisms, and its potential inhibition of sensitive species such as nitrifying microorganisms and potential formation of nitrous oxide. Its effect on sludge settling efficiency also should be explored.     

Physical parameter estimation from porcine ex vivo vocal fold dynamics in an inverse problem framework

This study presents a framework for a direct comparison of experimental vocal fold dynamics data to a numerical two-mass-model (2MM) by solving the corresponding inverse problem of which parameters lead to similar model behavior. The introduced 2MM features improvements such as a variable stiffness and a modified collision force. A set of physiologically sensible degrees of freedom is presented, and three optimization algorithms are compared on synthetic vocal fold trajectories. Finally, a total of 288 high-speed video recordings of six excised porcine larynges were optimized to validate the proposed framework. Particular focus lay on the subglottal pressure, as the experimental subglottal pressure is directly comparable to the model subglottal pressure. Fundamental frequency, amplitude and objective function values were also investigated. The employed 2MM is able to replicate the behavior of the porcine vocal folds very well. The model trajectories’ fundamental frequency matches the one of the experimental trajectories in \(98.6\%\) of the recordings. The relative error of the model trajectory amplitudes is on average \(9.5\%\). The experiments feature a mean subglottal pressure of 10.16 (SD \(= 2.31\)) \({\text {cmH}}_2{\text {O}}\); in the model, it was on average 7.61 (SD \(= 2.40\)) \({\text {cmH}}_2{\text {O}}\). A tendency of the model to underestimate the subglottal pressure is found, but the model is capable of inferring trends in the subglottal pressure. The average absolute error between the subglottal pressure in the model and the experiment is 2.90 (SD \(= 1.80\)) \({\text {cmH}}_2{\text {O}}\) or \(27.5\%\). A detailed analysis of the factors affecting the accuracy in matching the subglottal pressure is presented.     

Uptake of ammonium and soluble reactive phosphorus in forested streams: influence of dissolved organic matter composition

Many microbes responsible for inorganic nutrient uptake and transformation utilize dissolved organic matter (DOM) as a nutrient or energy source, but little is known about whether DOM composition is an important driver of nutrient uptake in streams. Our goal was to determine whether incorporating DOM composition metrics with other more commonly considered biological, physical, and chemical variables improved our ability to explain patterns of ammonium (\({\text{NH}}_{4}^{ + }\)–N) and soluble reactive phosphorus (SRP) uptake across 11 Lake Superior tributaries. Nutrient uptake velocities (V_f) ranged from undetectable to 14.6 mm min^?1 for \({\text{NH}}_{4}^{ + }\)–N and undetectable to 7.2 mm min^?1 for SRP. Logistic regressions suggested that DOM composition was a useful predictor of where SRP uptake occurred (4/11 sites) and \({\text{NH}}_{4}^{ + }\)–N concentration was a useful predictor of where \({\text{NH}}_{4}^{ + }\)–N uptake occurred (9/11 sites). Multiple regression analysis revealed that the best models included temperature, specific discharge, and canopy cover, and DOM composition as significant predictors of \({\text{NH}}_{4}^{ + }\)–N V_f. Partial least squares revealed fluorescence index (describing the source of aquatic fulvic acids), specific ultraviolet absorbance at 254 nm (an indicator of DOM aromaticity), temperature, and conductivity were highly influential predictors of \({\text{NH}}_{4}^{ + }\)–N V_f. Therefore, streams with higher temperatures, lower solute concentrations, more terrestrial DOM signal and greater aromaticity had greater \({\text{NH}}_{4}^{ + }\)–N V_f. Our results suggest that DOM composition may be an important, yet often overlooked, predictor of \({\text{NH}}_{4}^{ + }\)–N and SRP uptake in deciduous forest streams that should be considered along with commonly measured predictors.     

<Superscript>1</Superscript>H, <Superscript>13</Superscript>C, <Superscript>15</Superscript>N resonance assignments of the extracellular loop 1 domain (ECL1) of <Emphasis Type="Italic">Streptococcus pneumoniae</Emphasis> D39 FtsX,an essential cell division protein

FtsX is an integral membrane protein from Streptococcus pneumoniae (pneumococcus) that harbors an extracellular loop 1 domain (\({\text{FtsX}}^{\text{ECL1}}_{Spn}\)) that interacts with PcsB, an peptidoglycan hydrolase that is essential for cell growth and division. Here, we report nearly complete backbone and side chain resonance assignments and a secondary structural analysis of \({\text{FtsX}}^{\text{ECL1}}_{Spn}\) (residues 47–168 of FtsX) as first steps toward structure determination of \({\text{FtsX}}^{\text{ECL1}}_{Spn}\).     

Cell Volume Regulation in the Proximal Tubule of Rat Kidney

We developed a dynamic model of a rat proximal convoluted tubule cell in order to investigate cell volume regulation mechanisms in this nephron segment. We examined whether regulatory volume decrease (RVD), which follows exposure to a hyposmotic peritubular solution, can be achieved solely via stimulation of basolateral K\(^+\) and \(\hbox {Cl}^-\) channels and \(\hbox {Na}^+\)–\(\hbox {HCO}_3^-\) cotransporters. We also determined whether regulatory volume increase (RVI), which follows exposure to a hyperosmotic peritubular solution under certain conditions, may be accomplished by activating basolateral \(\hbox {Na}^+\)/H\(^+\) exchangers. Model predictions were in good agreement with experimental observations in mouse proximal tubule cells assuming that a 10% increase in cell volume induces a fourfold increase in the expression of basolateral K\(^+\) and \(\hbox {Cl}^-\) channels and \(\hbox {Na}^+\)–\(\hbox {HCO}_3^-\) cotransporters. Our results also suggest that in response to a hyposmotic challenge and subsequent cell swelling, \(\hbox {Na}^+\)–\(\hbox {HCO}^-_3\) cotransporters are more efficient than basolateral K\(^+\) and \(\hbox {Cl}^-\) channels at lowering intracellular osmolality and reducing cell volume. Moreover, both RVD and RVI are predicted to stabilize net transcellular \(\hbox {Na}^+\) reabsorption, that is, to limit the net \(\hbox {Na}^+\) flux decrease during a hyposmotic challenge or the net \(\hbox {Na}^+\) flux increase during a hyperosmotic challenge.     

Modeling of chromosome intermingling by partially overlapping uniform random polygons

During the early phase of the cell cycle the eukaryotic genome is organized into chromosome territories. The geometry of the interface between any two chromosomes remains a matter of debate and may have important functional consequences. The Interchromosomal Network model (introduced by Branco and Pombo) proposes that territories intermingle along their periphery. In order to partially quantify this concept we here investigate the probability that two chromosomes form an unsplittable link. We use the uniform random polygon as a crude model for chromosome territories and we model the interchromosomal network as the common spatial region of two overlapping uniform random polygons. This simple model allows us to derive some rigorous mathematical results as well as to perform computer simulations easily. We find that the probability that one uniform random polygon of length n that partially overlaps a fixed polygon is bounded below by \({1-O(\frac{1}{\sqrt n})}\). We use numerical simulations to estimate the dependence of the linking probability of two uniform random polygons (of lengths n and m, respectively) on the amount of overlapping. The degree of overlapping is parametrized by a parameter \({\epsilon\in [0,1]}\) such that \({\epsilon=0}\) indicates no overlapping and \({\epsilon=1}\) indicates total overlapping. We propose that this dependence relation may be modeled as \({f(\varepsilon, m, n) =1-{\frac{a(\epsilon)}{b(\epsilon)\sqrt{mn}+c(\epsilon)}}}\). Numerical evidence shows that this model works well when \({\epsilon}\) is relatively large \({(\varepsilon \ge 0.5)}\). We then use these results to model the data published by Branco and Pombo and observe that for the amount of overlapping observed experimentally the URPs have a non-zero probability of forming an unsplittable link.     

Mathematical Analysis of the Transmission Dynamics of HIV Syphilis Co-infection in the Presence of Treatment for Syphilis

The re-emergence of syphilis has become a global public health issue, and more persons are getting infected, especially in developing countries. This has also led to an increase in the incidence of human immunodeficiency virus (HIV) infections as some studies have shown in the recent decade. This paper investigates the synergistic interaction between HIV and syphilis using a mathematical model that assesses the impact of syphilis treatment on the dynamics of syphilis and HIV co-infection in a human population where HIV treatment is not readily available or accessible to HIV-infected individuals. In the absence of HIV, the syphilis-only model undergoes the phenomenon of backward bifurcation when the associated reproduction number (\({\mathcal {R}}_{T}\)) is less than unity, due to susceptibility to syphilis reinfection after recovery from a previous infection. The complete syphilis–HIV co-infection model also undergoes the phenomenon of backward bifurcation when the associated effective reproduction number (\({\mathcal {R}}_{C}\)) is less than unity for the same reason as the syphilis-only model. When susceptibility to syphilis reinfection after treatment is insignificant, the disease-free equilibrium of the syphilis-only model is shown to be globally asymptotically stable whenever the associated reproduction number (\({\mathcal {R}}_{T}\)) is less than unity. Sensitivity and uncertainty analysis show that the top three parameters that drive the syphilis infection (with respect to the associated response function, \({\mathcal {R}}_{T}\)) are the contact rate (\(\beta _S\)), modification parameter that accounts for the increased infectiousness of syphilis-infected individuals in the secondary stage of the infection (\(\theta _1\)) and treatment rate for syphilis-only infected individuals in the primary stage of the infection (\(r_1\)). The co-infection model was numerically simulated to investigate the impact of various treatment strategies for primary and secondary syphilis, in both singly and dually infected individuals, on the dynamics of the co-infection of syphilis and HIV. It is observed that if concerted effort is exerted in the treatment of primary and secondary syphilis (in both singly and dually infected individuals), especially with high treatment rates for primary syphilis, this will result in a reduction in the incidence of HIV (and its co-infection with syphilis) in the population.     

Study of the Population Dynamics of <Emphasis Type="Italic">Busseola fusca</Emphasis>, Maize Pest

Busseola fusca is a maize and sorghum pest that can cause significant damage to both crops. Given that maize is one of the main cereals grown in the worldwide, this pest is a major challenge for maize production and therefore for the economies of several countries . In this paper , based on the life cycle of B. fusca, we propose a mathematical model to study the population dynamics of this insect pest . A sensitivity analysis using the eFast method was performed to show the most important parameters of the model. We present the theoretical analysis of the model. More precisely, we derive a threshold parameter \({\mathcal {N}}_0\), called basic offspring number and show that the trivial equilibrium is globally asymptotically stable whenever \({\mathcal {N}}_0\le 1\), while if \({\mathcal {N}}_0>1\), the non trivial equilibrium is globally asymptotically stable. The theoretical results are supported by numerical simulations.     

Qualitative Analysis of an ODE Model of a Class of Enzymatic Reactions

The present paper analyzes an ODE model of a certain class of (open) enzymatic reactions. This type of model is used, for instance, to describe the interactions between messenger RNAs and microRNAs. It is shown that solutions defined by positive initial conditions are well defined and bounded on \([0, \infty )\) and that the positive octant of \({\mathbb {R}}^3\) is a positively invariant set. We prove further that in this positive octant there exists a unique equilibrium point, which is asymptotically stable and a global attractor for any initial state with positive components; a controllability property is emphasized. We also investigate the qualitative behavior of the QSSA system in the phase plane \({\mathbb {R}}^2\). For this planar system we obtain similar results regarding global stability by using Lyapunov theory, invariant regions and controllability.     

Performance and kinetics of ANAMMOX granular sludge with pH shock in a sequencing batch reactor

As an efficient and cost-effective nitrogen removal process, anaerobic ammonium oxidation (ANAMMOX) could be well operated at suitable pH condition. However, pH shock occurred in different kinds of wastewater and affected ANANNOX process greatly. The present research aimed at studying the performance and kinetics of ANAMMOX granular sludge with pH shock. When influent pH was below 7.5, effluent \({\text{NH}}_{4}^{ + }\)–N and \({\text{NO}}_{2}^{ - }\)–N increased with decreasing pH. At Ph 6.0, effluent \({\text{NO}}_{2}^{ - }\)–N approached 100 mg/L, and the ratios of \(\Delta {\text{NO}}_{2}^{ - } - {\text{N}}:\Delta {\text{NH}}_{4}^{ + } - {\text{N and }}\Delta {\text{NO}}_{3}^{ - } - {\text{N}}:\Delta {\text{NH}}_{4}^{ + } - {\text{N}}\) approached 2.2 and 1.3, respectively. Both greatly deviated from theoretical values. When influent pH was above 7.5, effluent \({\text{NH}}_{4}^{ + }\)–N and \({\text{NO}}_{2}^{ - }\)–N increased with increasing pH. At pH 9.0, ammonium removal rate (ARR) and nitrite removal rate (NRR) decreased to 0.011 ± 0.004 and 0.035 ± 0.004 kg/(m³·d), respectively. Besides, \(\Delta {\text{NO}}_{2}^{ - }\)–N:\(\Delta {\text{NH}}_{4}^{ + }\)–N deviated from theoretical value. Longer recovery time from pH 9.0 than from pH 6.0 indicated that alkaline surroundings inhibited anaerobic ammonium oxidizing bacteria (AAOB) greater. The sludge settling velocity was 2.15 cm/s at pH 7.5. However, it decreased to 2.02 cm/s when pH was 9.0. Acidic pH had little effect on sludge size, but disintegration of ANAMMOX granule was achieved with pH of 9.0. The Bell-shaped (A) model and the Ratkowsky model were more applicable to simulate the effect resulting from pH shock on ANAMMOX activity (R² > 0.95), and both could describe ANAMMOX activity well with pH shock. They indicated that q_max was 0.37 kg \(\Delta {\text{NH}}_{4}^{ + }\)–N/(kgMLSS·d) at the optimum pH value (7.47) in present study. The minimum pH during which ANAMMOX occurred was 5.68 while the maximum pH for ANAMMOX reaction was 9.26. Based on nitrogen removal performance with different pH, strongly acidic (pH ≤ 6.5) or alkaline (pH ≥ 8.5) inhibited ANAMMOX process. Besides, ANAMMOX appeared to be more susceptible to alkaline wastewater. Compared to extremely acidic condition (low pH), extremely alkaline condition (high pH) affected ANAMMOX granules much more.     

Comparison of the Functional Activities of Xanthine Oxidases Isolated from Microorganisms and from Cow’s Milk

The characteristics of the formation of the superoxide radical anion (\(\rm{O}_2^{\bullet-}\)) and hydrogen peroxide by xanthine oxidases isolated from microorganisms and from cow’s milk were investigated. The increase in pH led to an increase in the rate of xanthine oxidation with oxygen by both xanthine oxidases. The functioning of xanthine oxidase from milk along with the two-electron reduction of O₂ to H₂O₂ carries through the one-electron reduction of O₂ to \(\rm{O}_2^{\bullet-}\), and the rate and the fraction of generation of \(\rm{O}_2^{\bullet-}\) increased with increasing pH. Under operation of the microbial xanthine oxidase, the \(\rm{O}_2^{\bullet-}\) radical was not detected in the medium. The results suggest a difference in the operation of active centers of enzyme from different sources.     

When is a Phylogenetic Network Simply an Amalgamation of Two Trees?

Phylogenetic networks generalise phylogenetic (evolutionary) trees by allowing for the representation of reticulation (non-treelike) events. The structure of such networks is often viewed by the phylogenetic trees they embed. In this paper, we determine when a phylogenetic network \({\mathcal {N}}\) has two phylogenetic tree embeddings which collectively contain all of the edges of \({\mathcal {N}}\). This determination leads to a polynomial-time algorithm for recognising such networks and an unexpected characterisation of the class of reticulation-visible networks.     

Measuring the signs of the methyl <Superscript>1</Superscript>H chemical shift differences between major and ‘invisible’ minor protein conformational states using methyl <Superscript>1</Superscript>H multi-quantum spectroscopy

Carr–Purcell–Meiboom–Gill (CPMG) type relaxation dispersion experiments are now routinely used to characterise protein conformational dynamics that occurs on the μs to millisecond (ms) timescale between a visible major state and ‘invisible’ minor states. The exchange rate(s) (\( k_{{{\text{ex}}}} \)), population(s) of the minor state(s) and the absolute value of the chemical shift difference \(|{\Delta \varpi }|\) (ppm) between different exchanging states can be extracted from the CPMG data. However the sign of \({\Delta \varpi }\) that is required to reconstruct the spectrum of the ‘invisible’ minor state(s) cannot be obtained from CPMG data alone. Building upon the recently developed triple quantum (TQ) methyl \( ^{1} {\text{H}} \) CPMG experiment (Yuwen in Angew Chem 55:11490–11494, 2016) we have developed pulse sequences that use carbon detection to generate and evolve single quantum (SQ), double quantum (DQ) and TQ coherences from methyl protons in the indirect dimension to measure the chemical exchange-induced shifts of the SQ, DQ and TQ coherences from which the sign of \({\Delta \varpi }\) is readily obtained for two state exchange. Further a combined analysis of the CPMG data and the difference in exchange induced shifts between the SQ and DQ resonances and between the SQ and TQ resonances improves the estimates of exchange parameters like the population of the minor state. We demonstrate the use of these experiments on two proteins undergoing exchange: (1) the ~ 18 kDa cavity mutant of T4 Lysozyme (\( k_{{{\text{ex}}}} \sim\,3500{\text{ s}}^{{ - 1}} \)) and (2) the \(\sim\,4.7\) kDa Peripheral Sub-unit Binding Domain (PSBD) from the acetyl transferase of Bacillus stearothermophilus (\(k_{ex} \sim\,13,000\hbox { s}^{-1}\)).     

A Mathematical Model for the Macrophage Response to Respiratory Viral Infection in Normal and Asthmatic Conditions

Respiratory viral infections are common in the general population and one of the most important causes of asthma aggravation and exacerbation. Despite many studies, it is not well understood how viral infections cause more severe symptoms and exacerbations in asthmatics. We develop a mathematical model of two types of macrophages that play complementary roles in fighting viral infection: classically \((\hbox {CA}\)-\(\hbox {M}\Phi )\) and alternatively activated macrophages \((\hbox {AA}\)-\(\hbox {M}\Phi )\). \(\hbox {CA}\)-\(\hbox {M}\Phi \) destroy infected cells and tissues to remove viruses, while \(\hbox {AA}\)-\(\hbox {M}\Phi \) repair damaged tissues. We show that a higher viral load or longer duration of infection provokes a stronger immune response from the macrophage system. By adjusting the parameters, we model the differences in response to respiratory viral infection in normal and asthmatic subjects and show how this skews the system toward a response that generates more severe symptoms in asthmatic patients.     

Cellular mechanosensitivity to substrate stiffness decreases with increasing dissimilarity to cell stiffness

Computational modelling has received increasing attention to investigate multi-scale coupled problems in micro-heterogeneous biological structures such as cells. In the current study, we investigated for a single cell the effects of (1) different cell-substrate attachment (2) and different substrate modulus \(\textit{E}_\mathrm{s}\) on intracellular deformations. A fibroblast was geometrically reconstructed from confocal micrographs. Finite element models of the cell on a planar substrate were developed. Intracellular deformations due to substrate stretch of \(\lambda =1.1\), were assessed for: (1) cell-substrate attachment implemented as full basal contact (FC) and 124 focal adhesions (FA), respectively, and \(\textit{E}_\mathrm{s}\,=\,\)140 KPa and (2) \(\textit{E}_\mathrm{s}\,=\,10\), 140, 1000, and 10,000 KPa, respectively, and FA attachment. The largest strains in cytosol, nucleus and cell membrane were higher for FC (1.35\(\text {e}^{-2}\), 0.235\(\text {e}^{-2}\) and 0.6\(\text {e}^{-2}\)) than for FA attachment (0.0952\(\text {e}^{-2}\), 0.0472\(\text {e}^{-2}\) and 0.05\(\text {e}^{-2}\)). For increasing \(\textit{E}_\mathrm{s}\), the largest maximum principal strain was 4.4\(\text {e}^{-4}\), 5\(\text {e}^{-4}\), 5.3\(\text {e}^{-4}\) and 5.3\(\text {e}^{-4}\) in the membrane, 9.5\(\text {e}^{-4}\), 1.1\(\text {e}^{-4}\), 1.2\(\text {e}^{-3}\) and 1.2\(\text {e}^{-3}\) in the cytosol, and 4.5\(\text {e}^{-4}\), 5.3\(\text {e}^{-4}\), 5.7\(\text {e}^{-4}\) and 5.7\(\text {e}^{-4}\) in the nucleus. The results show (1) the importance of representing FA in cell models and (2) higher cellular mechanical sensitivity for substrate stiffness changes in the range of cell stiffness. The latter indicates that matching substrate stiffness to cell stiffness, and moderate variation of the former is very effective for controlled variation of cell deformation. The developed methodology is useful for parametric studies on cellular mechanics to obtain quantitative data of subcellular strains and stresses that cannot easily be measured experimentally.     

Phylogenetics beyond biology

Evolutionary processes have been described not only in biology but also for a wide range of human cultural activities including languages and law. In contrast to the evolution of DNA or protein sequences, the detailed mechanisms giving rise to the observed evolution-like processes are not or only partially known. The absence of a mechanistic model of evolution implies that it remains unknown how the distances between different taxa have to be quantified. Considering distortions of metric distances, we first show that poor choices of the distance measure can lead to incorrect phylogenetic trees. Based on the well-known fact that phylogenetic inference requires additive metrics, we then show that the correct phylogeny can be computed from a distance matrix \({\mathbf {D}}\) if there is a monotonic, subadditive function \(\zeta\) such that \(\zeta ^{-1}({\mathbf {D}})\) is additive. The required metric-preserving transformation \(\zeta\) can be computed as the solution of an optimization problem. This result shows that the problem of phylogeny reconstruction is well defined even if a detailed mechanistic model of the evolutionary process remains elusive.     

Glutathione <Emphasis Type="BoldItalic">S</Emphasis>-transferase P1 gene polymorphisms and susceptibility to coronary artery disease in a subgroup of north Indian population

The present study aimed to investigate the association of \(\hbox {g}.313\hbox {A}{>}\hbox {G}\) and \(\hbox {g}.341\hbox {C}{>}\hbox {T}\) polymorphisms of GSTP1 with coronary artery disease (CAD) in a subgroup of north Indian population. In the present case–control study, CAD patients (\(n = 200\)) and age-matched, sex-matched and ethnicity-matched healthy controls (\(n = 200\)) were genotyped for polymorphisms in GSTP1 using polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) method. Genotype distribution of \(\hbox {g}.313\hbox {A}{>}\hbox {G}\) and \(\hbox {g}.341\hbox {C}{>}\hbox {T}\) polymorphisms of GSTP1 gene was significantly different between cases and controls (\(P = 0.005\) and 0.024, respectively). Binary logistic regression analysis showed significant association of A/G (odds ratio (OR): 1.6, 95% CI: 1.08–2.49, \(P = 0.020\)) and G/G (OR: 3.1, 95% CI: 1.41–6.71, P \(=\) 0.005) genotypes of GSTP1 \(\hbox {g}.313\hbox {A}{\!>\!}\hbox {G}\), and C/T (OR: 5.8, 95% CI: 1.26–26.34, \(P = 0.024\)) genotype of GSTP1 \(\hbox {g}.341\hbox {C}{>}\hbox {T}\) with CAD. The A/G and G/G genotypes of \(\hbox {g}.313\hbox {A}{>}\hbox {G}\) and C/T genotype of \(\hbox {g}.341\hbox {C}{>}\hbox {T}\) conferred 6.5-fold increased risk for CAD (OR: 6.5, 95% CI: 1.37–31.27, \(P = 0.018\)). Moreover, the recessive model of GSTP1 \(\hbox {g}.313\hbox {A}{>}\hbox {G}\) is the best fit inheritance model to predict the susceptible gene effect (OR: 2.3, 95% CI: 1.11–4.92, \(P = 0.020\)). In conclusion, statistically significant associations of GSTP1 \(\hbox {g}.313\hbox {A}{>}\hbox {G}\) (A/G, G/G) and \(\hbox {g}.341\hbox {C}{>}\hbox {T}\) (C/T) genotypes with CAD were observed.