Integrated Proteogenomic Characterization of Clear Cell Renal Cell Carcinoma.

To elucidate the deregulated functional modules that drive clear cell renal cell carcinoma (ccRCC), we performed comprehensive genomic, epigenomic, transcriptomic, proteomic, and phosphoproteomic characterization of treatment-naive ccRCC and paired normal adjacent tissue samples. Genomic analyses identified a distinct molecular subgroup associated with genomic instability. Integration of proteogenomic measurements uniquely identified protein dysregulation of cellular mechanisms impacted by genomic alterations, including oxidative phosphorylation-related metabolism, protein translation processes, and phospho-signaling modules. To assess the degree of immune infiltration in individual tumors, we identified microenvironment cell signatures that delineated four immune-based ccRCC subtypes characterized by distinct cellular pathways. This study reports a large-scale proteogenomic analysis of ccRCC to discern the functional impact of genomic alterations and provides evidence for rational treatment selection stemming from ccRCC pathobiology.

Differentiation route determines the functional outputs of adult megakaryopoiesis.

Emerging evidence has revealed a direct differentiation route from hematopoietic stem cells to megakaryocytes (direct route), in addition to the classical differentiation route through a series of restricted hematopoietic progenitors (stepwise route). This raises the question of the importance of two alternative routes for megakaryopoiesis. Here, we developed fate-mapping systems to distinguish the two routes, comparing their quantitative and functional outputs. We found that megakaryocytes were produced through the two routes with comparable kinetics and quantity under homeostasis. Single-cell RNA sequencing of the fate-mapped megakaryocytes revealed that the direct and stepwise routes contributed to the niche-supporting and immune megakaryocytes, respectively, but contributed to the platelet-producing megakaryocytes together. Megakaryocytes derived from the two routes displayed different activities and were differentially regulated by chemotherapy and inflammation. Our work links differentiation route to the heterogeneity of megakaryocytes. Alternative differentiation routes result in variable combinations of functionally distinct megakaryocyte subpopulations poised for different physiological demands.

The mitochondrial intermembrane space protein mitofissin drives mitochondrial fission required for mitophagy.

Mitophagy plays an important role in mitochondrial homeostasis by selective degradation of mitochondria. During mitophagy, mitochondria should be fragmented to allow engulfment within autophagosomes, whose capacity is exceeded by the typical mitochondria mass. However, the known mitochondrial fission factors, dynamin-related proteins Dnm1 in yeasts and DNM1L/Drp1 in mammals, are dispensable for mitophagy. Here, we identify Atg44 as a mitochondrial fission factor that is essential for mitophagy in yeasts, and we therefore term Atg44 and its orthologous proteins mitofissin. In mitofissin-deficient cells, a part of the mitochondria is recognized by the mitophagy machinery as cargo but cannot be enwrapped by the autophagosome precursor, the phagophore, due to a lack of mitochondrial fission. Furthermore, we show that mitofissin directly binds to lipid membranes and brings about lipid membrane fragility to facilitate membrane fission. Taken together, we propose that mitofissin acts directly on lipid membranes to drive mitochondrial fission required for mitophagy.

Two-dimensional Perovskitoids Enhance Stability in Perovskite Solar Cells.

Two-dimensional (2D)/three-dimensional (3D) perovskite heterostructures have played a key role in advancing the performance of perovskite solar cells (PSCs)1,2. However, the migration of cations between 2D and 3D layers results in the disruption of octahedral networks that leads to degradation in performance over time3,4. We hypothesized that perovskitoids, with robust organic-inorganic networks enabled by edge- and face-sharing, could impede ion migration. We explored a set of perovskitoids of varying dimensionality, and found that cation migration within perovskitoid/perovskite heterostructures was suppressed compared to the 2D/3D perovskite case. Increasing the dimensionality of perovskitoids improves charge transport when they are interfaced with 3D perovskite surfaces - this the result of enhanced octahedral connectivity and out-of-plane orientation. The 2D perovskitoid (A6BfP)8Pb7I22 (A6BfP: N-aminohexyl-benz[f]-phthalimide) provides efficient passivation of perovskite surfaces and enables uniform large-area perovskite films. Devices based on perovskitoid/perovskite heterostructures achieve a certified quasi-steady-state power conversion efficiency of 24.6% for centimeter-area PSCs. We removed the fragile hole transport layers and showed stable operation of the underlying perovskitoid/perovskite heterostructure at 85°C for 1,250 hours for encapsulated large-area devices in an air ambient.

