The Monarch Initiative in 2024: an analytic platform integrating phenotypes, genes and diseases across species.

Bridging the gap between genetic variations, environmental determinants, and phenotypic outcomes is critical for supporting clinical diagnosis and understanding mechanisms of diseases. It requires integrating open data at a global scale. The Monarch Initiative advances these goals by developing open ontologies, semantic data models, and knowledge graphs for translational research. The Monarch App is an integrated platform combining data about genes, phenotypes, and diseases across species. Monarch's APIs enable access to carefully curated datasets and advanced analysis tools that support the understanding and diagnosis of disease for diverse applications such as variant prioritization, deep phenotyping, and patient profile-matching. We have migrated our system into a scalable, cloud-based infrastructure; simplified Monarch's data ingestion and knowledge graph integration systems; enhanced data mapping and integration standards; and developed a new user interface with novel search and graph navigation features. Furthermore, we advanced Monarch's analytic tools by developing a customized plugin for OpenAI's ChatGPT to increase the reliability of its responses about phenotypic data, allowing us to interrogate the knowledge in the Monarch graph using state-of-the-art Large Language Models. The resources of the Monarch Initiative can be found at monarchinitiative.org and its corresponding code repository at github.com/monarch-initiative/monarch-app.

Characterizing Intra-Tumor and Inter-Tumor Variability of Immune Cell Infiltrates in Murine Syngeneic Tumors.

Murine tumors are indispensable model systems in preclinical immuno-oncology research. While immunologic heterogeneity is well-known to be an important factor that can influence treatment outcome, there is a severe paucity of data concerning the nature of this heterogeneity in murine tumor models. Using serial sectioning methodology combined with IHC analysis and whole-slide image analysis, the depth-dependent variation in immune-cell abundance in tumor specimens was investigated at single-cell resolution. Specifically, intra- and intertumor variability in cell density of nine immune-cell biomarkers was quantified in multiple murine tumor models. The analysis showed that intertumor variability was typically the dominant source of variation in measurements of immune-cell densities. Statistical power analysis revealed the effect of group size and variance in immune-cell density on the predictive power of detecting a statistically meaningful fold-change in immune-cell density. Intertumor variability in the ratio of immune-cell densities showed distinct patterns in select tumor models and revealed the existence of strong correlations between select biomarker pairs. Furthermore, the relative proportion of immune cells at different depths across tumor samples was preserved in some but not all tumor models, thereby revealing the existence of compositional heterogeneity. Taken together, these results reveal novel insights into the nature of immunologic heterogeneity, which is not accessible through typical omics approaches.

Targeting Sagebrush (Artemisia Spp.) Restoration Following Wildfire with Greater Sage-Grouse (Centrocercus Urophasianus) Nest Selection and Survival Models.

Unprecedented conservation efforts for sagebrush (Artemisia spp.) ecosystems across the western United States have been catalyzed by risks from escalated wildfire activity that reduces habitat for sagebrush-obligate species such as Greater Sage-Grouse (Centrocercus urophasianus). However, post-fire restoration is challenged by spatial variation in ecosystem processes influencing resilience to disturbance and resistance to non-native invasive species, and spatial and temporal lags between slower sagebrush recovery processes and faster demographic responses of sage-grouse to loss of important habitat. Decision-support frameworks that account for these factors can help users strategically apply restoration efforts by predicting short and long-term ecological benefits of actions. Here, we developed a framework that strategically targets burned areas for restoration actions (e.g., seeding or planting sagebrush) that have the greatest potential to positively benefit sage-grouse populations through time. Specifically, we estimated sagebrush recovery following wildfire and risk of non-native annual grass invasion under four scenarios: passive recovery, grazing exclusion, active restoration with seeding, and active restoration with seedling transplants. We then applied spatial predictions of integrated nest site selection and survival models before wildfire, immediately following wildfire, and at 30 and 50 years post-wildfire based on each restoration scenario and measured changes in habitat. Application of this framework coupled with strategic planting designs aimed at developing patches of nesting habitat may help increase operational resilience for fire-impacted sagebrush ecosystems.

An orthogonal methods assessment of topical drug concentrations in skin and the impact for risk assessment in the viable epidermis.

