Impacts of Neanderthal-Introgressed Sequences on the Landscape of Human Gene Expression.

Regulatory variation influencing gene expression is a key contributor to phenotypic diversity, both within and between species. Unfortunately, RNA degrades too rapidly to be recovered from fossil remains, limiting functional genomic insights about our extinct hominin relatives. Many Neanderthal sequences survive in modern humans due to ancient hybridization, providing an opportunity to assess their contributions to transcriptional variation and to test hypotheses about regulatory evolution. We developed a flexible Bayesian statistical approach to quantify allele-specific expression (ASE) in complex RNA-seq datasets. We identified widespread expression differences between Neanderthal and modern human alleles, indicating pervasive cis-regulatory impacts of introgression. Brain regions and testes exhibited significant downregulation of Neanderthal alleles relative to other tissues, consistent with natural selection influencing the tissue-specific regulatory landscape. Our study demonstrates that Neanderthal-inherited sequences are not silent remnants of ancient interbreeding but have measurable impacts on gene expression that contribute to variation in modern human phenotypes.

The WHO estimates of excess mortality associated with the COVID-19 pandemic.

The World Health Organization has a mandate to compile and disseminate statistics on mortality, and we have been tracking the progression of the COVID-19 pandemic since the beginning of 20201. Reported statistics on COVID-19 mortality are problematic for many countries owing to variations in testing access, differential diagnostic capacity and inconsistent certification of COVID-19 as cause of death. Beyond what is directly attributable to it, the pandemic has caused extensive collateral damage that has led to losses of lives and livelihoods. Here we report a comprehensive and consistent measurement of the impact of the COVID-19 pandemic by estimating excess deaths, by month, for 2020 and 2021. We predict the pandemic period all-cause deaths in locations lacking complete reported data using an overdispersed Poisson count framework that applies Bayesian inference techniques to quantify uncertainty. We estimate 14.83 million excess deaths globally, 2.74 times more deaths than the 5.42 million reported as due to COVID-19 for the period. There are wide variations in the excess death estimates across the six World Health Organization regions. We describe the data and methods used to generate these estimates and highlight the need for better reporting where gaps persist. We discuss various summary measures, and the hazards of ranking countries' epidemic responses.

Rejoinder to the discussion on "The central role of the identifying assumption in population size estimation".

We organize the discussants' major comments into the following categories: sensitivity analyses, zero counts, model selection, the marginal no-highest-order interaction (NHOI) assumption, and the usefulness of our proposed framework.

The central role of the identifying assumption in population size estimation.

The problem of estimating the size of a population based on a subset of individuals observed across multiple data sources is often referred to as capture-recapture or multiple-systems estimation. This is fundamentally a missing data problem, where the number of unobserved individuals represents the missing data. As with any missing data problem, multiple-systems estimation requires users to make an untestable identifying assumption in order to estimate the population size from the observed data. If an appropriate identifying assumption cannot be found for a data set, no estimate of the population size should be produced based on that data set, as models with different identifying assumptions can produce arbitrarily different population size estimates-even with identical observed data fits. Approaches to multiple-systems estimation often do not explicitly specify identifying assumptions. This makes it difficult to decouple the specification of the model for the observed data from the identifying assumption and to provide justification for the identifying assumption. We present a re-framing of the multiple-systems estimation problem that leads to an approach that decouples the specification of the observed-data model from the identifying assumption, and discuss how common models fit into this framing. This approach takes advantage of existing software and facilitates various sensitivity analyses. We demonstrate our approach in a case study estimating the number of civilian casualties in the Kosovo war.

Estimating seroprevalence of SARS-CoV-2 in Ohio: A Bayesian multilevel poststratification approach with multiple diagnostic tests.

