Optimal Rates for Phylogenetic Inference and Experimental Design in the Era of Genome-Scale Data Sets

With the rise of genome-scale data sets, there has been a call for increased data scrutiny and careful selection of loci that are appropriate to use in an attempt to resolve a phylogenetic problem. Such loci should maximize phylogenetic information content while minimizing the risk of homoplasy. Theory posits the existence of characters that evolve at an optimum rate, and efforts to determine optimal rates of inference have been a cornerstone of phylogenetic experimental design for over two decades. However, both theoretical and empirical investigations of optimal rates have varied dramatically in their conclusions: spanning no relationship to a tight relationship between the rate of change and phylogenetic utility. Herein, we synthesize these apparently contradictory views, demonstrating both empirical and theoretical conditions under which each is correct. We find that optimal rates of characters-not genes-are generally robust to most experimental design decisions. Moreover, consideration of site rate heterogeneity within a given locus is critical to accurate predictions of utility. Factors such as taxon sampling or the targeted number of characters providing support for a topology are additionally critical to the predictions of phylogenetic utility based on the rate of character change. Further, optimality of rates and predictions of phylogenetic utility are not equivalent, demonstrating the need for further development of comprehensive theory of phylogenetic experimental design. [Divergence time; GC bias; homoplasy; incongruence; information content; internode length; optimal rates; phylogenetic informativeness; phylogenetic theory; phylogenetic utility; phylogenomics; signal and noise; subtending branch length; state space; taxon and character sampling.].

Genomic Comparison Among Global Isolates of L. interrogans Serovars Copenhageni and Icterohaemorrhagiae Identified Natural Genetic Variation Caused by an Indel

Leptospirosis is a worldwide zoonosis, responsible for more than 1 million cases and 60,000 deaths every year. Among the 13 pathogenic species of the genus Leptospira, serovars belonging to L. interrogans serogroup Icterohaemorrhagiae are considered to be the most virulent strains, and responsible for majority of the reported severe cases. Serovars Copenhageni and Icterohaemorrhagiae are major representatives of this serogroup and despite their public health relevance, little is known regarding the genetic differences between these two serovars. In this study, we analyzed the genome sequences of 67 isolates belonging to L. interrogans serovars Copenhageni and Icterohaemorrhagiae to investigate the influence of spatial and temporal variations on DNA sequence diversity. Out of the 1072 SNPs identified, 276 were in non-coding regions and 796 in coding regions. Indel analyses identified 258 indels, out of which 191 were found in coding regions and 67 in non-coding regions. Our phylogenetic analyses based on SNP dataset revealed that both serovars are closely related but showed distinct spatial clustering. However, likelihood ratio test of the indel data statistically confirmed the presence of a frameshift mutation within a homopolymeric tract of lic12008 gene (related to LPS biosynthesis) in all the L. interrogans serovar Icterohaemorrhagiae strains but not in the Copenhageni strains. Therefore, this internal indel identified can genetically distinguish L. interrogans serovar Copenhageni from serovar Icterohaemorrhagiae with high discriminatory power. To our knowledge, this is the first study to identify global sequence variations (SNPs and Indels) in L. interrogans serovars Copenhageni and Icterohaemorrhagiae.

Potential effectiveness of long-acting injectable pre-exposure prophylaxis for HIV prevention in men who have sex with men: a modelling study

BACKGROUND

Oral pre-exposure prophylaxis (PrEP) prevents HIV infection in men who have sex with men (MSM); however, adherence is an ongoing concern. Long-acting injectable PrEP is being tested in phase 3 trials and could address challenges associated with adherence. We examined the potential effectiveness of long-acting injectable PrEP compared with oral PrEP in MSM.

METHODS

We used an agent-based model to simulate HIV transmission in a dynamic network of 11 245 MSM in Atlanta, GA, USA. We used raw data from studies in macaque models and pharmacokinetic data from safety trials to estimate the time-varying efficacy of long-acting injectable PrEP. The effect of long-acting injectable PrEP on the cumulative number of new HIV infections over 10 years (2015-24) was compared with no PrEP and daily oral PrEP across a range of coverage levels. Sensitivity analyses were done with varying maximum efficacy and drug half-life values.

FINDINGS

