Sub-communities of the vaginal microbiota in pregnant and non-pregnant women.

Diverse and non-Lactobacillus-dominated vaginal microbial communities are associated with adverse health outcomes such as preterm birth and the acquisition of sexually transmitted infections. Despite the importance of recognizing and understanding the key risk-associated features of these communities, their heterogeneous structure and properties remain ill-defined. Clustering approaches are commonly used to characterize vaginal communities, but they lack sensitivity and robustness in resolving substructures and revealing transitions between potential sub-communities. Here, we address this need with an approach based on mixed membership topic models. Using longitudinal data from cohorts of pregnant and non-pregnant study participants, we show that topic models more accurately describe sample composition, longitudinal changes, and better predict the loss of Lactobacillus dominance. We identify several non-Lactobacillus-dominated sub-communities common to both cohorts and independent of reproductive status. In non-pregnant individuals, we find that the menstrual cycle modulates transitions between and within sub-communities, as well as the concentrations of half of the cytokines and 18% of metabolites. Overall, our analyses based on mixed membership models reveal substructures of vaginal ecosystems which may have important clinical and biological associations.

Food & You: A digital cohort on personalized nutrition.

Nutrition is a key contributor to health. Recently, several studies have identified associations between factors such as microbiota composition and health-related responses to dietary intake, raising the potential of personalized nutritional recommendations. To further our understanding of personalized nutrition, detailed individual data must be collected from participants in their day-to-day lives. However, this is challenging in conventional studies that require clinical measurements and site visits. So-called digital or remote cohorts allow in situ data collection on a daily basis through mobile applications, online services, and wearable sensors, but they raise questions about study retention and data quality. "Food & You" is a personalized nutrition study implemented as a digital cohort in which participants track food intake, physical activity, gut microbiota, glycemia, and other data for two to four weeks. Here, we describe the study protocol, report on study completion rates, and describe the collected data, focusing on assessing their quality and reliability. Overall, the study collected data from over 1000 participants, including high-resolution data of nutritional intake of more than 46 million kcal collected from 315,126 dishes over 23,335 participant days, 1,470,030 blood glucose measurements, 49,110 survey responses, and 1,024 stool samples for gut microbiota analysis. Retention was high, with over 60% of the enrolled participants completing the study. Various data quality assessment efforts suggest the captured high-resolution nutritional data accurately reflect individual diet patterns, paving the way for digital cohorts as a typical study design for personalized nutrition.

Abrupt perturbation and delayed recovery of the vaginal ecosystem following childbirth.

The vaginal ecosystem is closely tied to human health and reproductive outcomes, yet its dynamics in the wake of childbirth remain poorly characterized. Here, we profile the vaginal microbiota and cytokine milieu of participants sampled longitudinally throughout pregnancy and for at least one year postpartum. We show that delivery, regardless of mode, is associated with a vaginal pro-inflammatory cytokine response and the loss of Lactobacillus dominance. By contrast, neither the progression of gestation nor the approach of labor strongly altered the vaginal ecosystem. At 9.5-months postpartum-the latest timepoint at which cytokines were assessed-elevated inflammation coincided with vaginal bacterial communities that had remained perturbed (highly diverse) from the time of delivery. Time-to-event analysis indicated a one-year postpartum probability of transitioning to Lactobacillus dominance of 49.4%. As diversity and inflammation declined during the postpartum period, dominance by L. crispatus, the quintessential health-associated commensal, failed to return: its prevalence before, immediately after, and one year after delivery was 41%, 4%, and 9%, respectively. Revisiting our pre-delivery data, we found that a prior live birth was associated with a lower odds of L. crispatus dominance in pregnant participants-an outcome modestly tempered by a longer ( > 18-month) interpregnancy interval. Our results suggest that reproductive history and childbirth in particular remodel the vaginal ecosystem and that the timing and degree of recovery from delivery may help determine the subsequent health of the woman and of future pregnancies.

Multiscale analysis of count data through topic alignment.

Topic modeling is a popular method used to describe biological count data. With topic models, the user must specify the number of topics $K$. Since there is no definitive way to choose $K$ and since a true value might not exist, we develop a method, which we call topic alignment, to study the relationships across models with different $K$. In addition, we present three diagnostics based on the alignment. These techniques can show how many topics are consistently present across different models, if a topic is only transiently present, or if a topic splits into more topics when $K$ increases. This strategy gives more insight into the process of generating the data than choosing a single value of $K$ would. We design a visual representation of these cross-model relationships, show the effectiveness of these tools for interpreting the topics on simulated and real data, and release an accompanying R package, alto.

