Guidelines and definitions for research on epithelial-mesenchymal transition.

Epithelial-mesenchymal transition (EMT) encompasses dynamic changes in cellular organization from epithelial to mesenchymal phenotypes, which leads to functional changes in cell migration and invasion. EMT occurs in a diverse range of physiological and pathological conditions and is driven by a conserved set of inducing signals, transcriptional regulators and downstream effectors. With over 5,700 publications indexed by Web of Science in 2019 alone, research on EMT is expanding rapidly. This growing interest warrants the need for a consensus among researchers when referring to and undertaking research on EMT. This Consensus Statement, mediated by 'the EMT International Association' (TEMTIA), is the outcome of a 2-year-long discussion among EMT researchers and aims to both clarify the nomenclature and provide definitions and guidelines for EMT research in future publications. We trust that these guidelines will help to reduce misunderstanding and misinterpretation of research data generated in various experimental models and to promote cross-disciplinary collaboration to identify and address key open questions in this research field. While recognizing the importance of maintaining diversity in experimental approaches and conceptual frameworks, we emphasize that lasting contributions of EMT research to increasing our understanding of developmental processes and combatting cancer and other diseases depend on the adoption of a unified terminology to describe EMT.

Axon Degeneration Gated by Retrograde Activation of Somatic Pro-apoptotic Signaling.

During development, sensory axons compete for limiting neurotrophic support, and local neurotrophin insufficiency triggers caspase-dependent axon degeneration. The signaling driving axon degeneration upon local deprivation is proposed to reside within axons. Our results instead support a model in which, despite the apoptotic machinery being present in axons, the cell body is an active participant in gating axonal caspase activation and axon degeneration. Loss of trophic support in axons initiates retrograde activation of a somatic pro-apoptotic pathway, which, in turn, is required for distal axon degeneration via an anterograde pro-degenerative factor. At a molecular level, the cell body is the convergence point of two signaling pathways whose integrated action drives upregulation of pro-apoptotic Puma, which, unexpectedly, is confined to the cell body. Puma then overcomes inhibition by pro-survival Bcl-xL and Bcl-w and initiates the anterograde pro-degenerative program, highlighting the role of the cell body as an arbiter of large-scale axon removal.

Pathological axonal death through a MAPK cascade that triggers a local energy deficit.

Axonal death disrupts functional connectivity of neural circuits and is a critical feature of many neurodegenerative disorders. Pathological axon degeneration often occurs independently of known programmed death pathways, but the underlying molecular mechanisms remain largely unknown. Using traumatic injury as a model, we systematically investigate mitogen-activated protein kinase (MAPK) families and delineate a MAPK cascade that represents the early degenerative response to axonal injury. The adaptor protein Sarm1 is required for activation of this MAPK cascade, and this Sarm1-MAPK pathway disrupts axonal energy homeostasis, leading to ATP depletion before physical breakdown of damaged axons. The protective cytoNmnat1/Wld(s) protein inhibits activation of this MAPK cascade. Further, MKK4, a key component in the Sarm1-MAPK pathway, is antagonized by AKT signaling, which modulates the degenerative response by limiting activation of downstream JNK signaling. Our results reveal a regulatory mechanism that integrates distinct signals to instruct pathological axon degeneration.

An anti-CRISPR that pulls apart a CRISPR-Cas complex.

In biological systems, the activities of macromolecular complexes must sometimes be turned off. Thus, a wide variety of protein inhibitors has evolved for this purpose. These inhibitors function through diverse mechanisms, including steric blocking of crucial interactions, enzymatic modification of key residues or substrates, and perturbation of post-translational modifications1. Anti-CRISPRs-proteins that block the activity of CRISPR-Cas systems-are one of the largest groups of inhibitors described, with more than 90 families that function through diverse mechanisms2-4. Here, we characterize the anti-CRISPR AcrIF25, and we show that it inhibits the type I-F CRISPR-Cas system by pulling apart the fully assembled effector complex. AcrIF25 binds to the predominant CRISPR RNA-binding components of this complex, comprising six Cas7 subunits, and strips them from the RNA. Structural and biochemical studies indicate that AcrIF25 removes one Cas7 subunit at a time, starting at one end of the complex. Notably, this feat is achieved with no apparent enzymatic activity. To our knowledge, AcrIF25 is the first example of a protein that disassembles a large and stable macromolecular complex in the absence of an external energy source. As such, AcrIF25 establishes a paradigm for macromolecular complex inhibitors that may be used for biotechnological applications.