Interferon subverts an AHR-JUN axis to promote CXCL13+ T cells in lupus.

Systemic lupus erythematosus (SLE) is prototypical autoimmune disease driven by pathological T cell-B cell interactions1,2. Expansion of T follicular helper (TFH) and T peripheral helper (TPH) cells, two T cell populations that provide help to B cells, is a prominent feature of SLE3,4. Human TFH and TPH cells characteristically produce high levels of the B cell chemoattractant CXCL13 (refs. 5,6), yet regulation of T cell CXCL13 production and the relationship between CXCL13+ T cells and other T cell states remains unclear. Here, we identify an imbalance in CD4+ T cell phenotypes in patients with SLE, with expansion of PD-1+/ICOS+ CXCL13+ T cells and reduction of CD96hi IL-22+ T cells. Using CRISPR screens, we identify the aryl hydrocarbon receptor (AHR) as a potent negative regulator of CXCL13 production by human CD4+ T cells. Transcriptomic, epigenetic and functional studies demonstrate that AHR coordinates with AP-1 family member JUN to prevent CXCL13+ TPH/TFH cell differentiation and promote an IL-22+ phenotype. Type I interferon, a pathogenic driver of SLE7, opposes AHR and JUN to promote T cell production of CXCL13. These results place CXCL13+ TPH/TFH cells on a polarization axis opposite from T helper 22 (TH22) cells and reveal AHR, JUN and interferon as key regulators of these divergent T cell states.

Spin-vibronic coherence drives singlet-triplet conversion.

Design-specific control over the transitions between excited electronic states with different spin multiplicities is of the utmost importance in molecular and materials chemistry1-3. Previous studies have indicated that the combination of spin-orbit and vibronic effects, collectively termed the spin-vibronic effect, can accelerate quantum-mechanically forbidden transitions at non-adiabatic crossings4,5. However, it has been difficult to identify precise experimental manifestations of the spin-vibronic mechanism. Here we present coherence spectroscopy experiments that reveal the interplay between the spin, electronic and vibrational degrees of freedom that drive efficient singlet-triplet conversion in four structurally related dinuclear Pt(II) metal-metal-to-ligand charge-transfer (MMLCT) complexes. Photoexcitation activates the formation of a Pt-Pt bond, launching a stretching vibrational wavepacket. The molecular-structure-dependent decoherence and recoherence dynamics of this wavepacket resolve the spin-vibronic mechanism. We find that vectorial motion along the Pt-Pt stretching coordinates tunes the singlet and intermediate-state energy gap irreversibly towards the conical intersection and subsequently drives formation of the lowest stable triplet state in a ratcheting fashion. This work demonstrates the viability of using vibronic coherences as probes6-9 to clarify the interplay among spin, electronic and nuclear dynamics in spin-conversion processes, and this could inspire new modular designs to tailor the properties of excited states.

Inference and reconstruction of the heimdallarchaeial ancestry of eukaryotes.

In the ongoing debates about eukaryogenesis-the series of evolutionary events leading to the emergence of the eukaryotic cell from prokaryotic ancestors-members of the Asgard archaea play a key part as the closest archaeal relatives of eukaryotes1. However, the nature and phylogenetic identity of the last common ancestor of Asgard archaea and eukaryotes remain unresolved2-4. Here we analyse distinct phylogenetic marker datasets of an expanded genomic sampling of Asgard archaea and evaluate competing evolutionary scenarios using state-of-the-art phylogenomic approaches. We find that eukaryotes are placed, with high confidence, as a well-nested clade within Asgard archaea and as a sister lineage to Hodarchaeales, a newly proposed order within Heimdallarchaeia. Using sophisticated gene tree and species tree reconciliation approaches, we show that analogous to the evolution of eukaryotic genomes, genome evolution in Asgard archaea involved significantly more gene duplication and fewer gene loss events compared with other archaea. Finally, we infer that the last common ancestor of Asgard archaea was probably a thermophilic chemolithotroph and that the lineage from which eukaryotes evolved adapted to mesophilic conditions and acquired the genetic potential to support a heterotrophic lifestyle. Our work provides key insights into the prokaryote-to-eukaryote transition and a platform for better understanding the emergence of cellular complexity in eukaryotic cells.

