A multimodal atlas of hepatocellular carcinoma reveals convergent evolutionary paths and 'bad apple' effect on clinical trajectory.

BACKGROUND & AIMS: Hepatocellular Carcinoma (HCC) is a highly fatal cancer characterized by high intra-tumor heterogeneity (ITH). A panoramic understanding of its tumor evolution, in relation to its clinical trajectory, may provide novel prognostic and treatment strategies. METHODS: Through the Asia-Pacific Hepatocellular Carcinoma (AHCC) trials group (NCT03267641), we recruited one of the largest prospective cohorts of HCC with over 600 whole genome and transcriptome samples from 123 treatment-naïve patients. RESULTS: Using a multi-region sampling approach, we revealed seven convergent genetic evolutionary paths governed by the early driver mutations, late copy number variations and viral integrations, which stratify patient clinical trajectories after surgical resection. Furthermore, such evolutionary paths shaped the molecular profiles, leading to distinct transcriptomic subtypes. Most significantly, although we found the coexistence of multiple transcriptomic subtypes within certain tumors, patient prognosis was best predicted by the most aggressive cell fraction of the tumor, rather than by overall degree of transcriptomic ITH level - a phenomenon we termed the 'bad apple' effect. Finally, we found that characteristics throughout early and late tumor evolution provide significant and complementary prognostic power in predicting patient survival. CONCLUSIONS: Taken together, our study generated a comprehensive landscape of evolutionary history for HCC and provided a rich multi-omics resource for understanding tumor heterogeneity and clinical trajectories. CLINICAL TRIAL NUMBER: NCT03267641 (Observational cohort) IMPACT AND IMPLICATIONS: This prospective study, utilizing comprehensive multi-sector, multi-omics sequencing and clinical data from surgically resected HCC, reveals critical insights into the role of tumor evolution and intra-tumor heterogeneity (ITH) in determining the prognosis of Hepatocellular Carcinoma (HCC). These findings are invaluable for oncology researchers and clinicians, as they underscore the influence of distinct evolutionary paths and the 'bad apple' effect, where the most aggressive tumor fraction dictates disease progression. These insights not only enhance prognostic accuracy post-surgical resection but also pave the way for developing personalized therapies tailored to specific tumor evolutionary and transcriptomic profiles. The co-existence of multiple sub-types within the same tumor prompts a re-appraisal of the utilities of depending on single samples to represent the entire tumor and suggests the need for clinical molecular imaging. This research thus marks a significant step forward in the clinical understanding and management of HCC, underscoring the importance of integrating tumor evolutionary dynamics and multi-omics biomarkers into therapeutic decision-making.

Multifunctional sulfur doping in cobalt-based materials for high-energy density supercapacitors.

In this study, S-CCO@Co(OH)2('CCO' representing CuCo2O4/Cu2O; 'S-'representing sulfur doping) was synthesized by hydrothermal method followed by electrodeposition. The multiple effects of S doping were studied by S doping and constructing 3D core-shell structure. S doping induced the reduction of Cu2+and Co3+to Cu+and Co2+, respectively. Also, S partially replaces O and creates oxygen vacancies, which increases a number of active sites for the redox reaction enhancing the redox reaction activity. After the electrodeposition, S-Co bond is formed between the Co(OH)2shell and the S-CCO core, which suggests a synergistic effect between S doping and core-shell structure. The formation of S-Co bond is conducive to electron and ion transport, thus improving electrochemical performance. After modification, the specific capacitance of S-CCO@Co(OH)2is 4.28 times higher than CCO, up to 1730 Fg-1. Furthermore, the assembled S-CCO@Co(OH)2//activated carbon supercapacitor exhibits an energy density of 83.89 Whkg-1at 848.81 Wkg-1and a retention rate of 98.48% after 5000 charge and discharge cycles. Therefore, S doping and its mutual effect with the utilization of the core-shell structure considerably enhanced the electrochemical performance of the CCO-based electrodes, endowing its potential in further application.

