A Bayesian Approach for Modeling Cattle Movements in the United States: Scaling up a Partially Observed Network

Networks are rarely completely observed and prediction of unobserved edges is an important problem, especially in disease spread modeling where networks are used to represent the pattern of contacts. We focus on a partially observed cattle movement network in the U.S. and present a method for scaling up to a full network based on Bayesian inference, with the aim of informing epidemic disease spread models in the United States. The observed network is a 10% state stratified sample of Interstate Certificates of Veterinary Inspection that are required for interstate movement; describing approximately 20,000 movements from 47 of the contiguous states, with origins and destinations aggregated at the county level. We address how to scale up the 10% sample and predict unobserved intrastate movements based on observed movement distances. Edge prediction based on a distance kernel is not straightforward because the probability of movement does not always decline monotonically with distance due to underlying industry infrastructure. Hence, we propose a spatially explicit model where the probability of movement depends on distance, number of premises per county and historical imports of animals. Our model performs well in recapturing overall metrics of the observed network at the node level (U.S. counties), including degree centrality and betweenness; and performs better compared to randomized networks. Kernel generated movement networks also recapture observed global network metrics, including network size, transitivity, reciprocity, and assortativity better than randomized networks. In addition, predicted movements are similar to observed when aggregated at the state level (a broader geographic level relevant for policy) and are concentrated around states where key infrastructures, such as feedlots, are common. We conclude that the method generally performs well in predicting both coarse geographical patterns and network structure and is a promising method to generate full networks that incorporate the uncertainty of sampled and unobserved contacts.     

Risk of Cerebral Palsy and Childhood Epilepsy Related to Infections before or during Pregnancy

Background and Aim

Maternal infections during pregnancy have been associated with several neurological disorders in the offspring. However, given the lack of specificity for both the exposures and the outcomes, other factors related to infection such as impaired maternal immune function may be involved in the causal pathway. If impaired maternal immune function plays a role, we would expect infection before pregnancy to be associated with these neurological outcomes.

Methods/Principal Findings

The study population included all first-born singletons in Denmark between January 1 1982 and December 31 2004. We identified women who had hospital-recorded infections within the 5 year period before pregnancy, and women who had hospital-recorded infections during pregnancy. We grouped infections into either infections of the genitourinary system, or any other infections. Cox models were used to estimate adjusted hazard ratios (aHRs) with 95% confidence interval (CI). Maternal infection of the genitourinary system during pregnancy was associated with an increased risk of cerebral palsy (aHR = 1.63, 95% CI: 1.34–1.98) and epilepsy (aHR = 1.27, 95% CI: 1.13–1.42) in the children, compared to children of women without infections during pregnancy. Among women without hospital-recorded infections during pregnancy, maternal infection before pregnancy was associated with an increased risk of epilepsy (aHR = 1.35, 95% CI: 1.21–1.50 for infections of the genitourinary system, and HR = 1.12, 95% CI: 1.03–1.22 for any other infections) and a slightly higher risk of cerebral palsy (aHR = 1.20, 95% CI: 0.96–1.49 for infections of the genitourinary system, and HR = 1.23, 95% CI: 1.06–1.43 for any other infections) in the children, compared to children of women without infections before (and during) pregnancy.

Conclusions

These findings indicate that the maternal immune system, maternal infections, or factors related to maternal immune function play a role in the observed associations between maternal infections before pregnancy and cerebral diseases in the offspring.     

Quantitative Bioluminescent Imaging of Pre-Erythrocytic Malaria Parasite Infection Using Luciferase-Expressing Plasmodium yoelii