An AT-hook domain in MeCP2 determines the clinical course of Rett syndrome and related disorders.

Mutations in the X-linked MECP2 cause Rett syndrome, a devastating neurological disorder typified by a period of apparently normal development followed by loss of cognitive and psychomotor skills. Data from rare male patients suggest symptom onset and severity can be influenced by the location of the mutation, with amino acids 270 and 273 marking the difference between neonatal encephalopathy and death, on the one hand, and survival with deficits on the other. We therefore generated two mouse models expressing either MeCP2-R270X or MeCP2-G273X. The mice developed phenotypes at strikingly different rates and showed differential ATRX nuclear localization within the nervous system, over time, coinciding with phenotypic progression. We discovered that MeCP2 contains three AT-hook-like domains over a stretch of 250 amino acids, like HMGA DNA-bending proteins; one conserved AT-hook is disrupted in MeCP2-R270X, lending further support to the notion that one of MeCP2's key functions is to alter chromatin structure.

A signalling pathway for transcriptional regulation of sleep amount in mice.

In mice and humans, sleep quantity is governed by genetic factors and exhibits age-dependent variation1-3. However, the core molecular pathways and effector mechanisms that regulate sleep duration in mammals remain unclear. Here, we characterize a major signalling pathway for the transcriptional regulation of sleep in mice using adeno-associated virus-mediated somatic genetics analysis4. Chimeric knockout of LKB1 kinase-an activator of AMPK-related protein kinase SIK35-7-in adult mouse brain markedly reduces the amount and delta power-a measure of sleep depth-of non-rapid eye movement sleep (NREMS). Downstream of the LKB1-SIK3 pathway, gain or loss-of-function of the histone deacetylases HDAC4 and HDAC5 in adult brain neurons causes bidirectional changes of NREMS amount and delta power. Moreover, phosphorylation of HDAC4 and HDAC5 is associated with increased sleep need, and HDAC4 specifically regulates NREMS amount in posterior hypothalamus. Genetic and transcriptomic studies reveal that HDAC4 cooperates with CREB in both transcriptional and sleep regulation. These findings introduce the concept of signalling pathways targeting transcription modulators to regulate daily sleep amount and demonstrate the power of somatic genetics in mouse sleep research.

The association of cigarette smoking with DNA methylation and gene expression in human tissue samples.

Cigarette smoking adversely affects many aspects of human health, and epigenetic responses to smoking may reflect mechanisms that mediate or defend against these effects. Prior studies of smoking and DNA methylation (DNAm), typically measured in leukocytes, have identified numerous smoking-associated regions (e.g., AHRR). To identify smoking-associated DNAm features in typically inaccessible tissues, we generated array-based DNAm data for 916 tissue samples from the GTEx (Genotype-Tissue Expression) project representing 9 tissue types (lung, colon, ovary, prostate, blood, breast, testis, kidney, and muscle). We identified 6,350 smoking-associated CpGs in lung tissue (n = 212) and 2,735 in colon tissue (n = 210), most not reported previously. For all 7 other tissue types (sample sizes 38-153), no clear associations were observed (false discovery rate 0.05), but some tissues showed enrichment for smoking-associated CpGs reported previously. For 1,646 loci (in lung) and 22 (in colon), smoking was associated with both DNAm and local gene expression. For loci detected in both lung and colon (e.g., AHRR, CYP1B1, CYP1A1), top CpGs often differed between tissues, but similar clusters of hyper- or hypomethylated CpGs were observed, with hypomethylation at regulatory elements corresponding to increased expression. For lung tissue, 17 hallmark gene sets were enriched for smoking-associated CpGs, including xenobiotic- and cancer-related gene sets. At least four smoking-associated regions in lung were impacted by lung methylation quantitative trait loci (QTLs) that co-localize with genome-wide association study (GWAS) signals for lung function (FEV1/FVC), suggesting epigenetic alterations can mediate the effects of smoking on lung health. Our multi-tissue approach has identified smoking-associated regions in disease-relevant tissues, including effects that are shared across tissue types.

Three-dimensional electronic microfliers inspired by wind-dispersed seeds.

