Procrustes is a machine-learning approach that removes cross-platform batch effects from clinical RNA sequencing data.

With the increased use of gene expression profiling for personalized oncology, optimized RNA sequencing (RNA-seq) protocols and algorithms are necessary to provide comparable expression measurements between exome capture (EC)-based and poly-A RNA-seq. Here, we developed and optimized an EC-based protocol for processing formalin-fixed, paraffin-embedded samples and a machine-learning algorithm, Procrustes, to overcome batch effects across RNA-seq data obtained using different sample preparation protocols like EC-based or poly-A RNA-seq protocols. Applying Procrustes to samples processed using EC and poly-A RNA-seq protocols showed the expression of 61% of genes (N = 20,062) to correlate across both protocols (concordance correlation coefficient > 0.8, versus 26% before transformation by Procrustes), including 84% of cancer-specific and cancer microenvironment-related genes (versus 36% before applying Procrustes; N = 1,438). Benchmarking analyses also showed Procrustes to outperform other batch correction methods. Finally, we showed that Procrustes can project RNA-seq data for a single sample to a larger cohort of RNA-seq data. Future application of Procrustes will enable direct gene expression analysis for single tumor samples to support gene expression-based treatment decisions.

Evolutionary dynamics of tipifarnib in HRAS mutated head and neck squamous cell carcinoma.

Head and neck squamous cell carcinoma (HNSCC) is a highly prevalent malignancy worldwide, with a significant proportion of patients developing recurrent and/or metastatic (R/M) disease. Despite recent advances in therapy, the prognosis for patients with advanced HNSCC remains poor. Here, we present the case of a patient with recurrent metastatic HNSCC harboring an HRAS G12S mutation who achieved a durable response to treatment with tipifarnib, a selective inhibitor of farnesyltransferase. The patient was a 48-year-old woman who had previously received multiple lines of therapy with no significant clinical response. However, treatment with tipifarnib resulted in a durable partial response that lasted 8 months. Serial genomic and transcriptomic analyses demonstrated upregulation of YAP1 and AXL in metastatic lesions compared with the primary tumor, the evolution of the tumor microenvironment from an immune-enriched to a fibrotic subtype with increased angiogenesis, and activation of the PI3K/AKT/mTOR pathway in tipifarnib treatment. Lastly, in HRAS-mutated PDXs and in the syngeneic HRAS model, we demonstrated that tipifarnib efficacy is limited by activation of the AKT pathway, and dual treatment with tipifarnib and the PI3K inhibitor, BYL719, resulted in enhanced anti-tumor efficacy. Our case study highlights the potential of targeting HRAS mutations with tipifarnib in R/M HNSCC and identifies potential mechanisms of acquired resistance to tipifarnib, along with immuno-, chemo-, and radiation therapy. Preclinical results provide a firm foundation for further investigation of drug combinations of HRAS-and PI3K -targeting therapeutics in R/M HRAS-driven HNSCC.

Early-stage CCNG1+ HR+ HER2+ Invasive Breast Carcinoma in Older Women: Current Treatment and Future Perspectives for DeltaRex-G, a CCNG1 Inhibitor.

Women with HR+HER2+ early-stage breast cancer are disadvantaged by the lack of clinical trials focused on women ≥70 years of age. In the past years, there has been increasing controversy on the use of toxic chemotherapy as standard of care treatment for early- stage HR+ HER2+ breast carcinoma in older women. With precision medicine coming of age, molecular profiling of tumors and circulating tumor DNA has identified target oncogenes that could be used in designing an optimal treatment for this group of women. This article reviews the current treatment of early-stage triple receptor positive breast cancer, the risks of chemotherapy in older women, and CCNG1, a novel biomarker in development for the use of DeltaRex-G, a CCNG1 inhibitor. Further, future perspectives for DeltaRex-G in older women with early stage CCNG1+ HR+ HER2+ breast cancer are discussed.

Near-chromosomal de novo assembly of Bengal tiger genome reveals genetic hallmarks of apex predation.