Systemic toxicity assessments for oral or parenteral drugs often utilize the concentration of drug in plasma to enable safety margin calculations for human risk assessment. For topical drugs, there is no standard method for measuring drug concentrations in the stratum basale of the viable epidermis. This is particularly important since the superficial part of the epidermis, the stratum corneum (SC), is nonviable and where most of a topically applied drug remains, never penetrating deeper into the skin. We investigated the relative concentrations of a prototype kinase inhibitor using punch biopsy, laser capture microdissection, and imaging mass spectrometry methods in the SC, stratum basale, and dermis of minipig skin following topical application as a cream formulation. The results highlight the value of laser capture microdissection and mass spectrometry imaging in quantifying the large difference in drug concentration across the skin and even within the epidermis, and supports use of these methods for threshold-based toxicity risk assessments in specific anatomic locations of the skin, like of the stratum basale.

Wildfire and the ecological niche: Diminishing habitat suitability for an indicator species within semi-arid ecosystems.

Globally accelerating frequency and extent of wildfire threatens the persistence of specialist wildlife species through direct loss of habitat and indirect facilitation of exotic invasive species. Habitat specialists may be especially prone to rapidly changing environmental conditions because their ability to adapt lags behind the rate of habitat alteration. As a result, these populations may become increasingly susceptible to ecological traps by returning to suboptimal breeding habitats that were dramatically altered by disturbance. We demonstrate a multistage modeling approach that integrates habitat selection and survival during the key nesting life-stage of a bird species of high conservation concern, the greater sage-grouse (Centrocercus urophasianus; hereafter, sage-grouse). We applied these spatially explicit models to a spatiotemporally robust dataset of sage-grouse nest locations and fates across wildfire-altered sagebrush ecosystems of the Great Basin ecoregion, western United States. Female sage-grouse exhibited intricate habitat selection patterns that varied across regional gradients of ecological productivity among sagebrush communities, but often selected nest sites that disproportionately resulted in nest failure. For example, 23% of nests occurred in wildfire-affected habitats characterized by reduced sagebrush cover and greater composition of invasive annual grasses. We found survival of nests was negatively associated with wildfire-affected areas, but positively associated with higher elevations with increased ruggedness and overall shrub cover. Strong site fidelity likely drove sage-grouse to continue nesting in habitats degraded by wildfire. Hence, increasing frequency and extent of wildfire may contribute disproportionately to reduced reproductive success by creating ecological traps that act as population sinks. Identifying such habitat mismatches between selection and survival facilitates deeper understanding of the mechanisms driving reduced geographic niche space and population decline at broad spatiotemporal scales, while guiding management actions to areas that would be most beneficial to the species.

Territoriality drives preemptive habitat selection in recovering wolves: Implications for carnivore conservation.

According to the ideal-free distribution (IFD), individuals within a population are free to select habitats that maximize their chances of success. Assuming knowledge of habitat quality, the IFD predicts that average fitness will be approximately equal among individuals and between habitats, while density varies, implying that habitat selection will be density dependent. Populations are often assumed to follow an IFD, although this assumption is rarely tested with empirical data, and may be incorrect when territoriality indicates habitat selection tactics that deviate from the IFD (e.g. ideal-despotic distribution or ideal-preemptive distribution). When territoriality influences habitat selection, species' density will not directly reflect components of fitness such as reproductive success or survival. In such cases, assuming an IFD can lead to false conclusions about habitat quality. We tested theoretical models of density-dependent habitat selection on a species known to exhibit territorial behaviour in order to determine whether commonly applied habitat models are appropriate under these circumstances. We combined long-term radiotelemetry and census data from grey wolves Canis lupus in the Upper Peninsula of Michigan, USA to relate spatiotemporal variability in wolf density to underlying classifications of habitat within a hierarchical state-space modelling framework. We then iteratively applied isodar analysis to evaluate which distribution of habitat selection best described this recolonizing wolf population. The wolf population in our study expanded by >1,000% during our study (~50 to >600 individuals), and density-dependent habitat selection was most consistent with the ideal-preemptive distribution, as opposed to the ideal-free or ideal-despotic alternatives. Population density of terrestrial carnivores may not be positively correlated with the fitness value of their habitats, and density-dependent habitat selection patterns may help to explain complex predator-prey dynamics and cascading indirect effects. Source-sink population dynamics appear likely when species exhibit rapid growth and occupy interspersed habitats of contrasting quality. These conditions are likely and have implications for large carnivores in many systems, such as areas in North America and Europe where large predator species are currently recolonizing their former ranges.

