Large-scale snake genome analyses provide insights into vertebrate development.

Snakes are a remarkable squamate lineage with unique morphological adaptations, especially those related to the evolution of vertebrate skeletons, organs, and sensory systems. To clarify the genetic underpinnings of snake phenotypes, we assembled and analyzed 14 de novo genomes from 12 snake families. We also investigated the genetic basis of the morphological characteristics of snakes using functional experiments. We identified genes, regulatory elements, and structural variations that have potentially contributed to the evolution of limb loss, an elongated body plan, asymmetrical lungs, sensory systems, and digestive adaptations in snakes. We identified some of the genes and regulatory elements that might have shaped the evolution of vision, the skeletal system and diet in blind snakes, and thermoreception in infrared-sensitive snakes. Our study provides insights into the evolution and development of snakes and vertebrates.

Incomplete lineage sorting and phenotypic evolution in marsupials.

Incomplete lineage sorting (ILS) makes ancestral genetic polymorphisms persist during rapid speciation events, inducing incongruences between gene trees and species trees. ILS has complicated phylogenetic inference in many lineages, including hominids. However, we lack empirical evidence that ILS leads to incongruent phenotypic variation. Here, we performed phylogenomic analyses to show that the South American monito del monte is the sister lineage of all Australian marsupials, although over 31% of its genome is closer to the Diprotodontia than to other Australian groups due to ILS during ancient radiation. Pervasive conflicting phylogenetic signals across the whole genome are consistent with some of the morphological variation among extant marsupials. We detected hundreds of genes that experienced stochastic fixation during ILS, encoding the same amino acids in non-sister species. Using functional experiments, we confirm how ILS may have directly contributed to hemiplasy in morphological traits that were established during rapid marsupial speciation ca. 60 mya.

En masse organoid phenotyping informs metabolic-associated genetic susceptibility to NASH.

Genotype-phenotype associations for common diseases are often compounded by pleiotropy and metabolic state. Here, we devised a pooled human organoid-panel of steatohepatitis to investigate the impact of metabolic status on genotype-phenotype association. En masse population-based phenotypic analysis under insulin insensitive conditions predicted key non-alcoholic steatohepatitis (NASH)-genetic factors including the glucokinase regulatory protein (GCKR)-rs1260326:C>T. Analysis of NASH clinical cohorts revealed that GCKR-rs1260326-T allele elevates disease severity only under diabetic state but protects from fibrosis under non-diabetic states. Transcriptomic, metabolomic, and pharmacological analyses indicate significant mitochondrial dysfunction incurred by GCKR-rs1260326, which was not reversed with metformin. Uncoupling oxidative mechanisms mitigated mitochondrial dysfunction and permitted adaptation to increased fatty acid supply while protecting against oxidant stress, forming a basis for future therapeutic approaches for diabetic NASH. Thus, "in-a-dish" genotype-phenotype association strategies disentangle the opposing roles of metabolic-associated gene variant functions and offer a rich mechanistic, diagnostic, and therapeutic inference toolbox toward precision hepatology. VIDEO ABSTRACT.

Genomic architecture of autism from comprehensive whole-genome sequence annotation.

Fully understanding autism spectrum disorder (ASD) genetics requires whole-genome sequencing (WGS). We present the latest release of the Autism Speaks MSSNG resource, which includes WGS data from 5,100 individuals with ASD and 6,212 non-ASD parents and siblings (total n = 11,312). Examining a wide variety of genetic variants in MSSNG and the Simons Simplex Collection (SSC; n = 9,205), we identified ASD-associated rare variants in 718/5,100 individuals with ASD from MSSNG (14.1%) and 350/2,419 from SSC (14.5%). Considering genomic architecture, 52% were nuclear sequence-level variants, 46% were nuclear structural variants (including copy-number variants, inversions, large insertions, uniparental isodisomies, and tandem repeat expansions), and 2% were mitochondrial variants. Our study provides a guidebook for exploring genotype-phenotype correlations in families who carry ASD-associated rare variants and serves as an entry point to the expanded studies required to dissect the etiology in the ∼85% of the ASD population that remain idiopathic.

Mapping information-rich genotype-phenotype landscapes with genome-scale Perturb-seq.

