A simulation approach for power calculation in large cohort studies based on multistate models

Realistic power calculations for large cohort studies and nested case control studies are essential for successfully answering important and complex research questions in epidemiology and clinical medicine. For this, we provide a methodical framework for general realistic power calculations via simulations that we put into practice by means of an R‐based template. We consider staggered recruitment and individual hazard rates, competing risks, interaction effects, and the misclassification of covariates. The study cohort is assembled with respect to given age‐, gender‐, and community distributions. Nested case‐control analyses with a varying number of controls enable comparisons of power with a full cohort analysis. Time‐to‐event generation under competing risks, including delayed study‐entry times, is realized on the basis of a six‐state Markov model. Incidence rates, prevalence of risk factors and prefixed hazard ratios allow for the assignment of age‐dependent transition rates given in the form of Cox models. These provide the basis for a central simulation‐algorithm, which is used for the generation of sample paths of the underlying time‐inhomogeneous Markov processes. With the inclusion of frailty terms into the Cox models the Markov property is specifically biased. An “individual Markov process given frailty” creates some unobserved heterogeneity between individuals. Different left‐truncation‐ and right‐censoring patterns call for the use of Cox models for data analysis. p‐values are recorded over repeated simulation runs to allow for the desired power calculations. For illustration, we consider scenarios with a “testing” character as well as realistic scenarios. This enables the validation of a correct implementation of theoretical concepts and concrete sample size recommendations against an actual epidemiological background, here given with possible substudy designs within the German National Cohort.     

Effect of feeding boxes on the behavior of stereotyping amur tigers (Panthera tigris altaica) in the Zurich Zoo,Zurich, Switzerland

The stereotyped pacing shown by the two Amur tigers in the Zurich Zoo was hypothesized as being caused by permanently frustrated appetitive foraging behavior. Several electrically controlled feeding boxes were installed and access to each box was possible only twice a day for 15 min at semi‐random times. The boxes had to be opened actively by the tigers. Two trials were carried out: one with solitary confinement, and one with paired confinement. During box feeding, the female's stereotyped pacing was significantly reduced from 16% (solitary confinement, conventional feeding) and 7% (paired confinement, conventional feeding) to 1% (solitary confinement) and less than 0.01% (paired confinement) of the daily observed time. The female's sleeping increased significantly in both solitary and paired confinement. The male only showed a significant reduction in stereotyped pacing behavior when kept with the female (conventional feeding: 10%; box feeding: <0.01% of the daily observed time). On days with a box‐feeding regime in paired confinement, the male spent 25% (83 min) of the observed time with active behavior at the feeding boxes. The results support the hypothesis that permanently frustrated appetitive foraging behavior causes stereotyped pacing in adult tigers. Zoo Biol 21:573–584, 2002. © 2002 Wiley‐Liss, Inc.     

Climatic forcing and larval dispersal capabilities shape the replenishment of fishes and their habitat‐forming biota on a tropical coral reef

Fluctuations in marine populations often relate to the supply of recruits by oceanic currents. Variation in these currents is typically driven by large‐scale changes in climate, in particular ENSO (El Nino Southern Oscillation). The dependence on large‐scale climatic changes may, however, be modified by early life history traits of marine taxa. Based on eight years of annual surveys, along 150 km of coastline, we examined how ENSO influenced abundance of juvenile fish, coral spat, and canopy‐forming macroalgae. We then investigated what traits make populations of some fish families more reliant on the ENSO relationship than others. Abundance of juvenile fish and coral recruits was generally positively correlated with the Southern Oscillation Index (SOI), higher densities recorded during La Niña years, when the ENSO‐influenced Leeuwin Current is stronger and sea surface temperature higher. The relationship is typically positive and stronger among fish families with shorter pelagic larval durations and stronger swimming abilities. The relationship is also stronger at sites on the coral back reef, although the strongest of all relationships were among the lethrinids (r = .9), siganids (r = .9), and mullids (r = .8), which recruit to macroalgal meadows in the lagoon. ENSO effects on habitat seem to moderate SOI–juvenile abundance relationship. Macroalgal canopies are higher during La Niña years, providing more favorable habitat for juvenile fish and strengthening the SOI effect on juvenile abundance. Conversely, loss of coral following a La Niña‐related heat wave may have compromised postsettlement survival of coral dependent species, weakening the influence of SOI on their abundance. This assessment of ENSO effects on tropical fish and habitat‐forming biota and how it is mediated by functional ecology improves our ability to predict and manage changes in the replenishment of marine populations.     