In the absence of PrEP, the model predicted 2374 new HIV infections (95% simulation interval [SI] 2345-2412) between 2015 and 2024. The cumulative number of new HIV infections was reduced in all scenarios in which MSM received long-acting injectable PrEP compared with oral PrEP. At a coverage level of 35%, compared with no PrEP, long-acting injectable PrEP led to a 44% reduction in new HIV infections (1044 new infections averted [95% SI 1018-1077]) versus 33% (792 infections averted [763-821]) for oral PrEP. The relative benefit of long-acting injectable PrEP was sensitive to the assumed efficacy of injections received every 8 weeks, discontinuation rates, and terminal drug half-life.

INTERPRETATION

Long-acting injectable PrEP has the potential to produce larger reductions in HIV transmission in MSM than oral PrEP. However, the real-world, population-level impact of this approach will depend on uptake of this prevention method and its effectiveness, as well as retention of patients in clinical care.

FUNDING

National Institute on Drug Abuse and National Institute of Mental Health.

Neutral Theory and the Somatic Evolution of Cancer

Kimura's neutral theory argued that positive selection was not responsible for an appreciable fraction of molecular substitutions. Correspondingly, quantitative analysis reveals that the vast majority of substitutions in cancer genomes are not detectably under selection. Insights from the somatic evolution of cancer reveal that beneficial substitutions in cancer constitute a small but important fraction of the molecular variants. The molecular evolution of cancer community will benefit by incorporating the neutral theory of molecular evolution into their understanding and analysis of cancer evolution-and accepting the use of tractable, predictive models, even when there is some evidence that they are not perfect.

Mutation profiles in early-stage lung squamous cell carcinoma with clinical follow-up and correlation with markers of immune function

Background

Lung squamous cell carcinoma (LUSC) accounts for 20–30% of non-small cell lung cancers (NSCLCs). There are limited treatment strategies for LUSC in part due to our inadequate understanding of the molecular underpinnings of the disease. We performed whole-exome sequencing (WES) and comprehensive immune profiling of a unique set of clinically annotated early-stage LUSCs to increase our understanding of the pathobiology of this malignancy.

Methods

Matched pairs of surgically resected stage I-III LUSCs and normal lung tissues (n = 108) were analyzed by WES. Immunohistochemistry and image analysis-based profiling of 10 immune markers were done on a subset of LUSCs (n = 91). Associations among mutations, immune markers and clinicopathological variables were statistically examined using analysis of variance and Fisher’s exact test. Cox proportional hazards regression models were used for statistical analysis of clinical outcome.

Results

This early-stage LUSC cohort displayed an average of 209 exonic mutations per tumor. Fourteen genes exhibited significant enrichment for somatic mutation: TP53, MLL2, PIK3CA, NFE2L2, CDH8, KEAP1, PTEN, ADCY8, PTPRT, CALCR, GRM8, FBXW7, RB1 and CDKN2A. Among mutated genes associated with poor recurrence-free survival, MLL2 mutations predicted poor prognosis in both TP53 mutant and wild-type LUSCs. We also found that in treated patients, FBXW7 and KEAP1 mutations were associated with poor response to adjuvant therapy, particularly in TP53-mutant tumors. Analysis of mutations with immune markers revealed that ADCY8 and PIK3CA mutations were associated with markedly decreased tumoral PD-L1 expression, LUSCs with PIK3CA mutations exhibited elevated CD45ro levels and CDKN2A-mutant tumors displayed an up-regulated immune response.

Conclusion(s)

Our findings pinpoint mutated genes that may impact clinical outcome as well as personalized strategies for targeted immunotherapies in early-stage LUSC.

Analysis of mutation, selection, and epistasis: an informed approach to cancer clinical trials

Currently, drug development efforts and clinical trials to test them are often prioritized by targeting genes with high frequencies of somatic variants among tumors. However, differences in oncogenic mutation rate-not necessarily the effect the variant has on tumor growth-contribute enormously to somatic variant frequency. We argue that decoupling the contributions of mutation and cancer lineage selection to the frequency of somatic variants among tumors is critical to understanding-and predicting-the therapeutic potential of different interventions. To provide an indicator of that strength of selection and therapeutic potential, the frequency at which we observe a given variant across patients must be modulated by our expectation given the mutation rate and target size to provide an indicator of that strength of selection and therapeutic potential. Additionally, antagonistic and synergistic epistasis among mutations also impacts the potential therapeutic benefit of targeted drug development. Quantitative approaches should be fostered that use the known genetic architectures of cancer types, decouple mutation rate, and provide rigorous guidance regarding investment in targeted drug development. By integrating evolutionary principles and detailed mechanistic knowledge into those approaches, we can maximize our ability to identify those targeted therapies most likely to yield substantial clinical benefit.