Large, distributed collections of miniaturized, wireless electronic devices1,2 may form the basis of future systems for environmental monitoring3, population surveillance4, disease management5 and other applications that demand coverage over expansive spatial scales. Aerial schemes to distribute the components for such networks are required, and-inspired by wind-dispersed seeds6-we examined passive structures designed for controlled, unpowered flight across natural environments or city settings. Techniques in mechanically guided assembly of three-dimensional (3D) mesostructures7-9 provide access to miniature, 3D fliers optimized for such purposes, in processes that align with the most sophisticated production techniques for electronic, optoelectronic, microfluidic and microelectromechanical technologies. Here we demonstrate a range of 3D macro-, meso- and microscale fliers produced in this manner, including those that incorporate active electronic and colorimetric payloads. Analytical, computational and experimental studies of the aerodynamics of high-performance structures of this type establish a set of fundamental considerations in bio-inspired design, with a focus on 3D fliers that exhibit controlled rotational kinematics and low terminal velocities. An approach that represents these complex 3D structures as discrete numbers of blades captures the essential physics in simple, analytical scaling forms, validated by computational and experimental results. Battery-free, wireless devices and colorimetric sensors for environmental measurements provide simple examples of a wide spectrum of applications of these unusual concepts.

Biodiversity mitigates drought effects in the decomposer system across biomes.

Multiple facets of global change affect the earth system interactively, with complex consequences for ecosystem functioning and stability. Simultaneous climate and biodiversity change are of particular concern, because biodiversity may contribute to ecosystem resistance and resilience and may mitigate climate change impacts. Yet, the extent and generality of how climate and biodiversity change interact remain insufficiently understood, especially for the decomposition of organic matter, a major determinant of the biosphere-atmosphere carbon feedbacks. With an inter-biome field experiment using large rainfall exclusion facilities, we tested how drought, a common prediction of climate change models for many parts of the world, and biodiversity in the decomposer system drive decomposition in forest ecosystems interactively. Decomposing leaf litter lost less carbon (C) and especially nitrogen (N) in five different forest biomes following partial rainfall exclusion compared to conditions without rainfall exclusion. An increasing complexity of the decomposer community alleviated drought effects, with full compensation when large-bodied invertebrates were present. Leaf litter mixing increased diversity effects, with increasing litter species richness, which contributed to counteracting drought effects on C and N loss, although to a much smaller degree than decomposer community complexity. Our results show at a relevant spatial scale covering distinct climate zones that both, the diversity of decomposer communities and plant litter in forest floors have a strong potential to mitigate drought effects on C and N dynamics during decomposition. Preserving biodiversity at multiple trophic levels contributes to ecosystem resistance and appears critical to maintain ecosystem processes under ongoing climate change.

TFAM loss induces nuclear actin assembly upon mDia2 malonylation to promote liver cancer metastasis.

The mechanisms underlying cancer metastasis remain poorly understood. Here, we report that TFAM deficiency rapidly and stably induced spontaneous lung metastasis in mice with liver cancer. Interestingly, unexpected polymerization of nuclear actin was observed in TFAM-knockdown HCC cells when cytoskeleton was examined. Polymerization of nuclear actin is causally linked to the high-metastatic ability of HCC cells by modulating chromatin accessibility and coordinating the expression of genes associated with extracellular matrix remodeling, angiogenesis, and cell migration. Mechanistically, TFAM deficiency blocked the TCA cycle and increased the intracellular malonyl-CoA levels. Malonylation of mDia2, which drives actin assembly, promotes its nuclear translocation. Importantly, inhibition of malonyl-CoA production or nuclear actin polymerization significantly impeded the spread of HCC cells in mice. Moreover, TFAM was significantly downregulated in metastatic HCC tissues and was associated with overall survival and time to tumor recurrence of HCC patients. Taken together, our study connects mitochondria to the metastasis of human cancer via uncovered mitochondria-to-nucleus retrograde signaling, indicating that TFAM may serve as an effective target to block HCC metastasis.

Deciphering Photoinduced Catalytic Reaction Mechanisms in Natural and Artificial Photosynthetic Systems on Multiple Temporal and Spatial Scales Using X-ray Probes.