iDISCO: a simple, rapid method to immunolabel large tissue samples for volume imaging.

The visualization of molecularly labeled structures within large intact tissues in three dimensions is an area of intense focus. We describe a simple, rapid, and inexpensive method, iDISCO, that permits whole-mount immunolabeling with volume imaging of large cleared samples ranging from perinatal mouse embryos to adult organs, such as brains or kidneys. iDISCO is modeled on classical histology techniques, facilitating translation of section staining assays to intact tissues, as evidenced by compatibility with 28 antibodies to both endogenous antigens and transgenic reporters like GFP. When applied to degenerating neurons, iDISCO revealed unexpected variability in number of apoptotic neurons within individual sensory ganglia despite tight control of total number in all ganglia. It also permitted imaging of single degenerating axons in adult brain and the first visualization of cleaved Caspase-3 in degenerating embryonic sensory axons in vivo, even single axons. iDISCO enables facile volume imaging of immunolabeled structures in complex tissues. PAPERCLIP:

A tripartite rheostat controls self-regulated host plant resistance to insects.

Plants deploy receptor-like kinases and nucleotide-binding leucine-rich repeat receptors to confer host plant resistance (HPR) to herbivores1. These gene-for-gene interactions between insects and their hosts have been proposed for more than 50 years2. However, the molecular and cellular mechanisms that underlie HPR have been elusive, as the identity and sensing mechanisms of insect avirulence effectors have remained unknown. Here we identify an insect salivary protein perceived by a plant immune receptor. The BPH14-interacting salivary protein (BISP) from the brown planthopper (Nilaparvata lugens Stål) is secreted into rice (Oryza sativa) during feeding. In susceptible plants, BISP targets O. satvia RLCK185 (OsRLCK185; hereafter Os is used to denote O. satvia-related proteins or genes) to suppress basal defences. In resistant plants, the nucleotide-binding leucine-rich repeat receptor BPH14 directly binds BISP to activate HPR. Constitutive activation of Bph14-mediated immunity is detrimental to plant growth and productivity. The fine-tuning of Bph14-mediated HPR is achieved through direct binding of BISP and BPH14 to the selective autophagy cargo receptor OsNBR1, which delivers BISP to OsATG8 for degradation. Autophagy therefore controls BISP levels. In Bph14 plants, autophagy restores cellular homeostasis by downregulating HPR when feeding by brown planthoppers ceases. We identify an insect saliva protein sensed by a plant immune receptor and discover a three-way interaction system that offers opportunities for developing high-yield, insect-resistant crops.

A LaCl3-based lithium superionic conductor compatible with lithium metal.

Inorganic superionic conductors possess high ionic conductivity and excellent thermal stability but their poor interfacial compatibility with lithium metal electrodes precludes application in all-solid-state lithium metal batteries1,2. Here we report a LaCl3-based lithium superionic conductor possessing excellent interfacial compatibility with lithium metal electrodes. In contrast to a Li3MCl6 (M = Y, In, Sc and Ho) electrolyte lattice3-6, the UCl3-type LaCl3 lattice has large, one-dimensional channels for rapid Li+ conduction, interconnected by La vacancies via Ta doping and resulting in a three-dimensional Li+ migration network. The optimized Li0.388Ta0.238La0.475Cl3 electrolyte exhibits Li+ conductivity of 3.02 mS cm-1 at 30 °C and a low activation energy of 0.197 eV. It also generates a gradient interfacial passivation layer to stabilize the Li metal electrode for long-term cycling of a Li-Li symmetric cell (1 mAh cm-2) for more than 5,000 h. When directly coupled with an uncoated LiNi0.5Co0.2Mn0.3O2 cathode and bare Li metal anode, the Li0.388Ta0.238La0.475Cl3 electrolyte enables a solid battery to run for more than 100 cycles with a cutoff voltage of 4.35 V and areal capacity of more than 1 mAh cm-2. We also demonstrate rapid Li+ conduction in lanthanide metal chlorides (LnCl3; Ln = La, Ce, Nd, Sm and Gd), suggesting that the LnCl3 solid electrolyte system could provide further developments in conductivity and utility.

