How various factors, including demography, recombination or genome duplication, may impact the efficacy of natural selection and the burden of deleterious mutations,is a central question in evolutionary biology and genetics. In this study, we show that key evolutionary processes, including variations in i) effective population size (Ne) ii) recombination rates and iii) chromosome inheritance, have influenced the genetic load and efficacy of selection in Coho salmon (Oncorhynchus kisutch), a widely distributed salmonid species on the west coast of North America. Using whole genome resequencing data from 14 populations at different migratory distances from their southern glacial refugium, we found evidence supporting gene surfing, wherein reduced Ne at the postglacial recolonization front, leads to a decrease in the efficacy of selection and a surf of deleterious alleles in the northernmost populations. Furthermore, our results indicate that recombination rates play a prime role in shaping the load along the genome. Additionally, we identified variation in polyploidy as a contributing factor to within-genome variation of the load. Overall, our results align remarkably well with expectations under the nearly neutral theory of molecular evolution. We discuss the fundamental and applied implications of these findings for evolutionary and conservation genomics.

TAIMAN (TAI), the only insect ortholog of mammalian Steroid Receptor Coactivators (SRCs), is a critical modulator of ecdysone and juvenile hormone (JH) signaling pathways, which govern insect development and reproduction. The modulatory effect is mediated by JH-dependent TAI's heterodimerization with JH receptor Methoprene-tolerant and association with the Ecdysone Receptor complex. Insect hormones regulate insect physiology and development in concert with abiotic cues, such as photo- and thermoperiod. Here we tested the effects of JH and ecdysone signaling on the circadian clock by a combination of microsurgical operations, application of hormones and hormone mimics, and gene knockdowns in the linden bug Pyrrhocoris apterus males. Silencing taiman by each of three non-overlapping double-strand RNA fragments dramatically slowed the free-running period (FRP) to 27-29 hours, contrasting to 24 hours in controls. To further corroborate TAIMAN's clock modulatory function in the insect circadian clock, we performed taiman knockdown in the cockroach Blattella germanica. Although Blattella and Pyrrhocoris lineages separated ~380 mya, B. germanica taiman silencing slowed the FRP by more than 2 hours, suggesting a conserved TAI clock function in (at least) some insect groups. Interestingly, the pace of the linden bug circadian clock was neither changed by blocking JH and ecdysone synthesis, by application of the hormones or their mimics nor by the knockdown of corresponding hormone receptors. Our results promote TAI as a new circadian clock modulator, a role described for the first time in insects. We speculate that TAI participation in the clock is congruent with the mammalian SRC-2 role in orchestrating metabolism and circadian rhythms, and that TAI/SRCs might be conserved components of the circadian clock in animals.

Transcriptome-wide association studies (TWAS) aim to detect relationships between gene expression and a phenotype, and are commonly used for secondary analysis of genome-wide association study (GWAS) results. Results from TWAS analyses are often interpreted as indicating a genetic relationship between gene expression and a phenotype, but this interpretation is not consistent with the null hypothesis that is evaluated in the traditional TWAS framework. In this study we provide a mathematical outline of this TWAS framework, and elucidate what interpretations are warranted given the null hypothesis it actually tests. We then use both simulations and real data analysis to assess the implications of misinterpreting TWAS results as indicative of a genetic relationship between gene expression and the phenotype. Our simulation results show considerably inflated type 1 error rates for TWAS when interpreted this way, with 41% of significant TWAS associations detected in the real data analysis found to have insufficient statistical evidence to infer such a relationship. This demonstrates that in current implementations, TWAS cannot reliably be used to investigate genetic relationships between gene expression and a phenotype, but that local genetic correlation analysis can serve as a potential alternative.

Circular RNAs (circRNAs) have been recognized as critical regulators of skeletal muscle development. Myocyte enhancer factor 2A (MEF2A) is an evolutionarily conserved transcriptional factor that regulates myogenesis. However, it remains unclear whether MEF2A produces functional circRNAs. In this study, we identified two evolutionarily conserved circular MEF2A RNAs (circMEF2As), namely circMEF2A1 and circMEF2A2, in chicken and mouse muscle stem cells. Our findings revealed that circMEF2A1 promotes myogenesis by regulating the miR-30a-3p/PPP3CA/NFATC1 axis, whereas circMEF2A2 facilitates myogenic differentiation by targeting the miR-148a-5p/SLIT3/ROBO2/β-catenin signaling pathway. Furthermore, in vivo experiments demonstrated that circMEF2As both promote skeletal muscle growth. We also discovered that the linear MEF2A mRNA-derived MEF2A protein binds to its own promoter region, accelerating the transcription of MEF2A and upregulating the expression of both linear MEF2A and circMEF2As, forming a MEF2A autoregulated positive feedback loop. Moreover, circMEF2As positively regulate the expression of linear MEF2A by adsorbing miR-30a-3p and miR-148a-5p, which directly contribute to the MEF2A autoregulated feedback loop. Importantly, we found that mouse circMEF2As are essential for the myogenic differentiation of C2C12 cells. Collectively, our results demonstrated the evolution, function, and underlying mechanisms of circMEF2As in animal myogenesis, which may provide novel insight for both the farm animal meat industry and human medicine.