The tiger, a poster child for conservation, remains an endangered apex predator. Continued survival and recovery will require a comprehensive understanding of genetic diversity and the use of such information for population management. A high-quality tiger genome assembly will be an important tool for conservation genetics, especially for the Indian tiger, the most abundant subspecies in the wild. Here, we present high-quality near-chromosomal genome assemblies of a female and a male wild Indian tiger (Panthera tigris tigris). Our assemblies had a scaffold N50 of >140 Mb, with 19 scaffolds corresponding to the 19 numbered chromosomes, containing 95% of the genome. Our assemblies also enabled detection of longer stretches of runs of homozygosity compared to previous assemblies, which will help improve estimates of genomic inbreeding. Comprehensive genome annotation identified 26,068 protein-coding genes, including several gene families involved in key morphological features such as the teeth, claws, vision, olfaction, taste, and body stripes. We also identified 301 microRNAs, 365 small nucleolar RNAs, 632 transfer RNAs, and other noncoding RNA elements, several of which are predicted to regulate key biological pathways that likely contribute to the tiger's apex predatory traits. We identify signatures of positive selection in the tiger genome that are consistent with the Panthera lineage. Our high-quality genome will enable use of noninvasive samples for comprehensive assessment of genetic diversity, thus supporting effective conservation and management of wild tiger populations.

Multiomic analysis of papillary thyroid cancers identifies BAIAP2L1-BRAF fusion and requirement of TRIM25, PDE5A and PKCδ for tumorigenesis.

BACKGROUND: Papillary thyroid carcinoma (PTC) is one of the most common forms of thyroid cancer with a cure rate of over 90% after surgery. However, aggressive forms may still occur, and personalized therapeutic strategies are increasingly required. METHODS: We performed integrated genomic and proteomic analysis of PTC tumor samples from patients who did not harbor BRAF or RAS mutations. We validate the analysis and present in-depth molecular analysis of the identified genetic rearrangement by employing biochemical and cell biological assays. Finally, we employ 3D spheroid models, loss of function studies and chemical inhibitors to target the hitherto upregulated factors. The data are analysed with appropriate statistical tests which are mentioned in the legends section. RESULTS: In a 23-year-old patient with thyroiditis, we identified a novel rearrangement leading to a BAIAP2L1-BRAF fusion that transforms immortalized human thyroid cells in a kinase and CC-domain dependent manner. Moreover, quantitative proteomic analysis of the same patient samples revealed the upregulation of several proteins including the Ubiquitin E3 ligase TRIM25, PDE5A, and PKCδ. Further, in a cohort of PTC patients, we observed higher expression of TRIM25 and PKCδ in the tumor and metastatic lesions, when compared to the matched normal tissue. Inhibition of TRIM25, PDE5A and PKCδ with small molecules or RNA interference affected not only viability and proliferation of BAIAP2L1-BRAF transformed cells, but also the viability, growth and invasion of corresponding 3D spheroid cultures. CONCLUSIONS: Apart from unveiling a novel oncogenic BRAF fusion in PTCs, our data may open a novel avenue of therapeutic targeting in human PTCs.

Selected commensals educate the intestinal vascular and immune system for immunocompetence.

BACKGROUND: The intestinal microbiota fundamentally guides the development of a normal intestinal physiology, the education, and functioning of the mucosal immune system. The Citrobacter rodentium-carrier model in germ-free (GF) mice is suitable to study the influence of selected microbes on an otherwise blunted immune response in the absence of intestinal commensals. RESULTS: Here, we describe that colonization of adult carrier mice with 14 selected commensal microbes (OMM12 + MC2) was sufficient to reestablish the host immune response to enteric pathogens; this conversion was facilitated by maturation and activation of the intestinal blood vessel system and the step- and timewise stimulation of innate and adaptive immunity. While the immature colon of C. rodentium-infected GF mice did not allow sufficient extravasation of neutrophils into the gut lumen, colonization with OMM12 + MC2 commensals initiated the expansion and activation of the visceral vascular system enabling granulocyte transmigration into the gut lumen for effective pathogen elimination. CONCLUSIONS: Consortium modeling revealed that the addition of two facultative anaerobes to the OMM12 community was essential to further progress the intestinal development. Moreover, this study demonstrates the therapeutic value of a defined consortium to promote intestinal maturation and immunity even in adult organisms. Video Abstract.