Globally, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has infected more than 59 million people and killed more than 1.39 million. Designing and monitoring interventions to slow and stop the spread of the virus require knowledge of how many people have been and are currently infected, where they live, and how they interact. The first step is an accurate assessment of the population prevalence of past infections. There are very few population-representative prevalence studies of SARS-CoV-2 infections, and only two states in the United States-Indiana and Connecticut-have reported probability-based sample surveys that characterize statewide prevalence of SARS-CoV-2. One of the difficulties is the fact that tests to detect and characterize SARS-CoV-2 coronavirus antibodies are new, are not well characterized, and generally function poorly. During July 2020, a survey representing all adults in the state of Ohio in the United States collected serum samples and information on protective behavior related to SARS-CoV-2 and coronavirus disease 2019 (COVID-19). Several features of the survey make it difficult to estimate past prevalence: 1) a low response rate; 2) a very low number of positive cases; and 3) the fact that multiple poor-quality serological tests were used to detect SARS-CoV-2 antibodies. We describe a Bayesian approach for analyzing the biomarker data that simultaneously addresses these challenges and characterizes the potential effect of selective response. The model does not require survey sample weights; accounts for multiple imperfect antibody test results; and characterizes uncertainty related to the sample survey and the multiple imperfect, potentially correlated tests.

Global, regional, and national estimates and trends in stillbirths from 2000 to 2019: a systematic assessment.

BACKGROUND: Stillbirths are a major public health issue and a sensitive marker of the quality of care around pregnancy and birth. The UN Global Strategy for Women's, Children's and Adolescents' Health (2016-30) and the Every Newborn Action Plan (led by UNICEF and WHO) call for an end to preventable stillbirths. A first step to prevent stillbirths is obtaining standardised measurement of stillbirth rates across countries. We estimated stillbirth rates and their trends for 195 countries from 2000 to 2019 and assessed progress over time. METHODS: For a systematic assessment, we created a dataset of 2833 country-year datapoints from 171 countries relevant to stillbirth rates, including data from registration and health information systems, household-based surveys, and population-based studies. After data quality assessment and exclusions, we used 1531 datapoints to estimate country-specific stillbirth rates for 195 countries from 2000 to 2019 using a Bayesian hierarchical temporal sparse regression model, according to a definition of stillbirth of at least 28 weeks' gestational age. Our model combined covariates with a temporal smoothing process such that estimates were informed by data for country-periods with high quality data, while being based on covariates for country-periods with little or no data on stillbirth rates. Bias and additional uncertainty associated with observations based on alternative stillbirth definitions and source types, and observations that were subject to non-sampling errors, were included in the model. We compared the estimated stillbirth rates and trends to previously reported mortality estimates in children younger than 5 years. FINDINGS: Globally in 2019, an estimated 2·0 million babies (90% uncertainty interval [UI] 1·9-2·2) were stillborn at 28 weeks or more of gestation, with a global stillbirth rate of 13·9 stillbirths (90% UI 13·5-15·4) per 1000 total births. Stillbirth rates in 2019 varied widely across regions, from 22·8 stillbirths (19·8-27·7) per 1000 total births in west and central Africa to 2·9 (2·7-3·0) in western Europe. After west and central Africa, eastern and southern Africa and south Asia had the second and third highest stillbirth rates in 2019. The global annual rate of reduction in stillbirth rate was estimated at 2·3% (90% UI 1·7-2·7) from 2000 to 2019, which was lower than the 2·9% (2·5-3·2) annual rate of reduction in neonatal mortality rate (for neonates aged <28 days) and the 4·3% (3·8-4·7) annual rate of reduction in mortality rate among children aged 1-59 months during the same period. Based on the lower bound of the 90% UIs, 114 countries had an estimated decrease in stillbirth rate since 2000, with four countries having a decrease of at least 50·0%, 28 having a decrease of 25·0-49·9%, 50 having a decrease of 10·0-24·9%, and 32 having a decrease of less than 10·0%. For the remaining 81 countries, we found no decrease in stillbirth rate since 2000. Of these countries, 34 were in sub-Saharan Africa, 16 were in east Asia and the Pacific, and 15 were in Latin America and the Caribbean. INTERPRETATION: Progress in reducing the rate of stillbirths has been slow compared with decreases in the mortality rate of children younger than 5 years. Accelerated improvements are most needed in the regions and countries with high stillbirth rates, particularly in sub-Saharan Africa. Future prevention of stillbirths needs increased efforts to raise public awareness, improve data collection, assess progress, and understand public health priorities locally, all of which require investment. FUNDING: Bill & Melinda Gates Foundation and the UK Foreign, Commonwealth and Development Office.

