CyVerse: Cyberinfrastructure for open science.

CyVerse, the largest publicly-funded open-source research cyberinfrastructure for life sciences, has played a crucial role in advancing data-driven research since the 2010s. As the technology landscape evolved with the emergence of cloud computing platforms, machine learning and artificial intelligence (AI) applications, CyVerse has enabled access by providing interfaces, Software as a Service (SaaS), and cloud-native Infrastructure as Code (IaC) to leverage new technologies. CyVerse services enable researchers to integrate institutional and private computational resources, custom software, perform analyses, and publish data in accordance with open science principles. Over the past 13 years, CyVerse has registered more than 124,000 verified accounts from 160 countries and was used for over 1,600 peer-reviewed publications. Since 2011, 45,000 students and researchers have been trained to use CyVerse. The platform has been replicated and deployed in three countries outside the US, with additional private deployments on commercial clouds for US government agencies and multinational corporations. In this manuscript, we present a strategic blueprint for creating and managing SaaS cyberinfrastructure and IaC as free and open-source software.

Workflows for Rapid Functional Annotation of Diverse Arthropod Genomes.

Genome sequencing of a diverse array of arthropod genomes is already underway, and these genomes will be used to study human health, agriculture, biodiversity, and ecology. These new genomes are intended to serve as community resources and provide the foundational information required to apply 'omics technologies to a more diverse set of species. However, biologists require genome annotation to use these genomes and derive a better understanding of complex biological systems. Genome annotation incorporates two related, but distinct, processes: Demarcating genes and other elements present in genome sequences (structural annotation); and associating a function with genetic elements (functional annotation). While there are well-established and freely available workflows for structural annotation of gene identification in newly assembled genomes, workflows for providing the functional annotation required to support functional genomics studies are less well understood. Genome-scale functional annotation is required for functional modeling (enrichment, networks, etc.). A first-pass genome-wide functional annotation effort can rapidly identify under-represented gene sets for focused community annotation efforts. We present an open-source, open access, and containerized pipeline for genome-scale functional annotation of insect proteomes and apply it to various arthropod species. We show that the performance of the predictions is consistent across a set of arthropod genomes with varying assembly and annotation quality.

Chickspress: a resource for chicken gene expression.

