Upper bound estimators of the population size based on ordinal models for capture‐recapture experiments

Capture‐recapture studies have attracted a lot of attention over the past few decades, especially in applied disciplines where a direct estimate for the size of a population of interest is not available. Epidemiology, ecology, public health, and biodiversity are just a few examples. The estimation of the number of unseen units has been a challenge for theoretical statisticians, and considerable progress has been made in providing lower bound estimators for the population size. In fact, it is well known that consistent estimators for this cannot be provided in the very general case. Considering a case where capture‐recapture studies are summarized by a frequency of frequencies distribution, we derive a simple upper bound of the population size based on the cumulative distribution function. We introduce two estimators of this bound, without any specific parametric assumption on the distribution of the observed frequency counts. The behavior of the proposed estimators is investigated using several benchmark datasets and a large‐scale simulation experiment based on the scheme discussed by Pledger.     

陆地生态系统环境控制实验的研究方法及技术体系

生物及生态系统与环境变化间的反馈关系及其过程机制是生态学研究的重要内容。不同类型的生物环境因素控制实验以及大尺度的联网野外控制实验被认为是认识生态系统响应和适应环境变化过程机制、精细定量表达的有效手段及认知过程的加速器。近年来发展了大型野外物理模拟实验装置网络(如ECOTRON)、生态系统分析与实验平台(AnaEE)、国际干旱实验研究网络(Drought Network)、氮沉降联合实验网络(Nutrient Network),以及基于各区域性生态观测实验站的联网控制实验(如USA-ILTER)。发展大陆尺度联网实验研究平台事业正日益受到学术界的重视,将会在认知生态系统环境响应过程机制方面发挥更重要的作用。基于以上背景,本文综述了生态系统环境控制实验的研究方法和实验体系的发展,明确指出各种类型的生物环境控制实验需要形成联合协作体系,共同解决生态系统对环境变化的响应及适应的基本科学问题。目前的控制实验包括: 1) 实验室封闭装置内的生物生理生态学控制实验;2) 野外实验场的半开放部分环境要素控制实验;3) 近自然状态的野外环境控制实验;以及4) 基于野外生态站的联网控制实验。进而,本文还深入讨论了陆地生态系统的环境响应及适应过程机制实验系统设计的发展趋势,分析了基于大尺度自然环境梯度实验及生态站尺度的要素控制实验的优势,提出了整合两种实验技术、发展新一代的野外联网实验体系的科学设想,讨论了基于野外联网控制实验的研究体系,论证了研究生态系统对环境变化短期响应和长期适应的规律和机制、生态系统环境响应定量表达的技术途径。若本文提出的控制实验体系设计方案能够得以实施,必将大大促进我国乃至全球生态系统和环境变化科学的研究水平,对我国应对气候变化和生态环境建设具有重要的科学意义。     

Active and machine learning-based approaches to rapidly enhance microbial chemical production

In order to make renewable fuels and chemicals from microbes, new methods are required to engineer microbes more intelligently. Computational approaches, to engineer strains for enhanced chemical production typically rely on detailed mechanistic models (e.g., kinetic/stoichiometric models of metabolism)—requiring many experimental datasets for their parameterization—while experimental methods may require screening large mutant libraries to explore the design space for the few mutants with desired behaviors. To address these limitations, we developed an active and machine learning approach (ActiveOpt) to intelligently guide experiments to arrive at an optimal phenotype with minimal measured datasets. ActiveOpt was applied to two separate case studies to evaluate its potential to increase valine yields and neurosporene productivity in Escherichia coli. In both the cases, ActiveOpt identified the best performing strain in fewer experiments than the case studies used. This work demonstrates that machine and active learning approaches have the potential to greatly facilitate metabolic engineering efforts to rapidly achieve its objectives.     

Diatoms as indicators: The influences of experimental nitrogen enrichment on diatom assemblages in sub-Arctic streams

Indices using diatoms are widely used to assess water quality, but are usually constructed from field correlations and not tested through rigorous experimentation. We tested experimentally the performance of the Sørensen and the Shannon indices, and the trophic diatom index (TDI). Nitrogen was naturally limiting in the eight remote sub-Arctic streams used and we measured the effects of experimental nitrogen enrichment on diatom assemblages. Diatom densities increased significantly in the enriched reaches but there was no significant difference in invertebrate density between control and treatment reaches. Grazing effects were thus controlled for. Diversity within streams (Shannon index) was significantly reduced by nutrient addition but the Sørensen index did not change. The trophic diatom index (TDI), which is presumed to reflect nutrient concentration, was not influenced by nutrient addition and generally the values were low in both control and treatment reaches. Densities of the diatom genera Achnanthes and Gomphonema increased significantly with enrichment while those of Nitzschia and Fragilaria decreased significantly. Less abundant diatom species, which collectively constituted around 40% in relative abundance in the control reaches, were around 15–18% in treatment reaches. Growth forms were altered by the nutrients. Diatoms attached by mucilage pads were more abundant in treated reaches compared with control reaches. Motile diatoms became scarcer. The size of diatom species was unaffected by nutrient enrichment. This study showed that it is important to test experimentally indices that are developed for particular habitats before using them elsewhere.     

