First detection and complete genome sequence of a new potexvirus naturally infecting Adenium obesum.

Here, we report the detection and complete genome sequence of a novel potexvirus, tentatively named "Adenium obesum virus X" (AobVX), isolated from Adenium obesum, that was sent for virus screening at Australian Government post-entry quarantine (PEQ) facilities after being imported into Australia from China. The AobVX genome is 6781 nucleotides in length excluding the poly(A) tail and is predicted to encode conserved potexvirus proteins and sequence motifs across five open reading frames. The RNA-dependent RNA polymerase of this virus shares the highest amino acid sequence similarity with that of nerine potexvirus 1 (58.7% identity) and nerine virus X (58.58% identity). This is the first report of a positive-sense single-stranded RNA virus in A. obesum related to members of the genus Potexvirus in the family Alphaflexiviridae.

Prolactin and the evolution of male pregnancy.

Prolactin (PRL) is a multifunctional hormone of broad physiological importance, and is involved in many aspects of fish reproduction, including the regulation of live birth (viviparity) and both male and female parental care. Previous research suggests that PRL also plays an important reproductive role in syngnathid fishes (seahorses, pipefish and seadragons), a group with a highly derived reproductive strategy, male pregnancy - how the PRL axis has come to be co-opted for male pregnancy remains unclear. We investigated the molecular evolution and expression of the genes for prolactin and its receptor (PRLR) in an evolutionarily diverse sampling of syngnathid fishes to explore how the co-option of PRL for male pregnancy has impacted its evolution, and to clarify whether the PRL axis is also involved in regulating reproductive function in species with more rudimentary forms of male pregnancy. In contrast to the majority of teleost fishes, all syngnathid fishes tested carry single copies of PRL and PRLR that cluster genetically within the PRL1 and PRLRa lineages of teleosts, respectively. PRL1 gene expression in seahorses and pipefish is restricted to the pituitary, while PRLRa is expressed in all tissues, including the brood pouch of species with both rudimentary and complex brooding structures. Pituitary PRL1 expression remains stable throughout pregnancy, but PRLRa expression is specifically upregulated in the male brood pouch during pregnancy, consistent with the higher affinity of pouch tissues for PRL hormone during embryonic incubation. Finally, immunohistochemistry of brood pouch tissues reveals that both PRL1 protein and PRLRa and Na+/K+ ATPase-positive cells line the inner pouch epithelium, suggesting that pituitary-derived PRL1 may be involved in brood pouch osmoregulation during pregnancy. Our data provide a unique molecular perspective on the evolution and expression of prolactin and its receptor during male pregnancy, and provide the foundation for further manipulative experiments exploring the role of PRL in this unique form of reproduction.

Implementation of GA-VirReport, a Web-Based Bioinformatics Toolkit for Post-Entry Quarantine Screening of Virus and Viroids in Plants.

High-throughput sequencing (HTS) of host plant small RNA (sRNA) is a popular approach for plant virus and viroid detection. The major bottlenecks for implementing this approach in routine virus screening of plants in quarantine include lack of computational resources and/or expertise in command-line environments and limited availability of curated plant virus and viroid databases. We developed: (1) virus and viroid report web-based bioinformatics workflows on Galaxy Australia called GA-VirReport and GA-VirReport-Stats for detecting viruses and viroids from host plant sRNA extracts and (2) a curated higher plant virus and viroid database (PVirDB). We implemented sRNA sequencing with unique dual indexing on a set of plants with known viruses. Sequencing data were analyzed using GA-VirReport and PVirDB to validate these resources. We detected all known viruses in this pilot study with no cross-sample contamination. We then conducted a large-scale diagnosis of 105 imported plants processed at the post-entry quarantine facility (PEQ), Australia. We detected various pathogens in 14 imported plants and discovered that de novo assembly using 21-22 nt sRNA fraction and the megablast algorithm yielded better sensitivity and specificity. This study reports the successful, large-scale implementation of HTS and a user-friendly bioinformatics workflow for virus and viroid screening of imported plants at the PEQ.

Complete genome sequence of a novel potyvirus infecting Miscanthus sinensis (silver grass).

