As large-scale studies of gene expression with multiple sources of biological and technical variation become widely adopted, characterizing these drivers of variation becomes essential to ...understanding disease biology and regulatory genetics.
We describe a statistical and visualization framework, variancePartition, to prioritize drivers of variation based on a genome-wide summary, and identify genes that deviate from the genome-wide trend. Using a linear mixed model, variancePartition quantifies variation in each expression trait attributable to differences in disease status, sex, cell or tissue type, ancestry, genetic background, experimental stimulus, or technical variables. Analysis of four large-scale transcriptome profiling datasets illustrates that variancePartition recovers striking patterns of biological and technical variation that are reproducible across multiple datasets.
Our open source software, variancePartition, enables rapid interpretation of complex gene expression studies as well as other high-throughput genomics assays. variancePartition is available from Bioconductor: http://bioconductor.org/packages/variancePartition .
•We present various genotype to phenotype models to predict complex phenotypes with G×E as function of genotypic and environmental inputs.•We show how genotype-to-phenotype models can be generalized ...to incorporate additional phenotypic information to improve yield predictions.•We discuss how to evaluate the utility of information from new phenotyping techniques in the context of predictive genotype-to-phenotype models.
New types of phenotyping tools generate large amounts of data on many aspects of plant physiology and morphology with high spatial and temporal resolution. These new phenotyping data are potentially useful to improve understanding and prediction of complex traits, like yield, that are characterized by strong environmental context dependencies, i.e., genotype by environment interactions. For an evaluation of the utility of new phenotyping information, we will look at how this information can be incorporated in different classes of genotype-to-phenotype (G2P) models. G2P models predict phenotypic traits as functions of genotypic and environmental inputs. In the last decade, access to high-density single nucleotide polymorphism markers (SNPs) and sequence information has boosted the development of a class of G2P models called genomic prediction models that predict phenotypes from genome wide marker profiles. The challenge now is to build G2P models that incorporate simultaneously extensive genomic information alongside with new phenotypic information. Beyond the modification of existing G2P models, new G2P paradigms are required. We present candidate G2P models for the integration of genomic and new phenotyping information and illustrate their use in examples. Special attention will be given to the modelling of genotype by environment interactions. The G2P models provide a framework for model based phenotyping and the evaluation of the utility of phenotyping information in the context of breeding programs.
Heatstroke (HS), associated with the early activation of the coagulation system and frequently presenting with thrombocytopenia, poses a significant healthcare challenge. Understanding the ...relationship of nadir platelet count (PLT) within 24 h for adverse outcomes in HS patients is crucial for optimizing management strategies.
This retrospective cohort study, conducted in six tertiary care hospitals, involved patients diagnosed with HS and admitted to the emergency departments. The primary and secondary outcomes included in-hospital mortality and various acute complications, respectively, with logistic regression models utilized for assessing associations between nadir PLT and outcomes. The PLT count change curve was described using a generalized additive mixed model (GAMM), with additional analyses involving body temperature (BT) at 2 h also conducted.
Of the 152 patients included, 19 (12.5%) died in-hospital. The median nadir PLT within 24 h was 99.5 (58.8–145.0)*10^9/L. Notably, as a continuous variable (10*10^9/L), nadir PLT was significantly associated with in-hospital mortality (OR 0.76; 95% CI 0.64–0.91; P = 0.003) and other adverse outcomes like acute kidney and liver injury, even after adjustment for confounders. GAMM revealed a more rapid and significant PLT decline in the non-survival group over 24 h, with differential PLT dynamics also observed based on BT at 2 h.
Nadir PLT within 24 h were tied to in-hospital mortality and various adverse outcomes in HS patients. Early effective cooling measures demonstrated a positive impact on these associations, underscoring their importance in patient management.
The standard two-stage approach for estimating non-linear dose–response curves based on aggregated data typically excludes those studies with less than three exposure groups. We develop the one-stage ...method as a linear mixed model and present the main aspects of the methodology, including model specification, estimation, testing, prediction, goodness-of-fit, model comparison, and quantification of between-studies heterogeneity. Using both fictitious and real data from a published meta-analysis, we illustrated the main features of the proposed methodology and compared it to a traditional two-stage analysis. In a one-stage approach, the pooled curve and estimates of the between-studies heterogeneity are based on the whole set of studies without any exclusion. Thus, even complex curves (splines, spike at zero exposure) defined by several parameters can be estimated. We showed how the one-stage method may facilitate several applications, in particular quantification of heterogeneity over the exposure range, prediction of marginal and conditional curves, and comparison of alternative models. The one-stage method for meta-analysis of non-linear curves is implemented in the dosresmeta R package. It is particularly suited for dose–response meta-analyses of aggregated where the complexity of the research question is better addressed by including all the studies.
•Treatment by visit interactions exist in HAM-D scores and depression status in first month.•Baseline anxiety and major depressive disorder (MDD) associated with HAM-D scores and depression ...status.•Cocaine use disorder associated with HAM-D score and depression status in non-MDD individuals.•Cannabis use disorder associated with HAM-D scores and depression status in MDD individuals.•Amphetamine use disorder associated with HAM-D score just in non-MDD individuals.•Anxiety was associated with HAM-D scores regardless of MDD.
This study aimed to evaluate the longitudinal treatment effect on depression measured by Hamilton Depression Rating Scale (HAM-D) score in a randomized clinical trial for the treatment of opioid use disorder (OUD).
We conducted a secondary data analysis of data from the National Institute on Drug Abuse's Clinical Trials Network Protocol-0051. Patients with OUD (N = 570) were randomized to receive buprenorphine/naloxone (BUP-NX, n = 287) or extended-release naltrexone injection (XR-NTX, n = 283). The HAM-D score was completed at baseline and follow-up visit up to 36 weeks. A linear mixed model analysis was performed for log transformed HAM-D score and a generalized linear mixed model analysis was conducted for depression status.
Compared with BUP-NX, subjects randomized to XR-NTX had higher HAM-D scores at weeks 1 and 3 (p<0.05). There were significant interactions between treatment and visit on HAM-D score and depression status during the first four weeks of treatments in individuals without lifetime major depressive disorder (MDD). Past year cocaine use was associated with HAM-D score and depression status just in individuals without MDD, whereas past year cannabis use was associated with HAM-D score and depression status just in individuals with MDD. Past year amphetamine use was associated with HAM-D score just in individuals without MDD, however, lifetime anxiety was associated with HAM-D scores regardless of MDD.
When prescribing XR-NTX, particularly in the first month of treatment, it is essential to monitor for depressive symptoms. Screening for depression and multiple substance abuse may help clinicians identify appropriate treatment.
Meta-Analysis of Proportions Schwarzer, Guido; Rücker, Gerta
Methods in molecular biology (Clifton, N.J.),
2022, Volume:
2345
Journal Article
The meta-analysis of single proportions has become a popular application over the last two decades. Especially, systematic reviews of prevalence studies are conducted in various fields of science, ...including medicine, ecology, psychology, or social sciences. In this chapter, we illustrate meta-analysis methods to pool single proportions and to compare proportions from two groups. We introduce classic approaches based on the inverse variance method as well as generalized linear mixed models taking the binary structure of the data into account. The most common transformations of proportions and their back-transformations are described both for individual studies and in the meta-analysis setting.
Abstract
Background
Intensive care unit (ICU) length of stay (LOS) and the risk adjusted equivalent (RALOS) have been used as quality metrics. The latter measures entail either ratio or difference ...formulations or ICU random effects (RE), which have not been previously compared.
Methods
From calendar year 2016 data of an adult ICU registry-database (Australia & New Zealand Intensive Care Society (ANZICS) CORE), LOS predictive models were established using linear (LMM) and generalised linear (GLMM) mixed models. Model fixed effects quality-metric formulations were estimated as RALOSR for LMM (geometric mean derived from log(ICU LOS)) and GLMM (day) and observed minus expected ICU LOS (OMELOS from GLMM). Metric confidence intervals (95%CI) were estimated by bootstrapping; random effects (RE) were predicted for LMM and GLMM. Forest-plot displays of ranked quality-metric point-estimates (95%CI) were generated for ICU hospital classifications (metropolitan, private, rural/regional, and tertiary). Robust rank confidence sets (point estimate and 95%CI), both marginal (pertaining to a singular ICU) and simultaneous (pertaining to all ICU differences), were established.
Results
The ICU cohort was of 94,361 patients from 125 ICUs (metropolitan 16.9%, private 32.8%, rural/regional 6.4%, tertiary 43.8%). Age (mean, SD) was 61.7 (17.5) years; 58.3% were male; APACHE III severity-of-illness score 54.6 (25.7); ICU annual patient volume 1192 (702) and ICU LOS 3.2 (4.9). There was no concordance of ICU ranked model predictions, GLMM versus LMM, nor for the quality metrics used, RALOSR, OMELOS and site-specific RE for each of the ICU hospital classifications. Furthermore, there was no concordance between ICU ranking confidence sets, marginal and simultaneous for models or quality metrics.
Conclusions
Inference regarding adjusted ICU LOS was dependent upon the statistical estimator and the quality index used to quantify any LOS differences across ICUs. That is, there was no “one best model”; thus, ICU “performance” is determined by model choice and any rankings thereupon should be circumspect.
Display omitted
•Fire-excluded sandhills require density control of plant life forms for longleaf pine regeneration.•Midstory saplings do not affect longleaf pine regeneration.•Longleaf pine impedes ...its own regeneration when basal area exceeds threshold values.•Sand pine and evergreen hardwoods impede longleaf pine regeneration, but deciduous oaks facilitate it.
Longleaf pine researchers have quantified many aspects of regeneration, growth and yield, but have focused on productive sites such as pine flatwoods and upland pine. Sandhill longleaf pine woodlands represent >40 % of remnant stands of longleaf pine ecosystems; most have no history of agriculture due to poor site productivity, yet many have been degraded by resin collection, grazing and fire exclusion. Restoration of structure and function of sandhill longleaf pine depends on understanding early stand-development processes, which are governed by fire, propagule availability, and competition from a range of overstory tree species not limited to longleaf pine. We analyzed a 16-year monitoring dataset from Eglin Air Force Base in northwest Florida, USA, focusing on how four tree growth forms influenced longleaf pine juvenile stages of seedling and sapling numbers and ingrowth to the subadult size class. Analyses were done with generalized linear mixed models using Poisson, normal, or negative binomial error distributions according to the response variable.
Restoration activities were effective at reducing numbers of sapling evergreen hardwoods, and basal area of overstory sand pine and other non-longleaf conifers was reduced by 35 % over the 16-year study. We observed threshold values for adult longleaf pine effects on seedlings and ingrowth; seedling numbers were maximized at 7.8 m2/ha, and ingrowth was maximum at 5.1 m2/ha. Seedling density was strongly governed by mast year. Sand pine, an invasive native pine, was detrimental to longleaf pine at all juvenile life stages even at low basal area. Deciduous oaks, which have been observed to facilitate longleaf pine seedling establishment, instead facilitated ingrowth to the subadult size class. Evergreen hardwoods, which may dampen fire behavior, only decreased longleaf pine density at the sapling stage. Ongoing silvicultural restoration activities should continue to focus on frequent fire and sand pine removal; deciduous oaks should be retained for their facilitative value; and removal of evergreen hardwoods should not be a high priority as long as prescribed fire can propagate in these lower productivity sites.
Quantifying sex-specific additive genetic variance (VA) in fitness, and the cross-sex genetic correlation (rA), is prerequisite to predicting evolutionary dynamics and the magnitude of sexual ...conflict. Further, quantifying VA and rA in underlying fitness components, and genetic consequences of immigration and resulting gene flow, is required to identify mechanisms that maintain VA in fitness. However, these key parameters have rarely been estimated in wild populations experiencing natural environmental variation and immigration. We used comprehensive pedigree and life-history data from song sparrows (Melospiza melodia) to estimate VA and rA in sex-specific fitness and underlying fitness components, and to estimate additive genetic effects of immigrants alongside inbreeding depression. We found evidence of substantial VA in female and male fitness, with a moderate positive crosssex rA. There was also substantial VA in male but not female adult reproductive success, and moderate VA in juvenile survival but not adult annual survival. Immigrants introduced alleles with negative additive genetic effects on local fitness, potentially reducing population mean fitness through migration load, but alleviating expression of inbreeding depression. Our results show that VA for fitness can be maintained in the wild, and be broadly concordant between the sexes despite marked sex-specific VA in reproductive success.