Evidence for Extending Anomalous Miocene Volcanism at the Edge of the East Antarctic Craton

Using field observations followed by petrological, geochemical, geochronological, and geophysical data we infer the presence of a previously unknown Miocene subglacial volcanic center ~230 km from the South Pole. Evidence of volcanism is from boulders of olivine-bearing amygdaloidal/vesicular basalt and hyaloclastite deposited in a moraine in the southern Transantarctic Mountains. ⁴⁰Ar/³⁹Ar ages from five specimens plus U-Pb ages of detrital zircon from glacial till indicate igneous activity 25-17 Ma. The likely source of the volcanism is a circular -735 nT magnetic anomaly 60 km upflow from the sampling site. Subaqueous textures of the volcanics indicate eruption beneath ice or into water at the margin of an ice mass during the early Miocene. These rocks record the southernmost Cenozoic volcanism in Antarctica and expand the known extent of the oldest lavas associated with West Antarctic rift system. They may be an expression of lithospheric foundering beneath the southern Transantarctic Mountains.

CDKN2A Copy Number Loss Is an Independent Prognostic Factor in HPV-Negative Head and Neck Squamous Cell Carcinoma

Background

HPV infection is associated with high p16 expression and good prognosis in head and neck squamous cell carcinomas (HNSCCs). Analysis of CDKN2A, the gene encoding p16, may further elucidate the association between p16 expression and prognosis. We sought to determine whether CDKN2A copy number loss was associated with poor survival in HPV-negative HNSCCs.

Methods

The Cancer Genome Atlas HNSCC clinical and genomic data were obtained and integrated. Patients <80 years old with a primary tumor in the oral cavity, oropharynx, hypopharynx, or larynx were included. Stratifying by copy number loss status, CDKN2A mRNA and p16 protein expression levels were examined and overall survival (OS) and disease-free survival (DFS) were evaluated.

Results

401 patients with HPV-negative HNSCC were identified. 146 patients demonstrated CDKN2A copy number loss. The CDKN2A copy number loss group expressed significantly lower levels of CDKN2A mRNA and p16 protein than did the non-copy number loss group. Median OS for patients with and without CDKN2A copy number loss was 16.5 and 46.6 months, respectively (p = 0.007). Median DFS for both groups was 11.6 and 19.2 months, respectively (p = 0.03). In both univariate and multivariable analyses, stage IV designation, receipt of chemotherapy and CDKN2A copy number loss were predictive of OS.

Conclusion

CDKN2A copy number loss predicted poor survival independently of other patient and treatment factors and may be a clinically useful prognostic factor.

Detection of Regional Variation in Selection Intensity within Protein-Coding Genes Using DNA Sequence Polymorphism and Divergence

Numerous approaches have been developed to infer natural selection based on the comparison of polymorphism within species and divergence between species. These methods are especially powerful for the detection of uniform selection operating across a gene. However, empirical analyses have demonstrated that regions of protein-coding genes exhibiting clusters of amino acid substitutions are subject to different levels of selection relative to other regions of the same gene. To quantify this heterogeneity of selection within coding sequences, we developed Model Averaged Site Selection via Poisson Random Field (MASS-PRF). MASS-PRF identifies an ensemble of intragenic clustering models for polymorphic and divergent sites. This ensemble of models is used within the Poisson Random Field framework to estimate selection intensity on a site-by-site basis. Using simulations, we demonstrate that MASS-PRF has high power to detect clusters of amino acid variants in small genic regions, can reliably estimate the probability of a variant occurring at each nucleotide site in sequence data and is robust to historical demographic trends and recombination. We applied MASS-PRF to human gene polymorphism derived from the 1,000 Genomes Project and divergence data from the common chimpanzee. On the basis of this analysis, we discovered striking regional variation in selection intensity, indicative of positive or negative selection, in well-defined domains of genes that have previously been associated with neurological processing, immunity, and reproduction. We suggest that amino acid-altering substitutions within these regions likely are or have been selectively advantageous in the human lineage, playing important roles in protein function.

