Enabling precision medicine via standard communication of HTS provenance, analysis, and results.

A personalized approach based on a patient's or pathogen's unique genomic sequence is the foundation of precision medicine. Genomic findings must be robust and reproducible, and experimental data capture should adhere to findable, accessible, interoperable, and reusable (FAIR) guiding principles. Moreover, effective precision medicine requires standardized reporting that extends beyond wet-lab procedures to computational methods. The BioCompute framework (https://w3id.org/biocompute/1.3.0) enables standardized reporting of genomic sequence data provenance, including provenance domain, usability domain, execution domain, verification kit, and error domain. This framework facilitates communication and promotes interoperability. Bioinformatics computation instances that employ the BioCompute framework are easily relayed, repeated if needed, and compared by scientists, regulators, test developers, and clinicians. Easing the burden of performing the aforementioned tasks greatly extends the range of practical application. Large clinical trials, precision medicine, and regulatory submissions require a set of agreed upon standards that ensures efficient communication and documentation of genomic analyses. The BioCompute paradigm and the resulting BioCompute Objects (BCOs) offer that standard and are freely accessible as a GitHub organization (https://github.com/biocompute-objects) following the "Open-Stand.org principles for collaborative open standards development." With high-throughput sequencing (HTS) studies communicated using a BCO, regulatory agencies (e.g., Food and Drug Administration [FDA]), diagnostic test developers, researchers, and clinicians can expand collaboration to drive innovation in precision medicine, potentially decreasing the time and cost associated with next-generation sequencing workflow exchange, reporting, and regulatory reviews.

DUF1220 dosage is linearly associated with increasing severity of the three primary symptoms of autism.

One of the three most frequently documented copy number variations associated with autism spectrum disorder (ASD) is a 1q21.1 duplication that encompasses sequences encoding DUF1220 protein domains, the dosage of which we previously implicated in increased human brain size. Further, individuals with ASD frequently display accelerated brain growth and a larger brain size that is also associated with increased symptom severity. Given these findings, we investigated the relationship between DUF1220 copy number and ASD severity, and here show that in individuals with ASD (n = 170), the copy number (dosage) of DUF1220 subtype CON1 is highly variable, ranging from 56 to 88 copies following a Gaussian distribution. More remarkably, in individuals with ASD CON1 copy number is also linearly associated, in a dose-response manner, with increased severity of each of the three primary symptoms of ASD: social deficits (p = 0.021), communicative impairments (p = 0.030), and repetitive behaviors (p = 0.047). These data indicate that DUF1220 protein domain (CON1) dosage has an ASD-wide effect and, as such, is likely to be a key component of a major pathway underlying ASD severity. Finally, these findings, by implicating the dosage of a previously unexamined, copy number polymorphic and brain evolution-related gene coding sequence in ASD severity, provide an important new direction for further research into the genetic factors underlying ASD.

DUF1220 copy number is linearly associated with increased cognitive function as measured by total IQ and mathematical aptitude scores.

DUF1220 protein domains exhibit the greatest human lineage-specific copy number expansion of any protein-coding sequence in the genome, and variation in DUF1220 copy number has been linked to both brain size in humans and brain evolution among primates. Given these findings, we examined associations between DUF1220 subtypes CON1 and CON2 and cognitive aptitude. We identified a linear association between CON2 copy number and cognitive function in two independent populations of European descent. In North American males, an increase in CON2 copy number corresponded with an increase in WISC IQ (R (2) = 0.13, p = 0.02), which may be driven by males aged 6-11 (R (2) = 0.42, p = 0.003). We utilized ddPCR in a subset as a confirmatory measurement. This group had 26-33 copies of CON2 with a mean of 29, and each copy increase of CON2 was associated with a 3.3-point increase in WISC IQ (R (2) = 0.22, p = 0.045). In individuals from New Zealand, an increase in CON2 copy number was associated with an increase in math aptitude ability (R (2) = 0.10 p = 0.018). These were not confounded by brain size. To our knowledge, this is the first study to report a replicated association between copy number of a gene coding sequence and cognitive aptitude. Remarkably, dosage variations involving DUF1220 sequences have now been linked to human brain expansion, autism severity and cognitive aptitude, suggesting that such processes may be genetically and mechanistically inter-related. The findings presented here warrant expanded investigations in larger, well-characterized cohorts.

Communicating computational workflows in a regulatory environment.

The volume of nucleic acid sequence data has exploded recently, amplifying the challenge of transforming data into meaningful information. Processing data can require an increasingly complex ecosystem of customized tools, which increases difficulty in communicating analyses in an understandable way yet is of sufficient detail to enable informed decisions or repeats. This can be of particular interest to institutions and companies communicating computations in a regulatory environment. BioCompute Objects (BCOs; an instance of pipeline documentation that conforms to the IEEE 2791-2020 standard) were developed as a standardized mechanism for analysis reporting. A suite of BCOs is presented, representing interconnected elements of a computation modeled after those that might be found in a regulatory submission but are shared publicly - in this case a pipeline designed to identify viral contaminants in biological manufacturing, such as for vaccines.

Olduvai domain expression downregulates mitochondrial pathways: implications for human brain evolution and neoteny.

Encoded by the NBPF gene family, Olduvai (formerly DUF1220) protein domains have undergone the largest human lineage-specific copy number expansion of any coding region in the genome. Olduvai copy number shows a linear relationship with several brain size-related measures and cortical neuron number among primates and with normal and disease-associated (micro- and macrocephaly) variation in brain size in human populations. While Olduvai domains have been shown to promote proliferation of neural stem cells, the mechanism underlying such effects has remained unclear. Here, we investigate the function of Olduvai by transcriptome and proteome analyses of cells overexpressing NBPF1 , a gene encoding 7 Olduvai domains. Our results from both RNAseq and mass spectrometry approaches suggest a potential downregulation of mitochondria. In our proteomics study, a Gene Ontology (GO) enrichment analysis for the downregulated proteins revealed a striking overrepresentation of the biological process related to the mitochondrial electron transport chain ( p value: 1.81e-11) and identified deregulation of the NADH dehydrogenase activity ( p value: 2.43e-11) as the primary molecular function. We verify the reduction of apparent mitochondria via live-cell imaging experiments. Given these and previous Olduvai findings, we suggest that the Olduvai-mediated, dosage-dependent reduction in available energy via mitochondrial downregulation may have resulted in a developmental slowdown such that the neurogenic window among primates, and most extremely in humans, was expanded over a greater time interval, allowing for production of greater numbers of neurons and a larger brain. We further suggest that such a slowdown may extend to other developmental processes that also exhibit neotenic features.

Communicating regulatory high-throughput sequencing data using BioCompute Objects.

This project demonstrates the use of the IEEE 2791-2020 Standard (BioCompute Objects [BCO]) to enable the complete and concise communication of results from next generation sequencing (NGS) analysis. One arm of a clinical trial was replicated using synthetically generated data made to resemble real biological data and then two independent analyses were performed. The first simulated a pharmaceutical regulatory submission to the US Food and Drug Administration (FDA) including analysis of results and a BCO. The second simulated an FDA review that included an independent analysis of the submitted data. Of the 118 simulated patient samples generated, 117 (99.15%) were in agreement in the two analyses. This process exemplifies how a template BCO (tBCO), including a verification kit, facilitates transparency and reproducibility, thereby reinforcing confidence in the regulatory submission process.

Bioinformatics tools developed to support BioCompute Objects.

Developments in high-throughput sequencing (HTS) result in an exponential increase in the amount of data generated by sequencing experiments, an increase in the complexity of bioinformatics analysis reporting and an increase in the types of data generated. These increases in volume, diversity and complexity of the data generated and their analysis expose the necessity of a structured and standardized reporting template. BioCompute Objects (BCOs) provide the requisite support for communication of HTS data analysis that includes support for workflow, as well as data, curation, accessibility and reproducibility of communication. BCOs standardize how researchers report provenance and the established verification and validation protocols used in workflows while also being robust enough to convey content integration or curation in knowledge bases. BCOs that encapsulate tools, platforms, datasets and workflows are FAIR (findable, accessible, interoperable and reusable) compliant. Providing operational workflow and data information facilitates interoperability between platforms and incorporation of future dataset within an HTS analysis for use within industrial, academic and regulatory settings. Cloud-based platforms, including High-performance Integrated Virtual Environment (HIVE), Cancer Genomics Cloud (CGC) and Galaxy, support BCO generation for users. Given the 100K+ userbase between these platforms, BioCompute can be leveraged for workflow documentation. In this paper, we report the availability of platform-dependent and platform-independent BCO tools: HIVE BCO App, CGC BCO App, Galaxy BCO API Extension and BCO Portal. Community engagement was utilized to evaluate tool efficacy. We demonstrate that these tools further advance BCO creation from text editing approaches used in earlier releases of the standard. Moreover, we demonstrate that integrating BCO generation within existing analysis platforms greatly streamlines BCO creation while capturing granular workflow details. We also demonstrate that the BCO tools described in the paper provide an approach to solve the long-standing challenge of standardizing workflow descriptions that are both human and machine readable while accommodating manual and automated curation with evidence tagging. Database URL:  https://www.biocomputeobject.org/resources.

OncoMX: A Knowledgebase for Exploring Cancer Biomarkers in the Context of Related Cancer and Healthy Data.

PURPOSE: The purpose of OncoMX1 knowledgebase development was to integrate cancer biomarker and relevant data types into a meta-portal, enabling the research of cancer biomarkers side by side with other pertinent multidimensional data types. METHODS: Cancer mutation, cancer differential expression, cancer expression specificity, healthy gene expression from human and mouse, literature mining for cancer mutation and cancer expression, and biomarker data were integrated, unified by relevant biomedical ontologies, and subjected to rule-based automated quality control before ingestion into the database. RESULTS: OncoMX provides integrated data encompassing more than 1,000 unique biomarker entries (939 from the Early Detection Research Network [EDRN] and 96 from the US Food and Drug Administration) mapped to 20,576 genes that have either mutation or differential expression in cancer. Sentences reporting mutation or differential expression in cancer were extracted from more than 40,000 publications, and healthy gene expression data with samples mapped to organs are available for both human genes and their mouse orthologs. CONCLUSION: OncoMX has prioritized user feedback as a means of guiding development priorities. By mapping to and integrating data from several cancer genomics resources, it is hoped that OncoMX will foster a dynamic engagement between bioinformaticians and cancer biomarker researchers. This engagement should culminate in a community resource that substantially improves the ability and efficiency of exploring cancer biomarker data and related multidimensional data.

Recent Advances in Systems and Network Medicine: Meeting Report from the First International Conference in Systems and Network Medicine.

The First International Conference in Systems and Network Medicine gathered together 200 global thought leaders, scientists, clinicians, academicians, industry and government experts, medical and graduate students, postdoctoral scholars and policymakers. Held at Georgetown University Conference Center in Washington D.C. on September 11-13, 2019, the event featured a day of pre-conference lectures and hands-on bioinformatic computational workshops followed by two days of deep and diverse scientific talks, panel discussions with eminent thought leaders, and scientific poster presentations. Topics ranged from: Systems and Network Medicine in Clinical Practice; the role of -omics technologies in Health Care; the role of Education and Ethics in Clinical Practice, Systems Thinking, and Rare Diseases; and the role of Artificial Intelligence in Medicine. The conference served as a unique nexus for interdisciplinary discovery and dialogue and fostered formation of new insights and possibilities for health care systems advances.

The case for DUF1220 domain dosage as a primary contributor to anthropoid brain expansion.

Here we present the hypothesis that increasing copy number (dosage) of sequences encoding DUF1220 protein domains is a major contributor to the evolutionary increase in brain size, neuron number, and cognitive capacity that is associated with the primate order. We further propose that this relationship is restricted to the anthropoid sub-order of primates, with DUF1220 copy number markedly increasing in monkeys, further in apes, and most extremely in humans where the greatest number of copies (~272 haploid copies) is found. We show that this increase closely parallels the increase in brain size and neuron number that has occurred among anthropoid primate species. We also provide evidence linking DUF1220 copy number to brain size within the human species, both in normal populations and in individuals associated with brain size pathologies (1q21-associated microcephaly and macrocephaly). While we believe these and other findings presented here strongly suggest increase in DUF1220 copy number is a key contributor to anthropoid brain expansion, the data currently available rely largely on correlative measures that, though considerable, do not yet provide direct evidence for a causal connection. Nevertheless, we believe the evidence presented is sufficient to provide the basis for a testable model which proposes that DUF1220 protein domain dosage increase is a main contributor to the increase in brain size and neuron number found among the anthropoid primate species and that is at its most extreme in human.

Mode of genetic inheritance modifies the association of head circumference and autism-related symptoms: a cross-sectional study.

BACKGROUND: Frequently individuals with autism spectrum disorder (ASD) have been noted with a larger head circumference (HC) than their typical developing peers. Biologic hypotheses suggest that an overly rapid brain growth leads to the core symptoms of ASD by impairing connectivity. Literature is divided however where deleterious, protective and null associations of HC with ASD symptoms in individuals with ASD have been found. METHOD: Individuals (n = 1,416) from the Autism Genetic Resource Exchange with ASD were examined for associations of HC with ASD like symptoms. Mixed models controlling for sex, age, race/ethnicity, simplex/multiplex status and accounting for correlations between siblings were used. Interactions by simplex/multiplex were explored. Adjustments for height in a sub-population with available data were explored as well. RESULTS: A Significant interaction term (p = 0.03) suggested that the effect of HC was dependent on whether the individual was simplex or multiplex. In simplex individuals at mean age (8.9 years) 1 cm increase in head circumference was associated with a 24% increase in the odds of a high social diagnostic score from the Autism Diagnostic Interview-Revised (odds ratio  = 1.24, p = 0.01). There was no association in multiplex individuals. Additionally, individuals classified with a non-verbal IQ <70 were 90% simplex and had a significantly increased head circumference (0.7 cm p = 0.03) relative to a mid-range non-verbal IQ group. Interestingly, children classified with a >110 non-verbal IQ also had an increased HC (0.4 cm p = 0.04), relative to a mid-range non-verbal IQ group, and were 90% multiplex. HC effects do not appear to be confounded by height, however, larger samples with height information are needed. CONCLUSION: The potential link between brain growth and autism like symptoms is complex and could depend on specific etiologies. Further investigations accounting for a likely mode of inheritance will help identify an ASD subtype related to HC.

Impact of cigarette smoke exposure on innate immunity: a Caenorhabditis elegans model.

BACKGROUND: Cigarette smoking is the major cause of chronic obstructive pulmonary disease (COPD) and lung cancer. Respiratory bacterial infections have been shown to be involved in the development of COPD along with impaired airway innate immunity. METHODOLOGY/PRINCIPAL FINDINGS: To address the in vivo impact of cigarette smoke (CS) exclusively on host innate defense mechanisms, we took advantage of Caenorhabditis elegans (C. elegans), which has an innate immune system but lacks adaptive immune function. Pseudomonas aeruginosa (PA) clearance from intestines of C. elegans was dampened by CS. Microarray analysis identified 6 candidate genes with a 2-fold or greater reduction after CS exposure, that have a human orthologue, and that may participate in innate immunity. To confirm a role of CS-down-regulated genes in the innate immune response to PA, RNA interference (RNAi) by feeding was carried out in C. elegans to inhibit the gene of interest, followed by PA infection to determine if the gene affected innate immunity. Inhibition of lbp-7, which encodes a lipid binding protein, resulted in increased levels of intestinal PA. Primary human bronchial epithelial cells were shown to express mRNA of human Fatty Acid Binding Protein 5 (FABP-5), the human orthologue of lpb-7. Interestingly, FABP-5 mRNA levels from human smokers with COPD were significantly lower (p = 0.036) than those from smokers without COPD. Furthermore, FABP-5 mRNA levels were up-regulated (7-fold) after bacterial (i.e., Mycoplasma pneumoniae) infection in primary human bronchial epithelial cell culture (air-liquid interface culture). CONCLUSIONS: Our results suggest that the C. elegans model offers a novel in vivo approach to specifically study innate immune deficiencies resulting from exposure to cigarette smoke, and that results from the nematode may provide insight into human airway epithelial cell biology and cigarette smoke exposure.