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.     

<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}\).     

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.     

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.     

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.     

Trend analysis of variations in carbon stock using stock big data

Changes in land use affect the terrestrial carbon stock through changes in the land cover. Research on land use and analysis of variations in carbon stock have practical applications in the optimization of land use and the mitigation of climate change effects. This study was conducted in Baixiang and Julu counties in the Taihang Piedmont by employing the trend analysis method to characterize the variation in county land use and carbon stock. The findings show that in both counties, agricultural and unused land areas decreased while built-up land area increased, and the reduction in cropland was the main reason behind the agricultural land reduction. An inflection point appeared on the cropland curves of Julu, because the cropland area decreased by 1576.97 hm\(^{2}\) from 2004 to 2006. Cropland area in Baixiang decreased from 1996 to 1998 by a total of 129.89 hm\(^{2}\) and then remained relatively stable after 1998. The total carbon storage and variation in land use in the two counties displayed similar trends. Total carbon reserves in Julu increased by 2.76 \(\times \) 10\(^{4}\) tC (carbon equivalent), while those in Baixiang decreased by 0.63 \(\times \) 10\(^{4}\) tC. Carbon stock of built-up land in Julu and Baixiang increased by 2.44 \(\times \) 10\(^{4}\) and 1.22 \(\times \) 10\(^{4}\) tC, respectively.     

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.     

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.     

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}\)).     

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.     

Joint non-uniform sampling of all incremented time delays for quicker acquisition in protein relaxation studies

NMR relaxometry plays crucial role in studies of protein dynamics. The measurement of longitudinal and transverse relaxation rates of \(^{15}\)N is the main source of information on backbone motions. However, even the most basic approach exploiting a series of \(^{15}\)N HSQC spectra can require several hours of measurement time. Standard non-uniform sampling (NUS), i.e. random under-sampling of indirect time domain, typically cannot reduce this by more than 2–4\(\times\) due to relatively low “compressibility” of these spectra. In this paper we propose an extension of NUS to relaxation delays. The two-dimensional space of \(t_1\)/\(t_{relax}\) is sampled in a way similar to NUS of \(t_1\)/\(t_2\) domain in 3D spectra. The signal is also processed in a way similar to that known from 3D NUS spectra i.e. using one of the most popular compressed sensing algorithms, iterative soft thresholding. The 2D Fourier transform matrix is replaced with mixed inverse Laplace-Fourier transform matrix. The peak positions in resulting 3D spectrum are characterized by two frequency coordinates and relaxation rate and thus no additional fitting of exponential curves is required. The method is tested on three globular proteins, providing satisfactory results in a time corresponding to acquisition of two conventional \(^{15}\)N HSQC spectra.     

Understanding genetic diversity,spatial genetic structure,and mating system through microsatellite markers for the conservation and sustainable use of <Emphasis Type="Italic">Acrocomia aculeata</Emphasis> (Jacq.) Lodd. Ex Mart.

Acrocomia aculeata is a native palm distributed widely throughout Brazil that is used in a diverse array of products from the food industry to biodiesel oil production. This study uses nine microsatellite loci to assess the genetic diversity, spatial genetic structure (SGS), and mating system of A. aculeata. A total of 200 samples were collected from four populations (Fusquinha—FU, Padre Josimo—PJ, Gleba XV-GB, and Amparo—AM), in São Paulo State, Brazil. We also collected fruits from 20 randomly selected seed trees in the FU population to assess mating patterns, for a total of 246 genotyped embryos. We identified a total of 103 alleles and all loci were polymorphic. The average observed heterozygosity (\({H_o}\)) ranged from 0.410 (AM) to 0.531 (FU) and expected heterozygosity (\({H_e}\)) ranged from 0.547 (PJ) to 0.615 (GB). The average fixation index (\(F\)) ranged from 0.043 to 0.254 for FU and AM populations, respectively. The coancestry coefficient (\({\theta _{xy}}\)) was significant up to 38 m in PJ the population and 71 m in AM. Individual palm outcrossing rates were predominantly high (\({t_m}\)?=?0.985) and the paternity correlation (\({r_{{p_{(m)}}}}\)?=?0.02) was significantly low, indicating a high probability of the occurrence of half-sibs. Compared to other palm species, the studied populations show high levels of genetic diversity. Our results confirm that A. aculeata is primarily allogamous, with no significant paternity correlation, and seeds should be harvested from at least 40 trees to ensure high levels of genetic diversity in seed collection programs.     