Primary hepatic neuroendocrine carcinoma: report of two cases and literature review

Background

Primary hepatic neuroendocrine carcinoma (PHNEC) is extremely rare. The diagnosis of PHNEC remains challenging—partly due to its rarity, and partly due to its lack of unique clinical features. Available treatment options for PHNEC include surgical resection of the liver tumor(s), radiotherapy, liver transplant, transcatheter arterial chemoembolization (TACE), and administration of somatostatin analogues.

Case presentation

We report two male PHNEC cases and discuss the diagnosis and treatment options. Both cases presented with abdominal pain; case two also presented with symptoms of jaundice. The initial diagnosis for both cases was poorly differentiated grade 3 small-cell neuroendocrine carcinoma, based on imaging characteristics and the pathology of liver biopsies. Final diagnoses of PHNEC were arrived at by ruling out non-hepatic origins. Case one presented with a large tumor in the right liver lobe, and the patient was treated with TACE. Case two presented with tumors in both liver lobes, invasions into the left branch of hepatic portal vein, and metastasis in the hepatic hilar lymph node. This patient was ineligible for TACE and was allergic to the somatostatin analogue octreotide. This limited treatment options to supportive therapies such as albumin supplementation for liver protection. Patient one and two died at 61 and 109 days, respectively, following initial hospital admission.

Conclusions

We diagnosed both cases with poorly differentiated grade 3 small-cell PHNEC through imaging characteristics, immunohistochemical staining of liver biopsies, and examinations to eliminate non-hepatic origins. Neither TACE nor liver protection appeared to significantly extend survival time of the two patients, suggesting these treatments may be inadequate to improve survival of patients with poorly differentiated grade 3 small-cell PHNEC. The prognosis of poorly differentiated grade 3 small-cell PHNEC is poor due to limited and ineffective treatment options.

Electronic supplementary material

The online version of this article (10.1186/s12907-018-0070-7) contains supplementary material, which is available to authorized users.

Lvr, a Signaling System That Controls Global Gene Regulation and Virulence in Pathogenic Leptospira

Leptospirosis is an emerging zoonotic disease with more than 1 million cases annually. Currently there is lack of evidence for signaling pathways involved during the infection process of Leptospira. In our comprehensive genomic analysis of 20 Leptospira spp. we identified seven pathogen-specific Two-Component System (TCS) proteins. Disruption of two these TCS genes in pathogenic Leptospira strain resulted in loss-of-virulence in a hamster model of leptospirosis. Corresponding genes lvrA and lvrB (leptospira virulence regulator) are juxtaposed in an operon and are predicted to encode a hybrid histidine kinase and a hybrid response regulator, respectively. Transcriptome analysis of lvr mutant strains with disruption of one (lvrB) or both genes (lvrA/B) revealed global transcriptional regulation of 850 differentially expressed genes. Phosphotransfer assays demonstrated that LvrA phosphorylates LvrB and predicted further signaling downstream to one or more DNA-binding response regulators, suggesting that it is a branched pathway. Phylogenetic analyses indicated that lvrA and lvrB evolved independently within different ecological lineages in Leptospira via gene duplication. This study uncovers a novel-signaling pathway that regulates virulence in pathogenic Leptospira (Lvr), providing a framework to understand the molecular bases of regulation in this life-threatening bacterium.

Heterogeneity and mutation in KRAS and associated oncogenes: evaluating the potential for the evolution of resistance to targeting of KRAS G12C