Comparative nonclinical assessments of the biosimilar PF-06410293 and originator adalimumab.

Adalimumab, a recombinant fully human monoclonal antibody targeting tumor necrosis factor (TNF), is approved in the United States and Europe to treat various inflammatory and autoimmune indications. Biosimilars are approved biologics highly similar, but not identical, to approved biotherapeutics. To support clinical development of PF-06410293, an adalimumab biosimilar, nonclinical studies evaluated the structural, functional, toxicologic, and toxicokinetic similarity to originator adalimumab sourced from the United States (adalimumab-US) and European Union (adalimumab-EU). Structural similarity was assessed by peptide mapping. Biologic activity was measured via inhibition of TNF-induced apoptosis and Fc-based functionality assessments. In vivo nonclinical similarity was evaluated in a toxicity study in cynomolgus monkeys administered subcutaneous PF-06410293 or adalimumab-EU (0 or 157 mg/kg/week). Peptide mapping demonstrated PF-06410293, adalimumab-US, and adalimumab-EU had identical amino acid sequences. Comparative functional and binding assessments were similar. Effects of PF-06410293 and adalimumab-EU were similar and limited to pharmacologically mediated decreased cellularity of lymphoid follicles and germinal centers in spleen. Toxicokinetics were similar; maximum plasma concentration and area-under-the-concentration-time curve ratio of PF-06410293:adalimumab-EU ranged from 1.0 to 1.2. These studies supported PF-06410293 entry into clinical development. Many regulatory agencies now only request nonclinical in vivo testing if there is residual uncertainty regarding biosimilarity after in vitro analytical studies.

Patterns of divergence across the geographic and genomic landscape of a butterfly hybrid zone associated with a climatic gradient.

Hybrid zones are a valuable tool for studying the process of speciation and for identifying the genomic regions undergoing divergence and the ecological (extrinsic) and nonecological (intrinsic) factors involved. Here, we explored the genomic and geographic landscape of divergence in a hybrid zone between Papilio glaucus and Papilio canadensis. Using a genome scan of 28,417 ddRAD SNPs, we identified genomic regions under possible selection and examined their distribution in the context of previously identified candidate genes for ecological adaptations. We showed that differentiation was genomewide, including multiple candidate genes for ecological adaptations, particularly those involved in seasonal adaptation and host plant detoxification. The Z chromosome and four autosomes showed a disproportionate amount of differentiation, suggesting genes on these chromosomes play a potential role in reproductive isolation. Cline analyses of significantly differentiated genomic SNPs, and of species-diagnostic genetic markers, showed a high degree of geographic coincidence (81%) and concordance (80%) and were associated with the geographic distribution of a climate-mediated developmental threshold (length of the growing season). A relatively large proportion (1.3%) of the outliers for divergent selection were not associated with candidate genes for ecological adaptations and may reflect the presence of previously unrecognized intrinsic barriers between these species. These results suggest that exogenous (climate-mediated) and endogenous (unknown) clines may have become coupled and act together to reinforce reproductive isolation. This approach of assessing divergence across both the genomic and geographic landscape can provide insight about the interplay between the genetic architecture of reproductive isolation and endogenous and exogenous selection.

Tofacitinib attenuates pathologic immune pathways in patients with psoriasis: A randomized phase 2 study.