Simulated bone remodeling around tilted dental implants in the anterior maxilla

Dental implants have to be placed with the long axis in different angulations due to the change in bone morphology. The objective of this study was to investigate the different bone remodeling response induced by the tilted dental implants and to assess whether it could lead to bone loss and implant failure. In this study, bone remodeling due to palato-labially inclined dental implants placed in the anterior maxillary incisor region was simulated. CT-based finite element models of a maxillary bone with dental implants were created herein. Five dental implants were placed at \(+10^{\circ }\), \(+5^{\circ }\), \(0^{\circ }\), \(-5^{\circ }\) and \(-10^{\circ }\), respectively. The remodeling progression was recorded and compared. Model \(-10^{\circ }\) (palatal side) shows the highest bone density values, but the inclined implant at \(+10^{\circ }\) (labial side) leads to significant bone loss. From a biomechanical perspective, it is speculated that a palatally inclined implant is more likely to enhance the bone density in the maxillary anterior region, but labial inclination of implant could jeopardize its stability.     

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.     

Study of biological communities subject to imperfect detection: bias and precision of community <Emphasis Type="Italic">N</Emphasis>-mixture abundance models in small-sample situations

Community N-mixture abundance models for replicated counts provide a powerful and novel framework for drawing inferences related to species abundance within communities subject to imperfect detection. To assess the performance of these models, and to compare them to related community occupancy models in situations with marginal information, we used simulation to examine the effects of mean abundance \((\bar{\lambda }\): 0.1, 0.5, 1, 5), detection probability \((\bar{p}\): 0.1, 0.2, 0.5), and number of sampling sites (n _site : 10, 20, 40) and visits (n _visit : 2, 3, 4) on the bias and precision of species-level parameters (mean abundance and covariate effect) and a community-level parameter (species richness). Bias and imprecision of estimates decreased when any of the four variables \((\bar{\lambda }\), \(\bar{p}\), n _site , n _visit ) increased. Detection probability \(\bar{p}\) was most important for the estimates of mean abundance, while \(\bar{\lambda }\) was most influential for covariate effect and species richness estimates. For all parameters, increasing n _site was more beneficial than increasing n _visit . Minimal conditions for obtaining adequate performance of community abundance models were n _site  ≥ 20, \(\bar{p}\) ≥ 0.2, and \(\bar{\lambda }\) ≥ 0.5. At lower abundance, the performance of community abundance and community occupancy models as species richness estimators were comparable. We then used additive partitioning analysis to reveal that raw species counts can overestimate β diversity both of species richness and the Shannon index, while community abundance models yielded better estimates. Community N-mixture abundance models thus have great potential for use with community ecology or conservation applications provided that replicated counts are available.     

Plausibility and parameter sensitivity of micro-finite element-based joint load prediction at the proximal femur

A micro-finite element-based method to estimate the bone loading history based on bone architecture was recently presented in the literature. However, a thorough investigation of the parameter sensitivity and plausibility of this method to predict joint loads is still missing. The goals of this study were (1) to analyse the parameter sensitivity of the joint load predictions at one proximal femur and (2) to assess the plausibility of the results by comparing load predictions of ten proximal femora to in vivo hip joint forces measured with instrumented prostheses (available from www.orthoload.com). Joint loads were predicted by optimally scaling the magnitude of four unit loads (inclined \(-20^{\circ }\) to \(100^{\circ }\) with respect to the vertical axis) applied to micro-finite element models created from high-resolution computed tomography scans (\(30.3~\upmu \)m voxel size). Parameter sensitivity analysis was performed by varying a total of nine parameters and showed that predictions of the peak load directions (range 10\(^{\circ }\)–\(30^{\circ }\)) are more robust than the predicted peak load magnitudes (range 2344.8–4689.5 N). Comparing the results of all ten femora with the in vivo loading data of ten subjects showed that peak loads are plausible both in terms of the load direction (in vivo: \(18.2\pm 2.0^{\circ }\), predicted: \(20.0^{\circ }\)) and magnitude (in vivo: \(2707.6\pm 443.3~\hbox {N}\), predicted: \(3372.2\pm 597.9~\hbox {N}\)). Overall, this study suggests that micro-finite element-based joint load predictions are both plausible and robust in terms of the predicted peak load direction, but predicted load magnitudes should be interpreted with caution.     