Here, we describe the full-length genome sequence of a novel potyvirus, tentatively named "Miscanthus sinensis mosaic virus" (MsiMV), isolated from Miscanthus sinensis (silver grass) held in a post-entry quarantine facility after being imported into Western Australia, Australia. The MsiMV genome is 9604 nucleotides (nt) in length, encoding a 3071-amino-acid (aa) polyprotein with conserved sequence motifs. The MsiMV genome is most closely related to that of sorghum mosaic virus (SrMV), with 74% nt and 78.5% aa sequence identity to the SrMV polyprotein region. Phylogenetic analysis based on the polyprotein grouped MsiMV with SrMV, sugarcane mosaic virus (SCMV), and maize dwarf mosaic virus (MDMV). This is the first report of a novel monopartite ssRNA virus in Miscanthus sinensis related to members of the genus Potyvirus in the family Potyviridae.

Side-by-Side Comparison of Post-Entry Quarantine and High Throughput Sequencing Methods for Virus and Viroid Diagnosis.

Rapid and safe access to new plant genetic stocks is crucial for primary plant industries to remain profitable, sustainable, and internationally competitive. Imported plant species may spend several years in Post Entry Quarantine (PEQ) facilities, undergoing pathogen testing which can impact the ability of plant industries to quickly adapt to new global market opportunities by accessing new varieties. Advances in high throughput sequencing (HTS) technologies provide new opportunities for a broad range of fields, including phytosanitary diagnostics. In this study, we compare the performance of two HTS methods (RNA-Seq and sRNA-Seq) with that of existing PEQ molecular assays in detecting and identifying viruses and viroids from various plant commodities. To analyze the data, we tested several bioinformatics tools which rely on different approaches, including direct-read, de novo, and reference-guided assembly. We implemented VirusReport, a new portable, scalable, and reproducible nextflow pipeline that analyses sRNA datasets to detect and identify viruses and viroids. We raise awareness of the need to evaluate cross-sample contamination when analyzing HTS data routinely and of using methods to mitigate index cross-talk. Overall, our results suggest that sRNA analyzed using VirReport provides opportunities to improve quarantine testing at PEQ by detecting all regulated exotic viruses from imported plants in a single assay.

Comprehensive Mutational and Phenotypic Characterization of New Metastatic Cutaneous Squamous Cell Carcinoma Cell Lines Reveal Novel Drug Susceptibilities.

Cutaneous squamous cell carcinoma (cSCC) is a common skin cancer. Most patients who develop metastases (2-5%) present with advanced disease that requires a combination of radical surgery and adjuvant radiation therapy. There are few effective therapies for refractory disease. In this study, we describe novel patient-derived cell lines from cSCC metastases of the head and neck (designated UW-CSCC1 and UW-CSCC2). The cell lines genotypically and phenotypically resembled the original patient tumor and were tumorogenic in mice. Differences in cancer-related gene expression between the tumor and cell lines after various culturing conditions could be largely reversed by xenografting and reculturing. The novel drug susceptibilities of UW-CSCC1 and an irradiated subclone UW-CSCC1-R to drugs targeting cell cycle, PI3K/AKT/mTOR, and DNA damage pathways were observed using high-throughput anti-cancer and kinase-inhibitor compound libraries, which correlate with either copy number variations, targetable mutations and/or the upregulation of gene expression. A secondary screen of top hits in all three cell lines including PIK3CA-targeting drugs supports the utility of targeting the PI3K/AKT/mTOR pathway in this disease. UW-CSCC cell lines are thus useful preclinical models for determining targetable pathways and candidate therapeutics.

Whole genome, transcriptome and methylome profiling enhances actionable target discovery in high-risk pediatric cancer.

The Zero Childhood Cancer Program is a precision medicine program to benefit children with poor-outcome, rare, relapsed or refractory cancer. Using tumor and germline whole genome sequencing (WGS) and RNA sequencing (RNAseq) across 252 tumors from high-risk pediatric patients with cancer, we identified 968 reportable molecular aberrations (39.9% in WGS and RNAseq, 35.1% in WGS only and 25.0% in RNAseq only). Of these patients, 93.7% had at least one germline or somatic aberration, 71.4% had therapeutic targets and 5.2% had a change in diagnosis. WGS identified pathogenic cancer-predisposing variants in 16.2% of patients. In 76 central nervous system tumors, methylome analysis confirmed diagnosis in 71.1% of patients and contributed to a change of diagnosis in two patients (2.6%). To date, 43 patients have received a recommended therapy, 38 of whom could be evaluated, with 31% showing objective evidence of clinical benefit. Comprehensive molecular profiling resolved the molecular basis of virtually all high-risk cancers, leading to clinical benefit in some patients.