Genetic Variation in Reproductive Investment Across an Ephemerality Gradient in Daphnia pulex.

Species across the tree of life can switch between asexual and sexual reproduction. In facultatively sexual species, the ability to switch between reproductive modes is often environmentally dependent and subject to local adaptation. However, the ecological and evolutionary factors that influence the maintenance and turnover of polymorphism associated with facultative sex remain unclear. We studied the ecological and evolutionary dynamics of reproductive investment in the facultatively sexual model species, Daphnia pulex. We found that patterns of clonal diversity, but not genetic diversity varied among ponds consistent with the predicted relationship between ephemerality and clonal structure. Reconstruction of a multi-year pedigree demonstrated the coexistence of clones that differ in their investment into male production. Mapping of quantitative variation in male production using lab-generated and field-collected individuals identified multiple putative quantitative trait loci (QTL) underlying this trait, and we identified a plausible candidate gene. The evolutionary history of these QTL suggests that they are relatively young, and male limitation in this system is a rapidly evolving trait. Our work highlights the dynamic nature of the genetic structure and composition of facultative sex across space and time and suggests that quantitative genetic variation in reproductive strategy can undergo rapid evolutionary turnover.

Human ACE2 receptor polymorphisms and altered susceptibility to SARS-CoV-2.

COVID-19 is a respiratory illness caused by a novel coronavirus called SARS-CoV-2. The viral spike (S) protein engages the human angiotensin-converting enzyme 2 (ACE2) receptor to invade host cells with ~10-15-fold higher affinity compared to SARS-CoV S-protein, making it highly infectious. Here, we assessed if ACE2 polymorphisms can alter host susceptibility to SARS-CoV-2 by affecting this interaction. We analyzed over 290,000 samples representing >400 population groups from public genomic datasets and identified multiple ACE2 protein-altering variants. Using reported structural data, we identified natural ACE2 variants that could potentially affect virus-host interaction and thereby alter host susceptibility. These include variants S19P, I21V, E23K, K26R, T27A, N64K, T92I, Q102P and H378R that were predicted to increase susceptibility, while variants K31R, N33I, H34R, E35K, E37K, D38V, Y50F, N51S, M62V, K68E, F72V, Y83H, G326E, G352V, D355N, Q388L and D509Y were predicted to be protective variants that show decreased binding to S-protein. Using biochemical assays, we confirmed that K31R and E37K had decreased affinity, and K26R and T92I variants showed increased affinity for S-protein when compared to wildtype ACE2. Consistent with this, soluble ACE2 K26R and T92I were more effective in blocking entry of S-protein pseudotyped virus suggesting that ACE2 variants can modulate susceptibility to SARS-CoV-2.

Integrated genomic analysis reveals mutated ELF3 as a potential gallbladder cancer vaccine candidate.

Gallbladder cancer (GBC) is an aggressive gastrointestinal malignancy with no approved targeted therapy. Here, we analyze exomes (n = 160), transcriptomes (n = 115), and low pass whole genomes (n = 146) from 167 gallbladder cancers (GBCs) from patients in Korea, India and Chile. In addition, we also sequence samples from 39 GBC high-risk patients and detect evidence of early cancer-related genomic lesions. Among the several significantly mutated genes not previously linked to GBC are ETS domain genes ELF3 and EHF, CTNNB1, APC, NSD1, KAT8, STK11 and NFE2L2. A majority of ELF3 alterations are frame-shift mutations that result in several cancer-specific neoantigens that activate T-cells indicating that they are cancer vaccine candidates. In addition, we identify recurrent alterations in KEAP1/NFE2L2 and WNT pathway in GBC. Taken together, these define multiple targetable therapeutic interventions opportunities for GBC treatment and management.

The Indian cobra reference genome and transcriptome enables comprehensive identification of venom toxins.