Irreversible synthesis of an ultrastrong two-dimensional polymeric material.

Polymers that extend covalently in two dimensions have attracted recent attention1,2 as a means of combining the mechanical strength and in-plane energy conduction of conventional two-dimensional (2D) materials3,4 with the low densities, synthetic processability and organic composition of their one-dimensional counterparts. Efforts so far have proven successful in forms that do not allow full realization of these properties, such as polymerization at flat interfaces5,6 or fixation of monomers in immobilized lattices7-9. Another frequently employed synthetic approach is to introduce microscopic reversibility, at the cost of bond stability, to achieve 2D crystals after extensive error correction10,11. Here we demonstrate a homogenous 2D irreversible polycondensation that results in a covalently bonded 2D polymeric material that is chemically stable and highly processable. Further processing yields highly oriented, free-standing films that have a 2D elastic modulus and yield strength of 12.7 ± 3.8 gigapascals and 488 ± 57 megapascals, respectively. This synthetic route provides opportunities for 2D materials in applications ranging from composite structures to barrier coating materials.

Noncoding RNA circBtnl1 suppresses self-renewal of intestinal stem cells via disruption of Atf4 mRNA stability.

Intestinal stem cells (ISCs) at the crypt base are responsible for the regeneration of the intestinal epithelium. However, how ISC self-renewal is regulated still remains unclear. Here we identified a circular RNA, circBtnl1, that is highly expressed in ISCs. Loss of circBtnl1 in mice enhanced ISC self-renewal capacity and epithelial regeneration, without changes in mRNA and protein levels of its parental gene Btnl1. Mechanistically, circBtnl1 and Atf4 mRNA competitively bound the ATP-dependent RNA helicase Ddx3y to impair the stability of Atf4 mRNA in wild-type ISCs. Furthermore, ATF4 activated Sox9 transcription by binding to its promoter via a unique motif, to enhance the self-renewal capacity and epithelial regeneration of ISCs. In contrast, circBtnl1 knockout promoted Atf4 mRNA stability and enhanced ATF4 expression, which caused Sox9 transcription to potentiate ISC stemness. These data indicate that circBtnl1-mediated Atf4 mRNA decay suppresses Sox9 transcription that negatively modulates self-renewal maintenance of ISCs.

Large-scale across species transcriptomic analysis identifies genetic selection signatures associated with longevity in mammals.

Lifespan varies significantly among mammals, with more than 100-fold difference between the shortest and longest living species. This natural difference may uncover the evolutionary forces and molecular features that define longevity. To understand the relationship between gene expression variation and longevity, we conducted a comparative transcriptomics analysis of liver, kidney, and brain tissues of 103 mammalian species. We found that few genes exhibit common expression patterns with longevity in the three organs analyzed. However, pathways related to translation fidelity, such as nonsense-mediated decay and eukaryotic translation elongation, correlated with longevity across mammals. Analyses of selection pressure found that selection intensity related to the direction of longevity-correlated genes is inconsistent across organs. Furthermore, expression of methionine restriction-related genes correlated with longevity and was under strong selection in long-lived mammals, suggesting that a common strategy is utilized by natural selection and artificial intervention to control lifespan. Our results indicate that lifespan regulation via gene expression is driven through polygenic and indirect natural selection.

Manganese Increases the Sensitivity of the cGAS-STING Pathway for Double-Stranded DNA and Is Required for the Host Defense against DNA Viruses.

Manganese (Mn) is essential for many physiological processes, but its functions in innate immunity remain undefined. Here, we found that Mn2+ was required for the host defense against DNA viruses by increasing the sensitivity of the DNA sensor cGAS and its downstream adaptor protein STING. Mn2+ was released from membrane-enclosed organelles upon viral infection and accumulated in the cytosol where it bound directly to cGAS. Mn2+ enhanced the sensitivity of cGAS to double-stranded DNA (dsDNA) and its enzymatic activity, enabling cGAS to produce secondary messenger cGAMP in the presence of low concentrations of dsDNA that would otherwise be non-stimulatory. Mn2+ also enhanced STING activity by augmenting cGAMP-STING binding affinity. Mn-deficient mice showed diminished cytokine production and were more vulnerable to DNA viruses, and Mn-deficient STING-deficient mice showed no increased susceptibility. These findings indicate that Mn is critically involved and required for the host defense against DNA viruses.

