Lipidomic QTL in Diversity Outbred mice identifies a novel function for α/ß hydrolase domain 2 (Abhd2) as an enzyme that metabolizes phosphatidylcholine and cardiolipin.

We and others have previously shown that genetic association can be used to make causal connections between gene loci and small molecules measured by mass spectrometry in the bloodstream and in tissues. We identified a locus on mouse chromosome 7 where several phospholipids in liver showed strong genetic association to distinct gene loci. In this study, we integrated gene expression data with genetic association data to identify a single gene at the chromosome 7 locus as the driver of the phospholipid phenotypes. The gene encodes α/ß-hydrolase domain 2 (Abhd2), one of 23 members of the ABHD gene family. We validated this observation by measuring lipids in a mouse with a whole-body deletion of Abhd2. The Abhd2KO mice had a significant increase in liver levels of phosphatidylcholine and phosphatidylethanolamine. Unexpectedly, we also found a decrease in two key mitochondrial lipids, cardiolipin and phosphatidylglycerol, in male Abhd2KO mice. These data suggest that Abhd2 plays a role in the synthesis, turnover, or remodeling of liver phospholipids.

Identification of genetic drivers of plasma lipoprotein size in the Diversity Outbred mouse population.

Despite great progress in understanding lipoprotein physiology, there is still much to be learned about the genetic drivers of lipoprotein abundance, composition, and function. We used ion mobility spectrometry to survey 16 plasma lipoprotein subfractions in 500 Diversity Outbred mice maintained on a Western-style diet. We identified 21 quantitative trait loci (QTL) affecting lipoprotein abundance. To refine the QTL and link them to disease risk in humans, we asked if the human homologs of genes located at each QTL were associated with lipid traits in human genome-wide association studies. Integration of mouse QTL with human genome-wide association studies yielded candidate gene drivers for 18 of the 21 QTL. This approach enabled us to nominate the gene encoding the neutral ceramidase, Asah2, as a novel candidate driver at a QTL on chromosome 19 for large HDL particles (HDL-2b). To experimentally validate Asah2, we surveyed lipoproteins in Asah2-/- mice. Compared to wild-type mice, female Asah2-/- mice showed an increase in several lipoproteins, including HDL. Our results provide insights into the genetic regulation of circulating lipoproteins, as well as mechanisms by which lipoprotein subfractions may affect cardiovascular disease risk in humans.

Real-time estimation and forecasting of COVID-19 cases and hospitalizations in Wisconsin HERC regions for public health decision making processes.

BACKGROUND: The spread of the COVID-19 (SARS-CoV-2) and the surging number of cases across the United States have resulted in full hospitals and exhausted health care workers. Limited availability and questionable reliability of the data make outbreak prediction and resource planning difficult. Any estimates or forecasts are subject to high uncertainty and low accuracy to measure such components. The aim of this study is to apply, automate, and assess a Bayesian time series model for the real-time estimation and forecasting of COVID-19 cases and number of hospitalizations in Wisconsin healthcare emergency readiness coalition (HERC) regions. METHODS: This study makes use of the publicly available Wisconsin COVID-19 historical data by county. Cases and effective time-varying reproduction number [Formula: see text] by the HERC region over time are estimated using Bayesian latent variable models. Hospitalizations are estimated by the HERC region over time using a Bayesian regression model. Cases, effective Rt, and hospitalizations are forecasted over a 1-day, 3-day, and 7-day time horizon using the last 28 days of data, and the 20%, 50%, and 90% Bayesian credible intervals of the forecasts are calculated. The frequentist coverage probability is compared to the Bayesian credible level to evaluate performance. RESULTS: For cases and effective [Formula: see text], all three time horizons outperform the three credible levels of the forecast. For hospitalizations, all three time horizons outperform the 20% and 50% credible intervals of the forecast. On the contrary, the 1-day and 3-day periods underperform the 90% credible intervals. Questions about uncertainty quantification should be re-calculated using the frequentist coverage probability of the Bayesian credible interval based on observed data for all three metrics. CONCLUSIONS: We present an approach to automate the real-time estimation and forecasting of cases and hospitalizations and corresponding uncertainty using publicly available data. The models were able to infer short-term trends consistent with reported values at the HERC region level. Additionally, the models were able to accurately forecast and estimate the uncertainty of the measurements. This study can help identify the most affected regions and major outbreaks in the near future. The workflow can be adapted to other geographic regions, states, and even countries where decision-making processes are supported in real-time by the proposed modeling system.