Utilization of renewable energies for catalytically generating value-added chemicals is highly desirable in this era of rising energy demands and climate change impacts. Artificial photosynthetic systems or photocatalysts utilize light to convert abundant CO2, H2O, and O2 to fuels, such as carbohydrates and hydrogen, thus converting light energy to storable chemical resources. The emergence of intense X-ray pulses from synchrotrons, ultrafast X-ray pulses from X-ray free electron lasers, and table-top laser-driven sources over the past decades opens new frontiers in deciphering photoinduced catalytic reaction mechanisms on the multiple temporal and spatial scales. Operando X-ray spectroscopic methods offer a new set of electronic transitions in probing the oxidation states, coordinating geometry, and spin states of the metal catalytic center and photosensitizers with unprecedented energy and time resolution. Operando X-ray scattering methods enable previously elusive reaction steps to be characterized on different length scales and time scales. The methodological progress and their application examples collected in this review will offer a glimpse into the accomplishments and current state in deciphering reaction mechanisms for both natural and synthetic systems. Looking forward, there are still many challenges and opportunities at the frontier of catalytic research that will require further advancement of the characterization techniques.

Clades of huge phages from across Earth's ecosystems.

Bacteriophages typically have small genomes1 and depend on their bacterial hosts for replication2. Here we sequenced DNA from diverse ecosystems and found hundreds of phage genomes with lengths of more than 200 kilobases (kb), including a genome of 735 kb, which is-to our knowledge-the largest phage genome to be described to date. Thirty-five genomes were manually curated to completion (circular and no gaps). Expanded genetic repertoires include diverse and previously undescribed CRISPR-Cas systems, transfer RNAs (tRNAs), tRNA synthetases, tRNA-modification enzymes, translation-initiation and elongation factors, and ribosomal proteins. The CRISPR-Cas systems of phages have the capacity to silence host transcription factors and translational genes, potentially as part of a larger interaction network that intercepts translation to redirect biosynthesis to phage-encoded functions. In addition, some phages may repurpose bacterial CRISPR-Cas systems to eliminate competing phages. We phylogenetically define the major clades of huge phages from human and other animal microbiomes, as well as from oceans, lakes, sediments, soils and the built environment. We conclude that the large gene inventories of huge phages reflect a conserved biological strategy, and that the phages are distributed across a broad bacterial host range and across Earth's ecosystems.

Structural insights into the HDAC4-MEF2A-DNA complex and its implication in long-range transcriptional regulation.

Class IIa Histone deacetylases (HDACs), including HDAC4, 5, 7 and 9, play key roles in multiple important developmental and differentiation processes. Recent studies have shown that class IIa HDACs exert their transcriptional repressive function by interacting with tissue-specific transcription factors, such as members of the myocyte enhancer factor 2 (MEF2) family of transcription factors. However, the molecular mechanism is not well understood. In this study, we determined the crystal structure of an HDAC4-MEF2A-DNA complex. This complex adopts a dumbbell-shaped overall architecture, with a 2:4:2 stoichiometry of HDAC4, MEF2A and DNA molecules. In the complex, two HDAC4 molecules form a dimer through the interaction of their glutamine-rich domain (GRD) to form the stem of the 'dumbbell'; while two MEF2A dimers and their cognate DNA molecules are bridged by the HDAC4 dimer. Our structural observations were then validated using biochemical and mutagenesis assays. Further cell-based luciferase reporter gene assays revealed that the dimerization of HDAC4 is crucial in its ability to repress the transcriptional activities of MEF2 proteins. Taken together, our findings not only provide the structural basis for the assembly of the HDAC4-MEF2A-DNA complex but also shed light on the molecular mechanism of HDAC4-mediated long-range gene regulation.

WSB1/2 target chromatin-bound lysine-methylated RelA for proteasomal degradation and NF-κB termination.