A linear noise approximation for stochastic epidemic models fit to partially observed incidence counts.

Stochastic epidemic models (SEMs) fit to incidence data are critical to elucidating outbreak dynamics, shaping response strategies, and preparing for future epidemics. SEMs typically represent counts of individuals in discrete infection states using Markov jump processes (MJPs), but are computationally challenging as imperfect surveillance, lack of subject-level information, and temporal coarseness of the data obscure the true epidemic. Analytic integration over the latent epidemic process is impossible, and integration via Markov chain Monte Carlo (MCMC) is cumbersome due to the dimensionality and discreteness of the latent state space. Simulation-based computational approaches can address the intractability of the MJP likelihood, but are numerically fragile and prohibitively expensive for complex models. A linear noise approximation (LNA) that approximates the MJP transition density with a Gaussian density has been explored for analyzing prevalence data in large-population settings, but requires modification for analyzing incidence counts without assuming that the data are normally distributed. We demonstrate how to reparameterize SEMs to appropriately analyze incidence data, and fold the LNA into a data augmentation MCMC framework that outperforms deterministic methods, statistically, and simulation-based methods, computationally. Our framework is computationally robust when the model dynamics are complex and applies to a broad class of SEMs. We evaluate our method in simulations that reflect Ebola, influenza, and SARS-CoV-2 dynamics, and apply our method to national surveillance counts from the 2013-2015 West Africa Ebola outbreak.

Combining information to estimate adherence in studies of pre-exposure prophylaxis for HIV prevention: Application to HPTN 067.

In trials of oral HIV pre-exposure prophylaxis (PrEP), multiple approaches have been used to measure adherence, including self-report, pill counts, electronic dose monitoring devices, and biological measures such as drug levels in plasma, peripheral blood mononuclear cells, hair, and/or dried blood spots. No one of these measures is ideal and each has strengths and weaknesses. However, accurate estimates of adherence to oral PrEP are important as drug efficacy is closely tied to adherence, and secondary analyses of trial data within identified adherent/non-adherent subgroups may yield important insights into real-world drug effectiveness. We develop a statistical approach to combining multiple measures of adherence and show in simulated data that the proposed method provides a more accurate measure of true adherence than self-report. We then apply the method to estimate adherence in the ADAPT study (HPTN 067) in South African women.

Pointless spatial modeling.

The analysis of area-level aggregated summary data is common in many disciplines including epidemiology and the social sciences. Typically, Markov random field spatial models have been employed to acknowledge spatial dependence and allow data-driven smoothing. In the context of an irregular set of areas, these models always have an ad hoc element with respect to the definition of a neighborhood scheme. In this article, we exploit recent theoretical and computational advances to carry out modeling at the continuous spatial level, which induces a spatial model for the discrete areas. This approach also allows reconstruction of the continuous underlying surface, but the interpretation of such surfaces is delicate since it depends on the quality, extent and configuration of the observed data. We focus on models based on stochastic partial differential equations. We also consider the interesting case in which the aggregate data are supplemented with point data. We carry out Bayesian inference and, in the language of generalized linear mixed models, if the link is linear, an efficient implementation of the model is available via integrated nested Laplace approximations. For nonlinear links, we present two approaches: a fully Bayesian implementation using a Hamiltonian Monte Carlo algorithm and an empirical Bayes implementation, that is much faster and is based on Laplace approximations. We examine the properties of the approach using simulation, and then apply the model to the classic Scottish lip cancer data.

Child mortality estimation incorporating summary birth history data.

The United Nations' Sustainable Development Goal 3.2 aims to reduce under-five child mortality to 25 deaths per 1000 live births by 2030. Child mortality tends to be concentrated in developing regions where information needed to assess achievement of this goal often comes from surveys and censuses. In both, women are asked about their birth histories, but with varying degrees of detail. Full birth history (FBH) data contain the reported dates of births and deaths of every surveyed mother's children. In contrast, summary birth history (SBH) data contain only the total number of children born and total number of children who died for each mother. Specialized methods are needed to accommodate this type of data into analyses of child mortality trends. We develop a data augmentation scheme within a Bayesian framework where for SBH data, birth and death dates are introduced as auxiliary variables. Since we specify a full probability model for the data, many of the well-known biases that exist in this data can be accommodated, along with space-time smoothing on the underlying mortality rates. We illustrate our approach in a simulation, showing robustness to model misspecification and that uncertainty is reduced when incorporating SBH data over simply analyzing all available FBH data. We also apply our approach to data from the Central region of Malawi and compare with the well-known Brass method.