Labeling Self-Tracked Menstrual Health Records With Hidden Semi-Markov Models.

Globally, millions of women track their menstrual cycle and fertility via smartphone-based health apps, generating multivariate time series with frequent missing data. To leverage this type of data for studies of fertility or studies of the effect of the menstrual cycle on symptoms and diseases, it is critical to have methods for identifying reproductive events, such as ovulation, pregnancy losses or births. Here, we present a hierarchical approach relying on hidden semi-Markov models that adapts to changes in tracking behavior, explicitly captures variable- and state- dependent missingness, allows for variables of different type, and quantifies uncertainty. The accuracy on simulated data reaches 98% with no missing data and 90% with realistic missingness. On our partially labeled real-world time series, the accuracy reaches 93%. Our method also accurately predicts cycle length by learning user characteristics. Its implementation is publicly available (HiddenSemiMarkov R package) and transferable to any health time series, including self-reported symptoms and occasional tests.

Predicting pregnancy using large-scale data from a women's health tracking mobile application.

Predicting pregnancy has been a fundamental problem in women's health for more than 50 years. Previous datasets have been collected via carefully curated medical studies, but the recent growth of women's health tracking mobile apps offers potential for reaching a much broader population. However, the feasibility of predicting pregnancy from mobile health tracking data is unclear. Here we develop four models - a logistic regression model, and 3 LSTM models - to predict a woman's probability of becoming pregnant using data from a women's health tracking app, Clue by BioWink GmbH. Evaluating our models on a dataset of 79 million logs from 65,276 women with ground truth pregnancy test data, we show that our predicted pregnancy probabilities meaningfully stratify women: women in the top 10% of predicted probabilities have a 89% chance of becoming pregnant over 6 menstrual cycles, as compared to a 27% chance for women in the bottom 10%. We develop a technique for extracting interpretable time trends from our deep learning models, and show these trends are consistent with previous fertility research. Our findings illustrate the potential that women's health tracking data offers for predicting pregnancy on a broader population; we conclude by discussing the steps needed to fulfill this potential.

Assessment of menstrual health status and evolution through mobile apps for fertility awareness.

For most women of reproductive age, assessing menstrual health and fertility typically involves regular visits to a gynecologist or another clinician. While these evaluations provide critical information on an individual's reproductive health status, they typically rely on memory-based self-reports, and the results are rarely, if ever, assessed at the population level. In recent years, mobile apps for menstrual tracking have become very popular, allowing us to evaluate the reliability and tracking frequency of millions of self-observations, thereby providing an unparalleled view, both in detail and scale, on menstrual health and its evolution for large populations. In particular, the primary aim of this study was to describe the tracking behavior of the app users and their overall observation patterns in an effort to understand if they were consistent with previous small-scale medical studies. The secondary aim was to investigate whether their precision allowed the detection and estimation of ovulation timing, which is critical for reproductive and menstrual health. Retrospective self-observation data were acquired from two mobile apps dedicated to the application of the sympto-thermal fertility awareness method, resulting in a dataset of more than 30 million days of observations from over 2.7 million cycles for two hundred thousand users. The analysis of the data showed that up to 40% of the cycles in which users were seeking pregnancy had recordings every single day. With a modeling approach using Hidden Markov Models to describe the collected data and estimate ovulation timing, it was found that follicular phases average duration and range were larger than previously reported, with only 24% of ovulations occurring at cycle days 14 to 15, while the luteal phase duration and range were in line with previous reports, although short luteal phases (10 days or less) were more frequently observed (in up to 20% of cycles). The digital epidemiology approach presented here can help to lead to a better understanding of menstrual health and its connection to women's health overall, which has historically been severely understudied.

Circadian clock-dependent and -independent posttranscriptional regulation underlies temporal mRNA accumulation in mouse liver.