Genetic determinants of gut microbiota composition and bile acid profiles in mice.

The microbial communities that inhabit the distal gut of humans and other mammals exhibit large inter-individual variation. While host genetics is a known factor that influences gut microbiota composition, the mechanisms underlying this variation remain largely unknown. Bile acids (BAs) are hormones that are produced by the host and chemically modified by gut bacteria. BAs serve as environmental cues and nutrients to microbes, but they can also have antibacterial effects. We hypothesized that host genetic variation in BA metabolism and homeostasis influence gut microbiota composition. To address this, we used the Diversity Outbred (DO) stock, a population of genetically distinct mice derived from eight founder strains. We characterized the fecal microbiota composition and plasma and cecal BA profiles from 400 DO mice maintained on a high-fat high-sucrose diet for ~22 weeks. Using quantitative trait locus (QTL) analysis, we identified several genomic regions associated with variations in both bacterial and BA profiles. Notably, we found overlapping QTL for Turicibacter sp. and plasma cholic acid, which mapped to a locus containing the gene for the ileal bile acid transporter, Slc10a2. Mediation analysis and subsequent follow-up validation experiments suggest that differences in Slc10a2 gene expression associated with the different strains influences levels of both traits and revealed novel interactions between Turicibacter and BAs. This work illustrates how systems genetics can be utilized to generate testable hypotheses and provide insight into host-microbe interactions.

Interventions to Disrupt Coronavirus Disease Transmission at a University, Wisconsin, USA, August-October 2020.

University settings have demonstrated potential for coronavirus disease (COVID-19) outbreaks; they combine congregate living, substantial social activity, and a young population predisposed to mild illness. Using genomic and epidemiologic data, we describe a COVID-19 outbreak at the University of Wisconsin-Madison, Madison, Wisconsin, USA. During August-October 2020, a total of 3,485 students, including 856/6,162 students living in dormitories, tested positive. Case counts began rising during move-in week, August 25-31, 2020, then rose rapidly during September 1-11, 2020. The university initiated multiple prevention efforts, including quarantining 2 dormitories; a subsequent decline in cases was observed. Genomic surveillance of cases from Dane County, in which the university is located, did not find evidence of transmission from a large cluster of cases in the 2 quarantined dorms during the outbreak. Coordinated implementation of prevention measures can reduce COVID-19 spread in university settings and may limit spillover to the surrounding community.

The Transcription Factor Nfatc2 Regulates ß-Cell Proliferation and Genes Associated with Type 2 Diabetes in Mouse and Human Islets.

Human genome-wide association studies (GWAS) have shown that genetic variation at >130 gene loci is associated with type 2 diabetes (T2D). We asked if the expression of the candidate T2D-associated genes within these loci is regulated by a common locus in pancreatic islets. Using an obese F2 mouse intercross segregating for T2D, we show that the expression of ~40% of the T2D-associated genes is linked to a broad region on mouse chromosome (Chr) 2. As all but 9 of these genes are not physically located on Chr 2, linkage to Chr 2 suggests a genomic factor(s) located on Chr 2 regulates their expression in trans. The transcription factor Nfatc2 is physically located on Chr 2 and its expression demonstrates cis linkage; i.e., its expression maps to itself. When conditioned on the expression of Nfatc2, linkage for the T2D-associated genes was greatly diminished, supporting Nfatc2 as a driver of their expression. Plasma insulin also showed linkage to the same broad region on Chr 2. Overexpression of a constitutively active (ca) form of Nfatc2 induced ß-cell proliferation in mouse and human islets, and transcriptionally regulated more than half of the T2D-associated genes. Overexpression of either ca-Nfatc2 or ca-Nfatc1 in mouse islets enhanced insulin secretion, whereas only ca-Nfatc2 was able to promote ß-cell proliferation, suggesting distinct molecular pathways mediating insulin secretion vs. ß-cell proliferation are regulated by NFAT. Our results suggest that many of the T2D-associated genes are downstream transcriptional targets of NFAT, and may act coordinately in a pathway through which NFAT regulates ß-cell proliferation in both mouse and human islets.