Snakebite envenoming is a serious and neglected tropical disease that kills ~100,000 people annually. High-quality, genome-enabled comprehensive characterization of toxin genes will facilitate development of effective humanized recombinant antivenom. We report a de novo near-chromosomal genome assembly of Naja naja, the Indian cobra, a highly venomous, medically important snake. Our assembly has a scaffold N50 of 223.35 Mb, with 19 scaffolds containing 95% of the genome. Of the 23,248 predicted protein-coding genes, 12,346 venom-gland-expressed genes constitute the 'venom-ome' and this included 139 genes from 33 toxin families. Among the 139 toxin genes were 19 'venom-ome-specific toxins' (VSTs) that showed venom-gland-specific expression, and these probably encode the minimal core venom effector proteins. Synthetic venom reconstituted through recombinant VST expression will aid in the rapid development of safe and effective synthetic antivenom. Additionally, our genome could serve as a reference for snake genomes, support evolutionary studies and enable venom-driven drug discovery.

Author Correction: Comprehensive analysis of single molecule sequencing-derived complete genome and whole transcriptome of Hyposidra talaca nuclear polyhedrosis virus.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has been fixed in the paper.

Comprehensive analysis of single molecule sequencing-derived complete genome and whole transcriptome of Hyposidra talaca nuclear polyhedrosis virus.

We sequenced the Hyposidra talaca NPV (HytaNPV) double stranded circular DNA genome using PacBio single molecule sequencing technology. We found that the HytaNPV genome is 139,089 bp long with a GC content of 39.6%. It encodes 141 open reading frames (ORFs) including the 37 baculovirus core genes, 25 genes conserved among lepidopteran baculoviruses, 72 genes known in baculovirus, and 7 genes unique to the HytaNPV genome. It is a group II alphabaculovirus that codes for the F protein and lacks the gp64 gene found in group I alphabaculovirus viruses. Using RNA-seq, we confirmed the expression of the ORFs identified in the HytaNPV genome. Phylogenetic analysis showed HytaNPV to be closest to BusuNPV, SujuNPV and EcobNPV that infect other tea pests, Buzura suppressaria, Sucra jujuba, and Ectropis oblique, respectively. We identified repeat elements and a conserved non-coding baculovirus element in the genome. Analysis of the putative promoter sequences identified motif consistent with the temporal expression of the genes observed in the RNA-seq data.

Enhancer identification and activity evaluation in the red flour beetle, Tribolium castaneum.

Evolution of cis-regulatory elements (such as enhancers) plays an important role in the production of diverse morphology. However, a mechanistic understanding is often limited by the absence of methods for studying enhancers in species other than established model systems. Here, we sought to establish methods to identify and test enhancer activity in the red flour beetle, Tribolium castaneum To identify possible enhancer regions, we first obtained genome-wide chromatin profiles from various tissues and stages of Tribolium using FAIRE (formaldehyde-assisted isolation of regulatory elements)-sequencing. Comparison of these profiles revealed a distinct set of open chromatin regions in each tissue and at each stage. In addition, comparison of the FAIRE data with sets of computationally predicted (i.e. supervised cis-regulatory module-predicted) enhancers revealed a very high overlap between the two datasets. Second, using nubbin in the wing and hunchback in the embryo as case studies, we established the first universal reporter assay system that works in various contexts in Tribolium, and in a cross-species context. Together, these advances will facilitate investigation of cis-evolution and morphological diversity in Tribolium and other insects.

Redeployment of a conserved gene regulatory network during Aedes aegypti development.

