Translational PK/PD Modeling of Tumor Growth Inhibition and Target Inhibition to Support Dose Range Selection of the LMP7 Inhibitor M3258 in Relapsed/Refractory Multiple Myeloma.

M3258 is an orally bioavailable, potent, selective, reversible inhibitor of the large multifunctional peptidase 7 (LMP7, ß5i, PSMB8) proteolytic subunit of the immunoproteasome, a component of the cellular protein degradation machinery, highly expressed in malignant hematopoietic cells including multiple myeloma. Here we describe the fit-for-purpose pharmacokinetic (PK)/pharmacodynamic (PD)/efficacy modeling of M3258 based on preclinical data from several species. The inhibition of LMP7 activity (PD) and tumor growth (efficacy) were tested in human multiple myeloma xenografts in mice. PK and efficacy data were correlated yielding a free M3258 concentration of 45 nM for half-maximal tumor growth inhibition (KC50). As M3258 only weakly inhibits LMP7 in mouse cells, both in vitro and in vivo bridging studies were performed in rats, monkeys, and dogs for translational modeling. These data indicated that the PD response in human xenograft models was closely reflected in dog PBMCs. A PK/PD model was established, predicting a free IC50 value of 9 nM for M3258 in dogs in vivo, in close agreement with in vitro measurements. In parallel, the human PK parameters of M3258 were predicted by various approaches including in vitro extrapolation and allometric scaling. Using PK/PD/efficacy simulations, the efficacious dose range and corresponding PD response in human were predicted. Taken together, these efforts supported the design of a phase Ia study of M3258 in multiple myeloma patients (NCT04075721). At the lowest tested dose level, the predicted exposure matched well with the observed exposure while the duration of LMP7 inhibition was underpredicted by the model. SIGNIFICANCE STATEMENT: M3258 is a novel inhibitor of the immunoproteasome subunit LMP7. The human PK and human efficacious dose range of M3258 were predicted using in vitro-in vivo extrapolation and allometric scaling methods together with a fit-for-purpose PK/PD and efficacy model based on data from several species. A comparison with data from the Phase Ia clinical study showed that the human PK was accurately predicted, while the extent and duration of PD response were more pronounced than estimated.

Accuracy in Near-Perfect Virus Phylogenies.

Phylogenetic trees from real-world data often include short edges with very few substitutions per site, which can lead to partially resolved trees and poor accuracy. Theory indicates that the number of sites needed to accurately reconstruct a fully resolved tree grows at a rate proportional to the inverse square of the length of the shortest edge. However, when inferred trees are partially resolved due to short edges, "accuracy" should be defined as the rate of discovering false splits (clades on a rooted tree) relative to the actual number found. Thus, accuracy can be high even if short edges are common. Specifically, in a "near-perfect" parameter space in which trees are large, the tree length $\xi$ (the sum of all edge lengths) is small, and rate variation is minimal, the expected false positive rate is less than $\xi/3$; the exact value depends on tree shape and sequence length. This expected false positive rate is far below the false negative rate for small $\xi$ and often well below 5% even when some assumptions are relaxed. We show this result analytically for maximum parsimony and explore its extension to maximum likelihood using theory and simulations. For hypothesis testing, we show that measures of split "support" that rely on bootstrap resampling consistently imply weaker support than that implied by the false positive rates in near-perfect trees. The near-perfect parameter space closely fits several empirical studies of human virus diversification during outbreaks and epidemics, including Ebolavirus, Zika virus, and SARS-CoV-2, reflecting low substitution rates relative to high transmission/sampling rates in these viruses.[Ebolavirus; epidemic; HIV; homoplasy; mumps virus; perfect phylogeny; SARS-CoV-2; virus; West Nile virus; Yule-Harding model; Zika virus.].

Origin and Diversification of the Saguaro Cactus (Carnegiea gigantea): A Within-Species Phylogenomic Analysis.