A growth-factor-activated lysosomal K+ channel regulates Parkinson's pathology.

Lysosomes have fundamental physiological roles and have previously been implicated in Parkinson's disease1-5. However, how extracellular growth factors communicate with intracellular organelles to control lysosomal function is not well understood. Here we report a lysosomal K+ channel complex that is activated by growth factors and gated by protein kinase B (AKT) that we term lysoKGF. LysoKGF consists of a pore-forming protein TMEM175 and AKT: TMEM175 is opened by conformational changes in, but not the catalytic activity of, AKT. The minor allele at rs34311866, a common variant in TMEM175, is associated with an increased risk of developing Parkinson's disease and reduces channel currents. Reduction in lysoKGF function predisposes neurons to stress-induced damage and accelerates the accumulation of pathological α-synuclein. By contrast, the minor allele at rs3488217-another common variant of TMEM175, which is associated with a decreased risk of developing Parkinson's disease-produces a gain-of-function in lysoKGF during cell starvation, and enables neuronal resistance to damage. Deficiency in TMEM175 leads to a loss of dopaminergic neurons and impairment in motor function in mice, and a TMEM175 loss-of-function variant is nominally associated with accelerated rates of cognitive and motor decline in humans with Parkinson's disease. Together, our studies uncover a pathway by which extracellular growth factors regulate intracellular organelle function, and establish a targetable mechanism by which common variants of TMEM175 confer risk for Parkinson's disease.

Structures of the glucocorticoid-bound adhesion receptor GPR97-Go complex.

Adhesion G-protein-coupled receptors (GPCRs) are a major family of GPCRs, but limited knowledge of their ligand regulation or structure is available1-3. Here we report that glucocorticoid stress hormones activate adhesion G-protein-coupled receptor G3 (ADGRG3; also known as GPR97)4-6, a prototypical adhesion GPCR. The cryo-electron microscopy structures of GPR97-Go complexes bound to the anti-inflammatory drug beclomethasone or the steroid hormone cortisol revealed that glucocorticoids bind to a pocket within the transmembrane domain. The steroidal core of glucocorticoids is packed against the 'toggle switch' residue W6.53, which senses the binding of a ligand and induces activation of the receptor. Active GPR97 uses a quaternary core and HLY motif to fasten the seven-transmembrane bundle and to mediate G protein coupling. The cytoplasmic side of GPR97 has an open cavity, where all three intracellular loops interact with the Go protein, contributing to the high basal activity of GRP97. Palmitoylation at the cytosolic tail of the Go protein was found to be essential for efficient engagement with GPR97 but is not observed in other solved GPCR complex structures. Our work provides a structural basis for ligand binding to the seven-transmembrane domain of an adhesion GPCR and subsequent G protein coupling.

Structures of Neisseria meningitidis Cas9 Complexes in Catalytically Poised and Anti-CRISPR-Inhibited States.

High-resolution Cas9 structures have yet to reveal catalytic conformations due to HNH nuclease domain positioning away from the cleavage site. Nme1Cas9 and Nme2Cas9 are compact nucleases for in vivo genome editing. Here, we report structures of meningococcal Cas9 homologs in complex with sgRNA, dsDNA, or the AcrIIC3 anti-CRISPR protein. DNA-bound structures represent an early step of target recognition, a later HNH pre-catalytic state, the HNH catalytic state, and a cleaved-target-DNA-bound state. In the HNH catalytic state of Nme1Cas9, the active site is seen poised at the scissile phosphodiester linkage of the target strand, providing a high-resolution view of the active conformation. The HNH active conformation activates the RuvC domain. Our structures explain how Nme1Cas9 and Nme2Cas9 read distinct PAM sequences and how AcrIIC3 inhibits Nme1Cas9 activity. These structures provide insights into Cas9 domain rearrangements, guide-target engagement, cleavage mechanism, and anti-CRISPR inhibition, facilitating the optimization of these genome-editing platforms.

Spatial spillovers of violent conflict amplify the impacts of climate variability on malaria risk in sub-Saharan Africa.