The nematode Caenorhabditis elegans memorizes various external chemicals, such as ions and odorants, during feeding. Here we find that C. elegans is attracted to the monosaccharides glucose and fructose after exposure to these monosaccharides in the presence of food; however, it avoids them without conditioning. The attraction to glucose requires a gustatory neuron called ASEL. ASEL activity increases when glucose concentration decreases. Optogenetic ASEL stimulation promotes forward movements; however, after glucose conditioning, it promotes turning, suggesting that after glucose conditioning, the behavioral output of ASEL activation switches toward glucose. We previously reported that chemotaxis toward sodium ion (Na+), which is sensed by ASEL, increases after Na+ conditioning in the presence of food. Interestingly, glucose conditioning decreases Na+ chemotaxis, and conversely, Na+ conditioning decreases glucose chemotaxis, suggesting the reciprocal inhibition of learned chemotaxis to distinct chemicals. The activation of PKC-1, an nPKC ε/η ortholog, in ASEL promotes glucose chemotaxis and decreases Na+ chemotaxis after glucose conditioning. Furthermore, genetic screening identified ENSA-1, an ortholog of the protein phosphatase inhibitor ARPP-16/19, which functions in parallel with PKC-1 in glucose-induced chemotactic learning toward distinct chemicals. These findings suggest that kinase-phosphatase signaling regulates the balance between learned behaviors based on glucose conditioning in ASEL, which might contribute to migration toward chemical compositions where the animals were previously fed.

Cas12g is an endonuclease belonging to the type V RNA-guided CRISPR-Cas family. It is known for its ability to cleave RNA substrates using a conserved endonuclease active site located in the RuvC domain. In this study, we determined the crystal structure of apo-Cas12g, the cryo-EM structure of the Cas12g-sgRNA binary complex and investigated conformational changes that occur during the transition from the apo state to the Cas12g-sgRNA binary complex. The conserved zinc finger motifs in Cas12g undergo an ordered-to-disordered transition from the apo to the sgRNA-bound state and their mutations negatively impact on target RNA cleavage. Moreover, we identified a lid motif in the RuvC domain that undergoes transformation from a helix to loop to regulate the access to the RuvC active site and subsequent cleavage of the RNA substrate. Overall, our study provides valuable insights into the mechanisms by which Cas12g recognizes sgRNA and the conformational changes it undergoes from sgRNA binding to the activation of the RNase active site, thereby laying a foundation for the potential repurposing of Cas12g as a tool for RNA-editing.

Transcription factors (TFs) play critical roles in specifying many aspects of neuronal cell fate including dendritic morphology. How TFs are accurately regulated during neuronal morphogenesis is not fully understood. Here, we show that LIM homeodomain protein MEC-3, the key TF for C. elegans PVD dendrite morphogenesis, is regulated by both alternative splicing and an E3 ubiquitin ligase. The mec-3 gene generates several transcripts by alternative splicing. We find that mbl-1, the orthologue of the muscular dystrophy disease gene muscleblind-like (MBNL), is required for PVD dendrite arbor formation. Our data suggest mbl-1 regulates the alternative splicing of mec-3 to produce its long isoform. Deleting the long isoform of mec-3(deExon2) causes reduction of dendrite complexity. Through a genetic modifier screen, we find that mutation in the E3 ubiquitin ligase EEL-1 suppresses mbl-1 phenotype. eel-1 mutants also suppress mec-3(deExon2) mutant but not the mec-3 null phenotype. Loss of EEL-1 alone leads to excessive dendrite branches. Together, these results indicate that MEC-3 is fine-tuned by alternative splicing and the ubiquitin system to produce the optimal level of dendrite branches.