Activating mutations in RAS genes are associated with approximately 20% of all human cancers. New targeted therapies show preclinical promise in inhibiting the KRAS G12C variant. However, concerns exist regarding the effectiveness of such therapies in vivo given the possibilities of existing intratumor heterogeneity or de novo mutation leading to treatment resistance. We performed deep sequencing of 27 KRAS G12-positive lung tumors to determine the prevalence of other oncogenic mutations within KRAS or within commonly mutated downstream genes that could confer resistance at the time of treatment. We also passaged patient-derived xenografts to assess the potential for novel KRAS mutation to arise during subsequent tumor evolution. Furthermore, we estimate the de novo mutation rate in KRAS position 12 and in genes downstream of KRAS. Finally, we present an approach for estimation of the selection intensity for these point mutations that explains their high prevalence in tumors. We find no evidence of heterogeneity that may compromise KRAS G12C targeted therapy within sequenced lung tumors or passaged xenografts. We find that mutations that confer resistance are even less likely to occur downstream of KRAS than to occur within KRAS. Our approach predicts that BRAF V600E would provide the highest fitness advantage for de novo-resistant subclones. Overall, our findings suggest that resistance to targeted therapy of KRAS G12C-positive tumors is unlikely to be present at the time of treatment and, among the de novo mutations likely to confer resistance, mutations in BRAF, a gene with targeted inhibitors presently available, result in subclones with the highest fitness advantage.

Somatic evolutionary timings of driver mutations

BACKGROUND

A unified analysis of DNA sequences from hundreds of tumors concluded that the driver mutations primarily occur in the earliest stages of cancer formation, with relatively few driver mutation events detected in the late-arising subclones. However, emerging evidence from the sequencing of multiple tumors and tumor regions per individual suggests that late-arising subclones with additional driver mutations are underestimated in single-sample analyses.

METHODS

To test whether driver mutations generally map to early tumor development, we examined multi-regional tumor sequencing data from 101 individuals reported in 11 published studies. Following previous studies, we annotated mutations as early-arising when all tumors/regions had those mutations (ubiquitous). We then inferred the fraction of mutations occurring early and compared it with late-arising mutations that were found in only single tumors/regions.

RESULTS

While a large fraction of driver mutations in tumors occurred relatively early in cancers, later driver mutations occurred at least as frequently as the early drivers in a substantial number of patients. This result was robust to many different approaches to annotate driver mutations. The relative frequency of early and late driver mutations varied among patients of the same cancer type and in different cancer types. We found that previous reports of the preponderance of early driver mutations were primarily informed by analysis of single tumor variant allele profiles, with which it is challenging to clearly distinguish between early and late drivers.

CONCLUSIONS

The origin and preponderance of new driver mutations are not limited to early stages of tumor evolution, with different tumors and regions showing distinct driver mutations and, consequently, distinct characteristics. Therefore, tumors with extensive intratumor heterogeneity appear to have many newly acquired drivers.

Eye patches: Protein assembly of index-gradient squid lenses

A parabolic relationship between lens radius and refractive index allows spherical lenses to avoid spherical aberration. We show that in squid, patchy colloidal physics resulted from an evolutionary radiation of globular S-crystallin proteins. Small-angle x-ray scattering experiments on lens tissue show colloidal gels of S-crystallins at all radial positions. Sparse lens materials form via low-valence linkages between disordered loops protruding from the protein surface. The loops are polydisperse and bind via a set of hydrogen bonds between disordered side chains. Peripheral lens regions with low particle valence form stable, volume-spanning gels at low density, whereas central regions with higher average valence gel at higher densities. The proteins demonstrate an evolved set of linkers for self-assembly of nanoparticles into volumetric materials.

Light sensing by opsins and fungal ecology: NOP-1 modulates entry into sexual reproduction in response to environmental cues

Understanding the genetic basis of the switch from asexual to sexual lifestyles in response to sometimes rapid environmental changes is one of the major challenges in fungal ecology. Light appears to play a critical role in the asexual-sexual switch-but fungal genomes harbour diverse light sensors. Fungal opsins are homologous to bacterial green-light-sensory rhodopsins, and their organismal functions in fungi have not been well understood. Three of these opsin-like proteins were widely distributed across fungal genomes, but homologs of the Fusarium opsin-like protein CarO were present only in plant-associated fungi. Key amino acids, including potential retinal binding sites, functionally diverged on the phylogeny of opsins. This diversification of opsin-like proteins could be correlated with life history-associated differences among fungi in their expression and function during morphological development. In Neurospora crassa and related species, knockout of the opsin NOP-1 led to a phenotype in the regulation of the asexual-sexual switch, modulating response to both light and oxygen conditions. Sexual development commenced early in ∆nop-1 strains cultured in unsealed plates under constant blue and white light. Furthermore, comparative transcriptomics showed that the expression of nop-1 is light-dependent and that the ∆nop-1 strain abundantly expresses genes involved in oxidative stress response, genes enriched in NAD/NADP binding sites, genes with functions in proton transmembrane movement and catalase activity, and genes involved in the homeostasis of protons. Based on these observations, we contend that light and oxidative stress regulate the switch via light-responsive and ROS pathways in model fungus N. crassa and other fungi.

