GenomeMUSter mouse genetic variation service enables multitrait, multipopulation data integration and analysis.

Hundreds of inbred mouse strains and intercross populations have been used to characterize the function of genetic variants that contribute to disease. Thousands of disease-relevant traits have been characterized in mice and made publicly available. New strains and populations including consomics, the collaborative cross, expanded BXD, and inbred wild-derived strains add to existing complex disease mouse models, mapping populations, and sensitized backgrounds for engineered mutations. The genome sequences of inbred strains, along with dense genotypes from others, enable integrated analysis of trait-variant associations across populations, but these analyses are hampered by the sparsity of genotypes available. Moreover, the data are not readily interoperable with other resources. To address these limitations, we created a uniformly dense variant resource by harmonizing multiple data sets. Missing genotypes were imputed using the Viterbi algorithm with a data-driven technique that incorporates local phylogenetic information, an approach that is extendable to other model organisms. The result is a web- and programmatically accessible data service called GenomeMUSter, comprising single-nucleotide variants covering 657 strains at 106.8 million segregating sites. Interoperation with phenotype databases, analytic tools, and other resources enable a wealth of applications, including multitrait, multipopulation meta-analysis. We show this in cross-species comparisons of type 2 diabetes and substance use disorder meta-analyses, leveraging mouse data to characterize the likely role of human variant effects in disease. Other applications include refinement of mapped loci and prioritization of strain backgrounds for disease modeling to further unlock extant mouse diversity for genetic and genomic studies in health and disease.

Antenatal care is associated with adherence to iron supplementation among pregnant women in selected low-middle-income-countries of Asia, Africa, and Latin America & the Caribbean regions: Insights from Demographic and Health Surveys.

Anaemia is a global public health problem affecting 800 million women and children globally. Anaemia is associated with perinatal mortality, child morbidity and mortality, mental development, immune competence, susceptibility to lead poisoning and performance at work. The objective of this article is to identify whether antenatal care-seeking was associated with the uptake of iron supplementation among pregnant women, adjusting for a range of covariates. This article used data from the cross-sectional recent Demographic and Health Surveys (DHS) of 12 countries in Asia, Africa and Latin America & the Caribbean regions. The individual-level data from 273,144 women of reproductive age (15-49 years) were analysed from multi-country DHS. Multiple Logistic regression analyses were conducted using Predictive Analytics Software for Windows (PASW), Release 18.0. Receiving at least four antenatal care visits was significantly associated with the consumption of 90 or more iron-containing supplements in 12 low and middle income countries across three regions after adjusting for different household and respondent characteristics, while mass media exposure was found to be a significant predictor in India and Indonesia. Antenatal care seems to be the most important predictor of adherence to iron intake in the selected countries across Africa, Asia, Latin America and Caribbean regions.

Statistical analysis and modeling of the geometry and topology of plant roots.

The root is an important organ of a plant since it is responsible for water and nutrient uptake. Analyzing and modeling variabilities in the geometry and topology of roots can help in assessing the plant's health, understanding its growth patterns, and modeling relations between plant species and between plants and their environment. In this article, we develop a framework for the statistical analysis and modeling of the geometry and topology of plant roots. We represent root structures as points in a tree-shape space equipped with a metric that quantifies geometric and topological differences between pairs of roots. We then use these building blocks to compute geodesics, i.e., optimal deformations under the metric between root structures, and to perform statistical analysis on root populations. We demonstrate the utility of the proposed framework through an application to a dataset of wheat roots grown in different environmental conditions. We also show that the framework can be used in various applications including classification and regression.

Mapping population-based structural connectomes.