Space-time modeling of child mortality at the Admin-2 level in a low and middle income countries context.

The Sustainable Development Goals call for a total reduction of preventable child mortality before 2030. Further, the goals state the desirability to have subnational mortality estimates. Estimates at this level are required for health interventions at the subnational level. In a low and middle income countries context, the data on mortality typically consist of household surveys, which are carried out with a stratified, cluster design, and census microsamples. Most household surveys collect full birth history (FBH) data on birth and death dates of a mother's children, but censuses collect summary birth history (SBH) data which consist only of the number of children born and the number that died. In previous work, direct (survey-weighted) estimates with associated variances were derived from FBH data and smoothed in space and time. Unfortunately, the FBH data from household surveys are usually not sufficiently abundant to obtain yearly estimates at the Admin-2 level (at which interventions are often made). In this paper we describe four extensions to previous work: (i) combining SBH data with FBH data, (ii) modeling on a yearly scale, to combine data on a yearly scale with data at coarser time scales, (iii) adjusting direct estimates in Admin-2 areas where we do not observe any deaths due to small sample sizes, (iv) acknowledge differences in data sources by modeling potential bias arising from the various data sources. The methods are illustrated using household survey and census data from Kenya and Malawi, to produce mortality estimates from 1980 to the time of the most recent survey, and predictions to 2020.

Ecological inference for infectious disease data, with application to vaccination strategies.

Disease surveillance systems provide a rich source of data regarding infectious diseases, aggregated across geographical regions. The analysis of such ecological data is fraught with difficulties, and, unless care and suitable data summaries are available, will lead to biased estimates of individual-level parameters. We consider using surveillance data to study the impacts of vaccination. To catalog the problems of ecological inference, we start with an individual-level model, which contains familiar parameters, and derive an ecologically consistent model for infectious diseases in partially vaccinated populations. We compare with other popular model classes and highlight deficiencies. We explore the properties of the new model through simulation and demonstrate that, under standard assumptions, the ecological model provides less biased estimates. We then fit the new model to data collected on measles outbreaks in Germany from 2005-2007.

Small Area Estimation for Disease Prevalence Mapping.

Small area estimation (SAE) entails estimating characteristics of interest for domains, often geographical areas, in which there may be few or no samples available. SAE has a long history and a wide variety of methods have been suggested, from a bewildering range of philosophical standpoints. We describe design-based and model-based approaches and models that are specified at the area-level and at the unit-level, focusing on health applications and fully Bayesian spatial models. The use of auxiliary information is a key ingredient for successful inference when response data are sparse and we discuss a number of approaches that allow the inclusion of covariate data. SAE for HIV prevalence, using data collected from a Demographic Health Survey in Malawi in 2015-2016, is used to illustrate a number of techniques. The potential use of SAE techniques for outcomes related to COVID-19 is discussed.

Geogenomic Segregation and Temporal Trends of Human Pathogenic Escherichia coli O157:H7, Washington, USA, 2005-20141.

The often-noted and persistent increased incidence of Escherichia coli O157:H7 infections in rural areas is not well understood. We used a cohort of E. coli O157:H7 cases reported in Washington, USA, during 2005-2014, along with phylogenomic characterization of the infecting isolates, to identify geographic segregation of and temporal trends in specific phylogenetic lineages of E. coli O157:H7. Kernel estimation and generalized additive models demonstrated that pathogen lineages were spatially segregated during the period of analysis and identified a focus of segregation spanning multiple, predominantly rural, counties for each of the main clinical lineages, Ib, IIa, and IIb. These results suggest the existence of local reservoirs from which humans are infected. We also noted a secular increase in the proportion of lineage IIa and IIb isolates. Spatial segregation by phylogenetic lineage offers the potential to identify local reservoirs and intervene to prevent continued transmission.