BACKGROUND: Tofacitinib is an oral Janus kinase inhibitor being investigated for psoriasis. OBJECTIVE: We sought to elucidate the molecular mechanisms underlying the clinical efficacy of tofacitinib in patients with psoriasis. METHODS: Twelve patients with plaque psoriasis were randomized (3:1) to receive 10 mg of tofacitinib or placebo twice daily for 12 weeks. Biopsy specimens were taken from nonlesional (baseline) and lesional (baseline, days 1 and 3, and weeks 1, 2, 4, and 12) skin. Biopsy specimens were examined for psoriatic epidermal features (thickness, Ki67(+) keratinocytes and keratin 16 [KRT16] mRNA expression, and phosphorylated signal transducer and activator of transcription [pSTAT](+) nuclei) and T-cell and dendritic cell (DC) subsets by using immunohistochemistry, and mRNA transcripts were quantified by using a microarray. RESULTS: In lesional skin keratinocyte pSTAT1 and pSTAT3 staining was increased at baseline but reduced after 1 day of tofacitinib (baseline, median of 1290 pSTAT1(+) cells/µm(2); day 1, median of 332 pSTAT1(+) cells/µm(2); and nonlesional, median of 155 pSTAT1(+) cells/µm(2)). Epidermal thickness and KRT16 mRNA expression were significantly and progressively reduced after days 1 and 3 of tofacitinib administration, respectively (eg, KRT16 decreased 2.74-fold, day 3 vs baseline, P = .016). Decreases in DC and T-cell numbers were observed after weeks 1 and 2, respectively. At week 4, significant decreases in IL-23/TH17 pathways were observed that persisted through week 12. Improvements in clinical and histologic features were strongly associated with changes in expression of psoriasis-related genes and reduction in IL-17 gene expression. CONCLUSIONS: Tofacitinib has a multitiered response in patients with psoriasis: (1) rapid attenuation of keratinocyte Janus kinase/STAT signaling; (2) removal of keratinocyte-induced cytokine signaling, leading to reductions in pathologic DC and T-cell numbers to nonlesional levels; and (3) inhibition of the IL-23/TH17 pathway.

Two Theobroma cacao genotypes with contrasting pathogen tolerance show aberrant transcriptional and ROS responses after salicylic acid treatment.

Understanding the genetic basis of pathogen susceptibility in various crop plants is crucial to increasing the stability of food, feed, and fuel production. Varietal differences in defence responses provide insights into the mechanisms of resistance and are a key resource for plant breeders. To explore the role of salicylic acid in the regulation of defence in cacao, we demonstrated that SA treatment decreased susceptibility to a pod rot pathogen, Phytophthora tropicalis in two genotypes, Scavina 6 and Imperial College Selection 1, which differ in their resistance to several agriculturally important pathogens. Transient overexpression of TcNPR1, a major transcriptional regulator of the SA-dependent plant immune system, also increased pathogen tolerance in cacao leaves. To explore further the genetic basis of resistance in cacao, we used microarrays to measure gene expression profiles after salicylic acid (SA) treatment in these two cacao genotypes. The two genotypes displayed distinct transcriptional responses to SA. Unexpectedly, the expression profile of the susceptible genotype ICS1 included a larger number of pathogenesis-related genes that were induced by SA at 24h after treatment, whereas genes encoding many chloroplast and mitochondrial proteins implicated in reactive oxygen species production were up-regulated in the resistant genotype, Sca6. Sca6 accumulated significantly more superoxide at 24h after treatment of leaves with SA. These experiments revealed critical insights regarding the molecular differences between cacao varieties, which will allow a better understanding of defence mechanisms to help guide breeding programmes.

Identification of tubular injury microRNA biomarkers in urine: comparison of next-generation sequencing and qPCR-based profiling platforms.

BACKGROUND: MicroRNAs (miRNAs) are small, non-coding RNAs that regulate protein levels post-transcriptionally. miRNAs play important regulatory roles in many cellular processes and have been implicated in several diseases. Recent studies have reported significant levels of miRNAs in a variety of body fluids, raising the possibility that miRNAs could serve as useful biomarkers. Next-generation sequencing (NGS) is increasingly employed in biomedical investigations. Although concordance between this platform and qRT-PCR based assays has been reported in high quality specimens, information is lacking on comparisons in biofluids especially urine. Here we describe the changes in miRNA expression patterns in a rodent model of renal tubular injury (gentamicin). Our aim is to compare RNA sequencing and qPCR based miRNA profiling in urine specimen from control and rats with confirmed tubular injury. RESULTS: Our preliminary examination of the concordance between miRNA-seq and qRT-PCR in urine specimen suggests minimal agreement between platforms probably due to the differences in sensitivity. Our results suggest that although miRNA-seq has superior specificity, it may not detect low abundant miRNAs in urine samples. Specifically, miRNA-seq did not detect some sequences which were identified by qRT-PCR. On the other hand, the qRT-PCR analysis was not able to detect the miRNA isoforms, which made up the majority of miRNA changes detected by NGS. CONCLUSIONS: To our knowledge, this is the first time that miRNA profiling platforms including NGS have been compared in urine specimen. miRNAs identified by both platforms, let-7d, miR-203, and miR-320, may potentially serve as promising novel urinary biomarkers for drug induced renal tubular epithelial injury.