The liver stages of Plasmodium parasites are important targets for the development of anti-malarial vaccine candidates and chemoprophylaxis approaches that aim to prevent clinical infection. Analyzing the impact of interventions on liver stages in the murine malaria model system Plasmodium yoelii has been cumbersome and requires terminal procedures. In vivo imaging of bioluminescent parasites has previously been shown to be an effective and non-invasive alternative to monitoring liver stage burden in the Plasmodium berghei model. Here we report the generation and characterization of a transgenic P. yoelii parasite expressing the reporter protein luciferase throughout the parasite life cycle. In vivo bioluminescent imaging of these parasites allows for quantitative analysis of P. yoelii liver stage burden and parasite development, which is comparable to quantitative RT-PCR analysis of liver infection. Using this system, we show that both BALB/cJ and C57BL/6 mice show comparable susceptibility to P. yoelii infection with sporozoites and that bioluminescent imaging can be used to monitor protective efficacy of attenuated parasite immunizations. Thus, this rapid, simple and noninvasive method for monitoring P. yoelii infection in the liver provides an efficient system to screen and evaluate the effects of anti-malarial interventions in vivo and in real-time.     

Expression of Fungal Cutinase and Swollenin in Tobacco Chloroplasts Reveals Novel Enzyme Functions and/or Substrates

In order to produce low-cost biomass hydrolyzing enzymes, transplastomic lines were generated that expressed cutinase or swollenin within chloroplasts. While swollenin expressing plants were homoplasmic, cutinase transplastomic lines remained heteroplasmic. Both transplastomic lines showed interesting modifications in their phenotype, chloroplast structure, and functions. Ultrastructural analysis of chloroplasts from cutinase- and swollenin-expressing plants did not show typical lens shape and granal stacks. But, their thylakoid membranes showed unique scroll like structures and chloroplast envelope displayed protrusions, stretching into the cytoplasm. Unusual honeycomb structures typically observed in etioplasts were observed in mature chloroplasts expressing swollenin. Treatment of cotton fiber with chloroplast-derived swollenin showed enlarged segments and the intertwined inner fibers were irreversibly unwound and fully opened up due to expansin activity of swollenin, causing disruption of hydrogen bonds in cellulose fibers. Cutinase transplastomic plants showed esterase and lipase activity, while swollenin transplastomic lines lacked such enzyme activities. Higher plants contain two major galactolipids, monogalactosyldiacylglycerol (MGDG) and digalactosyldiacylglycerol (DGDG), in their chloroplast thylakoid membranes that play distinct roles in their structural organization. Surprisingly, purified cutinase effectively hydrolyzed DGDG to MGDG, showing alpha galactosidase activity. Such hydrolysis resulted in unstacking of granal thylakoids in chloroplasts and other structural changes. These results demonstrate DGDG as novel substrate and function for cutinase. Both MGDG and DGDG were reduced up to 47.7% and 39.7% in cutinase and 68.5% and 67.5% in swollenin expressing plants. Novel properties and functions of both enzymes reported here for the first time should lead to better understanding and enhanced biomass hydrolysis.     

Paclitaxel-Induced Apoptosis Is BAK-Dependent,but BAX and BIM-Independent in Breast Tumor

Paclitaxel (Taxol)-induced cell death requires the intrinsic cell death pathway, but the specific participants and the precise mechanisms are poorly understood. Previous studies indicate that a BH3-only protein BIM (BCL-2 Interacting Mediator of cell death) plays a role in paclitaxel-induced apoptosis. We show here that BIM is dispensable in apoptosis with paclitaxel treatment using bim^−/− MEFs (mouse embryonic fibroblasts), the bim^−/− mouse breast tumor model, and shRNA-mediated down-regulation of BIM in human breast cancer cells. In contrast, both bak ^−/− MEFs and human breast cancer cells in which BAK was down-regulated by shRNA were more resistant to paclitaxel. However, paclitaxel sensitivity was not affected in bax^−/− MEFs or in human breast cancer cells in which BAX was down-regulated, suggesting that paclitaxel-induced apoptosis is BAK-dependent, but BAX-independent. In human breast cancer cells, paclitaxel treatment resulted in MCL-1 degradation which was prevented by a proteasome inhibitor, MG132. A Cdk inhibitor, roscovitine, blocked paclitaxel-induced MCL-1 degradation and apoptosis, suggesting that Cdk activation at mitotic arrest could induce subsequent MCL-1 degradation in a proteasome-dependent manner. BAK was associated with MCL-1 in untreated cells and became activated in concert with loss of MCL-1 expression and its release from the complex. Our data suggest that BAK is the mediator of paclitaxel-induced apoptosis and could be an alternative target for overcoming paclitaxel resistance.     