Exogenous Methyl Jasmonate Mediated MiRNA-mRNA Network Improves Heat Tolerance of Perennial Ryegrass.

Heat stress can hinder the growth of perennial ryegrass (Lolium perenne L.). Methyl jasmonate (MeJA) applied exogenously can increase heat stress tolerance in plants; however, the regulatory mechanisms involved in heat tolerance mediated by MeJA are poorly understood in perennial ryegrass. Here, the microRNA (miRNA) expression profiles of perennial ryegrass were assessed to elucidate the regulatory pathways associated with heat tolerance induced by MeJA. Plants were subjected to four treatments, namely, control (CK), MeJA pre-treatment (T), heat stress treatment (H), and MeJA pre-treatment + heat stress (TH). According to the results, 102 miRNAs were up-regulated in all treatments, with 20, 27, and 33 miRNAs being up-regulated in the T, H, and TH treatment groups, respectively. The co-expression network analysis between the deferentially expressed miRNAs and their corresponding target genes showed that 20 miRNAs modulated 51 potential target genes. Notably, the miRNAs that targeted genes related to with regards to heat tolerance were driven by MeJA, and they were involved in four pathways: novel-m0258-5p mediated signal transduction, novel-m0350-5p mediated protein homeostasis, miR397-z, miR5658-z, and novel-m0008-5p involved in cell wall component, and miR1144-z and miR5185-z dominated chlorophyll degradation. Overall, the findings of this research paved the way for more research into the heat tolerance mechanism in perennial ryegrass and provided a theoretical foundation for developing cultivars with enhanced heat tolerance.

Population genomic analysis reveals distinct demographics and recent adaptation in the black flying fox (Pteropus alecto).

As the only mammalian group capable of powered flight, bats have many unique biological traits. Previous comparative genomic studies in bats have focused on long-term evolution. However, the micro-evolutionary processes driving recent evolution are largely under-explored. Using resequencing data from 50 black flying foxes (Pteropus alecto), one of the model species for bats, we find that black flying fox has much higher genetic diversity and lower levels of linkage disequilibrium than most of the mammalian species. Demographic inference reveals strong population fluctuations (>100 fold) coinciding with multiple historical events including the last glacial change and Toba super eruption, suggesting that the black flying fox is a very resilient species with strong recovery abilities. While long-term adaptation in the black flying fox is enriched in metabolic genes, recent adaptation in the black flying fox has a unique landscape where recently selected genes are not strongly enriched in any functional category. The demographic history and mode of adaptation suggest that black flying fox might be a well-adapted species with strong evolutionary resilience. Taken together, this study unravels a vibrant landscape of recent evolution for the black flying fox and sheds light on several unique evolutionary processes for bats comparing to other mammalian groups.

Genome instability is associated with ethnic differences between Asians and Europeans in hepatocellular carcinoma.

Hepatocellular carcinoma (HCC) is one of the deadliest cancer types with diverse etiological factors across the world. Although large scale genomic studies have been conducted in different countries, integrative analysis of HCC genomes and ethnic comparison across cohorts are lacking. Methods: We first integrated genomes of 1,349 HCC patients from five large cohorts across the world and applied multiple statistical methods in identifying driver genes. Subsequently, we systematically compared HCC genomes and transcriptomes between Asians and Europeans using the TCGA cohort. Results: We identified 29 novel candidate driver genes, many of which are infrequent tumor suppressors driving late-stage tumor progression. When we systematically compared ethnic differences in the genomic landscape between Asian and European HCCs using the TCGA cohort (n = 348), we found little differences in driver frequencies. Through multi-modal integrative analysis, we found higher genomic instability in Asians together with a collection of molecular events ranging from tumor mutation burden (TMB), copy number alterations as well as transcriptomic subtypes segregating distinctively between two ethnic backgrounds. Strikingly, we identified an Asian specific transcriptomic subtype with multiple ethnically enriched genomic alterations, in particular chromosome 16 deletion, leading to a clinically aggressive RNA subgroup unique to Asians. Integrating multi-modal information, we found that survival models predict patient prognosis much better in Asians than in Europeans, demonstrating a higher potential for precision medicine applications in Asia. Conclusion: For the first time, we have uncovered an unprecedented amount of genomic differences segregating distinctively across ethnicities in HCC and highlighted the importance of differential disease biology and management in HCC across ethnic backgrounds.