Maximizing Power in Phylogenetics and Phylogenomics: A Perspective Illuminated by Fungal Big Data

Since its original inception over 150 years ago by Darwin, we have made tremendous progress toward the reconstruction of the Tree of Life. In particular, the transition from analyzing datasets comprised of small numbers of loci to those comprised of hundreds of loci, if not entire genomes, has aided in resolving some of the most vexing of evolutionary problems while giving us a new perspective on biodiversity. Correspondingly, phylogenetic trees have taken a central role in fields that span ecology, conservation, and medicine. However, the rise of big data has also presented phylogenomicists with a new set of challenges to experimental design, quantitative analyses, and computation. The sequencing of a number of very first genomes presented significant challenges to phylogenetic inference, leading fungal phylogenomicists to begin addressing pitfalls and postulating solutions to the issues that arise from genome-scale analyses relevant to any lineage across the Tree of Life. Here we highlight insights from fungal phylogenomics for topics including systematics and species delimitation, ecological and phenotypic diversification, and biogeography while providing an overview of progress made on the reconstruction of the fungal Tree of Life. Finally, we provide a review of considerations to phylogenomic experimental design for robust tree inference. We hope that this special issue of Advances in Genetics not only excites the continued progress of fungal evolutionary biology but also motivates the interdisciplinary development of new theory and methods designed to maximize the power of genomic scale data in phylogenetic analyses.

The ancestral levels of transcription and the evolution of sexual phenotypes in filamentous fungi

Changes in gene expression have been hypothesized to play an important role in the evolution of divergent morphologies. To test this hypothesis in a model system, we examined differences in fruiting body morphology of five filamentous fungi in the Sordariomycetes, culturing them in a common garden environment and profiling genome-wide gene expression at five developmental stages. We reconstructed ancestral gene expression phenotypes, identifying genes with the largest evolved increases in gene expression across development. Conducting knockouts and performing phenotypic analysis in two divergent species typically demonstrated altered fruiting body development in the species that had evolved increased expression. Our evolutionary approach to finding relevant genes proved far more efficient than other gene deletion studies targeting whole genomes or gene families. Combining gene expression measurements with knockout phenotypes facilitated the refinement of Bayesian networks of the genes underlying fruiting body development, regulation of which is one of the least understood processes of multicellular development.

CDKN2A copy number loss in HPV- and HPV+ head and neck cancer to indicate poor prognosis: An integrated genomic and clinical TCGA analysis

6060

Background: HPV infection is associated with high p16 expression and relatively good prognosis in head and neck cancers. Analysis of CDKN2A, the gene that encodes the p16 tumor suppressor protein, may further elucidate the association between HPV status and prognosis in head and neck squamous cell carcinomas (HNSCCs). We aimed to identify whether CDKN2A copy number loss was associated with poor survival in HNSCCs stratified by HPV status. Methods: We analyzed The Cancer Genome Atlas (TCGA) head and neck cancer data, integrating genomic measurements with clinical metadata. Patients 85 years old or younger with a primary tumor in the oral cavity, oropharynx, hypopharynx, or larynx were included. Defining CDKN2A copy number loss as a relative log2 copy number ratio < −0.6, CDKN2A mRNA and p16 protein expression levels were compared to confirm significant differences in gene transcription and translation between the copy number loss and non-copy number loss patient groups. Overall survival (OS) and disease-free survival (DFS) were evaluated to characterize prognostic differences between genomic groups. Results: 397 patients negative for HPV (HPV−) and 91 patients positive for HPV (HPV+) HNSCC were identified. 139 HPV− patients and 9 HPV+ patients demonstrated CDKN2A copy number loss. The CDKN2A copy number loss group expressed significantly lower levels of CDKN2A mRNA and p16 protein than did the non-copy number loss group in both HPV+ and HPV− disease. Median OS for HPV− patients with and without CDKN2A copy number loss was 21.8 months and 46.0 months (P = 0.02). Median DFS was 12.0 and 19.4 months respectively (P < 0.05). Median OS for HPV+ patients with and without CDKN2A copy number loss was 12.7 months and 57.4 months (P = 0.004) and median DFS was 7.0 and 36.6 months respectively (P = 0.02). Conclusions: CDKN2A copy number loss was associated with low CDKN2A mRNA and p16 protein expression, with poor prognosis in terms of disease-free and overall survival.