Integrative analysis of a cross-loci regulation network identifies App as a gene regulating insulin secretion from pancreatic islets.

Complex diseases result from molecular changes induced by multiple genetic factors and the environment. To derive a systems view of how genetic loci interact in the context of tissue-specific molecular networks, we constructed an F2 intercross comprised of >500 mice from diabetes-resistant (B6) and diabetes-susceptible (BTBR) mouse strains made genetically obese by the Leptin(ob/ob) mutation (Lep(ob)). High-density genotypes, diabetes-related clinical traits, and whole-transcriptome expression profiling in five tissues (white adipose, liver, pancreatic islets, hypothalamus, and gastrocnemius muscle) were determined for all mice. We performed an integrative analysis to investigate the inter-relationship among genetic factors, expression traits, and plasma insulin, a hallmark diabetes trait. Among five tissues under study, there are extensive protein-protein interactions between genes responding to different loci in adipose and pancreatic islets that potentially jointly participated in the regulation of plasma insulin. We developed a novel ranking scheme based on cross-loci protein-protein network topology and gene expression to assess each gene's potential to regulate plasma insulin. Unique candidate genes were identified in adipose tissue and islets. In islets, the Alzheimer's gene App was identified as a top candidate regulator. Islets from 17-week-old, but not 10-week-old, App knockout mice showed increased insulin secretion in response to glucose or a membrane-permeant cAMP analog, in agreement with the predictions of the network model. Our result provides a novel hypothesis on the mechanism for the connection between two aging-related diseases: Alzheimer's disease and type 2 diabetes.

Positional cloning of Sorcs1, a type 2 diabetes quantitative trait locus.

We previously mapped the type 2 diabetes mellitus-2 locus (T2dm2), which affects fasting insulin levels, to distal chromosome 19 in a leptin-deficient obese F2 intercross derived from C57BL/6 (B6) and BTBR T+ tf/J (BTBR) mice. Introgression of a 7-Mb segment of the B6 chromosome 19 into the BTBR background (strain 1339A) replicated the reduced insulin linked to T2dm2. The 1339A mice have markedly impaired insulin secretion in vivo and disrupted islet morphology. We used subcongenic strains derived from 1339A to localize the T2dm2 quantitative trait locus (QTL) to a 242-kb segment comprising the promoter, first exon and most of the first intron of the Sorcs1 gene. This was the only gene in the 1339A strain for which we detected amino acid substitutions and expression level differences between mice carrying B6 and BTBR alleles of this insert, thereby identifying variation within the Sorcs1 gene as underlying the phenotype associated with the T2dm2 locus. SorCS1 binds platelet-derived growth factor, a growth factor crucial for pericyte recruitment to the microvasculature, and may thus have a role in expanding or maintaining the islet vasculature. Our identification of the Sorcs1 gene provides insight into the pathway underlying the pathophysiology of obesity-induced type 2 diabetes mellitus.

Positional cloning of a type 2 diabetes quantitative trait locus; tomosyn-2, a negative regulator of insulin secretion.

We previously mapped a type 2 diabetes (T2D) locus on chromosome 16 (Chr 16) in an F2 intercross from the BTBR T (+) tf (BTBR) Lep(ob/ob) and C57BL/6 (B6) Lep(ob/ob) mouse strains. Introgression of BTBR Chr 16 into B6 mice resulted in a consomic mouse with reduced fasting plasma insulin and elevated glucose levels. We derived a panel of sub-congenic mice and narrowed the diabetes susceptibility locus to a 1.6 Mb region. Introgression of this 1.6 Mb fragment of the BTBR Chr 16 into lean B6 mice (B6.16(BT36-38)) replicated the phenotypes of the consomic mice. Pancreatic islets from the B6.16(BT36-38) mice were defective in the second phase of the insulin secretion, suggesting that the 1.6 Mb region encodes a regulator of insulin secretion. Within this region, syntaxin-binding protein 5-like (Stxbp5l) or tomosyn-2 was the only gene with an expression difference and a non-synonymous coding single nucleotide polymorphism (SNP) between the B6 and BTBR alleles. Overexpression of the b-tomosyn-2 isoform in the pancreatic ß-cell line, INS1 (832/13), resulted in an inhibition of insulin secretion in response to 3 mM 8-bromo cAMP at 7 mM glucose. In vitro binding experiments showed that tomosyn-2 binds recombinant syntaxin-1A and syntaxin-4, key proteins that are involved in insulin secretion via formation of the SNARE complex. The B6 form of tomosyn-2 is more susceptible to proteasomal degradation than the BTBR form, establishing a functional role for the coding SNP in tomosyn-2. We conclude that tomosyn-2 is the major gene responsible for the T2D Chr 16 quantitative trait locus (QTL) we mapped in our mouse cross. Our findings suggest that tomosyn-2 is a key negative regulator of insulin secretion.