Advances in understanding the structural connectomes of human brain require improved approaches for the construction, comparison and integration of high-dimensional whole-brain tractography data from a large number of individuals. This article develops a population-based structural connectome (PSC) mapping framework to address these challenges. PSC simultaneously characterizes a large number of white matter bundles within and across different subjects by registering different subjects' brains based on coarse cortical parcellations, compressing the bundles of each connection, and extracting novel connection weights. A robust tractography algorithm and streamline post-processing techniques, including dilation of gray matter regions, streamline cutting, and outlier streamline removal are applied to improve the robustness of the extracted structural connectomes. The developed PSC framework can be used to reproducibly extract binary networks, weighted networks and streamline-based brain connectomes. We apply the PSC to Human Connectome Project data to illustrate its application in characterizing normal variations and heritability of structural connectomes in healthy subjects.

Exome sequencing reveals pathogenic mutations in 91 strains of mice with Mendelian disorders.

Spontaneously arising mouse mutations have served as the foundation for understanding gene function for more than 100 years. We have used exome sequencing in an effort to identify the causative mutations for 172 distinct, spontaneously arising mouse models of Mendelian disorders, including a broad range of clinically relevant phenotypes. To analyze the resulting data, we developed an analytics pipeline that is optimized for mouse exome data and a variation database that allows for reproducible, user-defined data mining as well as nomination of mutation candidates through knowledge-based integration of sample and variant data. Using these new tools, putative pathogenic mutations were identified for 91 (53%) of the strains in our study. Despite the increased power offered by potentially unlimited pedigrees and controlled breeding, about half of our exome cases remained unsolved. Using a combination of manual analyses of exome alignments and whole-genome sequencing, we provide evidence that a large fraction of unsolved exome cases have underlying structural mutations. This result directly informs efforts to investigate the similar proportion of apparently Mendelian human phenotypes that are recalcitrant to exome sequencing.

Coincidental loss of DOCK8 function in NLRP10-deficient and C3H/HeJ mice results in defective dendritic cell migration.

Dendritic cells (DCs) are the primary leukocytes responsible for priming T cells. To find and activate naïve T cells, DCs must migrate to lymph nodes, yet the cellular programs responsible for this key step remain unclear. DC migration to lymph nodes and the subsequent T-cell response are disrupted in a mouse we recently described lacking the NOD-like receptor NLRP10 (NLR family, pyrin domain containing 10); however, the mechanism by which this pattern recognition receptor governs DC migration remained unknown. Using a proteomic approach, we discovered that DCs from Nlrp10 knockout mice lack the guanine nucleotide exchange factor DOCK8 (dedicator of cytokinesis 8), which regulates cytoskeleton dynamics in multiple leukocyte populations; in humans, loss-of-function mutations in Dock8 result in severe immunodeficiency. Surprisingly, Nlrp10 knockout mice crossed to other backgrounds had normal DOCK8 expression. This suggested that the original Nlrp10 knockout strain harbored an unexpected mutation in Dock8, which was confirmed using whole-exome sequencing. Consistent with our original report, NLRP3 inflammasome activation remained unaltered in NLRP10-deficient DCs even after restoring DOCK8 function; however, these DCs recovered the ability to migrate. Isolated loss of DOCK8 via targeted deletion confirmed its absolute requirement for DC migration. Because mutations in Dock genes have been discovered in other mouse lines, we analyzed the diversity of Dock8 across different murine strains and found that C3H/HeJ mice also harbor a Dock8 mutation that partially impairs DC migration. We conclude that DOCK8 is an important regulator of DC migration during an immune response and is prone to mutations that disrupt its crucial function.

RASS: a web server for RNA alignment in the joint sequence-structure space.

Comparison of ribonucleic acid (RNA) molecules is important for revealing their evolutionary relationships, predicting their functions and predicting their structures. Many methods have been developed for comparing RNAs using either sequence or three-dimensional (3D) structure (backbone geometry) information. Sequences and 3D structures contain non-overlapping sets of information that both determine RNA functions. When comparing RNA 3D structures, both types of information need to be taken into account. However, few methods compare RNA structures using both sequence and 3D structure information. Recently, we have developed a new method based on elastic shape analysis (ESA) that compares RNA molecules by combining both sequence and 3D structure information. ESA treats RNA structures as 3D curves with sequence information encoded on additional coordinates so that the alignment can be performed in the joint sequence-structure space. The similarity between two RNA molecules is quantified by a formal distance, geodesic distance. In this study, we implement a web server for the method, called RASS, to make it publicly available to research community. The web server is located at http://cloud.stat.fsu.edu/RASS/.