High-throughput sequencing and proteomics technologies are markedly increasing the amount of RNA and peptide data that are available to researchers, which are typically made publicly available via data repositories such as the NCBI Sequence Read Archive and proteome archives, respectively. These data sets contain valuable information about when and where gene products are expressed, but this information is not readily obtainable from archived data sets. Here we report Chickspress (http://geneatlas.arl.arizona.edu), the first publicly available gene expression resource for chicken tissues. Since there is no single source of chicken gene models, Chickspress incorporates both NCBI and Ensembl gene models and links these gene sets with experimental gene expression data and QTL information. By linking gene models from both NCBI and Ensembl gene prediction pipelines, researchers can, for the first time, easily compare gene models from each of these prediction workflows to available experimental data for these products. We use Chickspress data to show the differences between these gene annotation pipelines. Chickspress also provides rapid search, visualization and download capacity for chicken gene sets based upon tissue type, developmental stage and experiment type. This first Chickspress release contains 161 gene expression data sets, including expression of mRNAs, miRNAs, proteins and peptides. We provide several examples demonstrating how researchers may use this resource.

Oxygen Perfusion (Persufflation) of Human Pancreata Enhances Insulin Secretion and Attenuates Islet Proinflammatory Signaling.

BACKGROUND: All human islets used in research and for the clinical treatment of diabetes are subject to ischemic damage during pancreas procurement, preservation, and islet isolation. A major factor influencing islet function is exposure of pancreata to cold ischemia during unavoidable windows of preservation by static cold storage (SCS). Improved preservation methods may prevent this functional deterioration. In the present study, we investigated whether pancreas preservation by gaseous oxygen perfusion (persufflation) better preserved islet function versus SCS. METHODS: Human pancreata were preserved by SCS or by persufflation in combination with SCS. Islets were subsequently isolated, and preparations in each group matched for SCS or total preservation time were compared using dynamic glucose-stimulated insulin secretion as a measure of ß-cell function and RNA sequencing to elucidate transcriptomic changes. RESULTS: Persufflated pancreata had reduced SCS time, which resulted in islets with higher glucose-stimulated insulin secretion compared to islets from SCS only pancreata. RNA sequencing of islets from persufflated pancreata identified reduced inflammatory and greater metabolic gene expression, consistent with expectations of reducing cold ischemic exposure. Portions of these transcriptional responses were not associated with time spent in SCS and were attributable to pancreatic reoxygenation. Furthermore, persufflation extended the total preservation time by 50% without any detectable decline in islet function or viability. CONCLUSIONS: These data demonstrate that pancreas preservation by persufflation rather than SCS before islet isolation reduces inflammatory responses and promotes metabolic pathways in human islets, which results in improved ß cell function.

In vitro characterization of neonatal, juvenile, and adult porcine islet oxygen demand, ß-cell function, and transcriptomes.

BACKGROUND: There is currently a shortage of human donor pancreata which limits the broad application of islet transplantation as a treatment for type 1 diabetes. Porcine islets have demonstrated potential as an alternative source, but a study evaluating islets from different donor ages under unified protocols has yet to be conducted. METHODS: Neonatal porcine islets (NPI; 1-3 days), juvenile porcine islets (JPI; 18-21 days), and adult porcine islets (API; 2+ years) were compared in vitro, including assessments of oxygen consumption rate, membrane integrity determined by FDA/PI staining, ß-cell proliferation, dynamic glucose-stimulated insulin secretion, and RNA sequencing. RESULTS: Oxygen consumption rate normalized to DNA was not significantly different between ages. Membrane integrity was age dependent, and API had the highest percentage of intact cells. API also had the highest glucose-stimulated insulin secretion response during a dynamic insulin secretion assay and had 50-fold higher total insulin content compared to NPI and JPI. NPI and JPI had similar glucose responsiveness, ß-cell percentage, and ß-cell proliferation rate. Transcriptome analysis was consistent with physiological assessments. API transcriptomes were enriched for cellular metabolic and insulin secretory pathways, while NPI exhibited higher expression of genes associated with proliferation. CONCLUSIONS: The oxygen demand, membrane integrity, ß-cell function and proliferation, and transcriptomes of islets from API, JPI, and NPI provide a comprehensive physiological comparison for future studies. These assessments will inform the optimal application of each age of porcine islet to expand the availability of islet transplantation.

Transcriptome response to heat stress in a chicken hepatocellular carcinoma cell line.

Heat stress triggers an evolutionarily conserved set of responses in cells. The transcriptome responds to hyperthermia by altering expression of genes to adapt the cell or organism to survive the heat challenge. RNA-seq technology allows rapid identification of environmentally responsive genes on a large scale. In this study, we have used RNA-seq to identify heat stress responsive genes in the chicken male white leghorn hepatocellular (LMH) cell line. The transcripts of 812 genes were responsive to heat stress (p < 0.01) with 235 genes upregulated and 577 downregulated following 2.5 h of heat stress. Among the upregulated were genes whose products function as chaperones, along with genes affecting collagen synthesis and deposition, transcription factors, chromatin remodelers, and genes modulating the WNT and TGF-beta pathways. Predominant among the downregulated genes were ones that affect DNA replication and repair along with chromosomal segregation. Many of the genes identified in this study have not been previously implicated in the heat stress response. These data extend our understanding of the transcriptome response to heat stress with many of the identified biological processes and pathways likely to function in adapting cells and organisms to hyperthermic stress. Furthermore, this study should provide important insight to future efforts attempting to improve species abilities to withstand heat stress through genome-wide association studies and breeding.

Transgenic soya bean seeds accumulating ß-carotene exhibit the collateral enhancements of oleate and protein content traits.

Transgenic soya bean (Glycine max) plants overexpressing a seed-specific bacterial phytoene synthase gene from Pantoea ananatis modified to target to plastids accumulated 845 µg ß carotene g(-1) dry seed weight with a desirable 12:1 ratio of ß to α. The ß carotene accumulating seeds exhibited a shift in oil composition increasing oleic acid with a concomitant decrease in linoleic acid and an increase in seed protein content by at least 4% (w/w). Elevated ß-carotene accumulating soya bean cotyledons contain 40% the amount of abscisic acid compared to nontransgenic cotyledons. Proteomic and nontargeted metabolomic analysis of the mid-maturation ß-carotene cotyledons compared to the nontransgenic did not reveal any significant differences that would account for the altered phenotypes of both elevated oleate and protein content. Transcriptomic analysis, confirmed by RT-PCR, revealed a number of significant differences in ABA-responsive transcripton factor gene expression in the crtB transgenics compared to nontransgenic cotyledons of the same maturation stage. The altered seed composition traits seem to be attributed to altered ABA hormone levels varying transcription factor expression. The elevated ß-carotene, oleic acid and protein traits in the ß-carotene soya beans confer a substantial additive nutritional quality to soya beans.

Impact of the genetic background on the composition of the chicken plasma MiRNome in response to a stress.

Circulating extra-cellular microRNAs (miRNAs) have emerged as promising minimally invasive markers in human medicine. We evaluated miRNAs isolated from total plasma as biomarker candidates of a response to an abiotic stress (feed deprivation) in a livestock species. Two chicken lines selected for high (R+) and low (R-) residual feed intake were chosen as an experimental model because of their extreme divergence in feed intake and energy metabolism. Adult R+ and R- cocks were sampled after 16 hours of feed deprivation and again four hours after re-feeding. More than 292 million sequence reads were generated by small RNA-seq of total plasma RNA. A total of 649 mature miRNAs were identified; after quality filtering, 148 miRNAs were retained for further analyses. We identified 23 and 19 differentially abundant miRNAs between feeding conditions and between lines respectively, with only two miRNAs identified in both comparisons. We validated a panel of six differentially abundant miRNAs by RT-qPCR on a larger number of plasma samples and checked their response to feed deprivation in liver. Finally, we evaluated the conservation and tissue distribution of differentially abundant miRNAs in plasma across a variety of red jungle fowl tissues. We show that the chicken plasma miRNome reacts promptly to the alteration of the animal physiological condition driven by a feed deprivation stress. The plasma content of stress-responsive miRNAs is strongly influenced by the genetic background, with differences reflecting the phenotypic divergence acquired through long-term selection, as evidenced by the profiles of conserved miRNAs with a regulatory role in energy metabolism (gga-miR-204, gga-miR-let-7f-5p and gga-miR-122-5p). These results reinforce the emerging view in human medicine that even small genetic differences can have a considerable impact on the resolution of biomarker studies, and provide support for the emerging interest in miRNAs as potential novel and minimally invasive biomarkers for livestock species.

Sequencing three crocodilian genomes to illuminate the evolution of archosaurs and amniotes.

The International Crocodilian Genomes Working Group (ICGWG) will sequence and assemble the American alligator (Alligator mississippiensis), saltwater crocodile (Crocodylus porosus) and Indian gharial (Gavialis gangeticus) genomes. The status of these projects and our planned analyses are described.

RNA-seq based transcriptional map of bovine respiratory disease pathogen "Histophilus somni 2336".

Genome structural annotation, i.e., identification and demarcation of the boundaries for all the functional elements in a genome (e.g., genes, non-coding RNAs, proteins and regulatory elements), is a prerequisite for systems level analysis. Current genome annotation programs do not identify all of the functional elements of the genome, especially small non-coding RNAs (sRNAs). Whole genome transcriptome analysis is a complementary method to identify "novel" genes, small RNAs, regulatory regions, and operon structures, thus improving the structural annotation in bacteria. In particular, the identification of non-coding RNAs has revealed their widespread occurrence and functional importance in gene regulation, stress and virulence. However, very little is known about non-coding transcripts in Histophilus somni, one of the causative agents of Bovine Respiratory Disease (BRD) as well as bovine infertility, abortion, septicemia, arthritis, myocarditis, and thrombotic meningoencephalitis. In this study, we report a single nucleotide resolution transcriptome map of H. somni strain 2336 using RNA-Seq method.The RNA-Seq based transcriptome map identified 94 sRNAs in the H. somni genome of which 82 sRNAs were never predicted or reported in earlier studies. We also identified 38 novel potential protein coding open reading frames that were absent in the current genome annotation. The transcriptome map allowed the identification of 278 operon (total 730 genes) structures in the genome. When compared with the genome sequence of a non-virulent strain 129Pt, a disproportionate number of sRNAs (∼30%) were located in genomic region unique to strain 2336 (∼18% of the total genome). This observation suggests that a number of the newly identified sRNAs in strain 2336 may be involved in strain-specific adaptations.

Transcriptome-based differentiation of closely-related Miscanthus lines.

BACKGROUND: Distinguishing between individuals is critical to those conducting animal/plant breeding, food safety/quality research, diagnostic and clinical testing, and evolutionary biology studies. Classical genetic identification studies are based on marker polymorphisms, but polymorphism-based techniques are time and labor intensive and often cannot distinguish between closely related individuals. Illumina sequencing technologies provide the detailed sequence data required for rapid and efficient differentiation of related species, lines/cultivars, and individuals in a cost-effective manner. Here we describe the use of Illumina high-throughput exome sequencing, coupled with SNP mapping, as a rapid means of distinguishing between related cultivars of the lignocellulosic bioenergy crop giant miscanthus (Miscanthus × giganteus). We provide the first exome sequence database for Miscanthus species complete with Gene Ontology (GO) functional annotations. RESULTS: A SNP comparative analysis of rhizome-derived cDNA sequences was successfully utilized to distinguish three Miscanthus × giganteus cultivars from each other and from other Miscanthus species. Moreover, the resulting phylogenetic tree generated from SNP frequency data parallels the known breeding history of the plants examined. Some of the giant miscanthus plants exhibit considerable sequence divergence. CONCLUSIONS: Here we describe an analysis of Miscanthus in which high-throughput exome sequencing was utilized to differentiate between closely related genotypes despite the current lack of a reference genome sequence. We functionally annotated the exome sequences and provide resources to support Miscanthus systems biology. In addition, we demonstrate the use of the commercial high-performance cloud computing to do computational GO annotation.

Identifying blood biomarkers and physiological processes that distinguish humans with superior performance under psychological stress.

BACKGROUND: Attrition of students from aviation training is a serious financial and operational concern for the U.S. Navy. Each late stage navy aviator training failure costs the taxpayer over $1,000,000 and ultimately results in decreased operational readiness of the fleet. Currently, potential aviators are selected based on the Aviation Selection Test Battery (ASTB), which is a series of multiple-choice tests that evaluate basic and aviation-related knowledge and ability. However, the ASTB does not evaluate a person's response to stress. This is important because operating sophisticated aircraft demands exceptional performance and causes high psychological stress. Some people are more resistant to this type of stress, and consequently better able to cope with the demands of naval aviation, than others. METHODOLOGY/PRINCIPAL FINDINGS: Although many psychological studies have examined psychological stress resistance none have taken advantage of the human genome sequence. Here we use high-throughput -omic biology methods and a novel statistical data normalization method to identify plasma proteins associated with human performance under psychological stress. We identified proteins involved in four basic physiological processes: innate immunity, cardiac function, coagulation and plasma lipid physiology. CONCLUSIONS/SIGNIFICANCE: The proteins identified here further elucidate the physiological response to psychological stress and suggest a hypothesis that stress-susceptible pilots may be more prone to shock. This work also provides potential biomarkers for screening humans for capability of superior performance under stress.

Sequential detergent extraction prior to mass spectrometry analysis.

Sequential detergent extraction of proteins from eukaryotic cells has been used to increase proteome coverage of 2D-PAGE. We have adapted sequential detergent extraction for use with the high-throughput non-electrophoretic proteomics method of liquid chromatography and electrospray ionisation tandem mass spectrometry. This method of extraction yields comprehensive proteomes that include up to twice as many membrane proteins as other published methods. Two thirds of these membrane proteins have more than one transmembrane domain and many of these have multiple transmembrane domains. Since sequential detergent extraction (SDE) separates proteins based upon their physicochemistry and sub-cellular localisation, this method also provides useful data about cellular localisation.

Localization of spindle checkpoint proteins in cells undergoing mitosis with unreplicated genomes.

CHO cells can be arrested with hydoxyurea at the beginning of the DNA synthesis phase of the cell cycle. Subsequent treatment with the xanthine, caffeine, induces cells to bypass the S-phase checkpoint and enter unscheduled mitosis [Schlegel and Pardee,1986, Science 232:1264-1266]. These treated cells build a normal spindle and distribute kinetochores, unattached to chromosomes, to their daughter cells [Brinkley et al.,1988, Nature 336:251-254; Zinkowski et al.,1991, J Cell Biol 113:1091-1110; Wise and Brinkley,1997, Cell Motil Cytoskeleton 36:291-302; Balczon et al.,2003, Chromosoma 112:96-102]. To investigate how these cells distribute kinetochores to daughter cells, we analyzed the spindle checkpoint components, Mad2, CENP-E, and the 3F3 phosphoepitope, using immunofluorescence and digital microscopy. Even though the kinetochores were unpaired and DNA was fragmented, the tension, alignment, and motor components of the checkpoint were found to be present and localized as predicted in prometaphase and metaphase. This unusual mitosis proves that a cell can successfully localize checkpoint proteins and divide even when kinetochores are unpaired and fragmented.

Modeling a whole organ using proteomics: the avian bursa of Fabricius.

While advances in proteomics have improved proteome coverage and enhanced biological modeling, modeling function in multicellular organisms requires understanding how cells interact. Here we used the chicken bursa of Fabricius, a common experimental system for B cell function, to model organ function from proteomics data. The bursa has two major functional cell types: B cells and the supporting stromal cells. We used differential detergent fractionation-multidimensional protein identification technology (DDF-MudPIT) to identify 5198 proteins from all cellular compartments. Of these, 1753 were B cell specific, 1972 were stroma specific and 1473 were shared between the two. By modeling programmed cell death (PCD), cell differentiation and proliferation, and transcriptional activation, we have improved functional annotation of chicken proteins and placed chicken-specific death receptors into the PCD process using phylogenetics. We have identified 114 transcription factors (TFs); 42 of the bursal B cell TFs have not been reported before in any B cells. We have also improved the structural annotation of a newly sequenced genome by confirming the in vivo expression of 4006 "predicted", and 6623 ab initio, ORFs. Finally, we have developed a novel method for facilitating structural annotation, "expressed peptide sequence tags" (ePSTs) and demonstrate its utility by identifying 521 potential novel proteins from the chicken "unassigned chromosome".