Dynamic phenotypic heterogeneity and the evolution of multiple RNA subtypes in hepatocellular carcinoma: the PLANET study.

Intra-tumor heterogeneity (ITH) is a key challenge in cancer treatment, but previous studies have focused mainly on the genomic alterations without exploring phenotypic (transcriptomic and immune) heterogeneity. Using one of the largest prospective surgical cohorts for hepatocellular carcinoma (HCC) with multi-region sampling, we sequenced whole genomes and paired transcriptomes from 67 HCC patients (331 samples). We found that while genomic ITH was rather constant across stages, phenotypic ITH had a very different trajectory and quickly diversified in stage II patients. Most strikingly, 30% of patients were found to contain more than one transcriptomic subtype within a single tumor. Such phenotypic ITH was found to be much more informative in predicting patient survival than genomic ITH and explains the poor efficacy of single-target systemic therapies in HCC. Taken together, we not only revealed an unprecedentedly dynamic landscape of phenotypic heterogeneity in HCC, but also highlighted the importance of studying phenotypic evolution across cancer types.

Genomic landscape of lung adenocarcinoma in East Asians.

Lung cancer is the world's leading cause of cancer death and shows strong ancestry disparities. By sequencing and assembling a large genomic and transcriptomic dataset of lung adenocarcinoma (LUAD) in individuals of East Asian ancestry (EAS; n = 305), we found that East Asian LUADs had more stable genomes characterized by fewer mutations and fewer copy number alterations than LUADs from individuals of European ancestry. This difference is much stronger in smokers as compared to nonsmokers. Transcriptomic clustering identified a new EAS-specific LUAD subgroup with a less complex genomic profile and upregulated immune-related genes, allowing the possibility of immunotherapy-based approaches. Integrative analysis across clinical and molecular features showed the importance of molecular phenotypes in patient prognostic stratification. EAS LUADs had better prediction accuracy than those of European ancestry, potentially due to their less complex genomic architecture. This study elucidated a comprehensive genomic landscape of EAS LUADs and highlighted important ancestry differences between the two cohorts.

PSiTE: a Phylogeny guided Simulator for Tumor Evolution.

SUMMARY: Simulating realistic clonal dynamics of tumors is an important topic in cancer genomics. Here, we present Phylogeny guided Simulator for Tumor Evolution, a tool that can simulate different types of tumor samples including single sector, multi-sector bulk tumor as well as single-cell tumor data under a wide range of evolutionary trajectories. Phylogeny guided Simulator for Tumor Evolution provides an efficient tool for understanding clonal evolution of cancer. AVAILABILITY AND IMPLEMENTATION: PSiTE is implemented in Python and is available at https://github.com/hchyang/PSiTE. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Selection and environmental adaptation along a path to speciation in the Tibetan frog Nanorana parkeri.

Tibetan frogs, Nanorana parkeri, are differentiated genetically but not morphologically along geographical and elevational gradients in a challenging environment, presenting a unique opportunity to investigate processes leading to speciation. Analyses of whole genomes of 63 frogs reveal population structuring and historical demography, characterized by highly restricted gene flow in a narrow geographic zone lying between matrilines West (W) and East (E). A population found only along a single tributary of the Yalu Zangbu River has the mitogenome only of E, whereas nuclear genes of W comprise 89-95% of the nuclear genome. Selection accounts for 579 broadly scattered, highly divergent regions (HDRs) of the genome, which involve 365 genes. These genes fall into 51 gene ontology (GO) functional classes, 14 of which are likely to be important in driving reproductive isolation. GO enrichment analyses of E reveal many overrepresented functional categories associated with adaptation to high elevations, including blood circulation, response to hypoxia, and UV radiation. Four genes, including DNAJC8 in the brain, TNNC1 and ADORA1 in the heart, and LAMB3 in the lung, differ in levels of expression between low- and high-elevation populations. High-altitude adaptation plays an important role in maintaining and driving continuing divergence and reproductive isolation. Use of total genomes enabled recognition of selection and adaptation in and between populations, as well as documentation of evolution along a stepped cline toward speciation.