A central goal of genetics is to define the relationships between genotypes and phenotypes. High-content phenotypic screens such as Perturb-seq (CRISPR-based screens with single-cell RNA-sequencing readouts) enable massively parallel functional genomic mapping but, to date, have been used at limited scales. Here, we perform genome-scale Perturb-seq targeting all expressed genes with CRISPR interference (CRISPRi) across >2.5 million human cells. We use transcriptional phenotypes to predict the function of poorly characterized genes, uncovering new regulators of ribosome biogenesis (including CCDC86, ZNF236, and SPATA5L1), transcription (C7orf26), and mitochondrial respiration (TMEM242). In addition to assigning gene function, single-cell transcriptional phenotypes allow for in-depth dissection of complex cellular phenomena-from RNA processing to differentiation. We leverage this ability to systematically identify genetic drivers and consequences of aneuploidy and to discover an unanticipated layer of stress-specific regulation of the mitochondrial genome. Our information-rich genotype-phenotype map reveals a multidimensional portrait of gene and cellular function.

A Pan-plant Protein Complex Map Reveals Deep Conservation and Novel Assemblies.

Plants are foundational for global ecological and economic systems, but most plant proteins remain uncharacterized. Protein interaction networks often suggest protein functions and open new avenues to characterize genes and proteins. We therefore systematically determined protein complexes from 13 plant species of scientific and agricultural importance, greatly expanding the known repertoire of stable protein complexes in plants. By using co-fractionation mass spectrometry, we recovered known complexes, confirmed complexes predicted to occur in plants, and identified previously unknown interactions conserved over 1.1 billion years of green plant evolution. Several novel complexes are involved in vernalization and pathogen defense, traits critical for agriculture. We also observed plant analogs of animal complexes with distinct molecular assemblies, including a megadalton-scale tRNA multi-synthetase complex. The resulting map offers a cross-species view of conserved, stable protein assemblies shared across plant cells and provides a mechanistic, biochemical framework for interpreting plant genetics and mutant phenotypes.

Senescence-Induced Vascular Remodeling Creates Therapeutic Vulnerabilities in Pancreas Cancer.

KRAS mutant pancreatic ductal adenocarcinoma (PDAC) is characterized by a desmoplastic response that promotes hypovascularity, immunosuppression, and resistance to chemo- and immunotherapies. We show that a combination of MEK and CDK4/6 inhibitors that target KRAS-directed oncogenic signaling can suppress PDAC proliferation through induction of retinoblastoma (RB) protein-mediated senescence. In preclinical mouse models of PDAC, this senescence-inducing therapy produces a senescence-associated secretory phenotype (SASP) that includes pro-angiogenic factors that promote tumor vascularization, which in turn enhances drug delivery and efficacy of cytotoxic gemcitabine chemotherapy. In addition, SASP-mediated endothelial cell activation stimulates the accumulation of CD8+ T cells into otherwise immunologically "cold" tumors, sensitizing tumors to PD-1 checkpoint blockade. Therefore, in PDAC models, therapy-induced senescence can establish emergent susceptibilities to otherwise ineffective chemo- and immunotherapies through SASP-dependent effects on the tumor vasculature and immune system.

Reconstructing Denisovan Anatomy Using DNA Methylation Maps.

Denisovans are an extinct group of humans whose morphology remains unknown. Here, we present a method for reconstructing skeletal morphology using DNA methylation patterns. Our method is based on linking unidirectional methylation changes to loss-of-function phenotypes. We tested performance by reconstructing Neanderthal and chimpanzee skeletal morphologies and obtained >85% precision in identifying divergent traits. We then applied this method to the Denisovan and offer a putative morphological profile. We suggest that Denisovans likely shared with Neanderthals traits such as an elongated face and a wide pelvis. We also identify Denisovan-derived changes, such as an increased dental arch and lateral cranial expansion. Our predictions match the only morphologically informative Denisovan bone to date, as well as the Xuchang skull, which was suggested by some to be a Denisovan. We conclude that DNA methylation can be used to reconstruct anatomical features, including some that do not survive in the fossil record.

Combinatorial Genetics Reveals a Scaling Law for the Effects of Mutations on Splicing.