Understanding the mechanism of detoxification initiation in arthropods after pesticide exposure is crucial. Although the identity of transcription factors that induce and regulate the expression of detoxification genes in response to pesticides is beginning to emerge, whether transcription factors directly interact with xenobiotics is unclear. The findings of this study revealed that a nuclear hormone receptor, Tetranychus cinnabarinus hormone receptor (HR) TcHR96h, regulates the overexpression of the detoxification gene TcGSTm02, which is involved in cyflumetofen resistance. The nuclear translocation of TcHR96h increased after cyflumetofen exposure, suggesting direct binding with cyflumetofen. The direct binding of TcHR96h and cyflumetofen was supported by several independent proteomic assays that quantify interactions with small molecules. Together, this study proposes a model for the initiation of xenobiotic detoxification in a polyphagous agricultural pest. These insights not only provide a better understanding of the mechanisms of xenobiotic detoxification and metabolism in arthropods, but also are crucial in understanding adaptation in polyphagous herbivores.

[This corrects the article DOI: 10.1371/journal.pgen.1010798.].

Fluctuating environments threaten fertility and viability. To better match the immediate, local environment, many organisms adopt alternative phenotypic states, a phenomenon called "phenotypic plasticity." Natural populations that predictably encounter fluctuating environments tend to be more plastic than conspecific populations that encounter a constant environment, suggesting that phenotypic plasticity can be adaptive. Despite pervasive evidence of such "adaptive phenotypic plasticity," gene regulatory mechanisms underlying plasticity remains poorly understood. Here we test the hypothesis that environment-dependent phenotypic plasticity is mediated by epigenetic factors. To test this hypothesis, we exploit the adaptive reproductive arrest of Drosophila melanogaster females, called diapause. Using an inbred line from a natural population with high diapause plasticity, we demonstrate that diapause is determined epigenetically: only a subset of genetically identical individuals enter diapause and this diapause plasticity is epigenetically transmitted for at least three generations. Upon screening a suite of epigenetic marks, we discovered that the active histone marks H3K4me3 and H3K36me1 are depleted in diapausing ovaries. Using ovary-specific knockdown of histone mark writers and erasers, we demonstrate that H3K4me3 and H3K36me1 depletion promotes diapause. Given that diapause is highly polygenic, that is, distinct suites of alleles mediate diapause plasticity across distinct genotypes, we also investigated the potential for genetic variation in diapause-determining epigenetic marks. Specifically, we asked if these histone marks were similarly depleted in diapause of a genotypically distinct line. We found evidence of divergence in both the gene expression program and histone mark abundance. This study reveals chromatin determinants of phenotypic plasticity and suggests that these determinants may be genotype-dependent, offering new insight into how organisms may exploit and evolve epigenetic mechanisms to persist in fluctuating environments.

Streptococcus pneumoniae (the pneumococcus) is well known for its ability to develop competence for natural DNA transformation. Competence development is regulated by an autocatalytic loop driven by variations in the basal level of transcription of the comCDE and comAB operons. These genes are part of the early gene regulon that controls expression of the late competence genes known to encode the apparatus of transformation. Several stressful conditions are known to promote competence development, although the induction pathways are remain poorly understood. Here we demonstrate that transient temperature elevation induces an immediate increase in the basal expression level of the comCDE operon and early genes that, in turn, stimulates its full induction, including that of the late competence regulon. This thermal regulation depends on the HtrA chaperone/protease and its proteolytic activity. We find that other competence induction stimulus, like norfloxacin, is not conveyed by the HtrA-dependent pathway. This finding strongly suggests that competence can be induced by at least two independent pathways and thus reinforces the view that competence is a general stress response system in the pneumococcus.

As species expand their geographic ranges, colonizing populations face novel ecological conditions, such as new environments and limited mates, and suffer from evolutionary consequences of demographic change through bottlenecks and mutation load accumulation. Self-fertilization is often observed at species range edges and, in addition to countering the lack of mates, is hypothesized as an evolutionary advantage against load accumulation through increased homozygosity and purging. We study how selfing impacts the accumulation of genetic load during range expansion via purging and/or speed of colonization. Using simulations, we disentangle inbreeding effects due to demography versus due to selfing and find that selfers expand faster, but still accumulate load, regardless of mating system. The severity of variants contributing to this load, however, differs across mating system: higher selfing rates purge large-effect recessive variants leaving a burden of smaller-effect alleles. We compare these predictions to the mixed-mating plant Arabis alpina, using whole-genome sequences from refugial outcrossing populations versus expanded selfing populations. Empirical results indicate accumulation of expansion load along with evidence of purging in selfing populations, concordant with our simulations, suggesting that while purging is a benefit of selfing evolving during range expansions, it is not sufficient to prevent load accumulation due to range expansion.