Long-term grazing affects relationships between nitrogen form uptake and biomass of alpine meadow plants

Trade-off between nutrient uptake rate and product accumulation has been found among species characterized as acquisitive and conservative strategies in resource utilization. However, long-term grazing causes changes in soil nutrient availability and plant species abundance by selective foraging and resource allocation between above- and belowground organs, which may cover up such trade-off. However, little is known whether the trade-off can be observed among species in community without grazing disturbance, and how grazing influences the trade-off. We conducted a ¹⁵N labelling experiment in winter grazing and grazing release alpine meadow communities on the Tibet Plateau. We examined changes in N form uptake of 11 common species and relationship of N chemical uptake rate with aboveground biomass. Grazing release increased soil \({\text{NH}}_{4}^{ + }\)–N and \({\text{NO}}_{3}^{ - }\)–N, and increased \({\text{NO}}_{3}^{ - }\)–N uptake rate in two species and \({\text{NH}}_{4}^{ + }\)–N uptake rate in three species. Meanwhile, grazing release decreased aboveground biomass of three species and two of them belong to those species’ increased N uptake rate. Contrarily, grazing release increased aboveground biomass of four species and none belongs to the changed N uptake rate. Thus, grazing release caused changes of plant nutrient uptake rate and aboveground production in different directions, which explains the negative relationship of N uptake rate with aboveground biomass in ungrazed community. Our results indicate that the increase in nutrient uptake is probably one of the mechanisms for acquisitive species to cope with the raising nutrient availability and/or competition from the conservative dominant grasses after grazing release.     

Isolation and characterization of microsatellite markers in the spiny lobster, <Emphasis Type="Italic">Panulirus echinatus</Emphasis> Smith, 1869 (Decapoda: Palinuridae) by Illumina MiSeq sequencing

Okra’s (Abelmoschus esculentus (L.) Moench) commercial cultivation is threatened in the tropics due to high incidence of yellow vein mosaic virus (YVMV) disease. Okra geneticists across the world tried to understand the inheritance pattern of YVMV disease tolerance without much success. Therefore, the inheritance pattern of YVMV disease in okra was revisited by employing six generations (\(\hbox {P}_{1}\), \(\hbox {P}_{2}\), \(\hbox {F}_{1}\), \(\hbox {F}_{2}\), \(\hbox {BC}_{1}\) and \(\hbox {BC}_{2}\)) of four selected crosses (one tolerant \(\times \) tolerant, two tolerant \(\times \) susceptible and one susceptible \(\times \) susceptible) using two tolerant (BCO-1 and Lal Bhendi) and two susceptible (Japanese Jhar Bhendi and PAN 2127) genotypes. Qualitative genetic analysis was done on the basis of segregation pattern of tolerant and susceptible plants in \(\hbox {F}_{2}\) and backcross generations of all the four crosses. It revealed that a single dominant gene along with some minor factors governed the disease tolerant trait in both the tolerant parents used. However, it was observed that genes governing disease tolerance identified in both the tolerant variety used was different. It could be concluded that the gene governing YVMV disease tolerance in okra was genotype specific. Further, duplicate gene action as evident from an approximate ratio of 15 : 1 (tolerant : susceptible) in the \(\hbox {F}_{2}\) population in the cross of two tolerant varieties gave a scope of increasing the tolerance level of the hybrid plants when both the tolerant genes are brought together. However, generation mean analysis revealed involvement of both additive and nonadditive effects in the inheritance of disease tolerance. Thus, the present study confirms that a complicated genetic inheritance pattern is involved in the disease tolerance against YVMV trait. The major tolerance genes could be transferred to other okra varieties, but the tolerance breaking virus strains might not allow them to achieve tolerance in stable condition. Therefore, accumulation of additional genes may be needed for a sustainable tolerance phenotype in okra.