An ENU-induced mutation in the mouse Rnf212 gene is associated with male meiotic failure and infertility.

The ENU-induced repro57 mutation was identified in an unbiased screen for the discovery of novel genes for fertility. Male repro57 homozygous mice are infertile and exhibit significantly reduced testis weight compared with WT mice. Histological examination of mutant testes revealed that spermatocytes degenerated during late prophase, and no mature spermatozoa were found in the seminiferous epithelium, suggesting that infertility is caused by the arrest of spermatogenesis at late meiotic prophase. Consistent with this hypothesis, the number of foci with MLH1, a protein essential for crossing over, is greatly reduced in repro57 mutant spermatocytes, which also lack chiasmata between homologs and exhibit premature dissociation of XY chromosomes. In repro57 mutant mice, we identified a mutation in the Rnf212 gene, encoding Ring finger protein 212. The overall phenotype of repro57 mice is consistent with the recently reported phenotype of the Rnf212 knockout mice; slight differences may be due to genetic background effects. Thus, the repro57 nonsense mutation provides a new allele of the mouse Rnf212 gene.

Development and validation of the JAX Cancer Treatment Profile™ for detection of clinically actionable mutations in solid tumors.

BACKGROUND: The continued development of targeted therapeutics for cancer treatment has required the concomitant development of more expansive methods for the molecular profiling of the patient's tumor. We describe the validation of the JAX Cancer Treatment Profile™ (JAX-CTP™), a next generation sequencing (NGS)-based molecular diagnostic assay that detects actionable mutations in solid tumors to inform the selection of targeted therapeutics for cancer treatment. METHODS: NGS libraries are generated from DNA extracted from formalin fixed paraffin embedded tumors. Using hybrid capture, the genes of interest are enriched and sequenced on the Illumina HiSeq 2500 or MiSeq sequencers followed by variant detection and functional and clinical annotation for the generation of a clinical report. RESULTS: The JAX-CTP™ detects actionable variants, in the form of single nucleotide variations and small insertions and deletions (≤50 bp) in 190 genes in specimens with a neoplastic cell content of ≥10%. The JAX-CTP™ is also validated for the detection of clinically actionable gene amplifications. CONCLUSIONS: There is a lack of consensus in the molecular diagnostics field on the best method for the validation of NGS-based assays in oncology, thus the importance of communicating methods, as contained in this report. The growing number of targeted therapeutics and the complexity of the tumor genome necessitate continued development and refinement of advanced assays for tumor profiling to enable precision cancer treatment.

RNA global alignment in the joint sequence-structure space using elastic shape analysis.

The functions of RNAs, like proteins, are determined by their structures, which, in turn, are determined by their sequences. Comparison/alignment of RNA molecules provides an effective means to predict their functions and understand their evolutionary relationships. For RNA sequence alignment, most methods developed for protein and DNA sequence alignment can be directly applied. RNA 3-dimensional structure alignment, on the other hand, tends to be more difficult than protein structure alignment due to the lack of regular secondary structures as observed in proteins. Most of the existing RNA 3D structure alignment methods use only the backbone geometry and ignore the sequence information. Using both the sequence and backbone geometry information in RNA alignment may not only produce more accurate classification, but also deepen our understanding of the sequence-structure-function relationship of RNA molecules. In this study, we developed a new RNA alignment method based on elastic shape analysis (ESA). ESA treats RNA structures as three dimensional curves with sequence information encoded on additional dimensions so that the alignment can be performed in the joint sequence-structure space. The similarity between two RNA molecules is quantified by a formal distance, geodesic distance. Based on ESA, a rigorous mathematical framework can be built for RNA structure comparison. Means and covariances of full structures can be defined and computed, and probability distributions on spaces of such structures can be constructed for a group of RNAs. Our method was further applied to predict functions of RNA molecules and showed superior performance compared with previous methods when tested on benchmark datasets. The programs are available at http://stat.fsu.edu/ ∼jinfeng/ESA.html.