Fluctuations in Species-Level Protein Expression Occur during Element and Nutrient Cycling in the Subsurface

While microbial activities in environmental systems play a key role in the utilization and cycling of essential elements and compounds, microbial activity and growth frequently fluctuates in response to environmental stimuli and perturbations. To investigate these fluctuations within a saturated aquifer system, we monitored a carbon-stimulated in situ Geobacter population while iron reduction was occurring, using 16S rRNA abundances and high-resolution tandem mass spectrometry proteome measurements. Following carbon amendment, 16S rRNA analysis of temporally separated samples revealed the rapid enrichment of Geobacter-like environmental strains with strong similarity to G. bemidjiensis. Tandem mass spectrometry proteomics measurements suggest high carbon flux through Geobacter respiratory pathways, and the synthesis of anapleurotic four carbon compounds from acetyl-CoA via pyruvate ferredoxin oxidoreductase activity. Across a 40-day period where Fe(III) reduction was occurring, fluctuations in protein expression reflected changes in anabolic versus catabolic reactions, with increased levels of biosynthesis occurring soon after acetate arrival in the aquifer. In addition, localized shifts in nutrient limitation were inferred based on expression of nitrogenase enzymes and phosphate uptake proteins. These temporal data offer the first example of differing microbial protein expression associated with changing geochemical conditions in a subsurface environment.     

A Dual Read-Out Assay to Evaluate the Potency of Compounds Active against Mycobacterium tuberculosis

Tuberculosis is a serious global health problem caused by the bacterium Mycobacterium tuberculosis. There is an urgent need for discovery and development of new treatments, but this can only be accomplished through rapid and reproducible M. tuberculosis assays designed to identify potent inhibitors. We developed an automated 96-well assay utilizing a recombinant strain of M. tuberculosis expressing a far-red fluorescent reporter to determine the activity of novel compounds; this allowed us to measure growth by monitoring both optical density and fluorescence. We determined that optical density and fluorescence were correlated with cell number during logarithmic phase growth. Fluorescence was stably maintained without antibiotic selection over 5 days, during which time cells remained actively growing. We optimized parameters for the assay, with the final format being 5 days’ growth in 96-well plates in the presence of 2% w/v DMSO. We confirmed reproducibility using rifampicin and other antibiotics. The dual detection method allows for a reproducible calculation of the minimum inhibitory concentration (MIC), at the same time detecting artefacts such as fluorescence quenching or compound precipitation. We used our assay to confirm anti-tubercular activity and establish the structure activity relationship (SAR) around the imidazo[1,2-a]pyridine-3-carboxamides, a promising series of M. tuberculosis inhibitors.     

Moving Domain Computational Fluid Dynamics to Interface with an Embryonic Model of Cardiac Morphogenesis