Liver and adipose expression associated SNPs are enriched for association to type 2 diabetes.

Genome-wide association studies (GWAS) have demonstrated the ability to identify the strongest causal common variants in complex human diseases. However, to date, the massive data generated from GWAS have not been maximally explored to identify true associations that fail to meet the stringent level of association required to achieve genome-wide significance. Genetics of gene expression (GGE) studies have shown promise towards identifying DNA variations associated with disease and providing a path to functionally characterize findings from GWAS. Here, we present the first empiric study to systematically characterize the set of single nucleotide polymorphisms associated with expression (eSNPs) in liver, subcutaneous fat, and omental fat tissues, demonstrating these eSNPs are significantly more enriched for SNPs that associate with type 2 diabetes (T2D) in three large-scale GWAS than a matched set of randomly selected SNPs. This enrichment for T2D association increases as we restrict to eSNPs that correspond to genes comprising gene networks constructed from adipose gene expression data isolated from a mouse population segregating a T2D phenotype. Finally, by restricting to eSNPs corresponding to genes comprising an adipose subnetwork strongly predicted as causal for T2D, we dramatically increased the enrichment for SNPs associated with T2D and were able to identify a functionally related set of diabetes susceptibility genes. We identified and validated malic enzyme 1 (Me1) as a key regulator of this T2D subnetwork in mouse and provided support for the association of this gene to T2D in humans. This integration of eSNPs and networks provides a novel approach to identify disease susceptibility networks rather than the single SNPs or genes traditionally identified through GWAS, thereby extracting additional value from the wealth of data currently being generated by GWAS.

Expression of tau and phosphorylated tau in the brain of normal and hemiparkinsonian rhesus macaques.

Tau is a neuronal protein involved in microtubule stabilization and intracellular vesicle transport in axons. In neurodegenerative disorders termed "tauopathies," like Alzheimer's and Parkinson's disease, tau becomes hyperphosphorylated and forms intracellular inclusions. Rhesus macaques are widely used for studying ageing processes and modeling neurodegenerative disorders, yet little is known about endogenous tau expression in their brains. In this study, immunohistochemical methods were used to map and characterize total tau, 3R- and 4R-tau isoforms, and phosphorylated tau (pThr231-tau and pSer202/Thr205-tau/AT8) expression bilaterally in 16 brain regions of normal and 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP)-induced hemiparkinsonian adult rhesus macaques. Tau-immunoreactivity (-ir), including both 3R and 4R isoforms, was observed throughout the brain, with varying regional intensities. The anterior cingulate cortex, entorhinal cortex, and hippocampus displayed the most robust tau-ir, while the subthalamic nucleus and white matter regions had minimal expression. Tau was present in neurons of gray matter regions; it was preferentially observed in fibers of the globus pallidus and substantia nigra and in cell bodies of the thalamus and subthalamic nucleus. In white matter regions, tau was abundantly present in oligodendrocytes. Additionally, neuronal pThr231-tau-ir was abundant in all brain regions, but not AT8-ir. Differences in regional and intracellular protein expression were not detected between control subjects and both brain hemispheres of MPTP-treated animals. Specifically, tau-ir in the substantia nigra of all subjects colocalized with GABAergic neurons. Overall, this report provides an in-depth characterization of tau expression in the rhesus macaque brain to facilitate future investigations for understanding and modeling tau pathology in this species.