The mammalian circadian clock coordinates physiology with environmental cycles through the regulation of daily oscillations of gene expression. Thousands of transcripts exhibit rhythmic accumulations across mouse tissues, as determined by the balance of their synthesis and degradation. While diurnally rhythmic transcription regulation is well studied and often thought to be the main factor generating rhythmic mRNA accumulation, the extent of rhythmic posttranscriptional regulation is debated, and the kinetic parameters (e.g., half-lives), as well as the underlying regulators (e.g., mRNA-binding proteins) are relatively unexplored. Here, we developed a quantitative model for cyclic accumulations of pre-mRNA and mRNA from total RNA-seq data, and applied it to mouse liver. This allowed us to identify that about 20% of mRNA rhythms were driven by rhythmic mRNA degradation, and another 15% of mRNAs regulated by both rhythmic transcription and mRNA degradation. The method could also estimate mRNA half-lives and processing times in intact mouse liver. We then showed that, depending on mRNA half-life, rhythmic mRNA degradation can either amplify or tune phases of mRNA rhythms. By comparing mRNA rhythms in wild-type and Bmal1-/- animals, we found that the rhythmic degradation of many transcripts did not depend on a functional BMAL1. Interestingly clock-dependent and -independent degradation rhythms peaked at distinct times of day. We further predicted mRNA-binding proteins (mRBPs) that were implicated in the posttranscriptional regulation of mRNAs, either through stabilizing or destabilizing activities. Together, our results demonstrate how posttranscriptional regulation temporally shapes rhythmic mRNA accumulation in mouse liver.

Non-circadian expression masking clock-driven weak transcription rhythms in U2OS cells.

U2OS cells harbor a circadian clock but express only a few rhythmic genes in constant conditions. We identified 3040 binding sites of the circadian regulators BMAL1, CLOCK and CRY1 in the U2OS genome. Most binding sites even in promoters do not correlate with detectable rhythmic transcript levels. Luciferase fusions reveal that the circadian clock supports robust but low amplitude transcription rhythms of representative promoters. However, rhythmic transcription of these potentially clock-controlled genes is masked by non-circadian transcription that overwrites the weaker contribution of the clock in constant conditions. Our data suggest that U2OS cells harbor an intrinsically rather weak circadian oscillator. The oscillator has the potential to regulate a large number of genes. The contribution of circadian versus non-circadian transcription is dependent on the metabolic state of the cell and may determine the apparent complexity of the circadian transcriptome.

The circadian clock coordinates ribosome biogenesis.

Biological rhythms play a fundamental role in the physiology and behavior of most living organisms. Rhythmic circadian expression of clock-controlled genes is orchestrated by a molecular clock that relies on interconnected negative feedback loops of transcription regulators. Here we show that the circadian clock exerts its function also through the regulation of mRNA translation. Namely, the circadian clock influences the temporal translation of a subset of mRNAs involved in ribosome biogenesis by controlling the transcription of translation initiation factors as well as the clock-dependent rhythmic activation of signaling pathways involved in their regulation. Moreover, the circadian oscillator directly regulates the transcription of ribosomal protein mRNAs and ribosomal RNAs. Thus the circadian clock exerts a major role in coordinating transcription and translation steps underlying ribosome biogenesis.

Genome-wide RNA polymerase II profiles and RNA accumulation reveal kinetics of transcription and associated epigenetic changes during diurnal cycles.

Interactions of cell-autonomous circadian oscillators with diurnal cycles govern the temporal compartmentalization of cell physiology in mammals. To understand the transcriptional and epigenetic basis of diurnal rhythms in mouse liver genome-wide, we generated temporal DNA occupancy profiles by RNA polymerase II (Pol II) as well as profiles of the histone modifications H3K4me3 and H3K36me3. We used these data to quantify the relationships of phases and amplitudes between different marks. We found that rhythmic Pol II recruitment at promoters rather than rhythmic transition from paused to productive elongation underlies diurnal gene transcription, a conclusion further supported by modeling. Moreover, Pol II occupancy preceded mRNA accumulation by 3 hours, consistent with mRNA half-lives. Both methylation marks showed that the epigenetic landscape is highly dynamic and globally remodeled during the 24-hour cycle. While promoters of transcribed genes had tri-methylated H3K4 even at their trough activity times, tri-methylation levels reached their peak, on average, 1 hour after Pol II. Meanwhile, rhythms in tri-methylation of H3K36 lagged transcription by 3 hours. Finally, modeling profiles of Pol II occupancy and mRNA accumulation identified three classes of genes: one showing rhythmicity both in transcriptional and mRNA accumulation, a second class with rhythmic transcription but flat mRNA levels, and a third with constant transcription but rhythmic mRNAs. The latter class emphasizes widespread temporally gated posttranscriptional regulation in the mouse liver.