Despite a wealth of molecular knowledge, quantitative laws for accurate prediction of biological phenomena remain rare. Alternative pre-mRNA splicing is an important regulated step in gene expression frequently perturbed in human disease. To understand the combined effects of mutations during evolution, we quantified the effects of all possible combinations of exonic mutations accumulated during the emergence of an alternatively spliced human exon. This revealed that mutation effects scale non-monotonically with the inclusion level of an exon, with each mutation having maximum effect at a predictable intermediate inclusion level. This scaling is observed genome-wide for cis and trans perturbations of splicing, including for natural and disease-associated variants. Mathematical modeling suggests that competition between alternative splice sites is sufficient to cause this non-linearity in the genotype-phenotype map. Combining the global scaling law with specific pairwise interactions between neighboring mutations allows accurate prediction of the effects of complex genotype changes involving >10 mutations.

Structural Alterations Driving Castration-Resistant Prostate Cancer Revealed by Linked-Read Genome Sequencing.

Nearly all prostate cancer deaths are from metastatic castration-resistant prostate cancer (mCRPC), but there have been few whole-genome sequencing (WGS) studies of this disease state. We performed linked-read WGS on 23 mCRPC biopsy specimens and analyzed cell-free DNA sequencing data from 86 patients with mCRPC. In addition to frequent rearrangements affecting known prostate cancer genes, we observed complex rearrangements of the AR locus in most cases. Unexpectedly, these rearrangements include highly recurrent tandem duplications involving an upstream enhancer of AR in 70%-87% of cases compared with <2% of primary prostate cancers. A subset of cases displayed AR or MYC enhancer duplication in the context of a genome-wide tandem duplicator phenotype associated with CDK12 inactivation. Our findings highlight the complex genomic structure of mCRPC, nominate alterations that may inform prostate cancer treatment, and suggest that additional recurrent events in the non-coding mCRPC genome remain to be discovered.

Comprehensive and Integrative Genomic Characterization of Hepatocellular Carcinoma.

Liver cancer has the second highest worldwide cancer mortality rate and has limited therapeutic options. We analyzed 363 hepatocellular carcinoma (HCC) cases by whole-exome sequencing and DNA copy number analyses, and we analyzed 196 HCC cases by DNA methylation, RNA, miRNA, and proteomic expression also. DNA sequencing and mutation analysis identified significantly mutated genes, including LZTR1, EEF1A1, SF3B1, and SMARCA4. Significant alterations by mutation or downregulation by hypermethylation in genes likely to result in HCC metabolic reprogramming (ALB, APOB, and CPS1) were observed. Integrative molecular HCC subtyping incorporating unsupervised clustering of five data platforms identified three subtypes, one of which was associated with poorer prognosis in three HCC cohorts. Integrated analyses enabled development of a p53 target gene expression signature correlating with poor survival. Potential therapeutic targets for which inhibitors exist include WNT signaling, MDM4, MET, VEGFA, MCL1, IDH1, TERT, and immune checkpoint proteins CTLA-4, PD-1, and PD-L1.

Membrane adhesion junctions regulate airway smooth muscle phenotype and function.

The local environment surrounding airway smooth muscle (ASM) cells has profound effects on the physiological and phenotypic properties of ASM tissues. ASM is continually subjected to the mechanical forces generated during breathing and to the constituents of its surrounding extracellular milieu. The smooth muscle cells within the airways continually modulate their properties to adapt to these changing environmental influences. Smooth muscle cells connect to the extracellular cell matrix (ECM) at membrane adhesion junctions that provide mechanical coupling between smooth muscle cells within the tissue. Membrane adhesion junctions also sense local environmental signals and transduce them to cytoplasmic and nuclear signaling pathways in the ASM cell. Adhesion junctions are composed of clusters of transmembrane integrin proteins that bind to ECM proteins outside the cell and to large multiprotein complexes in the submembranous cytoplasm. Physiological conditions and stimuli from the surrounding ECM are sensed by integrin proteins and transduced by submembranous adhesion complexes to signaling pathways to the cytoskeleton and nucleus. The transmission of information between the local environment of the cells and intracellular processes enables ASM cells to rapidly adapt their physiological properties to modulating influences in their extracellular environment: mechanical and physical forces that impinge on the cell, ECM constituents, local mediators, and metabolites. The structure and molecular organization of adhesion junction complexes and the actin cytoskeleton are dynamic and constantly changing in response to environmental influences. The ability of ASM to rapidly accommodate to the ever-changing conditions and fluctuating physical forces within its local environment is essential for its normal physiological function.