Lipidomic QTL in Diversity Outbred mice identifies a novel function for α/ß hydrolase domain 2 ( Abhd2 ) as an enzyme that metabolizes phosphatidylcholine and cardiolipin.

We and others have previously shown that genetic association can be used to make causal connections between gene loci and small molecules measured by mass spectrometry in the bloodstream and in tissues. We identified a locus on mouse chromosome 7 where several phospholipids in liver showed strong genetic association to distinct gene loci. In this study, we integrated gene expression data with genetic association data to identify a single gene at the chromosome 7 locus as the driver of the phospholipid phenotypes. The gene encodes α/ß-hydrolase domain 2 ( Abhd2 ), one of 23 members of the ABHD gene family. We validated this observation by measuring lipids in a mouse with a whole-body deletion of Abhd2 . The Abhd2 KO mice had a significant increase in liver levels of phosphatidylcholine and phosphatidylethanolamine. Unexpectedly, we also found a decrease in two key mitochondrial lipids, cardiolipin and phosphatidylglycerol, in male Abhd2 KO mice. These data suggest that Abhd2 plays a role in the synthesis, turnover, or remodeling of liver phospholipids.

Novel regulators of islet function identified from genetic variation in mouse islet Ca2+ oscillations.

Insufficient insulin secretion to meet metabolic demand results in diabetes. The intracellular flux of Ca2+ into ß-cells triggers insulin release. Since genetics strongly influences variation in islet secretory responses, we surveyed islet Ca2+ dynamics in eight genetically diverse mouse strains. We found high strain variation in response to four conditions: (1) 8 mM glucose; (2) 8 mM glucose plus amino acids; (3) 8 mM glucose, amino acids, plus 10 nM glucose-dependent insulinotropic polypeptide (GIP); and (4) 2 mM glucose. These stimuli interrogate ß-cell function, α- to ß-cell signaling, and incretin responses. We then correlated components of the Ca2+ waveforms to islet protein abundances in the same strains used for the Ca2+ measurements. To focus on proteins relevant to human islet function, we identified human orthologues of correlated mouse proteins that are proximal to glycemic-associated single-nucleotide polymorphisms in human genome-wide association studies. Several orthologues have previously been shown to regulate insulin secretion (e.g. ABCC8, PCSK1, and GCK), supporting our mouse-to-human integration as a discovery platform. By integrating these data, we nominate novel regulators of islet Ca2+ oscillations and insulin secretion with potential relevance for human islet function. We also provide a resource for identifying appropriate mouse strains in which to study these regulators.

Tsc2, a positional candidate gene underlying a quantitative trait locus for hepatic steatosis.

Nonalchoholic fatty liver disease (NAFLD) is the most common cause of liver dysfunction and is associated with metabolic diseases, including obesity, insulin resistance, and type 2 diabetes. We mapped a quantitative trait locus (QTL) for NAFLD to chromosome 17 in a cross between C57BL/6 (B6) and BTBR mouse strains made genetically obese with the Lep(ob/ob) mutation. We identified Tsc2 as a gene underlying the chromosome 17 NAFLD QTL. Tsc2 functions as an inhibitor of mammalian target of rapamycin, which is involved in many physiological processes, including cell growth, proliferation, and metabolism. We found that Tsc2(+/-) mice have increased lipogenic gene expression in the liver in an insulin-dependent manner. The coding single nucleotide polymorphism between the B6 and BTBR strains leads to a change in the ability to inhibit the expression of lipogenic genes and de novo lipogenesis in AML12 cells and to promote the proliferation of Ins1 cells. This difference is due to a different affinity of binding to Tsc1, which affects the stability of Tsc2.

Genetic mapping of paternal sorting of mitochondria in cucumber.

