VPBrowse: Genome-based representation of MS/MS spectra to quantify 10,000 bovine proteins.

SWATH is a data acquisition strategy acclaimed for generating quantitatively accurate and consistent measurements of proteins across multiple samples. Its utility for proteomics studies in nonlaboratory animals, however, is currently compromised by the lack of sufficiently comprehensive and reliable public libraries, either experimental or predicted, and relevant platforms that support their sharing and utilization in an intuitive manner. Here we describe the development of the Veterinary Proteome Browser, VPBrowse (http://browser.proteo.cloud/), an on-line platform for genome-based representation of the Bos taurus proteome, which is equipped with an interactive database and tools for searching, visualization, and building quantitative mass spectrometry assays. In its current version (VPBrowse 1.0), it contains high-quality fragmentation spectra acquired on QToF instrument for over 36,000 proteotypic peptides, the experimental evidence for over 10,000 proteins. Data can be downloaded in different formats to enable analysis using popular software packages for SWATH data processing whilst normalization to iRT scale ensures compatibility with diverse chromatography systems. When applied to published blood plasma dataset from the biomarker discovery study, the resource supported label-free quantification of additional proteins not reported by the authors previously including PSMA4, a tissue leakage protein and a promising candidate biomarker of animal's response to dehorning-related injury.

A cold case of hemolytic disease of the fetus and newborn resolved by genomic sequencing and population studies to define a new antigen in the Rh system.

BACKGROUND: We report an obstetric case involving an RhD-positive woman who had developed a red blood cell (RBC) antibody that was not detected until after delivery of a newborn, who presented with a positive direct antiglobulin test result. Immunohematology studies suggested that the maternal antibody was directed against a low-prevalence antigen on the paternal and newborn RBCs. RESULTS: Comprehensive blood group profiling by targeted exome sequencing revealed a novel nonsynonymous single nucleotide variant (SNV) RHCE c.486C>G (GenBank MZ326705) on the RHCE*Ce allele, for both the father and newborn. A subsequent genomic-based study to profile blood groups in an Indigenous Australian population revealed the same SNV in 2 of 247 individuals. Serology testing showed that the maternal antibody reacted specifically with RBCs from these two individuals. DISCUSSION: The maternal antibody was directed against a novel antigen in the Rh blood group system arising from an RHCE c.486C>G variant on the RHCE*Ce allele linked to RHD*01. The variant predicts a p.Asn162Lys change on the RhCE protein and has been registered as the 56th antigen in the Rh system, ISBT RH 004063. CONCLUSION: This antibody was of clinical significance, resulting in a mild to moderate hemolytic disease of the fetus and newborn (HDFN). In the past, the cause of such HDFN cases may have remained unresolved. Genomic sequencing combined with population studies now assists in resolving such cases. Further population studies have potential to inform the need to design population-specific red cell antibody typing panels for antibody screening in the Australian population.

Equitable Expanded Carrier Screening Needs Indigenous Clinical and Population Genomic Data.

Expanded carrier screening (ECS) for recessive monogenic diseases requires prior knowledge of genomic variation, including DNA variants that cause disease. The composition of pathogenic variants differs greatly among human populations, but historically, research about monogenic diseases has focused mainly on people with European ancestry. By comparison, less is known about pathogenic DNA variants in people from other parts of the world. Consequently, inclusion of currently underrepresented Indigenous and other minority population groups in genomic research is essential to enable equitable outcomes in ECS and other areas of genomic medicine. Here, we discuss this issue in relation to the implementation of ECS in Australia, which is currently being evaluated as part of the national Government's Genomics Health Futures Mission. We argue that significant effort is required to build an evidence base and genomic reference data so that ECS can bring significant clinical benefit for many Aboriginal and/or Torres Strait Islander Australians. These efforts are essential steps to achieving the Australian Government's objectives and its commitment "to leveraging the benefits of genomics in the health system for all Australians." They require culturally safe, community-led research and community involvement embedded within national health and medical genomics programs to ensure that new knowledge is integrated into medicine and health services in ways that address the specific and articulated cultural and health needs of Indigenous people. Until this occurs, people who do not have European ancestry are at risk of being, in relative terms, further disadvantaged.

Protein-coding potential of non-canonical open reading frames in human transcriptome.

In recent years, proteogenomics and ribosome profiling studies have identified a large number of proteins encoded by noncoding regions in the human genome. They are encoded by small open reading frames (sORFs) in the untranslated regions (UTRs) of mRNAs and long non-coding RNAs (lncRNAs). These sORF encoded proteins (SEPs) are often <150AA and show poor evolutionary conservation. A subset of them have been functionally characterized and shown to play an important role in fundamental biological processes including cardiac and muscle function, DNA repair, embryonic development and various human diseases. How many novel protein-coding regions exist in the human genome and what fraction of them are functionally important remains a mystery. In this review, we discuss current progress in unraveling SEPs, approaches used for their identification, their limitations and reliability of these identifications. We also discuss functionally characterized SEPs and their involvement in various biological processes and diseases. Lastly, we provide insights into their distinctive features compared to canonical proteins and challenges associated with annotating these in protein reference databases.

A comparative study of multi-omics integration tools for cancer driver gene identification and tumour subtyping.

Oncogenesis and cancer can arise as a consequence of a wide range of genomic aberrations including mutations, copy number alterations, expression changes and epigenetic modifications encompassing multiple omics layers. Integrating genomic, transcriptomic, proteomic and epigenomic datasets via multi-omics analysis provides the opportunity to derive a deeper and holistic understanding of the development and progression of cancer. There are two primary approaches to integrating multi-omics data: multi-staged (focused on identifying genes driving cancer) and meta-dimensional (focused on establishing clinically relevant tumour or sample classifications). A number of ready-to-use bioinformatics tools are available to perform both multi-staged and meta-dimensional integration of multi-omics data. In this study, we compared nine different integration tools using real and simulated cancer datasets. The performance of the multi-staged integration tools were assessed at the gene, function and pathway levels, while meta-dimensional integration tools were assessed based on the sample classification performance. Additionally, we discuss the influence of factors such as data representation, sample size, signal and noise on multi-omics data integration. Our results provide current and much needed guidance regarding selection and use of the most appropriate and best performing multi-omics integration tools.

Pharmacogenomic analysis of a genetically distinct Indigenous population.

Indigenous Australians face a disproportionately severe burden of chronic disease relative to other Australians, with elevated rates of morbidity and mortality. While genomics technologies are slowly gaining momentum in personalised treatments for many, a lack of pharmacogenomic research in Indigenous peoples could delay adoption. Appropriately implementing pharmacogenomics in clinical care necessitates an understanding of the frequencies of pharmacologically relevant genetic variants within Indigenous populations. We analysed whole-genome sequence data from 187 individuals from the Tiwi Islands and characterised the pharmacogenomic landscape of this population. Specifically, we compared variant profiles and allelic distributions of previously described pharmacologically significant genes and variants with other population groups. We identified 22 translationally relevant pharmacogenomic variants and 18 clinically actionable guidelines with implications for drug dosing and treatment of conditions including heart disease, diabetes and cancer. We specifically observed increased poor and intermediate metabolizer phenotypes in the CYP2C9 (PM:19%, IM:44%) and CYP2C19 (PM:18%, IM:44%) genes.

The genomic landscape of blood groups in Indigenous Australians in remote communities.

BACKGROUND: Red blood cell (RBC) membrane-associated blood group systems are clinically significant. Alloimmunisation is a persistent risk associated with blood transfusion owing to the antigen polymorphisms among these RBC-associated blood groups. Next-generation sequencing (NGS) offers an opportunity to characterize the blood group variant profile of a given individual. Australia comprises a large multiethnic population where most blood donors are Caucasian and blood group variants remain poorly studied among Indigenous Australians. In this study, we focused on the Tiwi Islanders, who have lived in relative isolation for thousands of years. METHODS AND MATERIALS: We predicted the blood group phenotype profiles in the Tiwi (457) and 1000 Genomes Phase 3 (1KGP3-2504) cohort individuals using RBCeq (https://www.rbceq.org/). The predicted phenotype prevalence was compared with the previous literature report. RESULTS: We report, for the first time, comprehensive blood group profiles corresponding to the 35 known blood group systems among the Indigenous Tiwi islander population and identify possible novel antigen variants therein. Our results demonstrate that the genetic makeup of the Tiwi participants is distinct from that of other populations, with a low prevalence of LU (Au[a-b+]) and ABO (A2) and D+C+c+E+e- phenotype, an absence of Diego blood group variants, and a unique RHD (DIII type4) variant. CONCLUSION: Our results may contribute to the development of a database of predicted phenotype donors among the Tiwi population and aid in improving transfusion safety for the ~2800 Tiwi people and the ~800,000 other Indigenous Australians throughout the nation.

BGvar: A comprehensive resource for blood group immunogenetics.

BACKGROUND: Blood groups form the basis of effective and safe blood transfusion. There are about 43 well-recognised human blood group systems presently known. Blood groups are molecularly determined by the presence of specific antigens on the red blood cells and are genetically determined and inherited following Mendelian principles. The lack of a comprehensive, relevant, manually compiled and genome-ready dataset of red cell antigens limited the widespread application of genomic technologies to characterise and interpret the blood group complement of an individual from genomic datasets. MATERIALS AND METHODS: A range of public datasets was used to systematically annotate the variation compendium for its functionality and allele frequencies across global populations. Details on phenotype or relevant clinical importance were collated from reported literature evidence. RESULTS: We have compiled the Blood Group Associated Genomic Variant Resource (BGvar), a manually curated online resource comprising all known human blood group related allelic variants including a total of 1700 International Society of Blood Transfusion approved alleles and 1706 alleles predicted and curated from literature reports. This repository includes 1682 single nucleotide variations (SNVs), 310 Insertions, Deletions (InDels) and Duplications (Copy Number Variations) and about 1360 combination mutations corresponding to 43 human blood group systems and 2 transcription factors. This compendium also encompasses gene fusion and rearrangement events occurring in human blood group genes. CONCLUSION: To the best of our knowledge, BGvar is a comprehensive and a user-friendly resource with most relevant collation of blood group alleles in humans. BGvar is accessible online at URL: http://clingen.igib.res.in/bgvar/.

Targeting epithelial-mesenchymal plasticity in cancer: clinical and preclinical advances in therapy and monitoring.

The concept of epithelial-mesenchymal plasticity (EMP), which describes the dynamic flux within the spectrum of phenotypic states that invasive carcinoma cells may reside, is being increasingly recognised for its role in cancer progression and therapy resistance. The myriad of events that are able to induce EMP, as well as the more recently characterised control loops, results in dynamic transitions of cancerous epithelial cells to more mesenchymal-like phenotypes through an epithelial-mesenchymal transition (EMT), as well as the reverse transition from mesenchymal phenotypes to an epithelial one. The significance of EMP, in its ability to drive local invasion, generate cancer stem cells and facilitate metastasis by the dissemination of circulating tumour cells (CTCs), highlights its importance as a targetable programme to combat cancer morbidity and mortality. The focus of this review is to consolidate the existing knowledge on the strategies currently in development to combat cancer progression via inhibition of specific facets of EMP. The prevalence of relapse due to therapy resistance and metastatic propensity that EMP endows should be considered when designing therapy regimes, and such therapies should synergise with existing chemotherapeutics to benefit efficacy. To further improve upon EMP-targeted therapies, it is imperative to devise monitoring strategies to assess the impact of such treatments on EMP-related phenomenon such as CTC burden, chemosensitivity/-resistance and micrometastasis in patients.

PGTools: A Software Suite for Proteogenomic Data Analysis and Visualization.

We describe PGTools, an open source software suite for analysis and visualization of proteogenomic data. PGTools comprises applications, libraries, customized databases, and visualization tools for analysis of mass-spectrometry data using combined proteomic and genomic backgrounds. A single command is sufficient to search databases, calculate false discovery rates, group and annotate proteins, generate peptide databases from RNA-Seq transcripts, identify altered proteins associated with cancer, and visualize genome scale peptide data sets using sophisticated visualization tools. We experimentally confirm a subset of proteogenomic peptides in human PANC-1 cells and demonstrate the utility of PGTools using a colorectal cancer data set that led to the identification of 203 novel protein coding regions missed by conventional proteomic approaches. PGTools should be equally useful for individual proteogenomic investigations as well as international initiatives such as chromosome-centric Human Proteome Project (C-HPP). PGTools is available at http://qcmg.org/bioinformatics/PGTools.

Proteogenomic analysis of Bradyrhizobium japonicum USDA110 using GenoSuite, an automated multi-algorithmic pipeline.

We present GenoSuite, an integrated proteogenomic pipeline to validate, refine and discover protein coding genes using high-throughput mass spectrometry (MS) data from prokaryotes. To demonstrate the effectiveness of GenoSuite, we analyzed proteomics data of Bradyrhizobium japonicum (USDA110), a model organism to study agriculturally important rhizobium-legume symbiosis. Our analysis confirmed 31% of known genes, refined 49 gene models for their translation initiation site (TIS) and discovered 59 novel protein coding genes. Notably, a novel protein which redefined the boundary of a crucial cytochrome P450 system related operon was discovered, known to be highly expressed in the anaerobic symbiotic bacteroids. A focused analysis on N-terminally acetylated peptides indicated downstream TIS for gene blr0594. Finally, ortho-proteogenomic analysis revealed three novel genes in recently sequenced B. japonicum USDA6(T) genome. The discovery of large number of missing genes and correction of gene models have expanded the proteomic landscape of B. japonicum and presents an unparalleled utility of proteogenomic analyses and versatility of GenoSuite for annotating prokaryotic genomes including pathogens.

Cloning and tissue distribution of novel splice variants of the ovine ghrelin gene.

BACKGROUND: The ghrelin axis is involved in the regulation of metabolism, energy balance, and the immune, cardiovascular and reproductive systems. The manipulation of this axis has potential for improving economically valuable traits in production animals, and polymorphisms in the ghrelin (GHRL) and ghrelin receptor (GHSR) genes have been associated with growth and carcass traits. Here we investigate the structure and expression of the ghrelin gene (GHRL) in sheep, Ovis aries. RESULTS: We identify two ghrelin mRNA isoforms, which we have designated Δex2 preproghrelin and Δex2,3 preproghrelin. Expression of Δex2,3 preproghrelin is likely to be restricted to ruminants, and would encode truncated ghrelin and a novel C-terminal peptide. Both Δex2 preproghrelin and canonical preproghrelin mRNA isoforms were expressed in a range of tissues. Expression of the Δex2,3 preproghrelin isoform, however, was restricted to white blood cells (WBC; where the wild-type preproghrelin isoform is not co-expressed), and gastrointestinal tissues. Expression of Δex2 preproghrelin and Δex2,3 preproghrelin mRNA was elevated in white blood cells in response to parasitic worm (helminth) infection in genetically susceptible sheep, but not in resistant sheep. CONCLUSIONS: The restricted expression of the novel preproghrelin variants and their distinct WBC expression pattern during parasite infection may indicate a novel link between the ghrelin axis and metabolic and immune function in ruminants.

Mapping the Pharmacogenetic Landscape in a Ugandan Population: Implications for Personalized Medicine in an Underrepresented Population.

Africans are extremely underrepresented in global genomic research. African populations face high burdens of communicable and non-communicable diseases and experience widespread polypharmacy. As population-specific genetic studies are crucial to understanding unique genetic profiles and optimizing treatments to reduce medication-related complications in this diverse population, the present study aims to characterize the pharmacogenomics profile of a rural Ugandan population. We analyzed low-pass whole genome sequencing data from 1998 Ugandans to investigate 18 clinically actionable pharmacogenes in this population. We utilized PyPGx to identify star alleles (haplotype patterns) and compared allele frequencies across populations using the Pharmacogenomics Knowledgebase PharmGKB. Clinical interpretations of the identified alleles were conducted following established dosing guidelines. Over 99% of participants displayed actionable phenotypes across the 18 pharmacogenes, averaging 3.5 actionable genotypes per individual. Several variant alleles known to affect drug metabolism (i.e., CYP3A5*1, CYP2B6*9, CYP3A5*6, CYP2D6*17, CYP2D6*29, and TMPT*3C)-which are generally more prevalent in African individuals-were notably enriched in the Ugandan cohort, beyond reported frequencies in other African peoples. More than half of the cohort exhibited a predicted impaired drug response associated with CFTR, IFNL3, CYP2B6, and CYP2C19, and approximately 31% predicted altered CYP2D6 metabolism. Potentially impaired CYP2C9, SLCO1B1, TPMT, and DPYD metabolic phenotypes were also enriched in Ugandans compared with other African populations. Ugandans exhibit distinct allele profiles that could impact drug efficacy and safety. Our findings have important implications for pharmacogenomics in Uganda, particularly with respect to the treatment of prevalent communicable and non-communicable diseases, and they emphasize the potential of pharmacogenomics-guided therapies to optimize healthcare outcomes and precision medicine in Uganda.

Novel genetic markers for chronic kidney disease in a geographically isolated population of Indigenous Australians: Individual and multiple phenotype genome-wide association study.

BACKGROUND: Chronic kidney disease (CKD) is highly prevalent among Indigenous Australians, especially those in remote regions. The Tiwi population has been isolated from mainland Australia for millennia and exhibits unique genetic characteristics that distinguish them from other Indigenous and non-Indigenous populations. Notably, the rate of end-stage renal disease is up to 20 times greater in this population compared to non-Indigenous populations. Despite the identification of numerous genetic loci associated with kidney disease through GWAS, the Indigenous population such as Tiwi remains severely underrepresented and the increased prevalence of CKD in this population may be due to unique disease-causing alleles/genes. METHODS: We used albumin-to-creatinine ratio (ACR) and estimated glomerular filtration rate (eGFR) to estimate the prevalence of kidney disease in the Tiwi population (N = 492) in comparison to the UK Biobank (UKBB) (N = 134,724) database. We then performed an exploratory factor analysis to identify correlations among 10 CKD-related phenotypes and identify new multi-phenotype factors. We subsequently conducted a genome-wide association study (GWAS) on all single and multiple phenotype factors using mixed linear regression models, adjusted for age, sex, population stratification, and genetic relatedness between individuals. RESULTS: Based on ACR, 20.3% of the population was at severely increased risk of CKD progression and showed elevated levels of ACR compared to the UKBB population independent of HbA1c. A GWAS of ACR revealed novel association loci in the genes MEG3 (chr14:100812018:T:A), RAB36 (rs11704318), and TIAM2 (rs9689640). Additionally, multiple phenotypes GWAS of ACR, eGFR, urine albumin, and serum creatinine identified a novel variant that mapped to the gene MEIS2 (chr15:37218869:A:G). Most of the identified variants were found to be either absent or rare in the UKBB population. CONCLUSIONS: Our study highlights the Tiwi population's predisposition towards elevated ACR, and the collection of novel genetic variants associated with kidney function. These associations may prove valuable in the early diagnosis and treatment of renal disease in this underrepresented population. Additionally, further research is needed to comprehensively validate the functions of the identified variants/genes.

Genomic characterisation of clinically significant blood group variants in Aboriginal Australians.

BACKGROUND: Hematological disorders are often treated with blood transfusions. Many blood group antigens and variants are population-specific, and for patients with rare blood types, extensive donor screening is required to find suitable matches for transfusion. There is a scarcity of knowledge regarding blood group variants in Aboriginal Australian populations, despite a higher need for transfusion due to the higher prevalence of renal diseases and anaemia. MATERIALS AND METHODS: In this study, we applied next-generation sequencing and analysis to 245 samples obtained from Aboriginal Australians from South-East Queensland, to predict antigen phenotypes for 36 blood group systems. RESULTS: We report potential weak antigens in blood group systems RH, FY and JR that have potential clinical implications in transfusion and pregnancy settings. These include partial DIII type 4, weak D type 33, and Del RHD (IVS2-2delA). The rare Rh phenotypes D+ C+ E+ c- e+ and D+ C+ E+ c+ e- were also detected. DISCUSSION: The comprehensive analyses of blood group genetic variant profiles identified in this study will provide insight and an opportunity to improve Aboriginal health by aiding in the identification of appropriate blood products for population-specific transfusion needs.

Pharmacological induction of chromatin remodeling drives chemosensitization in triple-negative breast cancer.

Targeted therapies have improved outcomes for certain cancer subtypes, but cytotoxic chemotherapy remains a mainstay for triple-negative breast cancer (TNBC). The epithelial-to-mesenchymal transition (EMT) is a developmental program co-opted by cancer cells that promotes metastasis and chemoresistance. There are no therapeutic strategies specifically targeting mesenchymal-like cancer cells. We report that the US Food and Drug Administration (FDA)-approved chemotherapeutic eribulin induces ZEB1-SWI/SNF-directed chromatin remodeling to reverse EMT that curtails the metastatic propensity of TNBC preclinical models. Eribulin induces mesenchymal-to-epithelial transition (MET) in primary TNBC in patients, but conventional chemotherapy does not. In the treatment-naive setting, but not after acquired resistance to other agents, eribulin sensitizes TNBC cells to subsequent treatment with other chemotherapeutics. These findings provide an epigenetic mechanism of action of eribulin, supporting its use early in the disease process for MET induction to prevent metastatic progression and chemoresistance. These findings warrant prospective clinical evaluation of the chemosensitizing effects of eribulin in the treatment-naive setting.

Association weight matrix for the genetic dissection of puberty in beef cattle.

We describe a systems biology approach for the genetic dissection of complex traits based on applying gene network theory to the results from genome-wide associations. The associations of single-nucleotide polymorphisms (SNP) that were individually associated with a primary phenotype of interest, age at puberty in our study, were explored across 22 related traits. Genomic regions were surveyed for genes harboring the selected SNP. As a result, an association weight matrix (AWM) was constructed with as many rows as genes and as many columns as traits. Each {i, j} cell value in the AWM corresponds to the z-score normalized additive effect of the ith gene (via its neighboring SNP) on the jth trait. Columnwise, the AWM recovered the genetic correlations estimated via pedigree-based restricted maximum-likelihood methods. Rowwise, a combination of hierarchical clustering, gene network, and pathway analyses identified genetic drivers that would have been missed by standard genome-wide association studies. Finally, the promoter regions of the AWM-predicted targets of three key transcription factors (TFs), estrogen-related receptor gamma (ESRRG), Pal3 motif, bound by a PPAR-gamma homodimer, IR3 sites (PPARG), and Prophet of Pit 1, PROP paired-like homeobox 1 (PROP1), were surveyed to identify binding sites corresponding to those TFs. Applied to our case, the AWM results recapitulate the known biology of puberty, captured experimentally validated binding sites, and identified candidate genes and gene-gene interactions for further investigation.

Machine Learning-Based Prediction of Distant Recurrence in Invasive Breast Carcinoma Using Clinicopathological Data: A Cross-Institutional Study.

Breast cancer is the most common type of cancer worldwide. Alarmingly, approximately 30% of breast cancer cases result in disease recurrence at distant organs after treatment. Distant recurrence is more common in some subtypes such as invasive breast carcinoma (IBC). While clinicians have utilized several clinicopathological measurements to predict distant recurrences in IBC, no studies have predicted distant recurrences by combining clinicopathological evaluations of IBC tumors pre- and post-therapy with machine learning (ML) models. The goal of our study was to determine whether classification-based ML techniques could predict distant recurrences in IBC patients using key clinicopathological measurements, including pathological staging of the tumor and surrounding lymph nodes assessed both pre- and post-neoadjuvant therapy, response to therapy via standard-of-care imaging, and binary status of adjuvant therapy administered to patients. We trained and tested four clinicopathological ML models using a dataset (144 and 17 patients for training and testing, respectively) from Duke University and validated the best-performing model using an external dataset (8 patients) from Dartmouth Hitchcock Medical Center. The random forest model performed better than the C-support vector classifier, multilayer perceptron, and logistic regression models, yielding AUC values of 1.0 in the testing set and 0.75 in the validation set (p < 0.002) across both institutions, thereby demonstrating the cross-institutional portability and validity of ML models in the field of clinical research in cancer. The top-ranking clinicopathological measurement impacting the prediction of distant recurrences in IBC were identified to be tumor response to neoadjuvant therapy as evaluated via SOC imaging and pathology, which included tumor as well as node staging.