Gene expression in closely related species mirrors local adaptation: consequences for responses to a warming world.

Local adaptation of populations could preclude or slow range expansions in response to changing climate, particularly when dispersal is limited. To investigate the differential responses of populations to changing climatic conditions, we exposed poleward peripheral and central populations of two Lepidoptera to reciprocal, common-garden climatic conditions and compared their whole-transcriptome expression. We found evidence of simple population differentiation in both species, and in the species with previously identified population structure and phenotypic local adaptation, we found several hundred genes that responded in a synchronized and localized fashion. These genes were primarily involved in energy metabolism and oxidative stress, and expression levels were most divergent between populations in the same environment in which we previously detected divergence for metabolism. We found no localized genes in the species with less population structure and for which no local adaptation was previously detected. These results challenge the assumption that species are functionally similar across their ranges and poleward peripheral populations are preadapted to warmer conditions. Rather, some taxa deserve population-level consideration when predicting the effects of climate change because they respond in genetically based, distinctive ways to changing conditions.

A comparative study of ripening among berries of the grape cluster reveals an altered transcriptional programme and enhanced ripening rate in delayed berries.

Transcriptional studies in relation to fruit ripening generally aim to identify the transcriptional states associated with physiological ripening stages and the transcriptional changes between stages within the ripening programme. In non-climacteric fruits such as grape, all ripening-related genes involved in this programme have not been identified, mainly due to the lack of mutants for comparative transcriptomic studies. A feature in grape cluster ripening (Vitis vinifera cv. Pinot noir), where all berries do not initiate the ripening at the same time, was exploited to study their shifted ripening programmes in parallel. Berries that showed marked ripening state differences in a véraison-stage cluster (ripening onset) ultimately reached similar ripeness states toward maturity, indicating the flexibility of the ripening programme. The expression variance between these véraison-stage berry classes, where 11% of the genes were found to be differentially expressed, was reduced significantly toward maturity, resulting in the synchronization of their transcriptional states. Defined quantitative expression changes (transcriptional distances) not only existed between the véraison transitional stages, but also between the véraison to maturity stages, regardless of the berry class. It was observed that lagging berries complete their transcriptional programme in a shorter time through altered gene expressions and ripening-related hormone dynamics, and enhance the rate of physiological ripening progression. Finally, the reduction in expression variance of genes can identify new genes directly associated with ripening and also assess the relevance of gene activity to the phase of the ripening programme.

Mechanistic investigations of test article-induced pancreatic toxicity at the endocrine-exocrine interface in the rat.

Pancreatic toxicity commonly affects the endocrine or exocrine pancreas. However, it can also occur at the endocrine-exocrine interface (EEI), where the capillary network of the islet merges with the capillaries of the surrounding acinar tissue, that is, the insulo-acinar portal system. The goal of this article is to describe a novel, test article-induced pancreatic toxicity that originated at the EEI and to summarize investigations into the mechanistic basis of the injury. This injury was initially characterized by light microscopy in 7/14 day-toxicity studies in Sprague-Dawley (Crl: CD®[SD]) rats with undisclosed test articles. Microvascular injury at the interface resulted in peri-islet serum exudation, fibrin deposition, hemorrhage, inflammation, and secondary degeneration/necrosis of surrounding exocrine tissue. More chronic injury presented as islet fibrosis and lobular atrophy. Direct cytotoxicity affecting the capillary endothelium at the EEI was confirmed ultrastructurally on day 4. Endothelial microparticle and blood flow studies further confirmed endothelial involvement. Similar lesions occurred less frequently in 2 other rat strains and not in the mouse, dog, or cynomolgus macaque. In summary, in vivo and investigative study data confirmed primary endothelial cytotoxicity in the pathogenesis of this lesion and suggested that the lesion may be rat/rat strain-specific and of uncertain relevance for human safety risk assessment.

Quantitative performance assessment of Ultivue multiplex panels in formalin-fixed, paraffin-embedded human and murine tumor specimens.