Reconstructing accurate historical relationships within a species poses numerous challenges, not least in many plant groups in which gene flow is high enough to extend well beyond species boundaries. Nonetheless, the extent of tree-like history within a species is an empirical question on which it is now possible to bring large amounts of genome sequence to bear. We assess phylogenetic structure across the geographic range of the saguaro cactus, an emblematic member of Cactaceae, a clade known for extensive hybridization and porous species boundaries. Using 200 Gb of whole genome resequencing data from 20 individuals sampled from 10 localities, we assembled two data sets comprising 150,000 biallelic single nucleotide polymorphisms (SNPs) from protein coding sequences. From these, we inferred within-species trees and evaluated their significance and robustness using five qualitatively different inference methods. Despite the low sequence diversity, large census population sizes, and presence of wide-ranging pollen and seed dispersal agents, phylogenetic trees were well resolved and highly consistent across both data sets and all methods. We inferred that the most likely root, based on marginal likelihood comparisons, is to the east and south of the region of highest genetic diversity, which lies along the coast of the Gulf of California in Sonora, Mexico. Together with striking decreases in marginal likelihood found to the north, this supports hypotheses that saguaro's current range reflects postglacial expansion from the refugia in the south of its range. We conclude with observations about practical and theoretical issues raised by phylogenomic data sets within species, in which SNP-based methods must be used rather than gene tree methods that are widely used when sequence divergence is higher. These include computational scalability, inference of gene flow, and proper assessment of statistical support in the presence of linkage effects. [Phylogenomics; phylogeography; rooting; Sonoran Desert.].

Prevalence of SARS-CoV-2 Antibodies in New York City Adults, June-October 2020: A Population-Based Survey.

BACKGROUND: Serosurveys help to ascertain burden of infection. Prior severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) serosurveys in New York City (NYC) used nonrandom samples. During June-October 2020, the NYC Health Department conducted a population-based survey estimating SARS-CoV-2 antibody prevalence in NYC adults. METHODS: Participants were recruited from the NYC 2020 Community Health Survey. We estimated citywide and stratified antibody prevalence using a hybrid design: serum tested with the DiaSorin LIAISON SARS-CoV-2 S1/S2 IgG assay and self-reported antibody test results were used together. We estimated univariate frequencies and 95% confidence intervals (CI), accounting for complex survey design. Two-sided P values ≤ .05 were statistically significant. RESULTS: There were 1074 respondents; 497 provided blood and 577 provided only a self-reported antibody test result. Weighted prevalence was 24.3% overall (95% CI, 20.7%-28.3%). Latino (30.7%; 95% CI, 24.1%-38.2%; P < .01) and black (30.7%; 95% CI, 21.9%-41.2%; P = .02) respondents had a higher weighted prevalence compared with white respondents (17.4%; 95% CI, 12.5%-23.7%). CONCLUSIONS: By October 2020, nearly 1 in 3 black and 1 in 3 Latino NYC adults had SARS-CoV-2 antibodies, highlighting unequal impacts of the coronavirus disease 2019 (COVID-19) pandemic on black and Latino NYC adults.

Population-Based Estimates of Coronavirus Disease 2019 (COVID-19)-like Illness, COVID-19 Illness, and Rates of Case Ascertainment, Hospitalizations, and Deaths-Noninstitutionalized New York City Residents, March-April 2020.

Using a population-based, representative telephone survey, ~930 000 New York City residents had COVID-19 illness beginning 20 March-30 April 2020, a period with limited testing. For every 1000 persons estimated with COVID-19 illness, 141.8 were tested and reported as cases, 36.8 were hospitalized, and 12.8 died, varying by demographic characteristics.

The Implications of Lineage-Specific Rates for Divergence Time Estimation.