Molecular patterns in salivary duct carcinoma identify prognostic subgroups.

Salivary duct carcinoma (SDCa) is a rare cancer with high rate of metastases and poor survival despite aggressive multimodality treatment. This study analyzes the genetic changes in SDCa, their impact on cancer pathways, and evaluates whether molecular patterns can identify subgroups with distinct clinical characteristics and outcome. Clinicopathologic details and tissue samples from 66 patients (48 males, 18 females) treated between 1995 and 2018 were obtained from multiple institutions. Androgen receptor (AR) was assessed by immunohistochemistry, and the Illumina TruSight 170 gene panel was used for DNA sequencing. Male gender, lympho-vascular invasion, lymph node metastasis, and smoking were significant predictors of disease-free survival. AR was present in 79%. Frequently encountered alterations were mutations in TP53 (51%), PIK3CA (32%) and HRAS (22%), as well as amplifications of CDK4/6 (22%), ERBB2 (21%), MYC (16%), and deletions of CDKN2A (13%). TP53 mutation and MYC amplifications were associated with decreased disease-free survival. Analysis of cancer pathways revealed that the PI3K pathway was most commonly affected. Alterations in the cell cycle pathway were associated with impaired disease-free survival (HR 2.6, P = 0.038). Three subgroups based on AR and ERBB2 status were identified, which featured distinct molecular patterns and outcome. Among AR positive SDCa, HRAS mutations were restricted to AR positive tumors without ERBB2 amplification and HRAS mutations featured high co-occurrence with PIK3CA alterations, which seems specific to SDCa. AR negative SDCa were associated with poor disease-free survival in multivariate analysis (HR 4.5, P = 0.010) and none of these tumors exhibited ERBB2 amplification or HRAS mutations. AR and ERBB2 status in SDCa thus classifies tumors with distinct molecular profiles relevant to future targeted therapy. Furthermore, clinical factors such as smoking and molecular features such as MYC amplification may serve as markers of poor prognosis of SDCa.

Development and validation of a targeted gene sequencing panel for application to disparate cancers.

Next generation sequencing has revolutionised genomic studies of cancer, having facilitated the development of precision oncology treatments based on a tumour's molecular profile. We aimed to develop a targeted gene sequencing panel for application to disparate cancer types with particular focus on tumours of the head and neck, plus test for utility in liquid biopsy. The final panel designed through Roche/Nimblegen combined 451 cancer-associated genes (2.01 Mb target region). 136 patient DNA samples were collected for performance and application testing. Panel sensitivity and precision were measured using well-characterised DNA controls (n = 47), and specificity by Sanger sequencing of the Aryl Hydrocarbon Receptor Interacting Protein (AIP) gene in 89 patients. Assessment of liquid biopsy application employed a pool of synthetic circulating tumour DNA (ctDNA). Library preparation and sequencing were conducted on Illumina-based platforms prior to analysis with our accredited (ISO15189) bioinformatics pipeline. We achieved a mean coverage of 395x, with sensitivity and specificity of >99% and precision of >97%. Liquid biopsy revealed detection to 1.25% variant allele frequency. Application to head and neck tumours/cancers resulted in detection of mutations aligned to published databases. In conclusion, we have developed an analytically-validated panel for application to cancers of disparate types with utility in liquid biopsy.

Genomic stratification and liquid biopsy in a rare adrenocortical carcinoma (ACC) case, with dual lung metastases.