Discovery of transgene insertion sites by high throughput sequencing of mate pair libraries.

BACKGROUND: Transgenesis by random integration of a transgene into the genome of a zygote has become a reliable and powerful method for the creation of new mouse strains that express exogenous genes, including human disease genes, tissue specific reporter genes or genes that allow for tissue specific recombination. Nearly 6,500 transgenic alleles have been created by random integration in embryos over the last 30 years, but for the vast majority of these strains, the transgene insertion sites remain uncharacterized. RESULTS: To obtain a complete understanding of how insertion sites might contribute to phenotypic outcomes, to more cost effectively manage transgenic strains, and to fully understand mechanisms of instability in transgene expression, we've developed methodology and a scoring scheme for transgene insertion site discovery using high throughput sequencing data. CONCLUSIONS: Similar to other molecular approaches to transgene insertion site discovery, high-throughput sequencing of standard paired-end libraries is hindered by low signal to noise ratios. This problem is exacerbated when the transgene consists of sequences that are also present in the host genome. We've found that high throughput sequencing data from mate-pair libraries are more informative when compared to data from standard paired end libraries. We also show examples of the genomic regions that harbor transgenes, which have in common a preponderance of repetitive sequences.

Assemblathon 1: a competitive assessment of de novo short read assembly methods.

Low-cost short read sequencing technology has revolutionized genomics, though it is only just becoming practical for the high-quality de novo assembly of a novel large genome. We describe the Assemblathon 1 competition, which aimed to comprehensively assess the state of the art in de novo assembly methods when applied to current sequencing technologies. In a collaborative effort, teams were asked to assemble a simulated Illumina HiSeq data set of an unknown, simulated diploid genome. A total of 41 assemblies from 17 different groups were received. Novel haplotype aware assessments of coverage, contiguity, structure, base calling, and copy number were made. We establish that within this benchmark: (1) It is possible to assemble the genome to a high level of coverage and accuracy, and that (2) large differences exist between the assemblies, suggesting room for further improvements in current methods. The simulated benchmark, including the correct answer, the assemblies, and the code that was used to evaluate the assemblies is now public and freely available from http://www.assemblathon.org/.

Landmark-free statistical analysis of the shape of plant leaves.

The shapes of plant leaves are important features to biologists, as they can help in distinguishing plant species, measuring their health, analyzing their growth patterns, and understanding relations between various species. Most of the methods that have been developed in the past focus on comparing the shape of individual leaves using either descriptors or finite sets of landmarks. However, descriptor-based representations are not invertible and thus it is often hard to map descriptor variability into shape variability. On the other hand, landmark-based techniques require automatic detection and registration of the landmarks, which is very challenging in the case of plant leaves that exhibit high variability within and across species. In this paper, we propose a statistical model based on the Squared Root Velocity Function (SRVF) representation and the Riemannian elastic metric of Srivastava et al. (2011) to model the observed continuous variability in the shape of plant leaves. We treat plant species as random variables on a non-linear shape manifold and thus statistical summaries, such as means and covariances, can be computed. One can then study the principal modes of variations and characterize the observed shapes using probability density models, such as Gaussians or Mixture of Gaussians. We demonstrate the usage of such statistical model for (1) efficient classification of individual leaves, (2) the exploration of the space of plant leaf shapes, which is important in the study of population-specific variations, and (3) comparing entire plant species, which is fundamental to the study of evolutionary relationships in plants. Our approach does not require descriptors or landmarks but automatically solves for the optimal registration that aligns a pair of shapes. We evaluate the performance of the proposed framework on publicly available benchmarks such as the Flavia, the Swedish, and the ImageCLEF2011 plant leaf datasets.

Elastic Shape Analysis of Tree-Like 3D Objects Using Extended SRVF Representation.