Whole-exome sequencing and immune profiling of early-stage lung adenocarcinoma with fully annotated clinical follow-up

Background

Lung adenocarcinomas (LUADs) lead to the majority of deaths attributable to lung cancer. We performed whole-exome sequencing (WES) and immune profiling analyses of a unique set of clinically annotated early-stage LUADs to better understand the pathogenesis of this disease and identify clinically relevant molecular markers.

Methods

We performed WES of 108 paired stage I-III LUADs and normal lung tissues using the Illumina HiSeq 2000 platform. Ten immune markers (PD-L1, PD-1, CD3, CD4, CD8, CD45ro, CD57, CD68, FOXP3 and Granzyme B) were profiled by imaging-based immunohistochemistry (IHC) in a subset of LUADs (n = 92). Associations among mutations, immune markers and clinicopathological variables were analyzed using ANOVA and Fisher's exact test. Cox proportional hazards regression models were used for multivariate analysis of clinical outcome.

Results

LUADs in this cohort exhibited an average of 243 coding mutations. We identified 28 genes with significant enrichment for mutation. SETD2-mutated LUADs exhibited relatively poor recurrence- free survival (RFS) and mutations in STK11 and ATM were associated with poor RFS among KRAS-mutant tumors. EGFR, KEAP1 and PIK3CA mutations were predictive of poor response to adjuvant therapy. Immune marker analysis revealed that LUADs in smokers and with relatively high mutation burdens exhibited increased levels of immune markers. Analysis of immunophenotypes revealed that LUADs with STK11 mutations exhibited relatively low levels of infiltrating CD4+/CD8+ T-cells indicative of a muted immune response. Tumoral PD-L1 was significantly elevated in TP53 mutant LUADs whereas PIK3CA mutant LUADs exhibited markedly down-regulated PD-L1 expression. LUADs with TP53 or KEAP1 mutations displayed relatively increased CD57 and Granzyme B levels indicative of augmented natural killer (NK) cell infiltration.

Conclusion(s)

Our study highlights molecular and immune phenotypes that warrant further analysis for their roles in clinical outcomes and personalized immune-based therapy of LUAD.

New insights on the sister lineage of percomorph fishes with an anchored hybrid enrichment dataset

Percomorph fishes represent over 17,100 species, including several model organisms and species of economic importance. Despite continuous advances in the resolution of the percomorph Tree of Life, resolution of the sister lineage to Percomorpha remains inconsistent but restricted to a small number of candidate lineages. Here we use an anchored hybrid enrichment (AHE) dataset of 132 loci with over 99,000 base pairs to identify the sister lineage of percomorph fishes. Initial analyses of this dataset failed to recover a strongly supported sister clade to Percomorpha, however, scrutiny of the AHE dataset revealed a bias towards high GC content at fast-evolving codon partitions (GC bias). By combining several existing approaches aimed at mitigating the impacts of convergence in GC bias, including RY coding and analyses of amino acids, we consistently recovered a strongly supported clade comprised of Holocentridae (squirrelfishes), Berycidae (Alfonsinos), Melamphaidae (bigscale fishes), Cetomimidae (flabby whalefishes), and Rondeletiidae (redmouth whalefishes) as the sister lineage to Percomorpha. Additionally, implementing phylogenetic informativeness (PI) based metrics as a filtration method yielded this same topology, suggesting PI based approaches will preferentially filter these fast-evolving regions and act in a manner consistent with other phylogenetic approaches aimed at mitigating GC bias. Our results provide a new perspective on a key issue for studies investigating the evolutionary history of more than one quarter of all living species of vertebrates.