Binding of methacycline to human serum albumin at subdomain IIA using multispectroscopic and molecular modeling methods

This study was designed to examine the interaction of methacyline (METC) with human serum albumin (HSA) by multispectroscopy and a molecular modeling method under simulative physiological conditions. The quenching mechanism was suggested to be static quenching based on fluorescence and ultraviolet–visible (UV–Vis) spectroscopy. According to the Vant' Hoff equation, the values of enthalpy (?H) and entropy change (?S) were calculated to be ?95.29 kJ/mol and ?218.13 J/mol/K, indicating that the main driving force of the interaction between HSA and METC were hydrogen bonds and van der Waals's forces. By performing displacement measurements, the specific binding of METC in the vicinity of Sudlow's site I of HSA was clarified. An apparent distance of 3.05 nm between Trp214 and METC was obtained via the fluorescence resonance energy transfer (FRET) method. Furthermore, the binding details between METC and HSA were further confirmed by molecular docking studies, which revealed that METC was bound at subdomain IIA through multiple interactions, such as hydrophobic effect, polar forces, hydrogen bonding, etc. The results of three‐dimensional fluorescence and Fourier transform infrared (FTIR) spectroscopy showed that METC caused conformational and some microenvironmental changes in HSA and reduced the α‐helix significantly in the range of 52.3?40.4% in HSA secondary structure. Moreover, the coexistence of metal ions such as Ca²⁺, Al³⁺, Fe³⁺, Zn²⁺, Cu²⁺, Cr³⁺ and Cd²⁺ can decrease the binding constants of METC–HSA. Copyright © 2012 John Wiley & Sons, Ltd.     

A design criterion for symmetric model discrimination based on flexible nominal sets

Experimental design applications for discriminating between models have been hampered by the assumption to know beforehand which model is the true one, which is counter to the very aim of the experiment. Previous approaches to alleviate this requirement were either symmetrizations of asymmetric techniques, or Bayesian, minimax, and sequential approaches. Here we present a genuinely symmetric criterion based on a linearized distance between mean-value surfaces and the newly introduced tool of flexible nominal sets. We demonstrate the computational efficiency of the approach using the proposed criterion and provide a Monte-Carlo evaluation of its discrimination performance on the basis of the likelihood ratio. An application for a pair of competing models in enzyme kinetics is given.     

Large scale maximum average power multiple inference on time-course count data with application to RNA-seq analysis

Experiments that longitudinally collect RNA sequencing (RNA-seq) data can provide transformative insights in biology research by revealing the dynamic patterns of genes. Such experiments create a great demand for new analytic approaches to identify differentially expressed (DE) genes based on large-scale time-course count data. Existing methods, however, are suboptimal with respect to power and may lack theoretical justification. Furthermore, most existing tests are designed to distinguish among conditions based on overall differential patterns across time, though in practice, a variety of composite hypotheses are of more scientific interest. Finally, some current methods may fail to control the false discovery rate. In this paper, we propose a new model and testing procedure to address the above issues simultaneously. Specifically, conditional on a latent Gaussian mixture with evolving means, we model the data by negative binomial distributions. Motivated by Storey (2007) and Hwang and Liu (2010), we introduce a general testing framework based on the proposed model and show that the proposed test enjoys the optimality property of maximum average power. The test allows not only identification of traditional DE genes but also testing of a variety of composite hypotheses of biological interest. We establish the identifiability of the proposed model, implement the proposed method via efficient algorithms, and demonstrate its good performance via simulation studies. The procedure reveals interesting biological insights, when applied to data from an experiment that examines the effect of varying light environments on the fundamental physiology of the marine diatom Phaeodactylum tricornutum.     

Media supplementation for targeted manipulation of monoclonal antibody galactosylation and fucosylation