We present a rigorous validation strategy to evaluate the performance of Ultivue multiplex immunofluorescence panels. We have quantified the accuracy and precision of four different multiplex panels (three human and one mouse) in tumor specimens with varying levels of T cell density. Our results show that Ultivue panels are typically accurate wherein the relative difference in cell proportion between a multiplex image and a 1-plex image is less than 20% for a given biomarker. Ultivue panels exhibited relatively high intra-run precision (CV ≤ 25%) and relatively low inter-run precision (CV >> 25%) which can be remedied by using local intensity thresholding to gate biomarker positivity. We also evaluated the reproducibility of cell-cell distance estimates measured from multiplex images which show high intra- and inter-run precision. We introduce a new metric, multiplex labeling efficiency, which can be used to benchmark the overall fidelity of the multiplex data across multiple batch runs. Taken together our results provide a comprehensive characterization of Ultivue panels and offer practical guidelines for analyzing multiplex images.

Harnessing Consumer Wearable Digital Biomarkers for Individualized Recognition of Postpartum Depression Using the All of Us Research Program Data Set: Cross-Sectional Study.

BACKGROUND: Postpartum depression (PPD) poses a significant maternal health challenge. The current approach to detecting PPD relies on in-person postpartum visits, which contributes to underdiagnosis. Furthermore, recognizing PPD symptoms can be challenging. Therefore, we explored the potential of using digital biomarkers from consumer wearables for PPD recognition. OBJECTIVE: The main goal of this study was to showcase the viability of using machine learning (ML) and digital biomarkers related to heart rate, physical activity, and energy expenditure derived from consumer-grade wearables for the recognition of PPD. METHODS: Using the All of Us Research Program Registered Tier v6 data set, we performed computational phenotyping of women with and without PPD following childbirth. Intraindividual ML models were developed using digital biomarkers from Fitbit to discern between prepregnancy, pregnancy, postpartum without depression, and postpartum with depression (ie, PPD diagnosis) periods. Models were built using generalized linear models, random forest, support vector machine, and k-nearest neighbor algorithms and evaluated using the κ statistic and multiclass area under the receiver operating characteristic curve (mAUC) to determine the algorithm with the best performance. The specificity of our individualized ML approach was confirmed in a cohort of women who gave birth and did not experience PPD. Moreover, we assessed the impact of a previous history of depression on model performance. We determined the variable importance for predicting the PPD period using Shapley additive explanations and confirmed the results using a permutation approach. Finally, we compared our individualized ML methodology against a traditional cohort-based ML model for PPD recognition and compared model performance using sensitivity, specificity, precision, recall, and F1-score. RESULTS: Patient cohorts of women with valid Fitbit data who gave birth included <20 with PPD and 39 without PPD. Our results demonstrated that intraindividual models using digital biomarkers discerned among prepregnancy, pregnancy, postpartum without depression, and postpartum with depression (ie, PPD diagnosis) periods, with random forest (mAUC=0.85; κ=0.80) models outperforming generalized linear models (mAUC=0.82; κ=0.74), support vector machine (mAUC=0.75; κ=0.72), and k-nearest neighbor (mAUC=0.74; κ=0.62). Model performance decreased in women without PPD, illustrating the method's specificity. Previous depression history did not impact the efficacy of the model for PPD recognition. Moreover, we found that the most predictive biomarker of PPD was calories burned during the basal metabolic rate. Finally, individualized models surpassed the performance of a conventional cohort-based model for PPD detection. CONCLUSIONS: This research establishes consumer wearables as a promising tool for PPD identification and highlights personalized ML approaches, which could transform early disease detection strategies.

Finding Long-COVID: Temporal Topic Modeling of Electronic Health Records from the N3C and RECOVER Programs.

Post-Acute Sequelae of SARS-CoV-2 infection (PASC), also known as Long-COVID, encompasses a variety of complex and varied outcomes following COVID-19 infection that are still poorly understood. We clustered over 600 million condition diagnoses from 14 million patients available through the National COVID Cohort Collaborative (N3C), generating hundreds of highly detailed clinical phenotypes. Assessing patient clinical trajectories using these clusters allowed us to identify individual conditions and phenotypes strongly increased after acute infection. We found many conditions increased in COVID-19 patients compared to controls, and using a novel method to associate patients with clusters over time, we additionally found phenotypes specific to patient sex, age, wave of infection, and PASC diagnosis status. While many of these results reflect known PASC symptoms, the resolution provided by this unprecedented data scale suggests avenues for improved diagnostics and mechanistic understanding of this multifaceted disease.