Rate variation adds considerable complexity to divergence time estimation in molecular phylogenies. Here, we evaluate the impact of lineage-specific rates-which we define as among-branch-rate-variation that acts consistently across the entire genome. We compare its impact to residual rates-defined as among-branch-rate-variation that shows a different pattern of rate variation at each sampled locus, and gene-specific rates-defined as variation in the average rate across all branches at each sampled locus. We show that lineage-specific rates lead to erroneous divergence time estimates, regardless of how many loci are sampled. Further, we show that stronger lineage-specific rates lead to increasing error. This contrasts to residual rates and gene-specific rates, where sampling more loci significantly reduces error. If divergence times are inferred in a Bayesian framework, we highlight that error caused by lineage-specific rates significantly reduces the probability that the 95% highest posterior density includes the correct value, and leads to sensitivity to the prior. Use of a more complex rate prior-which has recently been proposed to model rate variation more accurately-does not affect these conclusions. Finally, we show that the scale of lineage-specific rates used in our simulation experiments is comparable to that of an empirical data set for the angiosperm genus Ipomoea. Taken together, our findings demonstrate that lineage-specific rates cause error in divergence time estimates, and that this error is not overcome by analyzing genomic scale multilocus data sets. [Divergence time estimation; error; rate variation.].

Adapting Survey Data Collection to Respond to the COVID-19 Pandemic: Experiences From a Local Health Department.