Changes in gene regulatory networks (GRNs) underlie the evolution of morphological novelty and developmental system drift. The fruitfly Drosophila melanogaster and the dengue and Zika vector mosquito Aedes aegypti have substantially similar nervous system morphology. Nevertheless, they show significant divergence in a set of genes co-expressed in the midline of the Drosophila central nervous system, including the master regulator single minded and downstream genes including short gastrulation, Star, and NetrinA. In contrast to Drosophila, we find that midline expression of these genes is either absent or severely diminished in A. aegypti. Instead, they are co-expressed in the lateral nervous system. This suggests that in A. aegypti this "midline GRN" has been redeployed to a new location while lost from its previous site of activity. In order to characterize the relevant GRNs, we employed the SCRMshaw method we previously developed to identify transcriptional cis-regulatory modules in both species. Analysis of these regulatory sequences in transgenic Drosophila suggests that the altered gene expression observed in A. aegypti is the result of trans-dependent redeployment of the GRN, potentially stemming from cis-mediated changes in the expression of sim and other as-yet unidentified regulators. Our results illustrate a novel "repeal, replace, and redeploy" mode of evolution in which a conserved GRN acquires a different function at a new site while its original function is co-opted by a different GRN. This represents a striking example of developmental system drift in which the dramatic shift in gene expression does not result in gross morphological changes, but in more subtle differences in development and function of the late embryonic nervous system.

Identifying transcriptional cis-regulatory modules in animal genomes.

UNLABELLED: Gene expression is regulated through the activity of transcription factors (TFs) and chromatin-modifying proteins acting on specific DNA sequences, referred to as cis-regulatory elements. These include promoters, located at the transcription initiation sites of genes, and a variety of distal cis-regulatory modules (CRMs), the most common of which are transcriptional enhancers. Because regulated gene expression is fundamental to cell differentiation and acquisition of new cell fates, identifying, characterizing, and understanding the mechanisms of action of CRMs is critical for understanding development. CRM discovery has historically been challenging, as CRMs can be located far from the genes they regulate, have few readily identifiable sequence characteristics, and for many years were not amenable to high-throughput discovery methods. However, the recent availability of complete genome sequences and the development of next-generation sequencing methods have led to an explosion of both computational and empirical methods for CRM discovery in model and nonmodel organisms alike. Experimentally, CRMs can be identified through chromatin immunoprecipitation directed against TFs or histone post-translational modifications, identification of nucleosome-depleted 'open' chromatin regions, or sequencing-based high-throughput functional screening. Computational methods include comparative genomics, clustering of known or predicted TF-binding sites, and supervised machine-learning approaches trained on known CRMs. All of these methods have proven effective for CRM discovery, but each has its own considerations and limitations, and each is subject to a greater or lesser number of false-positive identifications. Experimental confirmation of predictions is essential, although shortcomings in current methods suggest that additional means of validation need to be developed. For further resources related to this article, please visit the WIREs website. CONFLICT OF INTEREST: The authors have declared no conflicts of interest for this article.

Evidence for deep regulatory similarities in early developmental programs across highly diverged insects.

Many genes familiar from Drosophila development, such as the so-called gap, pair-rule, and segment polarity genes, play important roles in the development of other insects and in many cases appear to be deployed in a similar fashion, despite the fact that Drosophila-like "long germband" development is highly derived and confined to a subset of insect families. Whether or not these similarities extend to the regulatory level is unknown. Identification of regulatory regions beyond the well-studied Drosophila has been challenging as even within the Diptera (flies, including mosquitoes) regulatory sequences have diverged past the point of recognition by standard alignment methods. Here, we demonstrate that methods we previously developed for computational cis-regulatory module (CRM) discovery in Drosophila can be used effectively in highly diverged (250-350 Myr) insect species including Anopheles gambiae, Tribolium castaneum, Apis mellifera, and Nasonia vitripennis. In Drosophila, we have successfully used small sets of known CRMs as "training data" to guide the search for other CRMs with related function. We show here that although species-specific CRM training data do not exist, training sets from Drosophila can facilitate CRM discovery in diverged insects. We validate in vivo over a dozen new CRMs, roughly doubling the number of known CRMs in the four non-Drosophila species. Given the growing wealth of Drosophila CRM annotation, these results suggest that extensive regulatory sequence annotation will be possible in newly sequenced insects without recourse to costly and labor-intensive genome-scale experiments. We develop a new method, Regulus, which computes a probabilistic score of similarity based on binding site composition (despite the absence of nucleotide-level sequence alignment), and demonstrate similarity between functionally related CRMs from orthologous loci. Our work represents an important step toward being able to trace the evolutionary history of gene regulatory networks and defining the mechanisms underlying insect evolution.