Monoclonal antibodies are critically important biologics as the largest class of molecules used to treat cancers, rheumatoid arthritis, and other chronic diseases. Antibody glycosylation is a critical quality attribute that has ramifications for patient safety and physiological efficacy—one that can be modified by such factors as media formulation and process conditions during production. Using a design-of-experiments approach, we examined the effect of 2-F-peracetyl fucose (2FP), uridine, and galactose on cell growth and metabolism, titer, and gene expression of key glycosylation-related proteins, and report how the glycoform distribution changed from Days 4 to 7 in a batch process used for IgG1 production from Chinese hamster ovary cells. We observed major glycosylation changes upon supplement addition, where the addition of 2FP decreased antibody fucosylation by up to 48%, galactose addition increased galactosylation by up to 21%, and uridine addition decreased fucosylation and increased galactosylation by 6% and 2%, respectively. Despite having major effects on glycosylation, neither galactose nor 2FP significantly affected cell culture growth, metabolism, or titer. Uridine improved peak cell densities by 23% but also reduced titer by ∼30%. The supplements caused significant changes in gene expression by Day 4 of the cultures where 2FP addition significantly reduced fucosyltransferase 8 and nucleotide sugar transporter gene expression (by ∼2-fold), and uridine addition significantly increased expression of UDP-GlcNAcT (SLC35A3) and B4GALT1–6 genes (by 1.5–3-fold). These gene expression data alongside glycosylation, metabolic, and growth data improve our understanding of the cellular mechanisms affected by media supplementation and suggest approaches for modifying antibody glycosylation in antibody production processes.     

Release from sheep-grazing appears to put some heart back into upland vegetation: A comparison of nutritional properties of plant species in long-term grazing experiments

Rewilding or wilding is a popularised means for enhancing the conservation value of marginal land. In the British uplands, it will involve a reduction, or complete removal, of livestock grazing (sheep), based on the belief that grazing has reduced plant species diversity, the ‘Wet Desert’ hypothesis. The hope is that if livestock is removed, diversity will recover. If true, we hypothesise that the species extirpated/reduced by grazing and then recover on its removal would more nutritious compared to those that persisted. We test this hypothesis at Moor House National Nature Reserve (North-Pennines), where seven sets of paired plots were established between 1953 and 1967 to compare ungrazed/sheep-grazed vegetation. Within these plot-pairs, we compared leaf properties of seven focal species that occurred only, or were present in much greater abundance, in the absence of grazing to those of 10 common species that were common in both grazed and ungrazed vegetation. Each sample was analysed for macro-nutrients, micro-nutrients, digestibility, palatability and decomposability. We ranked the species with respect to 22 variables based on effect size derived from Generalised Linear Modelling (GLM) and compared species using a Principal Components Analysis. We also assessed changes in abundance of the focal species through time using GLMs. Our results support the ‘Wet Desert’ hypothesis, that is, that long-term sheep grazing has selectively removed/reduced species like our focal ones and on recovery, they were more nutritious (macro-nutrients, some micro-nutrients) palatable, digestible and decomposable than common species. Measured changes in abundance of the focal species suggest that their recovery will take 10–20 years in blanket bog and 60 years in high-altitude grasslands. Collectively, these results suggest that sheep grazing has brought about biotic homogenization, and its removal in (re)wilding schemes will reverse this process eventually! The ‘white woolly maggots’ have eaten at least part of the heart out of the highlands/uplands, and it will take some time for recovery.     

REPEATED ASSESSMENT OF THE ENDOGENOUS 24-HOUR PROFILE OF BLOOD PRESSURE UNDER CONSTANT ROUTINE*

The impact of environmental and behavioral factors on the 24-h profile of blood pressure (BP) has been well established. Various attempts have been made to control these exogenous factors, in order to investigate a possible endogenous circadian variation of BP. Recently, we reported the results of the first environmentally and behaviorally controlled laboratory study with 24-h recordings of BP and heart rate (HR) during maintained wakefulness. In this constant-routine study, a pronounced endogenous circadian rhythm of HR was found, but circadian variation of BP was absent. This result suggested that the circadian rhythm of BP observed in earlier controlled studies, with sleep allowed, was evoked by the sleep–wake cycle as opposed to the endogenous circadian pacemaker. In order to verify our previous finding during maintained wakefulness, we repeated the experiment five times with six normotensive, healthy young subjects. Statistical analyses of the hourly measurements of BP and HR confirmed the replicable presence of an endogenous circadian rhythm of HR, as well as the consistent absence of an endogenous circadian variation of BP. Thus, this study provided additional evidence that the 24-h profile of BP—as observed under normal circumstances—is the sole result of environmental and behavioral factors such as the occurrence of sleep, and has no endogenous circadian component. (Chronobiology International, 18(1), 85–98, 2001)