Adrenocortical carcinoma is a rare malignancy with a poor prognosis and few treatment options. Molecular characterization of this cancer remains limited. We present a case of an adrenocortical carcinoma (ACC) in a 37-yr-old female, with dual lung metastases identified 1 yr following commencement of adjuvant mitotane therapy. As standard therapeutic regimens are often unsuccessful in ACC, we undertook a comprehensive genomic study into this case to identify treatment options and monitor disease progress. We performed targeted and whole-genome sequencing of germline, primary tumor, and both metastatic tumors from this patient and monitored recurrence over 2 years using liquid biopsy for ctDNA and steroid hormone measurements. Sequencing revealed the primary and metastatic tumors were hyperhaploid, with extensive loss of heterozygosity but few structural rearrangements. Loss-of-function mutations were identified in MSH2, TP53, RB1, and PTEN, resulting in tumors with mismatch repair signatures and microsatellite instability. At the cellular level, tumors were populated by mitochondria-rich oncocytes. Longitudinal ctDNA mutation and hormone profiles were unable to detect micrometastatic disease, consistent with clinical indicators of disease remission. The molecular signatures in our ACC case suggested immunotherapy in the event of disease progression; however, the patient remains free of cancer. The extensive molecular analysis presented here could be applied to other rare and/or poorly stratified cancers to identify novel or repurpose existing therapeutic options, thereby broadly improving diagnoses, treatments, and prognoses.

Mutational Patterns in Metastatic Cutaneous Squamous Cell Carcinoma.

Cutaneous squamous cell carcinoma from the head and neck typically metastasize to the lymph nodes of the neck and parotid glands. When a primary is not identified, they are difficult to distinguish from metastases of mucosal origin and primary salivary gland squamous cell carcinoma. UV radiation causes a mutation pattern that predominantly features cytosine to thymine transitions at dipyrimidine sites and has been associated with cutaneous squamous cell carcinoma. In this study, we used whole genome sequencing data from 15 cutaneous squamous cell carcinoma metastases and show that a UV mutation signature is pervasive across the cohort and distinct from mucosal squamous cell carcinoma. The mutational burden was exceptionally high and concentrated in some regions of the genome, especially insulator elements (mean 162 mutations/megabase). We therefore evaluated the likely impact of UV-induced mutations on the dipyrimidine-rich binding site of the main human insulator protein, CCCTC-binding factor, and the possible implications on CCCTC-binding factor function and the spatial organization of the genome. Our findings suggest that mutation signature analysis may be useful in determining the origin of metastases in the neck and the parotid gland. Furthermore, UV-induced DNA damage to insulator binding sites may play a role in the carcinogenesis and progression of cutaneous squamous cell carcinoma.

Platypus Induced Pluripotent Stem Cells: The Unique Pluripotency Signature of a Monotreme.

The platypus (Ornithorhynchus anatinus) is an egg-laying monotreme mammal whose ancestors diverged ∼166 million years ago from the evolutionary pathway that eventually gave rise to both marsupial and eutherian mammals. Consequently, its genome is an extraordinary amalgam of both ancestral reptilian and derived mammalian features. To gain insight into the evolution of mammalian pluripotency, we have generated induced pluripotent stem cells from the platypus (piPSCs). Deep sequencing of the piPSC transcriptome revealed that piPSCs robustly express the core eutherian pluripotency factors POU5F1/OCT4, SOX2, and NANOG. Given the more extensive role of SOX3 over SOX2 in avian pluripotency, our data indicate that between 315 and 166 million years ago, primitive mammals replaced the role of SOX3 in the vertebrate pluripotency network with SOX2. DAX1/NR0B1 is not expressed in piPSCs and an analysis of the platypus DAX1 promoter revealed the absence of a proximal SOX2-binding DNA motif known to be critical for DAX1 expression in eutherian pluripotent stem cells, suggesting that the acquisition of SOX2 responsiveness by DAX1 has facilitated its recruitment into the pluripotency network of eutherians. Using the RNAseq data, we were also able to demonstrate that in both fibroblasts and piPSCs, the expression ratio of X chromosomes to autosomes (X1-5 X1-5:AA) is approximately equal to 1, indicating that there is no upregulation of X-linked genes. Finally, the RNAseq data also allowed us to explore the process of X-linked gene inactivation in the platypus, where we determined that for any given gene, there is no preference for silencing of the maternal or paternal allele; that is, within a population of cells, the silencing of X-linked genes is not imprinted.