Mitochondria are organelles that have their own DNA; serve as the powerhouses of eukaryotic cells; play important roles in stress responses, programmed cell death, and ageing; and in the vast majority of eukaryotes, are maternally transmitted. Strict maternal transmission of mitochondria makes it difficult to select for better-performing mitochondria, or against deleterious mutations in the mitochondrial DNA. Cucumber is a useful plant for organellar genetics because its mitochondria are paternally transmitted and it possesses one of the largest mitochondrial genomes among all eukaryotes. Recombination among repetitive motifs in the cucumber mitochondrial DNA produces rearrangements associated with strongly mosaic (MSC) phenotypes. We previously reported nuclear control of sorting among paternally transmitted mitochondrial DNAs. The goal of this project was to map paternal sorting of mitochondria as a step towards its eventual cloning. We crossed single plants from plant introduction (PI) 401734 and Cucumis sativus var. hardwickii and produced an F(2) family. A total of 425 F(2) plants were genotyped for molecular markers and testcrossed as the female with MSC16. Testcross families were scored for frequencies of wild-type versus MSC progenies. Discrete segregations for percent wild-type progenies were not observed and paternal sorting of mitochondria was therefore analyzed as a quantitative trait. A major quantitative trait locus (QTL; LOD >23) was mapped between two simple sequence repeats encompassing a 459-kb region on chromosome 3. Nuclear genes previously shown to affect the prevalence of mitochondrial DNAs (MSH1, OSB1, and RECA homologs) were not located near this major QTL on chromosome 3. Sequencing of this region from PI 401734, together with improved annotation of the cucumber genome, should result in the eventual cloning of paternal sorting of mitochondria and provide insights about nuclear control of organellar-DNA sorting.

pH-Responsive Polymer Nanoparticles for Efficient Delivery of Cas9 Ribonucleoprotein With or Without Donor DNA.

Clustered regularly interspaced short palindromic repeat (CRISPR)-associated protein 9 (Cas9) may offer new therapeutics for genetic diseases through gene disruption via nonhomologous end joining (NHEJ) or gene correction via homology-directed repair (HDR). However, clinical translation of CRISPR technology is limited by the lack of safe and efficient delivery systems. Here, facilely fabricated pH-responsive polymer nanoparticles capable of safely and efficiently delivering Cas9 ribonucleoprotein alone (termed NHEJ-NP, diameter = 29.4 nm), or together with donor DNA (termed HDR-NP, diameter = 33.3 nm) are reported. Moreover, intravenously, intratracheally, and intramuscularly injected NHEJ-NP induces efficient gene editing in mouse liver, lung, and skeletal muscle, respectively. Intramuscularly injected HDR-NP also leads to muscle strength recovery in a Duchenne muscular dystrophy mouse model. NHEJ-NP and HDR-NP possess many desirable properties including high payload loading content, small and uniform sizes, high editing efficiency, good biocompatibility, low immunogenicity, and ease of production, storage, and transport, making them great interest for various genome editing applications with clinical potentials.

Genetic networks of liver metabolism revealed by integration of metabolic and transcriptional profiling.

Although numerous quantitative trait loci (QTL) influencing disease-related phenotypes have been detected through gene mapping and positional cloning, identification of the individual gene(s) and molecular pathways leading to those phenotypes is often elusive. One way to improve understanding of genetic architecture is to classify phenotypes in greater depth by including transcriptional and metabolic profiling. In the current study, we have generated and analyzed mRNA expression and metabolic profiles in liver samples obtained in an F2 intercross between the diabetes-resistant C57BL/6 leptin(ob/ob) and the diabetes-susceptible BTBR leptin(ob/ob) mouse strains. This cross, which segregates for genotype and physiological traits, was previously used to identify several diabetes-related QTL. Our current investigation includes microarray analysis of over 40,000 probe sets, plus quantitative mass spectrometry-based measurements of sixty-seven intermediary metabolites in three different classes (amino acids, organic acids, and acyl-carnitines). We show that liver metabolites map to distinct genetic regions, thereby indicating that tissue metabolites are heritable. We also demonstrate that genomic analysis can be integrated with liver mRNA expression and metabolite profiling data to construct causal networks for control of specific metabolic processes in liver. As a proof of principle of the practical significance of this integrative approach, we illustrate the construction of a specific causal network that links gene expression and metabolic changes in the context of glutamate metabolism, and demonstrate its validity by showing that genes in the network respond to changes in glutamine and glutamate availability. Thus, the methods described here have the potential to reveal regulatory networks that contribute to chronic, complex, and highly prevalent diseases and conditions such as obesity and diabetes.

Transcriptomic profiling of bovine IVF embryos revealed candidate genes and pathways involved in early embryonic development.