Africa carries a disproportionately high share of the global malaria burden, accounting for 94% of malaria cases and deaths worldwide in 2019. It is also a politically unstable region and the most vulnerable continent to climate change in recent decades. Knowledge about the modifying impacts of violent conflict on climate-malaria relationships remains limited. Here, we quantify the associations between violent conflict, climate variability, and malaria risk in sub-Saharan Africa using health surveys from 128,326 individuals, historical climate data, and 17,429 recorded violent conflicts from 2006 to 2017. We observe that spatial spillovers of violent conflict (SSVCs) have spatially distant effects on malaria risk. Malaria risk induced by SSVCs within 50 to 100 km from the households gradually increases from 0.1% (not significant, P>0.05) to 6.5% (95% CI: 0 to 13.0%). SSVCs significantly promote malaria risk within the average 20.1 to 26.9 °C range. At the 12-mo mean temperature of 22.5 °C, conflict deaths have the largest impact on malaria risk, with an approximately 5.8% increase (95% CI: 1.0 to 11.0%). Additionally, a pronounced association between SSVCs and malaria risk exists in the regions with 9.2 wet days per month. The results reveal that SSVCs increase population exposure to harsh environments, amplifying the effect of warm temperature and persistent precipitation on malaria transmission. Violent conflict therefore poses a substantial barrier to mosquito control and malaria elimination efforts in sub-Saharan Africa. Our findings support effective targeting of treatment programs and vector control activities in conflict-affected regions with a high malaria risk.

GPCR-MAPK signaling pathways underpin fitness trade-offs in whitefly.

Trade-offs between evolutionary gain and loss are prevalent in nature, yet their genetic basis is not well resolved. The evolution of insect resistance to insecticide is often associated with strong fitness costs; however, how the fitness trade-offs operates remains poorly understood. Here, we show that the mitogen-activated protein kinase (MAPK) pathway and its upstream and downstream actors underlie the fitness trade-offs associated with insecticide resistance in the whitefly Bemisia tabaci. Specifically, we find a key cytochrome P450 gene CYP6CM1, that confers neonicotinoids resistance to in B. tabaci, is regulated by the MAPKs p38 and ERK through their activation of the transcription factor cAMP-response element binding protein. However, phosphorylation of p38 and ERK also leads to the activation of the transcription repressor Cap "n" collar isoform C (CncC) that negatively regulates exuperantia (Ex), vasa (Va), and benign gonial cell neoplasm (Bg), key genes involved in oogenesis, leading to abnormal ovary growth and a reduction in female fecundity. We further demonstrate that the transmembrane G protein-coupled receptor (GPCR) neuropeptide FF receptor 2 (NPFF2) triggers the p38 and ERK pathways via phosphorylation. Additionally, a positive feedback loop between p38 and NPFF2 leads to the continuous activation of the MAPK pathways, thereby constitutively promoting neonicotinoids resistance but with a significant reproductive cost. Collectively, these findings provide fundamental insights into the role of cis-trans regulatory networks incurred by GPCR-MAPK signaling pathways in evolutionary trade-offs and applied knowledge that can inform the development of strategies for the sustainable pest control.

Dual mutations in the whitefly nicotinic acetylcholine receptor ß1 subunit confer target-site resistance to multiple neonicotinoid insecticides.

Neonicotinoid insecticides, which target insect nicotinic acetylcholine receptors (nAChRs), have been widely and intensively used to control the whitefly, Bemisia tabaci, a highly damaging, globally distributed, crop pest. This has inevitably led to the emergence of populations with resistance to neonicotinoids. However, to date, there have been no reports of target-site resistance involving mutation of B. tabaci nAChR genes. Here we characterize the nAChR subunit gene family of B. tabaci and identify dual mutations (A58T&R79E) in one of these genes (BTß1) that confer resistance to multiple neonicotinoids. Transgenic D. melanogaster, where the native nAChR Dß1 was replaced with BTß1A58T&R79E, were significantly more resistant to neonicotinoids than flies where Dß1 were replaced with the wildtype BTß1 sequence, demonstrating the causal role of the mutations in resistance. The two mutations identified in this study replace two amino acids that are highly conserved in >200 insect species. Three-dimensional modelling suggests a molecular mechanism for this resistance, whereby A58T forms a hydrogen bond with the R79E side chain, which positions its negatively-charged carboxylate group to electrostatically repulse a neonicotinoid at the orthosteric site. Together these findings describe the first case of target-site resistance to neonicotinoids in B. tabaci and provide insight into the molecular determinants of neonicotinoid binding and selectivity.

Perioperative Pembrolizumab for Early-Stage Non-Small-Cell Lung Cancer.