Distinct human Langerhans cell subsets orchestrate reciprocal functions and require different developmental regulation.

Langerhans cells (LCs) play a pivotal role in skin homeostasis, and the heterogeneity of LCs has long been considered. In this study, we have identified two steady-state (LC1 and LC2) and two activated LC subsets in the epidermis of human skin and in LCs derived from CD34+ hemopoietic stem cells (HSC-LCs) by utilizing single-cell RNA sequencing and mass cytometry. Analysis of HSC-LCs at multiple time-points during differentiation revealed that EGR1 and Notch signaling were among the top pathways regulating the bifurcation of LC1 and LC2. LC1 were characterized as classical LCs, mainly related to innate immunity and antigen processing. LC2 were similar to monocytes or myeloid dendritic cells, involving in immune responses and leukocyte activation. LC1 remained stable under inflammatory microenvironment, whereas LC2 were prone to being activated and demonstrated elevated expression of immuno-suppressive molecules. We revealed distinct human LC subsets that require different developmental regulation and orchestrate reciprocal functions.

Genotypes, Networks, Phenotypes: Moving Toward Plant Systems Genetics.

One of the central goals in biology is to understand how and how much of the phenotype of an organism is encoded in its genome. Although many genes that are crucial for organismal processes have been identified, much less is known about the genetic bases underlying quantitative phenotypic differences in natural populations. We discuss the fundamental gap between the large body of knowledge generated over the past decades by experimental genetics in the laboratory and what is needed to understand the genotype-to-phenotype problem on a broader scale. We argue that systems genetics, a combination of systems biology and the study of natural variation using quantitative genetics, will help to address this problem. We present major advances in these two mostly disconnected areas that have increased our understanding of the developmental processes of flowering time control and root growth. We conclude by illustrating and discussing the efforts that have been made toward systems genetics specifically in plants.

A FBXO7/EYA2-SCFFBXW7 axis promotes AXL-mediated maintenance of mesenchymal and immune evasion phenotypes of cancer cells.

A mesenchymal tumor phenotype associates with immunotherapy resistance, although the mechanism is unclear. Here, we identified FBXO7 as a maintenance regulator of mesenchymal and immune evasion phenotypes of cancer cells. FBXO7 bound and stabilized SIX1 co-transcriptional regulator EYA2, stimulating mesenchymal gene expression and suppressing IFNα/ß, chemokines CXCL9/10, and antigen presentation machinery, driven by AXL extracellular ligand GAS6. Ubiquitin ligase SCFFBXW7 antagonized this pathway by promoting EYA2 degradation. Targeting EYA2 Tyr phosphatase activity decreased mesenchymal phenotypes and enhanced cancer cell immunogenicity, resulting in attenuated tumor growth and metastasis, increased infiltration of cytotoxic T and NK cells, and enhanced anti-PD-1 therapy response in mouse tumor models. FBXO7 expression correlated with mesenchymal and immune-suppressive signatures in patients with cancer. An FBXO7-immune gene signature predicted immunotherapy responses. Collectively, the FBXO7/EYA2-SCFFBXW7 axis maintains mesenchymal and immune evasion phenotypes of cancer cells, providing rationale to evaluate FBXO7/EYA2 inhibitors in combination with immune-based therapies to enhance onco-immunotherapy responses.

Carbonic anhydrases in metazoan model organisms: molecules, mechanisms, and physiology.

During the past three decades, mice, zebrafish, fruit flies, and Caenorhabditis elegans have been the primary model organisms used for the study of various biological phenomena. These models have also been adopted and developed to investigate the physiological roles of carbonic anhydrases (CAs) and carbonic anhydrase-related proteins (CARPs). These proteins belong to eight CA families and are identified by Greek letters: α, ß, γ, δ, ζ, η, θ, and ι. Studies using model organisms have focused on two CA families, α-CAs and ß-CAs, which are expressed in both prokaryotic and eukaryotic organisms with species-specific distribution patterns and unique functions. This review covers the biological roles of CAs and CARPs in light of investigations performed in model organisms. Functional studies demonstrate that CAs are not only linked to the regulation of pH homeostasis, the classical role of CAs, but also contribute to a plethora of previously undescribed functions.

Genotype-Phenotype Relationships in the Context of Transcriptional Adaptation and Genetic Robustness.