How can one analyze detailed 3D biological objects, such as neuronal and botanical trees, that exhibit complex geometrical and topological variation? In this paper, we develop a novel mathematical framework for representing, comparing, and computing geodesic deformations between the shapes of such tree-like 3D objects. A hierarchical organization of subtrees characterizes these objects - each subtree has a main branch with some side branches attached - and one needs to match these structures across objects for meaningful comparisons. We propose a novel representation that extends the Square-Root Velocity Function (SRVF), initially developed for Euclidean curves, to tree-shaped 3D objects. We then define a new metric that quantifies the bending, stretching, and branch sliding needed to deform one tree-shaped object into the other. Compared to the current metrics such as the Quotient Euclidean Distance (QED) and the Tree Edit Distance (TED), the proposed representation and metric capture the full elasticity of the branches (i.e., bending and stretching) as well as the topological variations (i.e., branch death/birth and sliding). It completely avoids the shrinkage that results from the edge collapse and node split operations of the QED and TED metrics. We demonstrate the utility of this framework in comparing, matching, and computing geodesics between biological objects such as neuronal and botanical trees. We also demonstrate its application to various shape analysis tasks such as (i) symmetry analysis and symmetrization of tree-shaped 3D objects, (ii) computing summary statistics (means and modes of variations) of populations of tree-shaped 3D objects, (iii) fitting parametric probability distributions to such populations, and (iv) finally synthesizing novel tree-shaped 3D objects through random sampling from estimated probability distributions.

Statistical Analysis of Complex Shape Graphs.

This paper provides developments in statistical shape analysis of shape graphs, and demonstrates them using such complex objects as Retinal Blood Vessel (RBV) networks and neurons. The shape graphs are represented by sets of nodes and edges (articulated curves) connecting some nodes. The goals are to utilize nodes (locations, connectivity) and edges (edge weights and shapes) to: (1) characterize shapes, (2) quantify shape differences, and (3) model statistical variability. We develop a mathematical representation, elastic Riemannian metrics, and associated tools for shape graphs. Specifically, we derive tools for shape graph registration, geodesics, statistical summaries, shape modeling, and shape synthesis. Geodesics are convenient for visualizing optimal deformations, and PCA helps in dimension reduction and statistical modeling. One key challenge in comparing shape graphs with vastly different complexities (in number of nodes and edges). This paper introduces a novel multi-scale representation to handle this challenge. Using the notions of (1) "effective resistance" to cluster nodes and (2) elastic shape averaging of edge curves, it reduces graph complexity while retaining overall structures. This allows shape comparisons by bringing graphs to similar complexities. We demonstrate these ideas on 2D RBV networks from the STARE and DRIVE databases and 3D neurons from the NeuroMorpho database.

Nextflow pipeline for Visium and H&E data from patient-derived xenograft samples.

We designed a Nextflow DSL2-based pipeline, Spatial Transcriptomics Quantification (STQ), for simultaneous processing of 10x Genomics Visium spatial transcriptomics data and a matched hematoxylin and eosin (H&E)-stained whole-slide image (WSI), optimized for patient-derived xenograft (PDX) cancer specimens. Our pipeline enables the classification of sequenced transcripts for deconvolving the mouse and human species and mapping the transcripts to reference transcriptomes. We align the H&E WSI with the spatial layout of the Visium slide and generate imaging and quantitative morphology features for each Visium spot. The pipeline design enables multiple analysis workflows, including single or dual reference genome input and stand-alone image analysis. We show the utility of our pipeline on a dataset from Visium profiling of four melanoma PDX samples. The clustering of Visium spots and clustering of H&E imaging features reveal similar patterns arising from the two data modalities.

Assessment of Patient-Derived Xenograft Growth and Antitumor Activity: The NCI PDXNet Consensus Recommendations.