Assessing De Novo transcriptome assembly metrics for consistency and utility.

BACKGROUND: Transcriptome sequencing and assembly represent a great resource for the study of non-model species, and many metrics have been used to evaluate and compare these assemblies. Unfortunately, it is still unclear which of these metrics accurately reflect assembly quality. RESULTS: We simulated sequencing transcripts of Drosophila melanogaster. By assembling these simulated reads using both a "perfect" and a modern transcriptome assembler while varying read length and sequencing depth, we evaluated quality metrics to determine whether they 1) revealed perfect assemblies to be of higher quality, and 2) revealed perfect assemblies to be more complete as data quantity increased.Several commonly used metrics were not consistent with these expectations, including average contig coverage and length, though they became consistent when singletons were included in the analysis. We found several annotation-based metrics to be consistent and informative, including contig reciprocal best hit count and contig unique annotation count. Finally, we evaluated a number of novel metrics such as reverse annotation count, contig collapse factor, and the ortholog hit ratio, discovering that each assess assembly quality in unique ways. CONCLUSIONS: Although much attention has been given to transcriptome assembly, little research has focused on determining how best to evaluate assemblies, particularly in light of the variety of options available for read length and sequencing depth. Our results provide an important review of these metrics and give researchers tools to produce the highest quality transcriptome assemblies.

Diffusion tensor imaging detects treatment effects of FTY720 in experimental autoimmune encephalomyelitis mice.

Fingolimod (FTY720) is an orally available sphingosine-1-phosphate (S1P) receptor modulator reducing relapse frequency in patients with relapsing-remitting multiple sclerosis (RRMS). In addition to immunosuppression, neuronal protection by FTY720 has also been suggested, but remains controversial. Axial and radial diffusivities derived from in vivo diffusion tensor imaging (DTI) were employed as noninvasive biomarkers of axonal injury and demyelination to assess axonal protection by FTY720 in experimental autoimmune encephalomyelitis (EAE) mice. EAE was induced through active immunization of C57BL/6 mice using myelin oligodendrocyte glycoprotein peptide 35-55 (MOG(35-55)). We evaluated both the prophylactic and therapeutic treatment effect of FTY720 at doses of 3 and 10 mg/kg on EAE mice by daily clinical scoring and end-point in vivo DTI. Prophylactic administration of FTY720 suppressed the disease onset and prevented axon and myelin damage when compared with EAE mice without treatment. Therapeutic treatment by FTY720 did not prevent EAE onset, but reduced disease severity, improving axial and radial diffusivity towards the control values without statistical significance. Consistent with previous findings, in vivo DTI-derived axial and radial diffusivity correlated with clinical scores in EAE mice. The results support the use of in vivo DTI as an effective outcome measure for preclinical drug development.

Early induction of polyfunctional simian immunodeficiency virus (SIV)-specific T lymphocytes and rapid disappearance of SIV from lymph nodes of sooty mangabeys during primary infection.

Although the cellular immune response is essential for controlling SIV replication in Asian macaques, its role in maintaining nonpathogenic SIV infection in natural hosts such as sooty mangabeys (SM) remains to be defined. We have previously shown that similar to rhesus macaques (RM), SM are able to mount a T lymphocyte response against SIV infection. To investigate early control of SIV replication in natural hosts, we performed a detailed characterization of SIV-specific cellular immunity and viral control in the first 6 mo following SIV infection in SM. Detection of the initial SIV-specific IFN-γ ELISPOT response in SIVsmE041-infected SM coincided temporally with a decline in peak plasma viremia and was similar in magnitude, specificity, and breadth to SIVsmE041-infected and SIVmac239-infected RM. Despite these similarities, SM showed a greater reduction in postpeak plasma viremia and a more rapid disappearance of productively SIV-infected cells from the lymph node compared with SIVmac239-infected RM. The early Gag-specific CD8(+) T lymphocyte response was significantly more polyfunctional in SM compared with RM, and granzyme B-positive CD8(+) T lymphocytes were present at significantly higher frequencies in SM even prior to SIV infection. These findings suggest that the early SIV-specific T cell response may be an important determinant of lymphoid tissue viral clearance and absence of lymph node immunopathology in natural hosts of SIV infection.