Genetic manipulations with a robust and predictable outcome are critical to investigate gene function, as well as for therapeutic genome engineering. For many years, knockdown approaches and reagents including RNA interference and antisense oligonucleotides dominated functional studies; however, with the advent of precise genome editing technologies, CRISPR-based knockout systems have become the state-of-the-art tools for such studies. These technologies have helped decipher the role of thousands of genes in development and disease. Their use has also revealed how limited our understanding of genotype-phenotype relationships is. The recent discovery that certain mutations can trigger the transcriptional modulation of other genes, a phenomenon called transcriptional adaptation, has provided an additional explanation for the contradicting phenotypes observed in knockdown versus knockout models and increased awareness about the use of each of these approaches. In this review, we first cover the strengths and limitations of different gene perturbation strategies. Then we highlight the diverse ways in which the genotype-phenotype relationship can be discordant between these different strategies. Finally, we review the genetic robustness mechanisms that can lead to such discrepancies, paying special attention to the recently discovered phenomenon of transcriptional adaptation.

Adoptive T Cell Therapy Targeting Different Gene Products Reveals Diverse and Context-Dependent Immune Evasion in Melanoma.

Tumor immune escape limits durable responses to T cell therapy. Here, we examined how regulation and function of gene products that provide the target epitopes for CD8+ T cell anti-tumor immunity influence therapeutic efficacy and resistance. We used a CRISPR-Cas9-based method (CRISPitope) in syngeneic melanoma models to fuse the same model CD8+ T cell epitope to the C-termini of different endogenous gene products. Targeting melanosomal proteins or oncogenic CDK4R24C (Cyclin-dependent kinase 4) by adoptive cell transfer (ACT) of the same epitope-specific CD8+ T cells revealed diverse genetic and non-genetic immune escape mechanisms. ACT directed against melanosomal proteins, but not CDK4R24C, promoted melanoma dedifferentiation, and increased myeloid cell infiltration. CDK4R24C antigen persistence was associated with an interferon-high and T-cell-rich tumor microenvironment, allowing for immune checkpoint inhibition as salvage therapy. Thus, the choice of target antigen determines the phenotype and immune contexture of recurrent melanomas, with implications to the design of cancer immunotherapies.

K128 ubiquitination constrains RAS activity by expanding its binding interface with GAP proteins.

The RAS pathway is among the most frequently activated signaling nodes in cancer. However, the mechanisms that alter RAS activity in human pathologies are not entirely understood. The most prevalent post-translational modification within the GTPase core domain of NRAS and KRAS is ubiquitination at lysine 128 (K128), which is significantly decreased in cancer samples compared to normal tissue. Here, we found that K128 ubiquitination creates an additional binding interface for RAS GTPase-activating proteins (GAPs), NF1 and RASA1, thus increasing RAS binding to GAP proteins and promoting GAP-mediated GTP hydrolysis. Stimulation of cultured cancer cells with growth factors or cytokines transiently induces K128 ubiquitination and restricts the extent of wild-type RAS activation in a GAP-dependent manner. In KRAS mutant cells, K128 ubiquitination limits tumor growth by restricting RAL/ TBK1 signaling and negatively regulating the autocrine circuit induced by mutant KRAS. Reduction of K128 ubiquitination activates both wild-type and mutant RAS signaling and elicits a senescence-associated secretory phenotype, promoting RAS-driven pancreatic tumorigenesis.

What Has a Century of Quantitative Genetics Taught Us About Nature's Genetic Tool Kit?

The complexity of heredity has been appreciated for decades: Many traits are controlled not by a single genetic locus but instead by polymorphisms throughout the genome. The importance of complex traits in biology and medicine has motivated diverse approaches to understanding their detailed genetic bases. Here, we focus on recent systematic studies, many in budding yeast, which have revealed that large numbers of all kinds of molecular variation, from noncoding to synonymous variants, can make significant contributions to phenotype. Variants can affect different traits in opposing directions, and their contributions can be modified by both the environment and the epigenetic state of the cell. The integration of prospective (synthesizing and analyzing variants) and retrospective (examining standing variation) approaches promises to reveal how natural selection shapes quantitative traits. Only by comprehensively understanding nature's genetic tool kit can we predict how phenotypes arise from the complex ensembles of genetic variants in living organisms.