Balancing Detection and Eradication for Control of Epidemics: Sudden Oak Death in Mixed-Species Stands

Culling of infected individuals is a widely used measure for the control of several plant and animal pathogens but culling first requires detection of often cryptically-infected hosts. In this paper, we address the problem of how to allocate resources between detection and culling when the budget for disease management is limited. The results are generic but we motivate the problem for the control of a botanical epidemic in a natural ecosystem: sudden oak death in mixed evergreen forests in coastal California, in which species composition is generally dominated by a spreader species (bay laurel) and a second host species (coast live oak) that is an epidemiological dead-end in that it does not transmit infection but which is frequently a target for preservation. Using a combination of an epidemiological model for two host species with a common pathogen together with optimal control theory we address the problem of how to balance the allocation of resources for detection and epidemic control in order to preserve both host species in the ecosystem. Contrary to simple expectations our results show that an intermediate level of detection is optimal. Low levels of detection, characteristic of low effort expended on searching and detection of diseased trees, and high detection levels, exemplified by the deployment of large amounts of resources to identify diseased trees, fail to bring the epidemic under control. Importantly, we show that a slight change in the balance between the resources allocated to detection and those allocated to control may lead to drastic inefficiencies in control strategies. The results hold when quarantine is introduced to reduce the ingress of infected material into the region of interest.     

Human λ light chain locus: Organization and DNA sequences of three genomicJ regions

Evidence for the genomic organization of human lambda light chain joining (J) region gene segments is presented. A mouse J probe was used in Southern hybridizations to localize joining region sequences in a cosmid clone containing the genomic cluster of six human lambda constant (C) region gene segments. The results of these hybridizations suggest the presence of at least one J gene segment upstream from each constant region gene segment. The DNA sequences indicate that the human JI, J2, and J3 gene segments have consensus nonamer and heptamer sequences, proposed to be involved in V-J joining, are capable of encoding the known amino acid sequences for the respective J peptides, and have a sequence which could give functional RNA splice site at the end of their coding regions. Our data show that a single functional J is located 1.3 or 1.6 kb upstream of each of the C gene segments known to encode the Mcg, Kern^– Oz^–, and Kern^–Oz⁺ isotypes. Therefore, the gene organization of this region of the human lambda locus is J1 CI -J2C2-J3C3. The DNA sequences ofJ 1,J 2, andJ 3 presented in this paper establish that a singleJ gene segment precedes each expressed C gene segment, and support a model for the evolution of the human JC clusters where JICI andJ2C2-J3C3. arose from different ancestral JC units.     

Operculum closing as a defence against predatory leeches in four British freshwater prosobranch snails

The defence reaction of operculum closing in response to the presence of the molluscivorous leech Glossiphonia complanata (L.) and the non-molluscivorous Erpobdella octoculata (L.) was studied in four species of freshwater prosobranch gastropod. Bithynia tentaculata (L.) and Valvata piscinalis (Müller) can distinguish between the leeches, reacting only to G. complanata. V. piscinalis is capable of a greater degree of distance chemoreception of the leech ‘scent’. Valvata cristata Müller and Potamopyrgus jenkinsi (Smith) did not react to either leech. V. cristata may not be a potential prey item for G. complanata, while P. jenkinsi is fed on by the leech, but is a relative newcomer to the freshwater fauna. Animal Ecology Research Group, Department of Zoology, University of Oxford Commonwealth Forestry Institute, Department of Plant Sciences, University of Oxford     

Kinetic measurement of bacterial adherence to eukaryotic cells

We developed a kinetic assay using a monolayer of differentiated respiratory epithelium in culture to assess bacterial adherence. Mean residence time of bacteria in the tissue culture chamber was estimated from a model-independent (moment) analysis of the rate of bacterial washout from perfused Rose chambers. Results with this method compared favorably with visual assessment of adherence and double radiolabel method with H. influenzae. Adherence was assessed with low inoculae of H. influenzae, P. cepacia and P. aeruginosa avoiding cytotoxic effects seen when large inoculae are added to eukaryotic cells. This method will provide a means of assessing adherence of pathogenic respiratory bacteria to their cellular target at low inoculae.