BACKGROUND: Among patients with resectable early-stage non-small-cell lung cancer (NSCLC), a perioperative approach that includes both neoadjuvant and adjuvant immune checkpoint inhibition may provide benefit beyond either approach alone. METHODS: We conducted a randomized, double-blind, phase 3 trial to evaluate perioperative pembrolizumab in patients with early-stage NSCLC. Participants with resectable stage II, IIIA, or IIIB (N2 stage) NSCLC were assigned in a 1:1 ratio to receive neoadjuvant pembrolizumab (200 mg) or placebo once every 3 weeks, each of which was given with cisplatin-based chemotherapy for 4 cycles, followed by surgery and adjuvant pembrolizumab (200 mg) or placebo once every 3 weeks for up to 13 cycles. The dual primary end points were event-free survival (the time from randomization to the first occurrence of local progression that precluded the planned surgery, unresectable tumor, progression or recurrence, or death) and overall survival. Secondary end points included major pathological response, pathological complete response, and safety. RESULTS: A total of 397 participants were assigned to the pembrolizumab group, and 400 to the placebo group. At the prespecified first interim analysis, the median follow-up was 25.2 months. Event-free survival at 24 months was 62.4% in the pembrolizumab group and 40.6% in the placebo group (hazard ratio for progression, recurrence, or death, 0.58; 95% confidence interval [CI], 0.46 to 0.72; P<0.001). The estimated 24-month overall survival was 80.9% in the pembrolizumab group and 77.6% in the placebo group (P = 0.02, which did not meet the significance criterion). A major pathological response occurred in 30.2% of the participants in the pembrolizumab group and in 11.0% of those in the placebo group (difference, 19.2 percentage points; 95% CI, 13.9 to 24.7; P<0.0001; threshold, P = 0.0001), and a pathological complete response occurred in 18.1% and 4.0%, respectively (difference, 14.2 percentage points; 95% CI, 10.1 to 18.7; P<0.0001; threshold, P = 0.0001). Across all treatment phases, 44.9% of the participants in the pembrolizumab group and 37.3% of those in the placebo group had treatment-related adverse events of grade 3 or higher, including 1.0% and 0.8%, respectively, who had grade 5 events. CONCLUSIONS: Among patients with resectable, early-stage NSCLC, neoadjuvant pembrolizumab plus chemotherapy followed by resection and adjuvant pembrolizumab significantly improved event-free survival, major pathological response, and pathological complete response as compared with neoadjuvant chemotherapy alone followed by surgery. Overall survival did not differ significantly between the groups in this analysis. (Funded by Merck Sharp and Dohme; KEYNOTE-671 ClinicalTrials.gov number, NCT03425643.).

GC-biased gene conversion drives accelerated evolution of ultraconserved elements in mammalian and avian genomes.

Ultraconserved elements (UCEs) are the most conserved regions among the genomes of evolutionarily distant species and are thought to play critical biological functions. However, some UCEs rapidly evolved in specific lineages, and whether they contributed to adaptive evolution is still controversial. Here, using an increased number of sequenced genomes with high taxonomic coverage, we identified 2191 mammalian UCEs and 5938 avian UCEs from 95 mammal and 94 bird genomes, respectively. Our results show that these UCEs are functionally constrained and that their adjacent genes are prone to widespread expression with low expression diversity across tissues. Functional enrichment of mammalian and avian UCEs shows different trends indicating that UCEs may contribute to adaptive evolution of taxa. Focusing on lineage-specific accelerated evolution, we discover that the proportion of fast-evolving UCEs in nine mammalian and 10 avian test lineages range from 0.19% to 13.2%. Notably, up to 62.1% of fast-evolving UCEs in test lineages are much more likely to result from GC-biased gene conversion (gBGC). A single cervid-specific gBGC region embracing the uc.359 allele significantly alters the expression of Nova1 and other neural-related genes in the rat brain. Combined with the altered regulatory activity of ancient gBGC-induced fast-evolving UCEs in eutherians, our results provide evidence that synergy between gBGC and selection shaped lineage-specific substitution patterns, even in the most constrained regulatory elements. In summary, our results show that gBGC played an important role in facilitating lineage-specific accelerated evolution of UCEs, and further support the idea that a combination of multiple evolutionary forces shapes adaptive evolution.