Coherent modeling of longitudinal causal effects on binary outcomes

Analyses of biomedical studies often necessitate modeling longitudinal causal effects. The current focus on personalized medicine and effect heterogeneity makes this task even more challenging. Toward this end, structural nested mean models (SNMMs) are fundamental tools for studying heterogeneous treatment effects in longitudinal studies. However, when outcomes are binary, current methods for estimating multiplicative and additive SNMM parameters suffer from variation dependence between the causal parameters and the noncausal nuisance parameters. This leads to a series of difficulties in interpretation, estimation, and computation. These difficulties have hindered the uptake of SNMMs in biomedical practice, where binary outcomes are very common. We solve the variation dependence problem for the binary multiplicative SNMM via a reparameterization of the noncausal nuisance parameters. Our novel nuisance parameters are variation independent of the causal parameters, and hence allow for coherent modeling of heterogeneous effects from longitudinal studies with binary outcomes. Our parameterization also provides a key building block for flexible doubly robust estimation of the causal parameters. Along the way, we prove that an additive SNMM with binary outcomes does not admit a variation independent parameterization, thereby justifying the restriction to multiplicative SNMMs.     

Performance of a new family of modular,bed-supported,chromatography devices

Prepacked chromatography columns and cassette filtration units offer many advantages in bioprocessing. These include reduced labor costs and processing times, ease of storage, and enhanced process flexibility. Rectangular formats are particularly attractive as they can be easily stacked and multiplexed together for continuous processing. Cylindrical chromatography beds have dominated bioprocessing even though their bed support and pressure-flow performance vary with bed dimensions. This work presents the performance of novel, rhombohedral chromatography devices with internally supported beds. They are compatible with existing chromatography workstations and can be packed with any standard commercial resin. The devices offer pressure-flow characteristics independent of container-volume, simple multiplexing, and separation performance comparable to cylindrical columns. Their bi-planar, internal bed support allows mechanically less-rigid resins to be used at up to four times higher maximal linear velocities, and productivities approaching 200 g/L/h for affinity resins, compared to the 20 g/L/h typical of many column-based devices. Three 5 L devices should allow processing of up to 3 kg of monoclonal antibody per hour.     

Controlled release of adeno-associated virus from alginate hydrogel microbeads with enhanced sensitivity to ultrasound

Adeno-associated virus (AAV)-based gene therapy holds promise as a fundamental treatment for genetic disorders. For clinical applications, it is necessary to control AAV release timing to avoid an immune response to AAV. Here we propose an ultrasound (US)-triggered on-demand AAV release system using alginate hydrogel microbeads (AHMs) with a release enhancer. By using a centrifuge-based microdroplet shooting device, the AHMs encapsulating AAV with tungsten microparticles (W-MPs) are fabricated. Since W-MPs work as release enhancers, the AHMs have high sensitivity to the US with localized variation in acoustic impedance for improving the release of AAV. Furthermore, AHMs were coated with poly-l -lysine (PLL) to adjust the release of AAV. By applying US to the AAV encapsulating AHMs with W-MPs, the AAV was released on demand, and gene transfection to cells by AAV was confirmed without loss of AAV activity. This proposed US-triggered AAV release system expands methodological possibilities in gene therapy.     

Genetic inference of orchid population dynamics on different-aged lava flows in Costa Rica

Colonization is a fundamental ecological process that is important for the persistence of species, particularly when a changing environment necessitates range shifts. Vacant habitats available for colonization often arise from landscape disturbance. Colonization and population expansion processes can be inferred by examining the levels and spatial distribution of genetic variation of plant populations with known disturbance histories. Samples (N = 690) of the terrestrial orchid, Epidendrum radicans, were collected from five lava flow sites on the slopes of Volcán Arenal in Costa Rica that last experienced major eruptions in 1968 and 1992. Individuals were also sampled (N = 188) from four regional populations. Samples were characterized using 15 nuclear genetic markers and analyzed using population genetics statistics. Genetic diversity within sites was moderate (H_e = 0.092–0.192). Contrary to expectation, diversity tended to be lower on the older lava flows (0.131 vs. 0.172) which may reflect their more sheltered topography that restricted pollen/seed immigration, and/or greater intra- and interspecific competition. Genetic diversity measures indicate that the lava flows were colonized by numerous individuals that likely originated from multiple sources while spatial genetic structure (SGS) statistics indicate that most recruitment in the study sites subsequent to colonization resulted from in situ reproduction and localized seed deposition. Younger sites had significantly greater SGS over larger distances which reflects fewer reproductive events, and less spatial and temporal overlap of seed shadows relative to the older sites. Clones were also generally smaller on the younger sites (≤3 m vs. ≤8 m).     

