ABSTRACT
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.
Subject(s)
Genome , Snakes , Animals , Snakes/genetics , Adaptation, Physiological , Acclimatization , Evolution, Molecular , Phylogeny , Biological EvolutionABSTRACT
N7-methylguanosine (m7G) modification, routinely occurring at mRNA 5' cap or within tRNAs/rRNAs, also exists internally in messenger RNAs (mRNAs). Although m7G-cap is essential for pre-mRNA processing and protein synthesis, the exact role of mRNA internal m7G modification remains elusive. Here, we report that mRNA internal m7G is selectively recognized by Quaking proteins (QKIs). By transcriptome-wide profiling/mapping of internal m7G methylome and QKI-binding sites, we identified more than 1,000 high-confidence m7G-modified and QKI-bound mRNA targets with a conserved "GANGAN (N = A/C/U/G)" motif. Strikingly, QKI7 interacts (via C terminus) with the stress granule (SG) core protein G3BP1 and shuttles internal m7G-modified transcripts into SGs to regulate mRNA stability and translation under stress conditions. Specifically, QKI7 attenuates the translation efficiency of essential genes in Hippo signaling pathways to sensitize cancer cells to chemotherapy. Collectively, we characterized QKIs as mRNA internal m7G-binding proteins that modulate target mRNA metabolism and cellular drug resistance.
Subject(s)
DNA Helicases , RNA Helicases , DNA Helicases/metabolism , RNA Recognition Motif Proteins/genetics , RNA Recognition Motif Proteins/metabolism , RNA Helicases/metabolism , Stress Granules , Poly-ADP-Ribose Binding Proteins/genetics , Poly-ADP-Ribose Binding Proteins/metabolism , GTP-Binding Proteins/metabolism , RNA, Messenger/metabolism , Cytoplasmic Granules/metabolismABSTRACT
Interleukin-33 (IL-33), an epithelial cell-derived cytokine that responds rapidly to environmental insult, has a critical role in initiating airway inflammatory diseases. However, the molecular mechanism underlying IL-33 secretion following allergen exposure is not clear. Here, we found that two cell events were fundamental for IL-33 secretion after exposure to allergens. First, stress granule assembly activated by allergens licensed the nuclear-cytoplasmic transport of IL-33, but not the secretion of IL-33. Second, a neo-form murine amino-terminal p40 fragment gasdermin D (Gsdmd), whose generation was independent of inflammatory caspase-1 and caspase-11, dominated cytosolic secretion of IL-33 by forming pores in the cell membrane. Either the blockade of stress granule assembly or the abolishment of p40 production through amino acid mutation of residues 309-313 (ELRQQ) could efficiently prevent the release of IL-33 in murine epithelial cells. Our findings indicated that targeting stress granule disassembly and Gsdmd fragmentation could reduce IL-33-dependent allergic airway inflammation.
Subject(s)
Allergens , Interleukin-33 , Phosphate-Binding Proteins/metabolism , Pore Forming Cytotoxic Proteins/metabolism , Animals , Caspase 1/metabolism , Inflammation , Interleukin-1beta/metabolism , Interleukin-33/genetics , Interleukin-33/metabolism , Intracellular Signaling Peptides and Proteins/genetics , Intracellular Signaling Peptides and Proteins/metabolism , Mice , Peptide Hydrolases/metabolism , Stress GranulesABSTRACT
Photons with spin angular momentum possess intrinsic chirality, which underpins many phenomena including nonlinear optics1, quantum optics2, topological photonics3 and chiroptics4. Intrinsic chirality is weak in natural materials, and recent theoretical proposals5-7 aimed to enlarge circular dichroism by resonant metasurfaces supporting bound states in the continuum that enhance substantially chiral light-matter interactions. Those insightful works resort to three-dimensional sophisticated geometries, which are too challenging to be realized for optical frequencies8. Therefore, most of the experimental attempts9-11 showing strong circular dichroism rely on false/extrinsic chirality by using either oblique incidence9,10 or structural anisotropy11. Here we report on the experimental realization of true/intrinsic chiral response with resonant metasurfaces in which the engineered slant geometry breaks both in-plane and out-of-plane symmetries. Our result marks, to our knowledge, the first observation of intrinsic chiral bound states in the continuum with near-unity circular dichroism of 0.93 and a high quality factor exceeding 2,663 for visible frequencies. Our chiral metasurfaces may lead to a plethora of applications in chiral light sources and detectors, chiral sensing, valleytronics and asymmetric photocatalysis.
ABSTRACT
Although previous studies have identified human-specific accelerated regions as playing a key role in the recent evolution of the human brain, the characteristics and cellular functions of rapidly evolving conserved elements (RECEs) in ancestral primate lineages remain largely unexplored. Here, based on large-scale primate genome assemblies, we identify 888 RECEs that have been highly conserved in primates that exhibit significantly accelerated substitution rates in the ancestor of the Simiiformes. This primate lineage exhibits remarkable morphological innovations, including an expanded brain mass. Integrative multiomic analyses reveal that RECEs harbor sequences with potential cis-regulatory functions that are activated in the adult human brain. Importantly, genes linked to RECEs exhibit pronounced expression trajectories in the adult brain relative to the fetal stage. Furthermore, we observed an increase in the chromatin accessibility of RECEs in oligodendrocytes from individuals with Alzheimer's disease (AD) compared to that of a control group, indicating that these RECEs may contribute to brain aging and AD. Our findings serve to expand our knowledge of the genetic underpinnings of brain function during primate evolution.
Subject(s)
Alzheimer Disease , Animals , Humans , Alzheimer Disease/genetics , Evolution, Molecular , Primates/genetics , BrainABSTRACT
In recent years, the impact of age-related diseases on human health has become increasingly severe, and developing effective drugs to deal with these diseases has become an urgent task. Considering the essential regulatory role of hydrogen sulfide (H2S) in these diseases, it is regarded as a promising target for treatment. H2S is a novel gaseous transmitter involved in many critical physiological activities, including anti-oxidation, anti-inflammation, and angiogenesis. H2S also regulates cell activities such as cell proliferation, migration, invasion, apoptosis, and autophagy. These regulatory effects of H2S contribute to relieving and treating age-related diseases. In this review, we mainly focus on the pathogenesis and treatment prospects of H2S in regulating age-related diseases.
Subject(s)
Aging , Hydrogen Sulfide , Hydrogen Sulfide/metabolism , Hydrogen Sulfide/pharmacology , Humans , Aging/metabolism , Animals , Autophagy/drug effects , Apoptosis/drug effects , Cell Proliferation/drug effectsABSTRACT
Airborne bacteria are an influential component of the Earth's microbiomes, but their community structure and biogeographic distribution patterns have yet to be understood. We analyzed the bacterial communities of 370 air particulate samples collected from 63 sites around the world and constructed an airborne bacterial reference catalog with more than 27 million nonredundant 16S ribosomal RNA (rRNA) gene sequences. We present their biogeographic pattern and decipher the interlacing of the microbiome co-occurrence network with surface environments of the Earth. While the total abundance of global airborne bacteria in the troposphere (1.72 × 1024 cells) is 1 to 3 orders of magnitude lower than that of other habitats, the number of bacterial taxa (i.e., richness) in the atmosphere (4.71 × 108 to 3.08 × 109) is comparable to that in the hydrosphere, and its maximum occurs in midlatitude regions, as is also observed in other ecosystems. The airborne bacterial community harbors a unique set of dominant taxa (24 species); however, its structure appears to be more easily perturbed, due to the more prominent role of stochastic processes in shaping community assembly. This is corroborated by the major contribution of surface microbiomes to airborne bacteria (averaging 46.3%), while atmospheric conditions such as meteorological factors and air quality also play a role. Particularly in urban areas, human impacts weaken the relative importance of plant sources of airborne bacteria and elevate the occurrence of potential pathogens from anthropogenic sources. These findings serve as a key reference for predicting planetary microbiome responses and the health impacts of inhalable microbiomes with future changes in the environment.
Subject(s)
Air Microbiology , Microbiota , Anthropogenic Effects , Bacteria/genetics , Humans , Microbiota/genetics , RNA, Ribosomal, 16S/geneticsABSTRACT
Lorises are a group of globally threatened strepsirrhine primates that exhibit many unusual physiological and behavioral features, including a low metabolic rate, slow movement, and hibernation. Here, we assembled a chromosome-level genome sequence of the pygmy loris (Xanthonycticebus pygmaeus) and resequenced whole genomes from 50 pygmy lorises and 6 Bengal slow lorises (Nycticebus bengalensis). We found that many gene families involved in detoxification have been specifically expanded in the pygmy loris, including the GSTA gene family, with many newly derived copies functioning specifically in the liver. We detected many genes displaying evolutionary convergence between pygmy loris and koala, including PITRM1. Significant decreases in PITRM1 enzymatic activity in these two species may have contributed to their characteristic low rate of metabolism. We also detected many evolutionarily convergent genes and positively selected genes in the pygmy loris that are involved in muscle development. Functional assays demonstrated the decreased ability of one positively selected gene, MYOF, to up-regulate the fast-type muscle fiber, consistent with the lower proportion of fast-twitch muscle fibers in the pygmy loris. The protein product of another positively selected gene in the pygmy loris, PER2, exhibited weaker binding to the key circadian core protein CRY, a finding that may be related to this species' unusual circadian rhythm. Finally, population genomics analysis revealed that these two extant loris species, which coexist in the same habitat, have exhibited an inverse relationship in terms of their demography over the past 1 million years, implying strong interspecies competition after speciation.
Subject(s)
Adaptation, Biological , Biological Evolution , Lorisidae , Adaptation, Biological/genetics , Animals , Demography , Hibernation , Lorisidae/genetics , Metagenomics , Metalloendopeptidases/geneticsABSTRACT
Chiral nanostructures allow engineering of chiroptical responses; however, their design usually relies on empirical approaches and extensive numerical simulations. It remains unclear if a general strategy exists to enhance and maximize the intrinsic chirality of subwavelength photonic structures. Here, we suggest a microscopic theory and uncover the origin of strong chiral responses of resonant nanostructures. We reveal that the reactive helicity density is critically important for achieving maximum chirality at resonances. We demonstrate our general concept on the examples of planar photonic crystal slabs and metasurfaces, where out-of-plane mirror symmetry is broken by a bilayer design. Our findings provide a general recipe for designing photonic structures with maximum chirality, paving the way toward many applications, including chiral sensing, chiral emitters and detectors, and chiral quantum optics.
ABSTRACT
Controllable droplet manipulation has diverse applications; however, limited methods exist for externally manipulating droplets in confined spaces. Herein, we propose a portable triboelectric electrostatic tweezer (TET) by integrating electrostatic forces with a superhydrophobic surface that can even manipulate droplets in an enclosed space. Electrostatic induction causes the droplet to be subjected to an electrostatic force in an electrostatic field so that the droplet can be moved freely with the TET on a superhydrophobic platform. Characterized by its high precision, flexibility, and robust binding strength, TET can manipulate droplets under various conditions and achieve a wide range of representative fluid applications such as droplet microreactors, precise self-cleaning, cargo transportation, the targeted delivery of chemicals, liquid sorting, soft droplet robotics, and cell labeling. Specifically, TET demonstrated the ability to manipulate internal droplets from the outside of a closed system, such as performing cell labeling experiments within a sealed Petri dish without opening the culture system.
ABSTRACT
The emerging two-photon polymerization (TPP) technique enables high-resolution printing of complex 3D structures, revolutionizing micro/nano additive manufacturing. Various fast scanning and parallel processing strategies have been proposed to promote its efficiency. However, obtaining large numbers of uniform focal spots for parallel high-speed scanning remains challenging, which hampers the realization of higher throughput. We report a TPP printing platform that combines galvanometric mirrors and liquid crystal on silicon spatial light modulator (LCoS-SLM). By setting the target light field at LCoS-SLM's diffraction center, sufficient energy is acquired to support simultaneous polymerization of over 400 foci. With fast scanning, the maximum printing speed achieves 1.49 × 108 voxels s-1, surpassing the existing scanning-based TPP methods while maintaining high printing resolution and flexibility. To demonstrate the processing capability, functional 3D microstructure arrays are rapidly fabricated and applied in micro-optics and micro-object manipulation. Our method may expand the prospects of TPP in large-scale micro/nanomanufacturing.
ABSTRACT
Inspired by the reverse thrust generated by fuel injection, micromachines that are self-propelled by bubble ejection are developed, such as microrods, microtubes, and microspheres. However, controlling bubble ejection sites to build micromachines with programmable actuation and further enabling mechanical transmission remain challenging. Here, bubble-propelled mechanical microsystems are constructed by proposing a multimaterial femtosecond laser processing method, consisting of direct laser writing and selective laser metal reduction. The polymer frame of the microsystems is first printed, followed by the deposition of catalytic platinum into the desired local site of the microsystems by laser reduction. With this method, a variety of designable microrotors with selective bubble ejection sites are realized, which enable excellent mechanical transmission systems composed of single and multiple mechanical components, including a coupler, a crank slider, and a crank rocker system. We believe the presented bubble-propelled mechanical microsystems could be extended to applications in microrobotics, microfluidics, and microsensors.
ABSTRACT
BACKGROUND: Although many studies have been done to reveal artificial selection signatures in commercial and indigenous chickens, a limited number of genes have been linked to specific traits. To identify more trait-related artificial selection signatures and genes, we re-sequenced a total of 85 individuals of five indigenous chicken breeds with distinct traits from Yunnan Province, China. RESULTS: We found 30 million non-redundant single nucleotide variants and small indels (< 50 bp) in the indigenous chickens, of which 10 million were not seen in 60 broilers, 56 layers and 35 red jungle fowls (RJFs) that we compared with. The variants in each breed are enriched in non-coding regions, while those in coding regions are largely tolerant, suggesting that most variants might affect cis-regulatory sequences. Based on 27 million bi-allelic single nucleotide polymorphisms identified in the chickens, we found numerous selective sweeps and affected genes in each indigenous chicken breed and substantially larger numbers of selective sweeps and affected genes in the broilers and layers than previously reported using a rigorous statistical model. Consistent with the locations of the variants, the vast majority (~ 98.3%) of the identified selective sweeps overlap known quantitative trait loci (QTLs). Meanwhile, 74.2% known QTLs overlap our identified selective sweeps. We confirmed most of previously identified trait-related genes and identified many novel ones, some of which might be related to body size and high egg production traits. Using RT-qPCR, we validated differential expression of eight genes (GHR, GHRHR, IGF2BP1, OVALX, ELF2, MGARP, NOCT, SLC25A15) that might be related to body size and high egg production traits in relevant tissues of relevant breeds. CONCLUSION: We identify 30 million single nucleotide variants and small indels in the five indigenous chicken breeds, 10 million of which are novel. We predict substantially more selective sweeps and affected genes than previously reported in both indigenous and commercial breeds. These variants and affected genes are good candidates for further experimental investigations of genotype-phenotype relationships and practical applications in chicken breeding programs.
Subject(s)
Chickens , Polymorphism, Single Nucleotide , Quantitative Trait Loci , Selection, Genetic , Animals , Chickens/genetics , Genome , INDEL Mutation , Breeding , Phenotype , Genomics/methodsABSTRACT
BACKGROUND: Although multiple chicken genomes have been assembled and annotated, the numbers of protein-coding genes in chicken genomes and their variation among breeds are still uncertain due to the low quality of these genome assemblies and limited resources used in their gene annotations. To fill these gaps, we recently assembled genomes of four indigenous chicken breeds with distinct traits at chromosome-level. In this study, we annotated genes in each of these assembled genomes using a combination of RNA-seq- and homology-based approaches. RESULTS: We identified varying numbers (17,497-17,718) of protein-coding genes in the four indigenous chicken genomes, while recovering 51 of the 274 "missing" genes in birds in general, and 36 of the 174 "missing" genes in chickens in particular. Intriguingly, based on deeply sequenced RNA-seq data collected in multiple tissues in the four breeds, we found 571 ~ 627 protein-coding genes in each genome, which were missing in the annotations of the reference chicken genomes (GRCg6a and GRCg7b/w). After removing redundancy, we ended up with a total of 1,420 newly annotated genes (NAGs). The NAGs tend to be found in subtelomeric regions of macro-chromosomes (chr1 to chr5, plus chrZ) and middle chromosomes (chr6 to chr13, plus chrW), as well as in micro-chromosomes (chr14 to chr39) and unplaced contigs, where G/C contents are high. Moreover, the NAGs have elevated quadruplexes G frequencies, while both G/C contents and quadruplexes G frequencies in their surrounding regions are also high. The NAGs showed tissue-specific expression, and we were able to verify 39 (92.9%) of 42 randomly selected ones in various tissues of the four chicken breeds using RT-qPCR experiments. Most of the NAGs were also encoded in the reference chicken genomes, thus, these genomes might harbor more genes than previously thought. CONCLUSION: The NAGs are widely distributed in wild, indigenous and commercial chickens, and they might play critical roles in chicken physiology. Counting these new genes, chicken genomes harbor more genes than originally thought.
Subject(s)
Chickens , Genome , Molecular Sequence Annotation , Animals , Chickens/genetics , Base Composition , Telomere/genetics , Chromosomes/genetics , Genomics/methodsABSTRACT
Determining the functional consequences of karyotypic changes is invariably challenging because evolution tends to obscure many of its own footprints, such as accumulated mutations, recombination events, and demographic perturbations. Here, we describe the assembly of a chromosome-level reference genome of the gayal (Bos frontalis) thereby revealing the structure, at base-pair-level resolution, of a telo/acrocentric-to-telo/acrocentric Robertsonian translocation (2;28) (T/A-to-T/A rob[2;28]). The absence of any reduction in the recombination rate or genetic introgression within the fusion region of gayal served to challenge the long-standing view of a role for fusion-induced meiotic dysfunction in speciation. The disproportionate increase noted in the distant interactions across pro-chr2 and pro-chr28, and the change in open-chromatin accessibility following rob(2;28), may, however, have led to the various gene expression irregularities observed in the gayal. Indeed, we found that many muscle-related genes, located synthetically on pro-chr2 and pro-chr28, exhibited significant changes in expression. This, combined with genome-scale structural variants and expression alterations in genes involved in myofibril composition, may have driven the rapid sarcomere adaptation of gayal to its rugged mountain habitat. Our findings not only suggest that large-scale chromosomal changes can lead to alterations in genome-level expression, thereby promoting both adaptation and speciation, but also illuminate novel avenues for studying the relationship between karyotype evolution and speciation.
Subject(s)
Chromatin , Genome , Animals , CattleABSTRACT
Although the continual expansion of the brain during primate evolution accounts for our enhanced cognitive capabilities, the drivers of brain evolution have scarcely been explored in these ancestral nodes. Here, we performed large-scale comparative genomic, transcriptomic, and epigenomic analyses to investigate the evolutionary alterations acquired by brain genes and provide comprehensive listings of innovatory genetic elements along the evolutionary path from ancestral primates to human. The regulatory sequences associated with brain-expressed genes experienced rapid change, particularly in the ancestor of the Simiiformes. Extensive comparisons of single-cell and bulk transcriptomic data between primate and nonprimate brains revealed that these regulatory sequences may drive the high expression of certain genes in primate brains. Employing in utero electroporation into mouse embryonic cortex, we show that the primate-specific brain-biased gene BMP7 was recruited, probably in the ancestor of the Simiiformes, to regulate neuronal proliferation in the primate ventricular zone. Our study provides a comprehensive listing of genes and regulatory changes along the brain evolution lineage of ancestral primates leading to human. These data should be invaluable for future functional studies that will deepen our understanding not only of the genetic basis of human brain evolution but also of inherited disease.
Subject(s)
Brain , Primates , Mice , Humans , Animals , Primates/genetics , Brain/metabolism , Evolution, MolecularABSTRACT
HIV-1 is a highly host-specific retrovirus that infects humans but not most nonhuman primates. Thus, the lack of a suitable primate model that can be directly infected with HIV-1 hinders HIV-1/AIDS research. In the previous study, we have found that the northern pig-tailed macaques (NPMs) are susceptible to HIV-1 infection but show a nonpathogenic state. In this study, to understand this macaque-HIV-1 interaction, we assembled a de novo genome and longitudinal transcriptome for this species during the course of HIV-1 infection. Using comparative genomic analysis, a positively selected gene, Toll-like receptor 8, was identified with a weak ability to induce an inflammatory response in this macaque. In addition, an interferon-stimulated gene, interferon alpha inducible protein 27, was upregulated in acute HIV-1 infection and acquired an enhanced ability to inhibit HIV-1 replication compared with its human ortholog. These findings coincide with the observation of persistently downregulated immune activation and low viral replication and can partially explain the AIDS-free state in this macaque following HIV-1 infection. This study identified a number of unexplored host genes that may hamper HIV-1 replication and pathogenicity in NPMs and provided new insights into the host defense mechanisms in cross-species infection of HIV-1. This work will facilitate the adoption of NPM as a feasible animal model for HIV-1/AIDS research.
Subject(s)
HIV Infections , HIV-1 , Simian Immunodeficiency Virus , Animals , Humans , Macaca nemestrina , HIV-1/genetics , Genomics , Simian Immunodeficiency Virus/geneticsABSTRACT
Immune checkpoint inhibitors (ICIs) have shown promising efficacy in multiple cancers including biliary tract cancers (BTCs). However, the data focusing on the efficacy of ICIs in patients with gallbladder cancer (GBC) is still limited. In this study, we aim to assess the efficacy of ICIs in GBC and explore the clinicopathologic and molecular markers associated with ICI benefit. We retrospective analyzed 69 GBC patients who had received ICI therapy between January 2016 and December 2020. Tumor samples were obtained for genomic sequencing and immunohistochemical analysis. The median progression-free survival (PFS) and overall survival (OS) was 4.4 months and 8.5 months, respectively. Multivariate analysis indicated that alcohol intake history, carcinoma embryonic antigen (CEA) level ≥100 U/mL, and cutaneous immune-related adverse events (irAEs) were independent prognostic factors for PFS. CEA level ≥100 U/mL and cutaneous irAEs were independent prognostic factors for OS. The objective response rate and disease control rate (DCR) were 15.9% and 37.7%, respectively. Patients with cutaneous irAEs, high CD8+ T cell infiltrated or immune inflamed GBCs had higher DCR. Patients with high CD8+ T cell infiltrated or immune inflamed GBCs also had a notably improved prognosis. These results suggest that ICIs were effective in patients with GBC. High CEA level, cutaneous irAEs, high CD8+ T cell infiltration, and immune inflamed phenotype could be useful for predicting the efficacy of ICIs in GBC.
Subject(s)
Gallbladder Neoplasms , Immune Checkpoint Inhibitors , Humans , Gallbladder Neoplasms/drug therapy , Gallbladder Neoplasms/pathology , Gallbladder Neoplasms/immunology , Immune Checkpoint Inhibitors/therapeutic use , Immune Checkpoint Inhibitors/adverse effects , Male , Female , Middle Aged , Aged , Retrospective Studies , Prognosis , Adult , Aged, 80 and over , CD8-Positive T-Lymphocytes/immunology , CD8-Positive T-Lymphocytes/drug effects , Progression-Free Survival , Biomarkers, Tumor , Treatment OutcomeABSTRACT
Coexisting orders are key features of strongly correlated materials and underlie many intriguing phenomena from unconventional superconductivity to topological orders. Here, we report the coexistence of two interacting charge-density-wave (CDW) orders in EuTe_{4}, a layered crystal that has drawn considerable attention owing to its anomalous thermal hysteresis and a semiconducting CDW state despite the absence of perfect Fermi surface nesting. By accessing unoccupied conduction bands with time- and angle-resolved photoemission measurements, we find that monolayers and bilayers of Te in the unit cell host different CDWs that are associated with distinct energy gaps. The two gaps display dichotomous evolutions following photoexcitation, where the larger bilayer CDW gap exhibits less renormalization and faster recovery. Surprisingly, the CDW in the Te monolayer displays an additional momentum-dependent gap renormalization that cannot be captured by density-functional theory calculations. This phenomenon is attributed to interlayer interactions between the two CDW orders, which account for the semiconducting nature of the equilibrium state. Our findings not only offer microscopic insights into the correlated ground state of EuTe_{4} but also provide a general nonequilibrium approach to understand coexisting, layer-dependent orders in a complex system.