Development of a new patient-derived xenograft humanised mouse model to study human-specific tumour microenvironment and immunotherapy.

OBJECTIVE: As the current therapeutic strategies for human hepatocellular carcinoma (HCC) have been proven to have limited effectiveness, immunotherapy becomes a compelling way to tackle the disease. We aim to provide humanised mouse (humice) models for the understanding of the interaction between human cancer and immune system, particularly for human-specific drug testing. DESIGN: Patient-derived xenograft tumours are established with type I human leucocyte antigen matched human immune system in NOD-scid Il2rg-/- (NSG) mice. The longitudinal changes of the tumour and immune responses as well as the efficacy of immune checkpoint inhibitors are investigated. RESULTS: Similar to the clinical outcomes, the human immune system in our model is educated by the tumour and exhibits exhaustion phenotypes such as a significant declination of leucocyte numbers, upregulation of exhaustion markers and decreased the production of human proinflammatory cytokines. Notably, cytotoxic immune cells decreased more rapidly compared with other cell types. Tumour infiltrated T cells have much higher expression of exhaustion markers and lower cytokine production compared with peripheral T cells. In addition, tumour-associated macrophages and myeloid-derived suppressor cells are found to be highly enriched in the tumour microenvironment. Interestingly, the tumour also changes gene expression profiles in response to immune responses by upregulating immune checkpoint ligands. Most importantly, in contrast to the NSG model, our model demonstrates both therapeutic and side effects of immune checkpoint inhibitors pembrolizumab and ipilimumab. CONCLUSIONS: Our work provides a model for immune-oncology study and a useful parallel-to-human platform for anti-HCC drug testing, especially immunotherapy.

Whole-Genome Sequencing of African Dogs Provides Insights into Adaptations against Tropical Parasites.

Natural selection in domestic dogs is of great interest in evolutionary biology since dogs have migrated to every inhabited continent of the world alongside humans, and adapted to diverse environments. Here, we explored their demographic history and genetic basis of adaptation to the tropical African environment using whole genome analyses of 19 African indigenous dogs from Nigeria. Demographic analysis suggests that the ancestors of these dogs migrated into Africa from Eurasia 14,000 years ago and underwent a severe founder effect before population expansion. Admixture analysis further reveals that African dog genomes contain about 1.88-3.50% introgression from African golden wolves (Canis anthus). Population genetic analysis identifies 50 positively selected genes linked with immunity, angiogenesis, ultraviolet protection, as well as insulin secretion and sensitivity that may contribute to adaptation to tropical conditions. One of the positively selected genes, adhesion G protein-coupled receptor E1 (ADGRE1), has also been found to be association with severe malaria resistance in African human populations. Functional assessments showed that ADGRE1 provides protective host defense against Plasmodium infections. This result, together with the fact that the inflammatory response to canine babesiosis is similar to complicated falciparum malaria in humans, support the dogs as a model for the study of malaria control and treatment.

Characteristics of Organic Acid Secretion Associated with the Interaction between Burkholderia multivorans WS-FJ9 and Poplar Root System.