Heritable variation of mRNA decay rates in yeast.

Gene expression levels are determined by the balance between rates of mRNA transcription and decay, and genetic variation in either of these processes can result in heritable differences in transcript abundance. Although the genetics of gene expression has been a subject of intense interest, the contribution of heritable variation in mRNA decay rates to gene expression variation has received far less attention. To this end, we developed a novel statistical framework and measured allele-specific differences in mRNA decay rates in a diploid yeast hybrid created by mating two genetically diverse parental strains. We estimate that 31% of genes exhibit allelic differences in mRNA decay rates, of which 350 can be identified at a false discovery rate of 10%. Genes with significant allele-specific differences in mRNA decay rates have higher levels of polymorphism compared to other genes, with all gene regions contributing to allelic differences in mRNA decay rates. Strikingly, we find widespread evidence for compensatory evolution, such that variants influencing transcriptional initiation and decay have opposite effects, suggesting that steady-state gene expression levels are subject to pervasive stabilizing selection. Our results demonstrate that heritable differences in mRNA decay rates are widespread and are an important target for natural selection to maintain or fine-tune steady-state gene expression levels.

Time series modeling of pathogen-specific disease probabilities with subsampled data.

Many diseases arise due to exposure to one of multiple possible pathogens. We consider the situation in which disease counts are available over time from a study region, along with a measure of clinical disease severity, for example, mild or severe. In addition, we suppose a subset of the cases are lab tested in order to determine the pathogen responsible for disease. In such a context, we focus interest on modeling the probabilities of disease incidence given pathogen type. The time course of these probabilities is of great interest as is the association with time-varying covariates such as meteorological variables. In this set up, a natural Bayesian approach would be based on imputation of the unsampled pathogen information using Markov Chain Monte Carlo but this is computationally challenging. We describe a practical approach to inference that is easy to implement. We use an empirical Bayes procedure in a first step to estimate summary statistics. We then treat these summary statistics as the observed data and develop a Bayesian generalized additive model. We analyze data on hand, foot, and mouth disease (HFMD) in China in which there are two pathogens of primary interest, enterovirus 71 (EV71) and Coxackie A16 (CA16). We find that both EV71 and CA16 are associated with temperature, relative humidity, and wind speed, with reasonably similar functional forms for both pathogens. The important issue of confounding by time is modeled using a penalized B-spline model with a random effects representation. The level of smoothing is addressed by a careful choice of the prior on the tuning variance.

Seasonality of acquisition of respiratory bacterial pathogens in young children with cystic fibrosis.

BACKGROUND: Seasonal variations are often observed for respiratory tract infections; however, limited information is available regarding seasonal patterns of acquisition of common cystic fibrosis (CF)-related respiratory pathogens. We previously reported differential seasonal acquisition of Pseudomonas aeruginosa in young children with CF and no such variation for methicillin-susceptible Staphylococcus aureus acquisition. The purpose of this study was to describe and compare the seasonal incidence of acquisition of other respiratory bacterial pathogens in young children with CF. METHODS: We conducted a retrospective study to describe and compare the seasonal incidence of methicillin-resistant Staphylococcus aureus (MRSA), Stenotrophomonas maltophilia, Achromobacter xylosoxidans, and Haemophilus influenzae acquisition in young CF patients residing in the U.S. using the Cystic Fibrosis Foundation National Patient Registry, 2003-2009. Log-linear overdispersed Poisson regression was used to evaluate seasonal acquisition of each of these pathogens. RESULTS: A total of 4552 children met inclusion criteria. During follow-up 910 (20%), 1161 (26%), 228 (5%), and 2148 (47%) children acquired MRSA, S. maltophilia, A. xylosoxidans and H. influenzae, respectively. Compared to winter season, MRSA was less frequently acquired in spring (Incidence Rate Ratio [IRR]: 0.79; 95% Confidence Interval [CI]: 0.65, 0.96) and summer (IRR: 0.69; 95% CI: 0.57, 0.84) seasons. Similarly, a lower rate of A. xylosoxidans acquisition was observed in spring (IRR: 0.59; 95% CI: 0.39, 0.89). For H. influenzae, summer (IRR: 0.88; 95% CI: 0.78, 0.99) and autumn (IRR: 0.78; 95% CI: 0.69, 0.88) seasons were associated with lower acquisition rates compared to winter. No seasonal variation was observed for S. maltophilia acquisition. CONCLUSION: Acquisition of CF-related respiratory pathogens displays seasonal variation in young children with CF, with the highest rate of acquisition for most pathogens occurring in the winter. Investigation of factors underlying these observed associations may contribute to our understanding of the aetiology of these infections and guide future infection control strategies.