The New York City (NYC) Department of Health and Mental Hygiene ("Health Department") conducts routine surveys to describe the health of NYC residents. During the COVID-19 pandemic, the Health Department adjusted existing surveys and developed new ones to improve our understanding of the impact of the pandemic on physical health, mental health, and social determinants of health and to incorporate more explicit measures of racial inequities. The longstanding Community Health Survey was adapted in 2020 to ask questions about COVID-19 and recruit respondents for a population-based severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) serosurvey. A new survey panel, Healthy NYC, was launched in June 2020 and is being used to collect data on COVID-19, mental health, and social determinants of health. In addition, 7 Health Opinion Polls were conducted from March 2020 through March 2021 to learn about COVID-19-related knowledge, attitudes, and opinions, including vaccine intentions. We describe the contributions that survey data have made to the emergency response in NYC in ways that address COVID-19 and the profound inequities of the pandemic. (Am J Public Health. 2021;111(12):2176-2185. https://doi.org/10.2105/AJPH.2021.306515).

Evolution of l-DOPA 4,5-dioxygenase activity allows for recurrent specialisation to betalain pigmentation in Caryophyllales.

The evolution of l-DOPA 4,5-dioxygenase activity, encoded by the gene DODA, was a key step in the origin of betalain biosynthesis in Caryophyllales. We previously proposed that l-DOPA 4,5-dioxygenase activity evolved via a single Caryophyllales-specific neofunctionalisation event within the DODA gene lineage. However, this neofunctionalisation event has not been confirmed and the DODA gene lineage exhibits numerous gene duplication events, whose evolutionary significance is unclear. To address this, we functionally characterised 23 distinct DODA proteins for l-DOPA 4,5-dioxygenase activity, from four betalain-pigmented and five anthocyanin-pigmented species, representing key evolutionary transitions across Caryophyllales. By mapping these functional data to an updated DODA phylogeny, we then explored the evolution of l-DOPA 4,5-dioxygenase activity. We find that low l-DOPA 4,5-dioxygenase activity is distributed across the DODA gene lineage. In this context, repeated gene duplication events within the DODA gene lineage give rise to polyphyletic occurrences of elevated l-DOPA 4,5-dioxygenase activity, accompanied by convergent shifts in key functional residues and distinct genomic patterns of micro-synteny. In the context of an updated organismal phylogeny and newly inferred pigment reconstructions, we argue that repeated convergent acquisition of elevated l-DOPA 4,5-dioxygenase activity is consistent with recurrent specialisation to betalain synthesis in Caryophyllales.

A comparison of the four healthy days measures (HRQOL-4) with a single measure of self-rated general health in a population-based health survey in New York City.

BACKGROUND: Data on health-related quality of life (HRQOL) can be used to track health disparities, assess the impact of chronic diseases, and predict mortality. The Centers for Disease Control and Prevention's "Healthy Days Measures" (HRQOL-4) assesses four key domains: self-rated general health, physical health, mental health, and activity limitations. The domains are not easily combined to summarize overall HRQOL, and some evidence suggests that self-rated general health may be an adequate proxy indicator for overall HRQOL. This study compares self-rated general health as a solitary measure of HRQOL with two summary indices of the HRQOL-4 as a predictor of adverse health conditions in a representative sample of adult New York City residents. METHODS: The 2017 NYC Social Determinants of Health survey implemented by the New York City Department of Health and Mental Hygiene collected data from a representative sample of New Yorkers (n = 2335) via phone, mail, and web. We compared the information criteria and predictive power of self-rated general health with two alternative summary indices of the HRQOL-4 in predicting self-reported health conditions (hypertension, diabetes, obesity, non-specific psychological distress, and a summary indicator for at least one those four morbidities). RESULTS: Overall, 19.1% (95% CI: 16.9, 21.5) of respondents reported that they had fair or poor general health. Self-rated general health was significantly associated with days of poor physical health, poor mental health, and activity limitations (p < 0.001 for each). While the Akaike and Bayesian information criteria suggested that the summary indices of the HRQOL-4 produced marginally better models for predicting adverse health conditions, self-rated general health had slightly higher predictive power than did the summary indices in all models of physical health outcomes as measured by Tjur's pseudo-R2 and the area under the curve. CONCLUSION: We found very small differences between self-rated general health and the summary indices of the HRQOL-4 in predicting health conditions, suggesting self-rated general health is an appropriate proxy measure of overall HRQOL. Because it can be measured with a single question rather than four, it might be the most simple, efficient, and cost-effective method of summarizing HRQOL in large population-based surveys.

On the dynamical structure of calcium oscillations.

Oscillations in the concentration of free cytosolic Ca2+ are an important and ubiquitous control mechanism in many cell types. It is thus correspondingly important to understand the mechanisms that underlie the control of these oscillations and how their period is determined. We show that Class I Ca2+ oscillations (i.e., oscillations that can occur at a constant concentration of inositol trisphosphate) have a common dynamical structure, irrespective of the oscillation period. This commonality allows the construction of a simple canonical model that incorporates this underlying dynamical behavior. Predictions from the model are tested, and confirmed, in three different cell types, with oscillation periods ranging over an order of magnitude. The model also predicts that Ca2+ oscillation period can be controlled by modulation of the rate of activation by Ca2+ of the inositol trisphosphate receptor. Preliminary experimental evidence consistent with this hypothesis is presented. Our canonical model has a structure similar to, but not identical to, the classic FitzHugh-Nagumo model. The characterization of variables by speed of evolution, as either fast or slow variables, changes over the course of a typical oscillation, leading to a model without globally defined fast and slow variables.

Extensive gene tree discordance and hemiplasy shaped the genomes of North American columnar cacti.

Few clades of plants have proven as difficult to classify as cacti. One explanation may be an unusually high level of convergent and parallel evolution (homoplasy). To evaluate support for this phylogenetic hypothesis at the molecular level, we sequenced the genomes of four cacti in the especially problematic tribe Pachycereeae, which contains most of the large columnar cacti of Mexico and adjacent areas, including the iconic saguaro cactus (Carnegiea gigantea) of the Sonoran Desert. We assembled a high-coverage draft genome for saguaro and lower coverage genomes for three other genera of tribe Pachycereeae (Pachycereus, Lophocereus, and Stenocereus) and a more distant outgroup cactus, Pereskia We used these to construct 4,436 orthologous gene alignments. Species tree inference consistently returned the same phylogeny, but gene tree discordance was high: 37% of gene trees having at least 90% bootstrap support conflicted with the species tree. Evidently, discordance is a product of long generation times and moderately large effective population sizes, leading to extensive incomplete lineage sorting (ILS). In the best supported gene trees, 58% of apparent homoplasy at amino sites in the species tree is due to gene tree-species tree discordance rather than parallel substitutions in the gene trees themselves, a phenomenon termed "hemiplasy." The high rate of genomic hemiplasy may contribute to apparent parallelisms in phenotypic traits, which could confound understanding of species relationships and character evolution in cacti.

Associations Between Season, Processing Plant, and Hide Cleanliness Scores with Prevalence and Concentration of Major Shiga Toxin-Producing Escherichia coli on Beef Cattle Hides.

The objectives of this study were (1) to estimate the prevalence and concentration of the seven major Shiga toxin-producing Escherichia coli (STEC) serogroups (O26, O45, O103, O111, O121, O145, and O157), collectively called STEC-7, on cattle hides collected in different seasons and beef processing plants; and (2) to determine associations of season, plant, and hide cleanliness scores with the prevalence and concentration of STEC-7. A total of 720 hide surface samples (240/season) were collected over three seasons (summer and fall 2015 and spring 2016) from beef cattle carcasses in four commercial processing plants in the United States. Samples were subjected to selective culture and spiral plating methods. Overall model-adjusted mean prevalence (95% confidence interval) was 0.3% (0.03-2.3%) for STEC O26; 0.05% (<0.01-8.5%) for STEC O45; 0.2% (0.02-1.9%) for STEC O103; 0.05% (<0.01-8.5%) for STEC O145; and 3.1% (0.6-15.2%) for STEC O157. Four percent of hide samples were enumerable for STEC O157; mean concentration (standard deviation) = 2.1 (0.7) log10 colony-forming units (CFUs)/100 cm2. No samples were enumerable for non-O157 STEC. Hide-on prevalence of STEC O157 and STEC non-O157 (specifically of STEC O103) was higher in summer and spring, respectively. Across seasons and plants, the most common STEC non-O157 serogroups in this study (O26 and O103) were associated with a higher prevalence of STEC O157. Season and plant played a role in prevalence and concentration of STEC in beef cattle hides, varying by serogroup. Tailoring mitigation strategies at the plant can be challenging and processors would benefit from supplementary preharvest interventions to reduce overall contamination pressure at the plant, especially in fall and spring months when hide-on prevalence of STEC non-O157 is higher.

Intercellular Ca(2+) waves: mechanisms and function.

Intercellular calcium (Ca(2+)) waves (ICWs) represent the propagation of increases in intracellular Ca(2+) through a syncytium of cells and appear to be a fundamental mechanism for coordinating multicellular responses. ICWs occur in a wide diversity of cells and have been extensively studied in vitro. More recent studies focus on ICWs in vivo. ICWs are triggered by a variety of stimuli and involve the release of Ca(2+) from internal stores. The propagation of ICWs predominately involves cell communication with internal messengers moving via gap junctions or extracellular messengers mediating paracrine signaling. ICWs appear to be important in both normal physiology as well as pathophysiological processes in a variety of organs and tissues including brain, liver, retina, cochlea, and vascular tissue. We review here the mechanisms of initiation and propagation of ICWs, the key intra- and extracellular messengers (inositol 1,4,5-trisphosphate and ATP) mediating ICWs, and the proposed physiological functions of ICWs.

Quantification of Bacteria Indicative of Fecal and Environmental Contamination from Hides to Carcasses.

Fecal bacteria, which reside in the gastrointestinal tract of cattle, can contaminate beef carcasses during processing. In beef cattle slaughter plants, the presence and concentrations of generic Escherichia coli, coliforms, Enterobacteriaceae (EB), and total aerobic bacteria are monitored as indicator organisms of fecal and environmental contamination. The objectives of this study were as follows: (1) to determine the concentrations of generic E. coli, coliforms, EB, and aerobic bacteria on beef carcasses at different processing points in Midwestern commercial beef slaughter plants during the summer, spring, and fall seasons; and (2) to estimate bacterial transfer on carcasses during the hide removal and evisceration processes. Hide and carcass surface sample swabs were collected from slaughtered cattle at four large commercial processing plants. At each plant visit (3 visits to each of the 4 plants) and during 3 seasons, 20 samples were collected at 5 points: hide-on (hide of animal near exsanguination pit), hide-off carcass, pre-evisceration carcass, postevisceration carcass, and postintervention carcass, for a total of 3600 samples. Bacterial concentrations were determined using 3M™ Petrifilm™ plates. Associations between season and processing plant with concentrations of E. coli, coliforms, EB, and total aerobic bacteria, overall, between hide-on and hide-off, and between pre- and post-evisceration, were evaluated using multilevel mixed-effects linear regression models. Bacterial concentrations on beef carcasses significantly decreased throughout processing. Moreover, hide removal was an important source of carcass contamination, given bacterial concentrations detected on hide-off carcass samples were the highest, and bearing in mind that carcass muscle surfaces should be sterile. Results from this study indicate that the interventions applied by the processing plants were effective, as they probably contributed to the significant reduction of bacterial concentrations of carcasses.

Orosomucoid-like 3 (ORMDL3) upregulates airway smooth muscle proliferation, contraction, and Ca2+ oscillations in asthma.

BACKGROUND: Airway hyperresponsiveness is a major feature of asthma attributed predominantly to an extrinsic immune/inflammatory response increasing airway smooth muscle (ASM) contractility. OBJECTIVE: We investigated whether increased ASM expression of orosomucoid-like 3 (ORMDL3), a gene on chromosome 17q21 highly linked to asthma, induced increased ASM proliferation and contractility in vitro and influenced airway contractility and calcium flux in ASM in precision-cut lung slices (PCLSs) from wild-type and hORMDL3Zp3-Cre mice (which express increased levels of human ORMDL3 [hORMDL3]). METHODS: Levels of ASM proliferation and contraction were assessed in ASM cells transfected with ORMDL3 in vitro. In addition, airway contractility and calcium oscillations were quantitated in ASM cells in PCLSs derived from naive wild-type and naive hORMDL3Zp3-Cre mice, which do not have a blood supply. RESULTS: Increased ASM expression of ORMDL3 in vitro resulted in increased ASM proliferation and contractility. PCLSs derived from naive hORMDL3Zp3-Cre mice, which do not have airway inflammation, exhibit increased airway contractility with increased calcium oscillations in ASM cells. Increased ASM ORMDL3 expression increases levels of ASM sarcoplasmic reticulum Ca2+ ATPase 2b (SERCA2b), which increases ASM proliferation and contractility. CONCLUSION: Overall, these studies provide evidence that an intrinsic increase in ORMDL3 expression in ASM can induce increased ASM proliferation and contractility, which might contribute to increased airway hyperresponsiveness in the absence of airway inflammation in asthmatic patients.

The prevalence of terraced treescapes in analyses of phylogenetic data sets.

BACKGROUND: The pattern of data availability in a phylogenetic data set may lead to the formation of terraces, collections of equally optimal trees. Terraces can arise in tree space if trees are scored with parsimony or with partitioned, edge-unlinked maximum likelihood. Theory predicts that terraces can be large, but their prevalence in contemporary data sets has never been surveyed. We selected 26 data sets and phylogenetic trees reported in recent literature and investigated the terraces to which the trees would belong, under a common set of inference assumptions. We examined terrace size as a function of the sampling properties of the data sets, including taxon coverage density (the proportion of taxon-by-gene positions with any data present) and a measure of gene sampling "sufficiency". We evaluated each data set in relation to the theoretical minimum gene sampling depth needed to reduce terrace size to a single tree, and explored the impact of the terraces found in replicate trees in bootstrap methods. RESULTS: Terraces were identified in nearly all data sets with taxon coverage densities < 0.90. They were not found, however, in high-coverage-density (i.e., ≥ 0.94) transcriptomic and genomic data sets. The terraces could be very large, and size varied inversely with taxon coverage density and with gene sampling sufficiency. Few data sets achieved a theoretical minimum gene sampling depth needed to reduce terrace size to a single tree. Terraces found during bootstrap resampling reduced overall support. CONCLUSIONS: If certain inference assumptions apply, trees estimated from empirical data sets often belong to large terraces of equally optimal trees. Terrace size correlates to data set sampling properties. Data sets seldom include enough genes to reduce terrace size to one tree. When bootstrap replicate trees lie on a terrace, statistical support for phylogenetic hypotheses may be reduced. Although some of the published analyses surveyed were conducted with edge-linked inference models (which do not induce terraces), unlinked models have been used and advocated. The present study describes the potential impact of that inference assumption on phylogenetic inference in the context of the kinds of multigene data sets now widely assembled for large-scale tree construction.

The flagellar protein Enkurin is required for mouse sperm motility and for transport through the female reproductive tract.

Enkurin was identified initially in mouse sperm where it was suggested to act as an intracellular adaptor protein linking membrane calcium influx to intracellular signaling pathways. In order to examine the function of this protein, a targeted mutation was introduced into the mouse Enkurin gene. Males that were homozygous for this mutated allele were subfertile. This was associated with lower rates of sperm transport in the female reproductive tract, including reduced entry into the oviduct and slower migration to the site of fertilization in the distal oviduct, and with poor progressive motility in vitro. Flagella from wild-type animals exhibited symmetrical bending and progressive motility in culture medium, and demembranated flagella exhibited the "curlicue" response to Ca2+ in vitro. In contrast, flagella of mice homozygous for the mutated allele displayed only asymmetric bending, nonprogressive motility, and a loss of Ca2+-responsiveness following demembrantion. We propose that Enkurin is part of a flagellar Ca2+-sensor that regulates bending and that the motility defects following mutation of the locus are the proximate cause of subfertility.

Toward a Self-Updating Platform for Estimating Rates of Speciation and Migration, Ages, and Relationships of Taxa.

Rapidly growing biological data-including molecular sequences and fossils-hold an unprecedented potential to reveal how evolutionary processes generate and maintain biodiversity. However, researchers often have to develop their own idiosyncratic workflows to integrate and analyze these data for reconstructing time-calibrated phylogenies. In addition, divergence times estimated under different methods and assumptions, and based on data of various quality and reliability, should not be combined without proper correction. Here we introduce a modular framework termed SUPERSMART (Self-Updating Platform for Estimating Rates of Speciation and Migration, Ages, and Relationships of Taxa), and provide a proof of concept for dealing with the moving targets of evolutionary and biogeographical research. This framework assembles comprehensive data sets of molecular and fossil data for any taxa and infers dated phylogenies using robust species tree methods, also allowing for the inclusion of genomic data produced through next-generation sequencing techniques. We exemplify the application of our method by presenting phylogenetic and dating analyses for the mammal order Primates and for the plant family Arecaceae (palms). We believe that this framework will provide a valuable tool for a wide range of hypothesis-driven research questions in systematics, biogeography, and evolution. SUPERSMART will also accelerate the inference of a "Dated Tree of Life" where all node ages are directly comparable. [Bayesian phylogenetics; data mining; divide-and-conquer methods; GenBank; multilocus multispecies coalescent; next-generation sequencing; palms; primates; tree calibration.].

Store-operated calcium entry is required for sustained contraction and Ca2+ oscillations of airway smooth muscle.

KEY POINTS: Airway hyper-responsiveness in asthma is driven by excessive contraction of airway smooth muscle cells (ASMCs). Agonist-induced Ca2+ oscillations underlie this contraction of ASMCs and the magnitude of this contraction is proportional to the Ca2+ oscillation frequency. Sustained contraction and Ca2+ oscillations require an influx of extracellular Ca2+ , although the mechanisms and pathways mediating this Ca2+ influx during agonist-induced ASMC contraction are not well defined. By inhibiting store-operated calcium entry (SOCE) or voltage-gated Ca2+ channels (VGCCs), we show that SOCE, rather than Ca2+ influx via VGCCs, provides the major Ca2+ entry pathway into ASMCs to sustain ASMCs contraction and Ca2+ oscillations. SOCE may therefore serve as a potential target for new bronchodilators to reduce airway hyper-responsiveness in asthma. ABSTRACT: Asthma is characterized by airway hyper-responsiveness: the excessive contraction of airway smooth muscle. The extent of this airway contraction is proportional to the frequency of Ca2+ oscillations within airway smooth muscle cells (ASMCs). Sustained Ca2+ oscillations require a Ca2+ influx to replenish Ca2+ losses across the plasma membrane. Our previous studies implied store-operated calcium entry (SOCE) as the major pathway for this Ca2+ influx. In the present study, we explore this hypothesis, by examining the effects of SOCE inhibitors (GSK7975A and GSK5498A) as well as L-type voltage-gated Ca2+ channel inhibitors (nifedipine and nimodipine) on airway contraction and Ca2+ oscillations and SOCE-mediated Ca2+ influx in ASMCs within mouse precision-cut lung slices. We found that both GSK7975A and GSK5498A were able to fully relax methacholine-induced airway contraction by abolishing the Ca2+ oscillations, in a manner similar to that observed in zero extracellular Ca2+ ([Ca2+ ]e ). In addition, GSK7975A and GSK5498A inhibited increases in intracellular Ca2+ ([Ca2+ ]i ) in ASMCs with depleted Ca2+ -stores in response to increased [Ca2+ ]e , demonstrating a response consistent with the inhibition of SOCE. However, GSK7975A and GSK5498A did not reduce Ca2+ release via IP3 receptors stimulated with IP3 released from caged-IP3 . By contrast, nifedipine and nimodipine only partially reduced airway contraction, Ca2+ oscillation frequency and SOCE-mediated Ca2+ influx. These data suggest that SOCE is the major Ca2+ influx pathway for ASMCs with respect to sustaining agonist-induced airway contraction and the underlying Ca2+ oscillations. The mechanisms of SOCE may therefore form novel targets for new bronchodilators.

The relative contributions of store-operated and voltage-gated Ca2+ channels to the control of Ca2+ oscillations in airway smooth muscle.

KEY POINTS: Agonist-dependent oscillations in the concentration of free cytosolic calcium are a vital mechanism for the control of airway smooth muscle contraction and thus are a critical factor in airway hyper-responsiveness. Using a mathematical model, closely tied to experimental work, we show that the oscillations in membrane potential accompanying the calcium oscillations have no significant effect on the properties of the calcium oscillations. In addition, the model shows that calcium entry through store-operated calcium channels is critical for calcium oscillations, but calcium entry through voltage-gated channels has much less effect. The model predicts that voltage-gated channels are less important than store-operated channels in the control of airway smooth muscle tone. ABSTRACT: Airway smooth muscle contraction is typically the key mechanism underlying airway hyper-responsiveness, and the strength of muscle contraction is determined by the frequency of oscillations of intracellular calcium (Ca2+ ) concentration. In airway smooth muscle cells, these Ca2+ oscillations are caused by cyclic Ca2+ release from the sarcoplasmic reticulum, although Ca2+ influx via plasma membrane channels is also necessary to sustain the oscillations over longer times. To assess the relative contributions of store-operated and voltage-gated Ca2+ channels to this Ca2+ influx, we generated a comprehensive mathematical model, based on experimental Ca2+ measurements in mouse precision-cut lung slices, to simulate Ca2+ oscillations and changes in membrane potential. Agonist-induced Ca2+ oscillations are accompanied by oscillations in membrane potential, although the membrane potential oscillations are too small to generate large Ca2+ currents through voltage-gated Ca2+ channels, and thus have little effect on the Ca2+ oscillations. Ca2+ entry through voltage-gated channels only becomes important when the cell is depolarized (e.g. by a high external K+ concentration). As a result, agonist-induced Ca2+ oscillations are critically dependent on Ca2+ entry through store-operated channels but do not depend strongly on Ca2+ entry though voltage-gated channels.