The objective of this study was to investigate whether plant-bacteria interaction affects the secretion of organic acids by both organisms and to assess whether the production of IAA by the bacterium increases the secretion of organic acids by root exudates, and if the stress produced by low available phosphorus (P) affects the production of organic acids by bacteria, by roots, or by root exudates in presence of bacterial cultures. With this purpose, we used as a biological model poplar plants and one strain of Burkholderia multivorans able to solubilize P. High performance liquid chromatography was utilized to measure organic acids. The tests, the inductive effects of exogenous indole-3-acetic acid (IAA) on secretion of organic acids, the 2 × 4 × 2 factorial design experiment, and the ability of organic acids to solubilize tricalcium phosphate were performed to investigate the interactive effects. The results showed that, after B. multivorans WS-FJ9 interacted with the poplar root system, the key phosphate-solubilizing driving force was gluconic acid (GA) which was produced in three ways: (1) secreted by the root system in the presence of IAA produced by B. multivorans WS-FJ9; (2) secreted by B. multivorans WS-FJ9; and (3) secreted by the poplar root system in the presence of phosphorus stress. When phosphorus stress was absent, the GA was produced as outlined in (1) and (2) above. These results demonstrated that inoculating B. multivorans WS-FJ9 into the poplar root system could increase the amount of GA secretion and implied that the interaction between B. multivorans WS-FJ9 and the poplar root system could contribute to the increase of P available fraction for poplar plants.

The origin of chow chows in the light of the East Asian breeds.

BACKGROUND: East Asian dog breeds are one of the most ancient groups of dogs that radiated after the domestication of the dog and represent the most basal lineages of dog evolution. Among these, the Chow Chow is an ancient breed that embodies very distinct morphological and physiological features, such as sturdy build, dense coat, and blue/purple tongue. RESULTS: Using a Restricted site Associated DNA (RAD) sequencing approach, we sequenced the genomes of nine Chow Chows from China. Combined with a dataset of 37 canid whole genome sequencing (WGS) from several published works, we found that the Chow Chow is one of the most basal lineages, which originated together with other East Asian breeds, such as the Shar-Pei and Akita. Demographic analysis found that Chow Chows originated from the Chinese indigenous dog about 8300 years ago. The bottleneck leading to Chow Chows was not strong and genetic migration between Chow Chows and other populations is low. Two classes of genes show strong evidence of positive selection along the Chow Chow lineage, namely genes related to metabolism and digestion as well as muscle/heart development and differentiation. CONCLUSIONS: Dog breeds from East Asia, including the Chow Chow, originated from Chinese indigenous dogs very early in time. The genetic bottleneck leading to Chow Chows and migrations with other populations are found to be quite mild. Our current study represents an early endeavor to characterize the origin of East Asian dog breeds and establishes an important reference point for understanding the origin of ancient breeds in Asia.

Out of southern East Asia: the natural history of domestic dogs across the world.

The origin and evolution of the domestic dog remains a controversial question for the scientific community, with basic aspects such as the place and date of origin, and the number of times dogs were domesticated, open to dispute. Using whole genome sequences from a total of 58 canids (12 gray wolves, 27 primitive dogs from Asia and Africa, and a collection of 19 diverse breeds from across the world), we find that dogs from southern East Asia have significantly higher genetic diversity compared to other populations, and are the most basal group relating to gray wolves, indicating an ancient origin of domestic dogs in southern East Asia 33 000 years ago. Around 15 000 years ago, a subset of ancestral dogs started migrating to the Middle East, Africa and Europe, arriving in Europe at about 10 000 years ago. One of the out of Asia lineages also migrated back to the east, creating a series of admixed populations with the endemic Asian lineages in northern China before migrating to the New World. For the first time, our study unravels an extraordinary journey that the domestic dog has traveled on earth.

Re-evaluating data quality of dog mitochondrial, Y chromosomal, and autosomal SNPs genotyped by SNP array.