BACKGROUND: Early embryonic loss is a large contributor to infertility in cattle. Although genetic factors are known to affect early embryonic development, the discovery of such factors has been a serious challenge. The objective of this study was to identify genes differentially expressed between blastocysts and degenerative embryos at early stages of development. RESULTS: Using microarrays, genome-wide RNA expression was profiled and compared for in vitro fertilization (IVF) - derived blastocysts and embryos undergoing degenerative development up to the same time point. Surprisingly similar transcriptomic profiles were found in degenerative embryos and blastocysts. Nonetheless, we identified 67 transcripts that significantly differed between these two groups of embryos at a 15% false discovery rate, including 33 transcripts showing at least a two-fold difference. Several signaling and metabolic pathways were found to be associated with the developmental status of embryos, among which were previously known important steroid biosynthesis and cell communication pathways in early embryonic development. CONCLUSIONS: This study presents the first direct and comprehensive comparison of transcriptomes between IVF blastocysts and degenerative embryos, providing important information for potential genes and pathways associated with early embryonic development.

A model selection approach for the identification of quantitative trait loci in experimental crosses, allowing epistasis.

The identification of quantitative trait loci (QTL) and their interactions is a crucial step toward the discovery of genes responsible for variation in experimental crosses. The problem is best viewed as one of model selection, and the most important aspect of the problem is the comparison of models of different sizes. We present a penalized likelihood approach, with penalties on QTL and pairwise interactions chosen to control false positive rates. This extends the work of Broman and Speed to allow for pairwise interactions among QTL. A conservative version of our penalized LOD score provides strict control over the rate of extraneous QTL and interactions; a more liberal criterion is more lenient on interactions but seeks to maintain control over the rate of inclusion of false loci. The key advance is that one needs only to specify a target false positive rate rather than a prior on the number of QTL and interactions. We illustrate the use of our model selection criteria as exploratory tools; simulation studies demonstrate reasonable power to detect QTL. Our liberal criterion is comparable in power to two Bayesian approaches.

Association of Mobile Phone Location Data Indications of Travel and Stay-at-Home Mandates With COVID-19 Infection Rates in the US.

Importance: A stay-at-home social distancing mandate is a key nonpharmacological measure to reduce the transmission rate of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), but a high rate of adherence is needed. Objective: To examine the association between the rate of human mobility changes and the rate of confirmed cases of SARS-CoV-2 infection. Design, Setting, and Participants: This cross-sectional study used daily travel distance and home dwell time derived from millions of anonymous mobile phone location data from March 11 to April 10, 2020, provided by the Descartes Labs and SafeGraph to quantify the degree to which social distancing mandates were followed in the 50 US states and District of Columbia and the association of mobility changes with rates of coronavirus disease 2019 (COVID-19) cases. Exposure: State-level stay-at-home orders during the COVID-19 pandemic. Main Outcomes and Measures: The main outcome was the association of state-specific rates of COVID-19 confirmed cases with the change rates of median travel distance and median home dwell time of anonymous mobile phone users. The increase rates are measured by the exponent in curve fitting of the COVID-19 cumulative confirmed cases, while the mobility change (increase or decrease) rates were measured by the slope coefficient in curve fitting of median travel distance and median home dwell time for each state. Results: Data from more than 45 million anonymous mobile phone devices were analyzed. The correlation between the COVID-19 increase rate and travel distance decrease rate was -0.586 (95% CI, -0.742 to -0.370) and the correlation between COVID-19 increase rate and home dwell time increase rate was 0.526 (95% CI, 0.293 to 0.700). Increases in state-specific doubling time of total cases ranged from 1.0 to 6.9 days (median [interquartile range], 2.7 [2.3-3.3] days) before stay-at-home orders were enacted to 3.7 to 30.3 days (median [interquartile range], 6.0 [4.8-7.1] days) after stay-at-home social distancing orders were put in place, consistent with pandemic modeling results. Conclusions and Relevance: These findings suggest that stay-at-home social distancing mandates, when they were followed by measurable mobility changes, were associated with reduction in COVID-19 spread. These results come at a particularly critical period when US states are beginning to relax social distancing policies and reopen their economies. These findings support the efficacy of social distancing and could help inform future implementation of social distancing policies should they need to be reinstated during later periods of COVID-19 reemergence.