Peristaltic contraction of the embryonic heart tube produces time- and spatial-varying wall shear stress (WSS) and pressure gradients (∇P) across the atrioventricular (AV) canal. Zebrafish (Danio rerio) are a genetically tractable system to investigate cardiac morphogenesis. The use of Tg(fli1a:EGFP)^y1 transgenic embryos allowed for delineation and two-dimensional reconstruction of the endocardium. This time-varying wall motion was then prescribed in a two-dimensional moving domain computational fluid dynamics (CFD) model, providing new insights into spatial and temporal variations in WSS and ∇P during cardiac development. The CFD simulations were validated with particle image velocimetry (PIV) across the atrioventricular (AV) canal, revealing an increase in both velocities and heart rates, but a decrease in the duration of atrial systole from early to later stages. At 20-30 hours post fertilization (hpf), simulation results revealed bidirectional WSS across the AV canal in the heart tube in response to peristaltic motion of the wall. At 40-50 hpf, the tube structure undergoes cardiac looping, accompanied by a nearly 3-fold increase in WSS magnitude. At 110-120 hpf, distinct AV valve, atrium, ventricle, and bulbus arteriosus form, accompanied by incremental increases in both WSS magnitude and ∇P, but a decrease in bi-directional flow. Laminar flow develops across the AV canal at 20-30 hpf, and persists at 110-120 hpf. Reynolds numbers at the AV canal increase from 0.07±0.03 at 20-30 hpf to 0.23±0.07 at 110-120 hpf (p< 0.05, n=6), whereas Womersley numbers remain relatively unchanged from 0.11 to 0.13. Our moving domain simulations highlights hemodynamic changes in relation to cardiac morphogenesis; thereby, providing a 2-D quantitative approach to complement imaging analysis.     

Differences in morphological growth trajectories between greenside darter Etheostoma blennioides populations from two distinct drainages

The greenside darter Etheostoma blennioides is a small-bodied benthic stream fish that occurs in multiple drainages of the eastern USA. Previous studies have revealed some morphological differentiations among greensides from isolated drainage systems but growth relationships among functional morphological characteristics have not been examined within and between populations. Specifically, we tested for differences in functional morphology and allometric distinctions in morphological growth trajectories between greenside populations from two drainages, the Osage River and White River. Morphological differences between individuals of the two drainages included more dorso-ventrally compressed bodies, longer snouts and larger jaw features in White River darters. Furthermore, body depth, snout length and jaw width grow at significantly different trajectories in individuals of the two drainages. Individuals of the White River drainage (generally a higher gradient stream system) become comparatively more streamlined and develop proportionately larger feeding related traits during growth. These results suggest that individuals of both drainages share a similar morphology during early life but deviate during continued development and this may be partially due to environmentally induced plastic response. This study represents a novel approach for comparing morphological development among fish populations and presents environmental factors (e.g., stream gradient and elevation) that potentially influence functional characteristics.     

Will Lynx Lose Their Edge? Canada Lynx Occupancy in Washington

Populations of species located at southern range edges may be particularly vulnerable to the effects of climate change as warming temperatures and subsequent changes to ecosystems exceed species-specific tolerances. One such species is Canada lynx (Lynx canadensis), a cold-adapted mesocarnivore that maintains a large core population in Alaska, USA, and Canada but exists within several peripheral populations in the contiguous United States. Increases in temperature, declines in snow pack, and climate-influenced increases in fire frequency and intensity, could negatively affect lynx populations, threatening their long-term persistence in the continental United States. Despite these threats, our understanding of broad-scale effects on lynx occupancy and the extent of current lynx distribution in many of these peripheral populations is minimal. We conducted an occupancy survey of lynx in Washington, USA, using a spatially extensive camera-trapping array covering 7,000 km² of potential lynx habitat. We used the resulting database of detection data to develop single-season occupancy models to examine the abiotic and biotic effects on current lynx occupancy and predict future lynx distribution based on climate change forecasts. Our results show lynx occupancy across the Washington landscape is restricted and dictated largely by abiotic factors, disturbance regimes, and distance from source populations in Canada. Predictions of future distribution suggest lynx will be increasingly challenged by climatic changes, particularly at the southern and lower elevation portions of their range in Washington. Our results paint an alarming picture for lynx persistence in Washington that is relevant to current deliberations regarding lynx delisting from the Endangered Species Act. Our simple camera design was a highly effective method for surveying lynx across broad spatial scales, and could be a key monitoring tool for lynx that is easy to implement by researchers and government agencies. © 2020 The Wildlife Society.