Quality deficiencies in single nucleotide polymorphism (SNP) analyses have important implications. We used missingness rates to investigate the quality of a recently published dataset containing 424 mitochondrial, 211 Y chromosomal, and 160 432 autosomal SNPs generated by a semicustom Illumina SNP array from 5 392 dogs and 14 grey wolves. Overall, the individual missingness rate for mitochondrial SNPs was ~43.8%, with 980 (18.1%) individuals completely missing mitochondrial SNP genotyping (missingness rate=1). In males, the genotype missingness rate was ~28.8% for Y chromosomal SNPs, with 374 males recording rates above 0.96. These 374 males also exhibited completely failed mitochondrial SNPs genotyping, indicative of a batch effect. Individual missingness rates for autosomal markers were greater than zero, but less than 0.5. Neither mitochondrial nor Y chromosomal SNPs achieved complete genotyping (locus missingness rate=0), whereas 5.9% of autosomal SNPs had a locus missingness rate=1. The high missingness rates and possible batch effect show that caution and rigorous measures are vital when genotyping and analyzing SNP array data for domestic animals. Further improvements of these arrays will be helpful to future studies.

DoGSD: the dog and wolf genome SNP database.

The rapid advancement of next-generation sequencing technology has generated a deluge of genomic data from domesticated dogs and their wild ancestor, grey wolves, which have simultaneously broadened our understanding of domestication and diseases that are shared by humans and dogs. To address the scarcity of single nucleotide polymorphism (SNP) data provided by authorized databases and to make SNP data more easily/friendly usable and available, we propose DoGSD (http://dogsd.big.ac.cn), the first canidae-specific database which focuses on whole genome SNP data from domesticated dogs and grey wolves. The DoGSD is a web-based, open-access resource comprising ∼ 19 million high-quality whole-genome SNPs. In addition to the dbSNP data set (build 139), DoGSD incorporates a comprehensive collection of SNPs from two newly sequenced samples (1 wolf and 1 dog) and collected SNPs from three latest dog/wolf genetic studies (7 wolves and 68 dogs), which were taken together for analysis with the population genetic statistics, Fst. In addition, DoGSD integrates some closely related information including SNP annotation, summary lists of SNPs located in genes, synonymous and non-synonymous SNPs, sampling location and breed information. All these features make DoGSD a useful resource for in-depth analysis in dog-/wolf-related studies.

Genetic convergence in the adaptation of dogs and humans to the high-altitude environment of the tibetan plateau.

The high-altitude hypoxic environment represents one of the most extreme challenges for mammals. Previous studies of humans on the Tibetan plateau and in the Andes Mountains have identified statistical signatures of selection in different sets of loci. Here, we first measured the hemoglobin levels in village dogs from Tibet and those from Chinese lowlands. We found that the hemoglobin levels are very similar between the two groups, suggesting that Tibetan dogs might share similar adaptive strategies as the Tibetan people. Through a whole-genome sequencing approach, we have identified EPAS1 and HBB as candidate genes for the hypoxic adaptation on the Tibetan plateau. The population genetic analysis shows a significant convergence between humans and dogs in Tibet. The similarities in the sets of loci that exhibit putative signatures of selection and the hemoglobin levels between humans and dogs of the same environment, but not between human populations in different regions, suggests an extraordinary landscape of convergent evolution between human beings and their best friend on the Tibetan plateau.

The genomics of selection in dogs and the parallel evolution between dogs and humans.

The genetic bases of demographic changes and artificial selection underlying domestication are of great interest in evolutionary biology. Here we perform whole-genome sequencing of multiple grey wolves, Chinese indigenous dogs and dogs of diverse breeds. Demographic analysis show that the split between wolves and Chinese indigenous dogs occurred 32,000 years ago and that the subsequent bottlenecks were mild. Therefore, dogs may have been under human selection over a much longer time than previously concluded, based on molecular data, perhaps by initially scavenging with humans. Population genetic analysis identifies a list of genes under positive selection during domestication, which overlaps extensively with the corresponding list of positively selected genes in humans. Parallel evolution is most apparent in genes for digestion and metabolism, neurological process and cancer. Our study, for the first time, draws together humans and dogs in their recent genomic evolution.