Air pollution exposure is associated with MRSA acquisition in young U.S. children with cystic fibrosis.

BACKGROUND: The role of air pollution in increasing susceptibility to respiratory tract infections in the cystic fibrosis (CF) population has not been well described. We recently demonstrated that chronic PM2.5 exposure is associated with an increased risk of initial Pseudomonas aeruginosa acquisition in young children with CF. The purpose of this study was to determine whether PM2.5 exposure is a risk factor for acquisition of other respiratory pathogens in young children with CF. METHODS: We conducted a retrospective study of initial acquisition of methicillin susceptible and methicillin resistant Staphylococcus aureus (MSSA and MRSA), Stenotrophomonas maltophilia and Achromobacter xylosoxidans in U.S. children <6 years of age with CF using the CF Foundation Patient Registry, 2003-2009. Multivariable Weibull regression with interval-censored outcomes was used to evaluate the association of PM2.5 concentration in the year prior to birth and risk of acquisition of each organism. RESULTS: During follow-up 63%, 17%, 24%, and 5% of children acquired MSSA, MRSA, S. maltophilia, and A. xylosoxidans, respectively. A 10 µg/m3 increase in PM2.5 exposure was associated with a 68% increased risk of MRSA acquisition (Hazard Ratio: 1.68; 95% Confidence Interval: 1.24, 2.27). PM2.5 was not associated with acquisition of other respiratory pathogens. CONCLUSIONS: Fine particulate matter is an independent risk factor for initial MRSA acquisition in young children with CF. These results support the increasing evidence that air pollution contributes to pulmonary morbidities in the CF community.

Evolution and genetic architecture of chromatin accessibility and function in yeast.

Chromatin accessibility is an important functional genomics phenotype that influences transcription factor binding and gene expression. Genome-scale technologies allow chromatin accessibility to be mapped with high-resolution, facilitating detailed analyses into the genetic architecture and evolution of chromatin structure within and between species. We performed Formaldehyde-Assisted Isolation of Regulatory Elements sequencing (FAIRE-Seq) to map chromatin accessibility in two parental haploid yeast species, Saccharomyces cerevisiae and Saccharomyces paradoxus and their diploid hybrid. We show that although broad-scale characteristics of the chromatin landscape are well conserved between these species, accessibility is significantly different for 947 regions upstream of genes that are enriched for GO terms such as intracellular transport and protein localization exhibit. We also develop new statistical methods to investigate the genetic architecture of variation in chromatin accessibility between species, and find that cis effects are more common and of greater magnitude than trans effects. Interestingly, we find that cis and trans effects at individual genes are often negatively correlated, suggesting widespread compensatory evolution to stabilize levels of chromatin accessibility. Finally, we demonstrate that the relationship between chromatin accessibility and gene expression levels is complex, and a significant proportion of differences in chromatin accessibility might be functionally benign.

Integrative phenomics reveals insight into the structure of phenotypic diversity in budding yeast.

To better understand the quantitative characteristics and structure of phenotypic diversity, we measured over 14,000 transcript, protein, metabolite, and morphological traits in 22 genetically diverse strains of Saccharomyces cerevisiae. More than 50% of all measured traits varied significantly across strains [false discovery rate (FDR) = 5%]. The structure of phenotypic correlations is complex, with 85% of all traits significantly correlated with at least one other phenotype (median = 6, maximum = 328). We show how high-dimensional molecular phenomics data sets can be leveraged to accurately predict phenotypic variation between strains, often with greater precision than afforded by DNA sequence information alone. These results provide new insights into the spectrum and structure of phenotypic diversity and the characteristics influencing the ability to accurately predict phenotypes.