Recovery of the coral Montastrea annularis in the Florida Keys after the 1987 Caribbean “bleaching event”

Many reef-building corals and other cnidarians lost photosynthetic pigments and symbiotic algae (zooxanthellae) during the coral bleaching event in the Caribbean in 1987. The Florida Reef Tract included some of the first documented cases, with widespread bleaching of the massive coral Montastrea annularis beginning in late August. Phototransects at Carysfort Reef showed discoloration of >90% of colonies of this species in March 1988 compared to 0% in July 1986; however no mortality was observed between 1986 and 1988. Samples of corals collected in February and June 1988 had zooxanthellae densities ranging from 0.1 in the most lightly colored corals, to 1.6x10⁶ cells/cm² in the darker corals. Minimum densities increased to 0.5x10⁶ cells/cm² by August 1989. Chlorophyll-a content of zooxanthellae and zooxanthellar mitotic indices were significantly higher in corals with lower densities of zooxanthellae, suggesting that zooxanthellar at low densities may be more nutrientsufficient than those in unbleached corals. Ash-free dry weight of coral tissue was positively correlated with zooxanthellae density at all sample times and was significantly lower in June 1988 compared to August 1989. Proteins and lipids per cm² were significantly higher in August 1989 than in February or June, 1988. Although recovery of zooxanthellae density and coral pigmentation to normal levels may occur in less than one year, regrowth of tissue biomass and energy stores lost during the period of low symbiont densities may take significantly longer.     

Disposition kinetics of ML-1035 sulfoxide enantiomers and the prochiral sulfide in rats

ML-1035, 4-amino-5-chloro-2-[2-(methylsulfinyl)ethoxy]-N-[2-(diethylamino)ethyl]benzamide, is a sulfoxide compound and a racemic gastroprokinetic agent with a chiral center at the sulfur atom. We have investigated the disposition kinetics of (R)-ML-1035 sulfoxide (R) and (S)-ML-1035 sulfoxide (S) after the single enantiomers and the racemic mixture were administered to rats in separate experiments. There was no noticeable chiral inversion after either enantiomer dose. Both enantiomers were rapidly absorbed. After dosing with enantiomers or with the racemate, the resulting plasma concentration-time curve of R was closely parallel to that of S in both intravenous and oral experiments, suggesting that the two enantiomers have approximately the same disposition kinetics. After intravenous enantiomer doses, only S underwent conversion to sulfide, suggesting that sulfidation in the liver is enantioselective. However, the enantioselective sulfidation after intravenous dosing did not introduce a difference in the global plasma disposition profiles between R and S, since the reduction reaction is a minor metabolic process. Other metabolic reactions such as sulfonation and mono-N-desethylations were not enantioselective. After oral administration, conversion to sulfide was observed for both enantioners, implicating the existence of a nonhepatic pathway in sulfidation. Administration of a prochiral sulfide dose was associated with an enantioselective sulfoxidation, in which the R/S concentration ratios increased as a function of time. In addition, enantiomeric interaction causing changes in pharmacokinetic parameters was observed after the oral racemate dose, while the interaction is negligible after an intravenous racemate dose, indicating a route dependency in enantiomeric interaction. © 1993 Wiley-Liss, Inc.     

Role of plant volatiles in host location by the specialist parasitoidMicroplitis croceipes cresson (Braconidae: Hymenoptera)

A study was conducted to determine the primary source of volatile cues within the plant-host complex used by hostseeking freeflying female Microplitis cro-ceipesCresson in flight tunnel bioassays. In single-source and two-choice tests, using wasps given an oviposition experience on either cotton (Gossypium hirsutum)or cowpea (Vigna unguiculata)seedlings damaged by corn earworm (CEW; Helicoverpa zeaBoddie), the damaged seedlings were significantly more attractive than the CEW frass, which was in turn more attractive than the larvae themselves. In a series of two-choice wind-tunnel tests, the discriminatory ability of the wasps was examined, following various oviposition experiences. Significantly more wasps flew to plants with old damage than to plants with fresh damage, regardless of whether they had experience on fresh or old damage. In a comparison of plant species, wasps with only one experience on either hostdamaged cotton or host-damaged cowpea were unable to distinguish between them, and showed no preference for either plant, whereas wasps with multiple experiences on a particular plant preferentially flew to that plant in the choice test. In comparing hosts with nonhosts, wasps successfully learned to distinguish CEW from beet armyworm (BAW; Spodoptera exigua)on cotton but were unable to distinguish CEW from either BAW or cabbage looper (Trichoplusia ni)on cowpea. The results show the important role played by plant volatiles in the location of hosts by M. croceipesand indicate the wasps limitations in discriminating among the various odors. The ecological advantages and disadvantages of this behavior are discussed.     