Although patient-derived xenografts (PDX) are commonly used for preclinical modeling in cancer research, a standard approach to in vivo tumor growth analysis and assessment of antitumor activity is lacking, complicating the comparison of different studies and determination of whether a PDX experiment has produced evidence needed to consider a new therapy promising. We present consensus recommendations for assessment of PDX growth and antitumor activity, providing public access to a suite of tools for in vivo growth analyses. We expect that harmonizing PDX study design and analysis and assessing a suite of analytical tools will enhance information exchange and facilitate identification of promising novel therapies and biomarkers for guiding cancer therapy.

Estimating summary statistics in the spike-train space.

Estimating sample averages and sample variability is important in analyzing neural spike trains data in computational neuroscience. Current approaches have focused on advancing the use of parametric or semiparametric probability models of the underlying stochastic process, where the probabilistic distribution is characterized at each time point with basic statistics such as mean and variance. To directly capture and analyze the average and variability in the observation space of the spike trains, we focus on a data-driven approach where statistics are defined and computed in a function space in which the spike trains are viewed as individual points. Based on the definition of a "Euclidean" metric, a recent paper introduced the notion of the mean of a set of spike trains and developed an efficient algorithm to compute it under some restrictive conditions. Here we extend this study by: (1) developing a novel algorithm for mean computation that is quite general, and (2) introducing a notion of covariance of a set of spike trains. Specifically, we estimate the covariance matrix using the geometry of the warping functions that map the mean spike train to each of the spike trains in the dataset. Results from simulations as well as a neural recording in primate motor cortex indicate that the proposed mean and covariance successfully capture the observed variability in spike trains. In addition, a "Gaussian-type" probability model (defined using the estimated mean and covariance) reasonably characterizes the distribution of the spike trains and achieves a desirable performance in the classification of the spike trains.

Elastic statistical analysis of interval-valued time series.

We investigate the problem of statistical analysis of interval-valued time series data - two nonintersecting real-valued functions, representing lower and upper limits, over a period of time. Specifically, we pay attention to the two concepts of phase (or horizontal) variability and amplitude (or vertical) variability, and propose a phase-amplitude separation method. We view interval-valued time series as elements of a function (Hilbert) space and impose a Riemannian structure on it. We separate phase and amplitude variability in observed interval functions using a metric-based alignment solution. The key idea is to map an interval to a point in R 2 , view interval-valued time series as parameterized curves in R 2 , and borrow ideas from elastic shape analysis of planar curves, including PCA, to perform registration, summarization, analysis, and modeling of multiple series. The proposed phase-amplitude separation provides a new way of PCA and modeling for interval-valued time series, and enables shape clustering of interval-valued time series. We apply this framework to three different applications, including finance, meteorology and physiology, proves the effectiveness of proposed methods, and discovers some underlying patterns in the data. Experimental results on simulated data show that our method applies to the point-valued time series.

ACDA: implementation of an augmented drug synergy prediction algorithm.

Motivation: Drug synergy prediction is approached with machine learning techniques using molecular and pharmacological data. The published Cancer Drug Atlas (CDA) predicts a synergy outcome in cell-line models from drug target information, gene mutations and the models' monotherapy drug sensitivity. We observed low performance of the CDA, 0.339, measured by Pearson correlation of predicted versus measured sensitivity on DrugComb datasets. Results: We augmented the approach CDA by applying a random forest regression and optimization via cross-validation hyper-parameter tuning and named it Augmented CDA (ACDA). We benchmarked the ACDA's performance, which is 68% higher than that of the CDA when trained and validated on the same dataset spanning 10 tissues. We compared the performance of ACDA to one of the winning methods of the DREAM Drug Combination Prediction Challenge, the performance of which was lower than ACDA in 16 out of 19 cases. We further trained the ACDA on Novartis Institutes for BioMedical Research PDX encyclopedia data and generated sensitivity predictions for PDX models. Finally, we developed a novel approach to visualize synergy-prediction data. Availability and implementation: The source code is available at https://github.com/TheJacksonLaboratory/drug-synergy and the software package at PyPI. Supplementary information: Supplementary data are available at Bioinformatics Advances online.