Proteasome-mediated degradation of chromatin-bound NF-κB is critical in terminating the transcription of pro-inflammatory genes and can be triggered by Set9-mediated lysine methylation of the RelA subunit. However, the E3 ligase targeting methylated RelA remains unknown. Here, we find that two structurally similar substrate-recognizing components of Cullin-RING E3 ligases, WSB1 and WSB2, can recognize chromatin-bound methylated RelA for polyubiquitination and proteasomal degradation. We showed that WSB1/2 negatively regulated a subset of NF-κB target genes via associating with chromatin where they targeted methylated RelA for ubiquitination, facilitating the termination of NF-κB-dependent transcription. WSB1/2 specifically interacted with methylated lysines (K) 314 and 315 of RelA via their N-terminal WD-40 repeat (WDR) domains, thereby promoting ubiquitination of RelA. Computational modeling further revealed that a conserved aspartic acid (D) at position 158 within the WDR domain of WSB2 coordinates K314/K315 of RelA, with a higher affinity when either of the lysines is methylated. Mutation of D158 abolished WSB2's ability to bind to and promote ubiquitination of methylated RelA. Together, our study identifies a novel function and the underlying mechanism for WSB1/2 in degrading chromatin-bound methylated RelA and preventing sustained NF-κB activation, providing potential new targets for therapeutic intervention of NF-κB-mediated inflammatory diseases.

Sequencing-based fine-mapping and in silico functional characterization of the 10q24.32 arsenic metabolism efficiency locus across multiple arsenic-exposed populations.

Inorganic arsenic is highly toxic and carcinogenic to humans. Exposed individuals vary in their ability to metabolize arsenic, and variability in arsenic metabolism efficiency (AME) is associated with risks of arsenic-related toxicities. Inherited genetic variation in the 10q24.32 region, near the arsenic methyltransferase (AS3MT) gene, is associated with urine-based measures of AME in multiple arsenic-exposed populations. To identify potential causal variants in this region, we applied fine mapping approaches to targeted sequencing data generated for exposed individuals from Bangladeshi, American Indian, and European American populations (n = 2,357, 557, and 648 respectively). We identified three independent association signals for Bangladeshis, two for American Indians, and one for European Americans. The size of the confidence sets for each signal varied from 4 to 85 variants. There was one signal shared across all three populations, represented by the same SNP in American Indians and European Americans (rs191177668) and in strong linkage disequilibrium (LD) with a lead SNP in Bangladesh (rs145537350). Beyond this shared signal, differences in LD patterns, minor allele frequency (MAF) (e.g., rs12573221 ~13% in Bangladesh ~0.2% among American Indians), and/or heterogeneity in effect sizes across populations likely contributed to the apparent population specificity of the additional identified signals. One of our potential causal variants influences AS3MT expression and nearby DNA methylation in numerous GTEx tissue types (with rs4919690 as a likely causal variant). Several SNPs in our confidence sets overlap transcription factor binding sites and cis-regulatory elements (from ENCODE). Taken together, our analyses reveal multiple potential causal variants in the 10q24.32 region influencing AME, including a variant shared across populations, and elucidate potential biological mechanisms underlying the impact of genetic variation on AME.

Integrative multiomics profiling of passion fruit reveals the genetic basis for fruit color and aroma.

Passion fruit (Passiflora edulis) possesses a complex aroma and is widely grown in tropical and subtropical areas. Here, we conducted the de novo assembly, annotation, and comparison of PPF (P. edulis Sims) and YPF (P. edulis f. flavicarpa) reference genomes using PacBio, Illumina, and Hi-C technologies. Notably, we discovered evidence of recent whole-genome duplication events in P. edulis genomes. Comparative analysis revealed 7.6∼8.1 million single nucleotide polymorphisms, 1 million insertions/deletions, and over 142 Mb presence/absence variations among different P. edulis genomes. During the ripening of yellow passion fruit, metabolites related to flavor, aroma, and color were substantially accumulated or changed. Through joint analysis of genomic variations, differentially expressed genes, and accumulated metabolites, we explored candidate genes associated with flavor, aroma, and color distinctions. Flavonoid biosynthesis pathways, anthocyanin biosynthesis pathways, and related metabolites are pivotal factors affecting the coloration of passion fruit, and terpenoid metabolites accumulated more in PPF. Finally, by heterologous expression in yeast (Saccharomyces cerevisiae), we functionally characterized 12 terpene synthases. Our findings revealed that certain TPS homologs in both YPF and PPF varieties produce identical terpene products, while others yield distinct compounds or even lose their functionality. These discoveries revealed the genetic and metabolic basis of unique characteristics in aroma and flavor between the 2 passion fruit varieties. This study provides resources for better understanding the genome architecture and accelerating genetic improvement of passion fruits.