Piping plover chick ecology following landscape-level disturbance

Population declines of disturbance-dependent species due to suppression of natural disturbances are realized across ecosystems. The piping plover (Charadrius melodus; plover), a disturbance-dependent and conservation-reliant shorebird that nests on sandy beaches and barrier islands on the Atlantic Coast, was listed under the United States Endangered Species Act in 1986. In 2012, Hurricane Sandy landed on Fire and Westhampton islands, barrier island nesting sites for plovers in New York, USA. Hurricane Sandy was a natural disturbance in this system, creating abundant nesting habitat. The number of chicks produced by a pair, or a population, is a direct measure of reproductive output, and gaining a better understanding of productivity and chick behavior following large-scale habitat creation may improve plover habitat management and potentially species persistence. We evaluated the effects of landscape features on habitat selection, behavior, and survival of plover broods using logistic regression, generalized linear mixed effects models, and survival models. Plover broods selected flatter sites with less dense vegetation than available at random. Chick foraging rates were highest in moist substrates and were lower in areas of higher nesting plover density. Chick survival was greater for broods that hatched earlier in the breeding season and increased as chicks aged. Generally, providing access to sites with flatter, moist substrates will likely result in higher quality brood rearing habitat on the landscape. Ultimately, vegetation removal and habitat management may be needed to reduce plover nesting density and ensure sufficient habitat, which may in turn improve plover chick survival. Moreover, sustaining natural landscape disturbances such as those resulting from storms, and not taking actions to prevent hurricane-created overwash, will allow these landscape features to persist.     

Metropolis Monte Carlo calculations of DNA structure using internal coordinates and NMR distance restraints: An alternative method for generating a high-resolution solution structure

Summary A new method, a restrained Monte Carlo (rMC) calculation, is demonstrated for generating high-resolution structures of DNA oligonucleotides in solution from interproton distance restraints and bounds derived from complete relaxation matrix analysis of two-dimensional nuclear Overhauser effect (NOE) spectral peak intensities. As in the case of restrained molecular dynamics (rMD) refinement of structures, the experimental distance restraints and bounds are incorporated as a pseudo-energy term (or penalty function) into the mathematical expression for the molecular energy. However, the use of generalized helical parameters, rather than Cartesian coordinates, to define DNA conformation increases efficiency by decreasing by an order of magnitude the number of parameters needed to describe a conformation and by simplifying the potential energy profile. The Metropolis Monte Carlo method is employed to simulate an annealing process. The rMC method was applied to experimental 2D NOE data from the octamer duplex d(GTA-TAATG)·d(CATTATAC). Using starting structures from different locations in conformational space (e.g. A-DNA and B-DNA), the rMC calculations readily converged, with a root-mean-square deviation (RMSD) of <0.3 Å between structures generated using different protocols and starting structures. Theoretical 2D NOE peak intensities were calculated for the rMC-generated structures using the complete relaxation matrix program CORMA, enabling a comparison with experimental intensities via residual indices. Simulation of the vicinal proton coupling constants was carried out for the structures generated, enabling a comparison with the experimental deoxyribose ring coupling constants, which were not utilized in the structure determination in the case of the rMC simulations. Agreement with experimental 2D NOE and scalar coupling data was good in all cases. The rMC structures are quite similar to that refined by a traditional restrained MD approach (RMSD<0.5 Å) despite the different force fields used and despite the fact that MD refinement was conducted with additional restraints imposed on the endocyclic torsion angles of deoxyriboses. The computational time required for the rMC and rMD calculations is about the same. A comparison of structural parameters is made and some limitations of both methods are discussed with regard to the average nature of the experimental restraints used in the refinement.Abbreviations MC Monte Carlo - rMC restrained Monte Carlo - MD molecular dynamics - rMD restrained molecular dynamics - DG distance geometry - EM energy minimization - 2D NOE two-dimensional nuclear Overhauser effect - DQF-COSY double-quantum-filtered correlation spectroscopy - RMSD root-mean-square deviation To whom correspondence should be addressed.