Search

Background: Elucidating the dynamics of gene expression across developmental stages, including the genomic characteristics of brain expression during infancy, is pivotal in deciphering human psychiatric and neurological disorders and providing insights into developmental disorders. Results: Leveraging comprehensive human GWAS associations with temporal and spatial brain expression data, we discovered a distinctive co-expression cluster comprising 897 genes highly expressed specifically during infancy, enriched in functions related to the neuronal system. This gene cluster notably harbors the highest ratio of genes linked to psychiatric and neurological disorders. Through computational analysis, MYT1L emerged as a potential central transcription factor governing these genes. Remarkably, the infancy-specific expressed genes, including SYT1, exhibit prominent colocalization within human accelerated regions. Additionally, chromatin state analysis unveiled prevalent epigenetic markers associated with enhancer-specific modifications. In addition, this cluster of genes has demonstrated to be specifically highly expressed in cell-types including excitatory neurons, medial ganglionic eminence and caudal ganglionic eminence. Conclusions: This study comprehensively characterizes the genomics and epigenomics of genes specifically expressed during infancy, identifying crucial hub genes and transcription factors. These findings offer valuable insights into early detection strategies and interventions for psychiatric and neurological disorders. [ABSTRACT FROM AUTHOR]

Abstract Background Elucidating the dynamics of gene expression across developmental stages, including the genomic characteristics of brain expression during infancy, is pivotal in deciphering human psychiatric and neurological disorders and providing insights into developmental disorders. Results Leveraging comprehensive human GWAS associations with temporal and spatial brain expression data, we discovered a distinctive co-expression cluster comprising 897 genes highly expressed specifically during infancy, enriched in functions related to the neuronal system. This gene cluster notably harbors the highest ratio of genes linked to psychiatric and neurological disorders. Through computational analysis, MYT1L emerged as a potential central transcription factor governing these genes. Remarkably, the infancy-specific expressed genes, including SYT1, exhibit prominent colocalization within human accelerated regions. Additionally, chromatin state analysis unveiled prevalent epigenetic markers associated with enhancer-specific modifications. In addition, this cluster of genes has demonstrated to be specifically highly expressed in cell-types including excitatory neurons, medial ganglionic eminence and caudal ganglionic eminence. Conclusions This study comprehensively characterizes the genomics and epigenomics of genes specifically expressed during infancy, identifying crucial hub genes and transcription factors. These findings offer valuable insights into early detection strategies and interventions for psychiatric and neurological disorders.

This study systematically rejects the long‐standing notion of cospeciation as the dominant driver of codiversification between flowering plants and their specialist pollinators. Through cophylogenetic analysis of six classical specialized pollination systems, the research finds that cospeciation events are consistently outnumbered by non‐cospeciation events, such as host‐switch, duplication, and association losses. The findings support a more dynamic and diffuse codiversification paradigm, highlighting the importance of considering a broader range of evolutionary events in understanding plant–pollinator codiversification. This new understanding is robust across diverse pollination systems and has significant implications for conservation strategies in the face of environmental change. [ABSTRACT FROM AUTHOR]

In insects, the odorant receptor (OR) multigene family evolves by the birth‐and‐death evolutionary model, according to which the OR repertoire of each species has undergone specific gene gains and losses depending on their chemical environment, resulting in taxon‐specific OR lineage radiations with different sizes in the phylogenetic trees. Despite the general divergence in the gene family across different insect orders, the ORs in moths seem to be genetically conserved across species, clustered into 23 major clades containing multiple orthologous groups with single‐copy gene from each species. We hypothesized that ORs in these orthologous groups are tuned to ecologically important compounds and functionally conserved. cis‐Jasmone is one of the compounds that not only primes the plant defense of neighboring receiver plants, but also functions as a behavior regulator to various insects. To test our hypothesis, using Xenopus oocyte recordings, we functionally assayed the orthologues of BmorOR56, which has been characterized as a specific receptor for cis‐jasmone. Our results showed highly conserved response specificity of the BmorOR56 orthologues, with all receptors within this group exclusively responding to cis‐jasmone. This is supported by the dN/dS analysis, showing that strong purifying selection is acting on this group. Moreover, molecular docking showed that the ligand binding pockets of BmorOR56 orthologues to cis‐jasmone are similar. Taken together, our results suggest the high conservation of OR for ecologically important compounds across Heterocera. [ABSTRACT FROM AUTHOR]

Summary: The metabolism of massively accumulated chlorogenic acid is crucial for the successful germination of purple coneflower (Echinacea purpurea (L.) Menoch). A serine carboxypeptidase‐like (SCPL) acyltransferase (chicoric acid synthase, CAS) utilizes chlorogenic acid to produce chicoric acid during germination. However, it seems that the generation of chicoric acid lags behind the decrease in chlorogenic acid, suggesting an earlier route of chlorogenic acid metabolism.We discovered another chlorogenic acid metabolic product, 3,5‐dicaffeoylquinic acid, which is produced before chicoric acid, filling the lag phase. Then, we identified two additional typical clade IA SCPL acyltransferases, named chlorogenic acid condensing enzymes (CCEs), that catalyze the biosynthesis of 3,5‐dicaffeoylquinic acid from chlorogenic acid with different kinetic characteristics.Chlorogenic acid inhibits radicle elongation in a dose‐dependent manner, explaining the potential biological role of SCPL acyltransferases‐mediated continuous chlorogenic acid metabolism during germination. Both CCE1 and CCE2 are highly conserved among Echinacea species, supporting the observed metabolism of chlorogenic acid to 3,5‐dicaffeoylquinic acid in two Echinacea species without chicoric acid accumulation.The discovery of SCPL acyltransferase involved in the biosynthesis of 3,5‐dicaffeoylquinic acid suggests convergent evolution. Our research clarifies the metabolism strategy of chlorogenic acid in Echinacea species and provides more insight into plant metabolism. [ABSTRACT FROM AUTHOR]

SUMMARY: Chromatin is a dynamic network that regulates genome organization and gene expression. Different types of chromatin regulators are highly conserved among Archaeplastida, including unicellular algae, while some chromatin genes are only present in land plant genomes. Here, we review recent advances in understanding the function of conserved chromatin factors in basal land plants and algae. We focus on the role of Polycomb‐group genes which mediate H3K27me3‐based silencing and play a role in balancing gene dosage and regulating haploid‐to‐diploid transitions by tissue‐specific repression of the transcription factors KNOX and BELL in many representatives of the green lineage. Moreover, H3K27me3 predominantly occupies repetitive elements which can lead to their silencing in a unicellular alga and basal land plants, while it covers mostly protein‐coding genes in higher land plants. In addition, we discuss the role of nuclear matrix constituent proteins as putative functional lamin analogs that are highly conserved among land plants and might have an ancestral function in stress response regulation. In summary, our review highlights the importance of studying chromatin regulation in a wide range of organisms in the Archaeplastida. Significance Statement: Here we review recent findings on the role of chromatin regulators in basal Archaeplastida. These studies shed light on the ancestral function of chromatin factors in unicellular organisms and for land plant evolution. [ABSTRACT FROM AUTHOR]

Inositol-requiring enzyme 1 (IRE1), a mediator of the unfolded protein response, shows the highest degree of evolutionary conservation. Vertebrates express two IRE1 paralogs: IRE1α, which is universally expressed and IRE1β, which shows specific expression within mucus secreted cells in respiratory and gastrointestinal tracts. The biological properties and regulation of the two IRE1 duplicates show evolutionary differences. As recently suggested, IRE1β-deficient mice display impairment in secreted protein expression and mucosal homeostasis. Abnormal changes in IRE1β caused by external and internal factors can disrupt mucosal homeostasis and further lead to respiratory and gastrointestinal diseases. Here, we highlight the physiological functions of IRE1β in the respiratory and gastrointestinal tracts in response to environmental microbes, viruses, toxins, and food components.

During embryonic development, lymphatic endothelial cell (LEC) precursors are distinguished from blood endothelial cells by the expression of Prospero-related homeobox 1 (Prox1), which is essential for lymphatic vasculature formation in mouse and zebrafish. Prox1 expression initiation precedes LEC sprouting and migration, serving as the marker of specified LECs. Despite its crucial role in lymphatic development, Prox1 upstream regulation in LECs remains to be uncovered. SOX18 and COUP-TFII are thought to regulate Prox1 in mice by binding its promoter region. However, the specific regulation of Prox1 expression in LECs remains to be studied in detail. Here, we used evolutionary conservation and chromatin accessibility to identify enhancers located in the proximity of zebrafish prox1a active in developing LECs. We confirmed the functional role of the identified sequences through CRISPR/Cas9 mutagenesis of a lymphatic valve enhancer. The deletion of this region results in impaired valve morphology and function. Overall, our results reveal an intricate control of prox1a expression through a collection of enhancers. Ray-finned fish-specific distal enhancers drive pan-lymphatic expression, whereas vertebrate-conserved proximal enhancers refine expression in functionally distinct subsets of lymphatic endothelium. [ABSTRACT FROM AUTHOR]

Aging is accompanied by the gradual deregulation of the transcriptome. However, whether age-dependent changes in the transcriptome are evolutionarily conserved or diverged remains largely unexplored. Here, we performed a meta-analysis examining the age-dependent changes in the transcriptome using publicly available datasets of 11 representative metazoans, ranging from Caenorhabditis elegans to humans. To identify the transcriptomic changes associated with aging, we analyzed various aspects of the transcriptome, including genome composition, RNA processing, and functional consequences. The use of introns and novel splice sites tended to increase with age, particularly in the brain. In addition, our analysis suggests that the age-dependent accumulation of premature termination codon-containing transcripts is a common feature of aging across multiple animal species. Using C. elegans as a test model, we showed that several splicing factors that are evolutionarily conserved and age-dependently downregulated were required to maintain a normal lifespan. Thus, aberrant RNA processing appears to be associated with aging and a short lifespan in various species. [ABSTRACT FROM AUTHOR]

MiR399 plays an important role in plant growth and development. The objective of the present study was to elucidate the evolutionary characteristics of the MIR399 gene family in grapevine and investigate its role in stress response. To comprehensively investigate the functions of miR399 in grapevine, nine members of the Vvi-MIR399 family were identified based on the genome, using a miRBase database search, located on four chromosomes (Chr 2, Chr 10, Chr 15, and Chr 16). The lengths of the Vvi-miR399 precursor sequences ranged from 82 to 122 nt and they formed stable stem–loop structures, indicating that they could produce microRNAs (miRNAs). Furthermore, our results suggested that the 2 to 20 nt region of miR399 mature sequences were relatively conserved among family members. Phylogenetic analysis revealed that the Vvi-MIR399 members of dicots (Arabidopsis, tomato, and sweet orange) and monocots (rice and grapevine) could be divided into three clades, and most of the Vvi-MIR399s were closely related to sweet orange in dicots. Promoter analysis of Vvi-MIR399s showed that the majority of the predicted cis-elements were related to stress response. A total of 66.7% (6/9) of the Vvi-MIR399 promoters harbored drought, GA, and SA response elements, and 44.4% (4/9) of the Vvi-MIRR399 promoters also presented elements involved in ABA and MeJA response. The expression trend of Vvi-MIR399s was consistent in different tissues, with the lowest expression level in mature and young fruits and the highest expression level in stems and young leaves. However, nine Vvi-MIR399s and four target genes showed different expression patterns when exposed to low light, high light, heat, cold, drought, and salt stress. Interestingly, a putative target of Vvi-MIR399 targeted multiple genes; for example, seven Vvi-MIR399s simultaneously targeted VIT_213s0067g03280.1. Furthermore, overexpression of Vvi_MIR399e and Vvi_MIR399f in Arabidopsis enhanced tolerance to drought compared with wild-type (WT). In contrast, the survival rate of Vvi_MIR399d-overexpressed plants were zero after drought stress. In conclusion, Vvi-MIR399e and Vvi-MIR399f, which are related to drought tolerance in grapevine, provide candidate genes for future drought resistance breeding. [ABSTRACT FROM AUTHOR]

In recent years, several experimental evidences suggest that amino acid repeats are closely linked to many disease conditions, as they have a significant role in evolution of disordered regions of the polypeptide segments. Even though many algorithms and databases were developed for such analysis, each algorithm has some caveats, like limitation on the number of amino acids within the repeat patterns and number of query protein sequences. To this end, in the present work, a new method called the internal sequence repeats across multiple protein sequences (ISRMPS) is proposed for the first time to identify identical repeats across multiple protein sequences. It also identifies distantly located repeat patterns in various protein sequences. Our method can be applied to study evolutionary relationships, epitope mapping, CRISPR-Cas sequencing methods, and other comparative analytical assessments of protein sequences. [ABSTRACT FROM AUTHOR]

Initiation of meiosis in budding yeast does not commit the cells for meiosis. Thus, two distinct signaling cascades may differentially regulate meiosis initiation and commitment in budding yeast. To distinguish between the role of these signaling cascades, we reconstructed protein–protein interaction networks and gene regulatory networks with upregulated genes in meiosis initiation and commitment. Analyzing the integrated networks, we identified four master regulators (MRs) [Ume6p, Msn2p, Met31p, Ino2p], three transcription factors (TFs), and 279 target genes (TGs) unique for meiosis initiation, and three MRs [Ndt80p, Aro80p, Rds2p], 11 TFs, and 948 TGs unique for meiosis commitment. Functional enrichment analysis of these distinct members from the transcriptional cascades for meiosis initiation and commitment revealed that nutritional cues rewire gene expression for initiating meiosis and chromosomal recombination commits cells to meiosis. As meiotic chromosomal recombination is highly conserved in eukaryotes, we compared the evolutionary rate of unique members in the transcriptional cascade of two meiotic phases of Saccharomyces cerevisiae with members of the phylum Ascomycota, revealing that the transcriptional cascade governing chromosomal recombination during meiosis commitment has experienced greater purifying selection pressure (P value = 0.0013, 0.0382, 0.0448, 0.0369, 0.02967, 0.04937, 0.03046, 0.03357 and

SUMMARY: Understanding the underlying mechanisms of plant development is crucial to successfully steer or manipulate plant growth in a targeted manner. Leaves, the primary sites of photosynthesis, are vital organs for many plant species, and leaf growth is controlled by a tight temporal and spatial regulatory network. In this review, we focus on the genetic networks governing leaf cell proliferation, one major contributor to final leaf size. First, we provide an overview of six regulator families of leaf growth in Arabidopsis: DA1, PEAPODs, KLU, GRFs, the SWI/SNF complexes, and DELLAs, together with their surrounding genetic networks. Next, we discuss their evolutionary conservation to highlight similarities and differences among species, because knowledge transfer between species remains a big challenge. Finally, we focus on the increase in knowledge of the interconnectedness between these genetic pathways, the function of the cell cycle machinery as their central convergence point, and other internal and environmental cues. Significance Statement: To understand complex developmental processes such as leaf growth, information on individual gene functions as well as a broader system‐wide view of the process are crucial. Here, we summarize known information on several important leaf cell proliferation‐controlling pathways as well as their interconnections among each other in Arabidopsis thaliana and draw comparisons to other plant species to provide an overview of where current knowledge stands and how this field of study might move forward. [ABSTRACT FROM AUTHOR]

Oocyte maturation is accompanied by changes in abundances of thousands of mRNAs, many degraded and many preferentially stabilized. mRNA stability can be regulated by diverse features including GC content, codon bias, and motifs within the 30- untranslated region (UTR) interacting with RNA binding proteins (RBPs) and miRNAs. Many studies have identified factors participating in mRNA splicing, bulk mRNA storage, and translational recruitment in mammalian oocytes, but the roles of potentially hundreds of expressed factors, how they regulate cohorts of thousands of mRNAs, and to what extent their functions are conserved across species has not been determined. We performed an extensive in silico cross-species analysis of features associated with mRNAs of different stability classes during oocyte maturation (stable, moderately degraded, and highly degraded) for five mammalian species. Using publicly available RNA sequencing data for germinal vesicle (GV) and MII oocyte transcriptomes, we determined that 30-UTR length and synonymous codon usage are positively associated with stability, while greater GC content is negatively associated with stability. By applying machine learning and feature selection strategies, we identified RBPs and miRNAs that are predictive of mRNA stability, including some across multiple species and others more species-restricted. The results provide new insight into the mechanisms regulating maternal mRNA stabilization or degradation. NEW & NOTEWORTHY Conservation across species of mRNA features regulating maternal mRNA stability during mammalian oocyte maturation was analyzed. 30-Untranslated region length and synonymous codon usage are positively associated with stability, while GC content is negatively associated. Just three RNA binding protein motifs were predicted to regulate mRNA stability across all five species examined, but associated pathways and functions are shared, indicating oocytes of different species arrive at comparable physiological destinations via different routes. [ABSTRACT FROM AUTHOR]

The heterochronic genes of Caenorhabditis elegans comprise the best-studied pathway controlling the timing of tissue and organ formation in an animal. To begin to understand the evolution of this pathway and the significance of the relationships among its components, we characterized 11 Caenorhabditis briggsae orthologs of C. elegans heterochronic genes. Using CRISPR/Cas9, we made a variety of alleles and found that several mutant phenotypes differ in significant ways from those of C. elegans. Although most mutant orthologs displayed defects in developmental timing, their phenotypes could differ in which stages were affected, the penetrance and expressivity of the phenotypes, or by having additional pleiotropies that were not obviously connected to developmental timing. However, when examining pairwise epistasis and synergistic relationships, we found those paralleled the known relationships between their C. elegans orthologs, suggesting that the arrangements of these genes in functional modules are conserved, but the modules’ relationships to each other and/or to their targets has drifted since the time of the species’ last common ancestor. Furthermore, our investigation has revealed a relationship between this pathway to other aspects of the animal’s growth and development, including gonad development, which is relevant to both species. [ABSTRACT FROM AUTHOR]

Oocyte maturation is accompanied by changes in abundances of thousands of mRNAs, many degraded and many preferentially stabilized. mRNA stability can be regulated by diverse features including GC content, codon bias, and motifs within the 30-untranslated region (UTR) interacting with RNA binding proteins (RBPs) and miRNAs. Many studies have identified factors participating in mRNA splicing, bulk mRNA storage, and translational recruitment in mammalian oocytes, but the roles of potentially hundreds of expressed factors, how they regulate cohorts of thousands of mRNAs, and to what extent their functions are conserved across species has not been determined. We performed an extensive in silico cross-species analysis of features associated with mRNAs of different stability classes during oocyte maturation (stable, moderately degraded, and highly degraded) for five mammalian species. Using publicly available RNA sequencing data for germinal vesicle (GV) and MII oocyte transcriptomes, we determined that 30-UTR length and synonymous codon usage are positively associated with stability, while greater GC content is negatively associated with stability. By applying machine learning and feature selection strategies, we identified RBPs and miRNAs that are predictive of mRNA stability, including some across multiple species and others more species-restricted. The results provide new insight into the mechanisms regulating maternal mRNA stabilization or degradation. NEW & NOTEWORTHY Conservation across species of mRNA features regulating maternal mRNA stability during mammalian oocyte maturation was analyzed. 30-Untranslated region length and synonymous codon usage are positively associated with stability, while GC content is negatively associated. Just three RNA binding protein motifs were predicted to regulate mRNA stability across all five species examined, but associated pathways and functions are shared, indicating oocytes of different species arrive at comparable physiological destinations via different routes. [ABSTRACT FROM AUTHOR]

Initiation of meiosis in budding yeast does not commit the cells for meiosis. Thus, two distinct signaling cascades may differentially regulate meiosis initiation and commitment in budding yeast. To distinguish between the role of these signaling cascades, we reconstructed protein–protein interaction networks and gene regulatory networks with upregulated genes in meiosis initiation and commitment. Analyzing the integrated networks, we identified four master regulators (MRs) [Ume6p, Msn2p, Met31p, Ino2p], three transcription factors (TFs), and 279 target genes (TGs) unique for meiosis initiation, and three MRs [Ndt80p, Aro80p, Rds2p], 11 TFs, and 948 TGs unique for meiosis commitment. Functional enrichment analysis of these distinct members from the transcriptional cascades for meiosis initiation and commitment revealed that nutritional cues rewire gene expression for initiating meiosis and chromosomal recombination commits cells to meiosis. As meiotic chromosomal recombination is highly conserved in eukaryotes, we compared the evolutionary rate of unique members in the transcriptional cascade of two meiotic phases of Saccharomyces cerevisiae with members of the phylum Ascomycota, revealing that the transcriptional cascade governing chromosomal recombination during meiosis commitment has experienced greater purifying selection pressure (P value = 0.0013, 0.0382, 0.0448, 0.0369, 0.02967, 0.04937, 0.03046, 0.03357 and < 0.00001 for Ashbya gossypii, Yarrowia lipolytica, Debaryomyces hansenii, Aspergillus fumigatus, Neurospora crassa, Kluyveromyces lactis, Schizosaccharomyces pombe, Schizosaccharomyces cryophilus, and Schizosaccharomyces octosporus, respectively). This study demarcates crucial players driving meiosis initiation and commitment and demonstrates their differential rate of evolution in budding yeast. [ABSTRACT FROM AUTHOR]

NK2 genes (NKX2 gene cluster in humans) encode for homeodomain-containing transcription factors that are conserved along the phylogeny. According to the most detailed classifications, vertebrate NKX2 genes are classified into two distinct families, NK2.1 and NK2.2. The former is constituted by NKX2-1 and NKX2-4 genes, which are homologous to the Drosophila scro gene; the latter includes NKX2-2 and NKX2-8 genes, which are homologous to the Drosophila vnd gene. Conservation of these genes is not only related to molecular structure and expression, but also to biological functions. In Drosophila and vertebrates, NK2 genes share roles in the development of ventral regions of the central nervous system. In vertebrates, NKX2 genes have a relevant role in the development of several other organs such as the thyroid, lung, and pancreas. Loss-of-function mutations in NKX2-1 and NKX2-2 are the monogenic cause of the brain-lung-thyroid syndrome and neonatal diabetes, respectively. Alterations in NKX2-4 and NKX2-8 genes may play a role in multifactorial diseases, autism spectrum disorder, and neural tube defects, respectively. NKX2-1, NKX2-2, and NKX2-8 are expressed in various cancer types as either oncogenes or tumor suppressor genes. Several data indicate that evaluation of their expression in tumors has diagnostic and/or prognostic value.

Tumor necrosis factor alpha (TNFa) is an inflammatory cytokine that is involved in the pathogenesis of various inflammatory disorders including rheumatoid arthritis. TNF -alpha receptor I (TNFR1a) is one of the receptors TNFa binds with for its activation. Any variation in this receptor might affect the role of TNFa in successive events. Amino acid residue substitutions might happen in TNFR1a through non-synonymous single nucleotide polymorphisms (nsSNPs) which may alter the functioning of TNFa , hence, identifying any such substitutions is of paramount significance. In this study, six nsSNPs at five different evolutionary conserved regions are predicted to be detrimental to the structure and/or function of TNFR1a by using numerous computational tools. Their 3D models are also proposed in this study. Besides, they were found to reduce the stability and affect the molecular mechanisms of this protein. Two contrasting possibilities might happen because of these substitutions. One, they might reduce the production of TNFa which is overexpressed in inflammatory diseases, hence can play therapeutic role in such diseases. Second, they might possibly hinder the apoptosis to occur which can effectuate the uncontrolled division of cells, hence can be pathogenic in diseases like cancer. Further investigations on these nsSNPs using animal models and at cellular level will open doors to understand the underlying mechanisms behind various diseases. [ABSTRACT FROM AUTHOR]

SUMMARY: Enhancers are critical cis‐regulatory elements controlling gene expression during cell development and differentiation. However, genome‐wide enhancer characterization has been challenging due to the lack of a well‐defined relationship between enhancers and genes. Function‐based methods are the gold standard for determining the biological function of cis‐regulatory elements; however, these methods have not been widely applied to plants. Here, we applied a massively parallel reporter assay on Arabidopsis to measure enhancer activities across the genome. We identified 4327 enhancers with various combinations of epigenetic modifications distinctively different from animal enhancers. Furthermore, we showed that enhancers differ from promoters in their preference for transcription factors. Although some enhancers are not conserved and overlap with transposable elements forming clusters, enhancers are generally conserved across thousand Arabidopsis accessions, suggesting they are selected under evolution pressure and could play critical roles in the regulation of important genes. Moreover, comparison analysis reveals that enhancers identified by different strategies do not overlap, suggesting these methods are complementary in nature. In sum, we systematically investigated the features of enhancers identified by functional assay in A. thaliana, which lays the foundation for further investigation into enhancers' functional mechanisms in plants. Significance Statement: Genome‐wide enhancer characterization has been challenging due to the lack of a well‐defined relationship between enhancers and genes. In this study, we systematically investigated the features of enhancers identified by functional assay in A. thaliana, which lays the foundation for further investigation into enhancers' functional mechanisms in plants. [ABSTRACT FROM AUTHOR]

RNA-sequencing (RNA-seq) technology has led to a surge of neuroscience research using animal models to probe the complex molecular mechanisms underlying brain function and behavior, including substance use disorders. However, findings from rodent studies often fail to be translated into clinical treatments. Here, we developed a novel pipeline for narrowing candidate genes from preclinical studies by translational potential and demonstrated its utility in 2 RNA-seq studies of rodent self-administration. This pipeline uses evolutionary conservation and preferential expression of genes across brain tissues to prioritize candidate genes, increasing the translational utility of RNA-seq in model organisms. Initially, we demonstrate the utility of our prioritization pipeline using an uncorrected P-value. However, we found no differentially expressed genes in either dataset after correcting for multiple testing with false discovery rate (FDR < 0.05 or <0.1). This is likely due to low statistical power that is common across rodent behavioral studies, and, therefore, we additionally illustrate the use of our pipeline on a third dataset with differentially expressed genes corrected for multiple testing (FDR < 0.05). We also advocate for improved RNA-seq data collection, statistical testing, and metadata reporting that will bolster the field's ability to identify reliable candidate genes and improve the translational value of bioinformatics in rodent research. [ABSTRACT FROM AUTHOR]

Abstract Open reading frames (ORFs) with fewer than 100 codons are generally not annotated in genomes, although bona fide genes of that size are known. Newer biochemical studies have suggested that thousands of small protein-coding ORFs (smORFs) may exist in the human genome, but the true number and the biological significance of the micropeptides they encode remain uncertain. Here, we used a comparative genomics approach to identify high-confidence smORFs that are likely protein-coding. We identified 3,326 high-confidence smORFs using constraint within human populations and evolutionary conservation as additional lines of evidence. Next, we validated that, as a group, our high-confidence smORFs are conserved at the amino-acid level rather than merely residing in highly conserved non-coding regions. Finally, we found that high-confidence smORFs are enriched among disease-associated variants from GWAS. Overall, our results highlight that smORF-encoded peptides likely have important functional roles in human disease.

Abstract Background Phenotypic evolution is mainly explained by selection for phenotypic variation arising from factors including mutation and environmental noise. Recent theoretical and experimental studies have suggested that phenotypes with greater developmental stability tend to have a constant phenotype and gene expression level within a particular genetic and environmental condition, and this positively correlates with stronger evolutionary conservation, even after the accumulation of genetic changes. This could reflect a novel mechanism that contributes to evolutionary conservation; however, it remains unclear whether developmental stability is the cause, or whether at least it contributes to their evolutionary conservation. Here, using Japanese medaka lines, we tested experimentally whether developmental stages and gene expression levels with greater stability led to their evolutionary conservation. Results We first measured the stability of each gene expression level and developmental stage (defined here as the whole embryonic transcriptome) in the inbred F0 medaka population. We then measured their evolutionary conservation in the F3 generation by crossing the F0 line with the distantly related Japanese medaka line (Teradomori), followed by two rounds of intra-generational crossings. The results indicated that the genes and developmental stages that had smaller variations in the F0 generation showed lower diversity in the hybrid F3 generation, which implies a causal relationship between stability and evolutionary conservation. Conclusions These findings suggest that the stability in phenotypes, including the developmental stages and gene expression levels, leads to their evolutionary conservation; this most likely occurs due to their low potential to generate phenotypic variation. In addition, since the highly stable developmental stages match with the body-plan-establishment stage, it also implies that the developmental stability potentially contributed to the strict conservation of animal body plan.

Abstract Background Two main subclasses of macrophages are found in almost all solid tissues: embryo-derived resident tissue macrophages and bone marrow-derived infiltrated macrophages. These macrophage subtypes show transcriptional and functional divergence, and the programs that have shaped the evolution of renal macrophages and related signaling pathways remain poorly understood. To clarify these processes, we performed data analysis based on single-cell transcriptional profiling of renal tissue-resident and infiltrated macrophages in human, mouse and rat. Results In this study, we (i) characterized the transcriptional divergence among species and (ii) illustrated variability in expression among cells of each subtype and (iii) compared the gene regulation network and (iv) ligand-receptor pairs in human and mouse. Using single-cell transcriptomics, we mapped the promoter architecture during homeostasis. Conclusions Transcriptionally divergent genes, such as the differentially TF-encoding genes expressed in resident and infiltrated macrophages across the three species, vary among cells and include distinct promoter structures. The gene regulatory network in infiltrated macrophages shows comparatively better species-wide consistency than resident macrophages. The conserved transcriptional gene regulatory network in infiltrated macrophages among species is uniquely enriched in pathways related to kinases, and TFs associated with largely conserved regulons among species are uniquely enriched in kinase-related pathways.

Topologically associating domains (TADs) are fundamental units of three-dimensional (3D) nuclear organization. The regions bordering TADs-TAD boundaries-contribute to the regulation of gene expression by restricting interactions of cis-regulatory sequences to their target genes. TAD and TAD-boundary disruption have been implicated in rare-disease pathogenesis; however, we have a limited framework for integrating TADs and their variation across cell types into the interpretation of common-trait-associated variants. Here, we investigate an attribute of 3D genome architecture-the stability of TAD boundaries across cell types-and demonstrate its relevance to understanding how genetic variation in TADs contributes to complex disease. By synthesizing TAD maps across 37 diverse cell types with 41 genome-wide association studies (GWASs), we investigate the differences in disease association and evolutionary pressure on variation in TADs versus TAD boundaries. We demonstrate that genetic variation in TAD boundaries contributes more to complex-trait heritability, especially for immunologic, hematologic, and metabolic traits. We also show that TAD boundaries are more evolutionarily constrained than TADs. Next, stratifying boundaries by their stability across cell types, we find substantial variation. Compared to boundaries unique to a specific cell type, boundaries stable across cell types are further enriched for complex-trait heritability, evolutionary constraint, CTCF binding, and housekeeping genes. Thus, considering TAD boundary stability across cell types provides valuable context for understanding the genome's functional landscape and enabling variant interpretation that takes 3D structure into account.

MiR399 plays an important role in plant growth and development. The objective of the present study was to elucidate the evolutionary characteristics of the MIR399 gene family in grapevine and investigate its role in stress response. To comprehensively investigate the functions of miR399 in grapevine, nine members of the Vvi-MIR399 family were identified based on the genome, using a miRBase database search, located on four chromosomes (Chr 2, Chr 10, Chr 15, and Chr 16). The lengths of the Vvi-miR399 precursor sequences ranged from 82 to 122 nt and they formed stable stem–loop structures, indicating that they could produce microRNAs (miRNAs). Furthermore, our results suggested that the 2 to 20 nt region of miR399 mature sequences were relatively conserved among family members. Phylogenetic analysis revealed that the Vvi-MIR399 members of dicots (Arabidopsis, tomato, and sweet orange) and monocots (rice and grapevine) could be divided into three clades, and most of the Vvi-MIR399s were closely related to sweet orange in dicots. Promoter analysis of Vvi-MIR399s showed that the majority of the predicted cis-elements were related to stress response. A total of 66.7% (6/9) of the Vvi-MIR399 promoters harbored drought, GA, and SA response elements, and 44.4% (4/9) of the Vvi-MIRR399 promoters also presented elements involved in ABA and MeJA response. The expression trend of Vvi-MIR399s was consistent in different tissues, with the lowest expression level in mature and young fruits and the highest expression level in stems and young leaves. However, nine Vvi-MIR399s and four target genes showed different expression patterns when exposed to low light, high light, heat, cold, drought, and salt stress. Interestingly, a putative target of Vvi-MIR399 targeted multiple genes; for example, seven Vvi-MIR399s simultaneously targeted VIT_213s0067g03280.1. Furthermore, overexpression of Vvi_MIR399e and Vvi_MIR399f in Arabidopsis enhanced tolerance to drought compared with wild-type (WT). In contrast, the survival rate of Vvi_MIR399d-overexpressed plants were zero after drought stress. In conclusion, Vvi-MIR399e and Vvi-MIR399f, which are related to drought tolerance in grapevine, provide candidate genes for future drought resistance breeding.

In contrast to the other pentameric ligand-gated ion channels in the Cys-loop receptor superfamily, the ZACN gene encoding for the Zinc-Activated Channel (ZAC) is exclusively found in the mammalian genome. Human ZAC assembles into homomeric cation-selective channels gated by Zn2+, Cu2+ and H+, but the function of the receptor in human physiology is presently poorly understood. In this study, the degree of evolutionary conservation of a functional ZAC in mammals was probed by investigating the abilities of a selection of ZACs from 10 other mammalian species than human to be expressed at the protein level and assemble into cell surface-expressed functional receptors in mammalian cells and in Xenopus oocytes. In an enzyme-linked immunosorbent assay, transient transfections of tsA201 cells with cDNAs of hemagglutinin (HA)-epitope-tagged versions of these 10 ZACs resulted in robust total expression and cell surface expression levels of all proteins. Moreover, injection of cRNAs for 6 of these ZACs in oocytes resulted in the formation of functional receptors in two-electrode voltage-clamp recordings. The ZACs exhibited robust current amplitudes in response to Zn2+ (10 mM) and H+ (pH 4.0), and the concentration-response relationships displayed by Zn2+ at these channels were largely comparable to that at human ZAC. In conclusion, the findings suggest that the functionality of ZAC at the molecular level may be conserved throughout mammalian species, and that the channel thus may govern physiological functions in mammals, including humans.

Open reading frames (ORFs) with fewer than 100 codons are generally not annotated in genomes, although bona fide genes of that size are known. Newer biochemical studies have suggested that thousands of small protein-coding ORFs (smORFs) may exist in the human genome, but the true number and the biological significance of the micropeptides they encode remain uncertain. Here, we used a comparative genomics approach to identify high-confidence smORFs that are likely protein-coding. We identified 3,326 high-confidence smORFs using constraint within human populations and evolutionary conservation as additional lines of evidence. Next, we validated that, as a group, our high-confidence smORFs are conserved at the amino-acid level rather than merely residing in highly conserved non-coding regions. Finally, we found that high-confidence smORFs are enriched among disease-associated variants from GWAS. Overall, our results highlight that smORF-encoded peptides likely have important functional roles in human disease. [ABSTRACT FROM AUTHOR]

The genomes of RNA viruses encode the information required for replication in host cells both in their linear sequence and in complex higher-order structures. A subset of these RNA genome structures show clear sequence conservation, and have been extensively described for well-characterized viruses. However, the extent to which viral RNA genomes contain functional structural elements--unable to be detected by sequence alone--that nonetheless are critical to viral fitness is largely unknown. Here, we devise a structure-first experimental strategy and use it to identify 22 structure-similar motifs across the coding sequences of the RNA genomes for the four dengue virus serotypes. At least 10 of these motifs modulate viral fitness, revealing a significant unnoticed extent of RNA structure-mediated regulation within viral coding sequences. These viral RNA structures promote a compact global genome architecture, interact with proteins, and regulate the viral replication cycle. These motifs are also thus constrained at the levels of both RNA structure and protein sequence and are potential resistance-refractory targets for antivirals and live-attenuated vaccines. Structure-first identification of conserved RNA structure enables efficient discovery of pervasive RNA-mediated regulation in viral genomes and, likely, other cellular RNAs. [ABSTRACT FROM AUTHOR]

The question of whether amino acids critical to protein folding kinetics are evolutionarily conserved has been investigated intensively in the past, but no consensus has yet been reached. Recently, we have demonstrated that the transition state, dictating folding kinetics, is characterized as the state of maximum dynamic cooperativity, i.e., the state of maximum correlations between amino acid contact formations. Here, we investigate the evolutionary conservation of those amino acids contributing significantly to the dynamic cooperativity. We find a strong indication of a new kind of relationship—necessary but not sufficient causality—between the evolutionary conservation and the dynamic cooperativity: larger contributions to the dynamic cooperativity arise from more conserved residues, but not vice versa. This holds for all the protein systems for which long folding simulation trajectories are available. To our knowledge, this is the first systematic demonstration of any kind of evolutionary conservation of amino acids relevant to folding kinetics. [ABSTRACT FROM AUTHOR]

Background: Phenotypic evolution is mainly explained by selection for phenotypic variation arising from factors including mutation and environmental noise. Recent theoretical and experimental studies have suggested that phenotypes with greater developmental stability tend to have a constant phenotype and gene expression level within a particular genetic and environmental condition, and this positively correlates with stronger evolutionary conservation, even after the accumulation of genetic changes. This could reflect a novel mechanism that contributes to evolutionary conservation; however, it remains unclear whether developmental stability is the cause, or whether at least it contributes to their evolutionary conservation. Here, using Japanese medaka lines, we tested experimentally whether developmental stages and gene expression levels with greater stability led to their evolutionary conservation. Results: We first measured the stability of each gene expression level and developmental stage (defined here as the whole embryonic transcriptome) in the inbred F0 medaka population. We then measured their evolutionary conservation in the F3 generation by crossing the F0 line with the distantly related Japanese medaka line (Teradomori), followed by two rounds of intra-generational crossings. The results indicated that the genes and developmental stages that had smaller variations in the F0 generation showed lower diversity in the hybrid F3 generation, which implies a causal relationship between stability and evolutionary conservation. Conclusions: These findings suggest that the stability in phenotypes, including the developmental stages and gene expression levels, leads to their evolutionary conservation; this most likely occurs due to their low potential to generate phenotypic variation. In addition, since the highly stable developmental stages match with the body-plan-establishment stage, it also implies that the developmental stability potentially contributed to the strict conservation of animal body plan. [ABSTRACT FROM AUTHOR]

Background: Two main subclasses of macrophages are found in almost all solid tissues: embryo-derived resident tissue macrophages and bone marrow-derived infiltrated macrophages. These macrophage subtypes show transcriptional and functional divergence, and the programs that have shaped the evolution of renal macrophages and related signaling pathways remain poorly understood. To clarify these processes, we performed data analysis based on single-cell transcriptional profiling of renal tissue-resident and infiltrated macrophages in human, mouse and rat. Results: In this study, we (i) characterized the transcriptional divergence among species and (ii) illustrated variability in expression among cells of each subtype and (iii) compared the gene regulation network and (iv) ligand-receptor pairs in human and mouse. Using single-cell transcriptomics, we mapped the promoter architecture during homeostasis. Conclusions: Transcriptionally divergent genes, such as the differentially TF-encoding genes expressed in resident and infiltrated macrophages across the three species, vary among cells and include distinct promoter structures. The gene regulatory network in infiltrated macrophages shows comparatively better species-wide consistency than resident macrophages. The conserved transcriptional gene regulatory network in infiltrated macrophages among species is uniquely enriched in pathways related to kinases, and TFs associated with largely conserved regulons among species are uniquely enriched in kinase-related pathways. [ABSTRACT FROM AUTHOR]

The ConSurf web‐sever for the analysis of proteins, RNA, and DNA provides a quick and accurate estimate of the per‐site evolutionary rate among homologues. The analysis reveals functionally important regions, such as catalytic and ligand‐binding sites, which often evolve slowly. Since the last report in 2016, ConSurf has been improved in multiple ways. It now has a user‐friendly interface that makes it easier to perform the analysis and to visualize the results. Evolutionary rates are calculated based on a set of homologous sequences, collected using hidden Markov model‐based search tools, recently embedded in the pipeline. Using these, and following the removal of redundancy, ConSurf assembles a representative set of effective homologues for protein and nucleic acid queries to enable informative analysis of the evolutionary patterns. The analysis is particularly insightful when the evolutionary rates are mapped on the macromolecule structure. In this respect, the availability of AlphaFold model structures of essentially all UniProt proteins makes ConSurf particularly relevant to the research community. The UniProt ID of a query protein with an available AlphaFold model can now be used to start a calculation. Another important improvement is the Python re‐implementation of the entire computational pipeline, making it easier to maintain. This Python pipeline is now available for download as a standalone version. We demonstrate some of ConSurf's key capabilities by the analysis of caveolin‐1, the main protein of membrane invaginations called caveolae. [ABSTRACT FROM AUTHOR]

Abstract Background Elevated temperature and drought stress have substantial impacts on fruit quality, especially in terms of sugar metabolism and content. β-Amylase (BAM) plays a critical role in regulating jujube fruit sugar levels and abiotic stress response. Nevertheless, little is known about the regulatory functions of the BAM genes in jujube fruit. Results Nine jujube BAM genes were identified, clustered into four groups, and characterized to elucidate their structure, function, and distribution. Multiple sequence alignment and gene structure analysis showed that all ZjBAM genes contain Glu-186 and Glu-380 residues and are highly conserved. Phylogenetic and synteny analysis further indicated that the ZjBAM gene family is evolutionarily conserved and formed collinear pairs with the BAM genes of peach, apple, poplar, Arabidopsis thaliana, and cucumber. A single tandem gene pair was found within the ZjBAM gene family and is indicative of putative gene duplication events. We also explored the physicochemical properties, conserved motifs, and chromosomal and subcellular localization of ZjBAM genes as well as the interaction networks and 3D structures of ZjBAM proteins. A promoter cis-acting element analysis suggested that ZjBAM promoters comprise elements related to growth, development, phytohormones, and stress response. Furthermore, a metabolic pathways annotation analysis showed that ZjBAMs are significantly upregulated in the starch and sucrose metabolism, thereby controlling starch-maltose interconversion and hydrolyzing starch to maltose. Transcriptome and qRT-PCR analyses revealed that ZjBAMs respond positively to elevated temperature and drought stress. Specifically, ZjBAM1, ZjBAM2, ZjBAM5, and ZjBAM6 are significantly upregulated in response to severe drought. Bimolecular fluorescence complementation analysis demonstrated ZjBAM1-ZjAMY3, ZjBAM8-ZjDPE1, and ZjBAM7-ZjDPE1 protein interactions that were mainly present in the plasma membrane and nucleus. Conclusion The jujube BAM gene family exhibits high evolutionary conservation. The various expression patterns of ZjBAM gene family members indicate that they play key roles in jujube growth, development, and abiotic stress response. Additionally, ZjBAMs interact with α-amylase and glucanotransferase. Collectively, the present study provides novel insights into the structure, evolution, and functions of the jujube BAM gene family, thus laying a foundation for further exploration of ZjBAM functional mechanisms in response to elevated temperature and drought stress, while opening up avenues for the development of economic forests in arid areas.

The phylotranscriptomic analysis of development in several species revealed the expression of older and more conserved genes in midembryonic stages and younger and more divergent genes in early and late embryonic stages, which supported the hourglass mode of development. However, previous work only studied the transcriptome age of whole embryos or embryonic sublineages, leaving the cellular basis of the hourglass pattern and the variation of transcriptome ages among cell types unexplored. By analyzing both bulk and single-cell transcriptomic data, we studied the transcriptome age of the nematode Caenorhabditis elegans throughout development. Using the bulk RNA-seq data, we identified the morphogenesis phase in midembryonic development as the phylotypic stage with the oldest transcriptome and confirmed the results using whole-embryo transcriptome assembled from single-cell RNA-seq data. The variation in transcriptome ages among individual cell types remained small in early and midembryonic development and grew bigger in late embryonic and larval stages as cells and tissues differentiate. Lineages that give rise to certain tissues (e.g., hypodermis and some neurons) but not all recapitulated the hourglass pattern across development at the single-cell transcriptome level. Further analysis of the variation in transcriptome ages among the 128 neuron types in C. elegans nervous system found that a group of chemosensory neurons and their downstream interneurons expressed very young transcriptomes and may contribute to adaptation in recent evolution. Finally, the variation in transcriptome age among the neuron types, as well as the age of their cell fate regulators, led us to hypothesize the evolutionary history of some neuron types. [ABSTRACT FROM AUTHOR]

Introduction: Cassava (Manihot esculenta) is an annual root cropwhich provides themajorsourceof calories foroverhalf a billionpeoplearoundtheworld.Sinceits domestication ~10,000 years ago, cassava has been largely clonally propagated through stemcuttings. Minimal sexual recombination has led to an accumulation of deleteriousmutations made evident by heavy inbreeding depression. Methods: To locate and characterize these deleterious mutations, and to measure selection pressure across the cassava genome, we aligned 52 related Euphorbiaceae andother related species representingmillions of yearsof evolution. With singlebasepair resolution of genetic conservation, weused protein structuremodels, aminoacid impact, and evolutionary conservation across the Euphorbiaceae to estimate evolutionary constraint. With known deleterious mutations, we aimed to improve genomic evaluations of plant performance through genomic prediction. We first tested this hypothesis through simulation utilizing multi-kernel GBLUP to predict simulated phenotypes across separate populations of cassava. Results: Simulations showed a sizable increase of prediction accuracy when incorporating functional variants in the model when the trait was determined by<100 quantitative trait loci (QTL). Utilizing deleteriousmutations and functional weights informed through evolutionary conservation, we saw improvements in genomic prediction accuracy that were dependent on trait and prediction. Conclusion: Weshowed the potential for using evolutionary information to track functional variation across the genome, in order to improve whole genome trait prediction. We anticipate that continuedwork to improve genotype accuracy and deleterious mutation assessment will lead to improved genomic assessments of cassava clones. [ABSTRACT FROM AUTHOR]

Detailed knowledge of a protein's key residues may assist in understanding its function and designing inhibitors against it. Consequently, such knowledge of one of the severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2)'s proteins is advantageous since the virus is the etiological agent behind one of the biggest health crises of recent times. To that end, we constructed an exhaustive library of bacteria differing from each other by the mutated version of the virus's ORF3a viroporin they harbor. Since the protein is harmful to bacterial growth due to its channel activity, genetic selection followed by deep sequencing could readily identify mutations that abolish the protein's function. Our results have yielded numerous mutations dispersed throughout the sequence that counteract ORF3a's ability to slow bacterial growth. Comparing these data with the conservation pattern of ORF3a within the coronavirinae provided interesting insights: Deleterious mutations obtained in our study corresponded to conserved residues in the protein. However, despite the comprehensive nature of our mutagenesis coverage (108 average mutations per site), we could not reveal all of the protein's conserved residues. Therefore, it is tempting to speculate that our study unearthed positions in the protein pertinent to channel activity, while other conserved residues may correspond to different functionalities of ORF3a. In conclusion, our study provides important information on a key component of SARS‐CoV‐2 and establishes a procedure to analyze other viroporins comprehensively. [ABSTRACT FROM AUTHOR]

The core G protein signaling module, which consists of Gα and extra-large Gα (XLG) subunits coupled with the Gβγ dimer, is a master regulator of various stress responses. In this study, we compared the basal and salt stress-induced transcriptomic, metabolomic and phenotypic profiles in Gα, Gβ, and XLG-null mutants of two plant species, Arabidopsis thaliana and Marchantia polymorpha , and showed that G protein mediates the shift of transcriptional and metabolic homeostasis to stress readiness status. We demonstrated that such stress readiness serves as an intrinsic protection mechanism against further stressors through enhancing the phenylpropanoid pathway and abscisic acid responses. Furthermore, WRKY transcription factors were identified as key intermediates of G protein-mediated homeostatic shifts. Statistical and mathematical model comparisons between A. thaliana and M. polymorpha revealed evolutionary conservation of transcriptional and metabolic networks over land plant evolution, whereas divergence has occurred in the function of plant-specific atypical XLG subunit. Taken together, our results indicate that the shifts in transcriptional and metabolic homeostasis at least partially act as the mechanisms of G protein-coupled stress responses that are conserved between two distantly related plants. Heterotrimeric G protein mutants show altered plant responses to various abiotic and biotic stresses, but the underlying molecular mechanisms remain largely unclear. This study shows that G protein regulates cell stress readiness through modulating WRKY transcriptional pathways and phenylpropanoid metabolic homeostasis. The transcriptional and metabolic networks regulated by G proteins are largely conserved between Arabidopsis thaliana and Marchantia polymorpha. [ABSTRACT FROM AUTHOR]

The results of conventional gene-based analyses which combine epigenome and transcriptome data, including those conducted by the ENCODE/modENCODE projects, suggest various histone modifications performing regulatory functions in controlling mRNA expression (referred to as a histone code) in several model animals. While some histone codes were found to be universally adopted across organisms, “species-specific” histone codes have also been defined. We found that the characterization of these histone codes was confounded by factors (e.g. gene essentiality, expression breadth) that are independent of, but correlated with, gene expression levels. Hence, we attempted to decode histone marks in mouse (Mus musculus), fly (Drosophila melanogaster), and worm (Caenorhabditis elegans) genomes by examining ratios of RNA sequencing (and chromatin immunoprecipitation sequencing) intensities between paralog genes to remove confounding effects that would otherwise be present in a gene-based approach. With this paralog-based approach, associations between four histone modifications (H3K4me3, H3K27ac, H3K9ac, and H3K36me3) and gene expression are substantially revised. For example, we demonstrate that H3K27ac and H3K9ac represent universal active marks in promoters, rather than worm-specific marks as previously reported. Second, acting regions of the studied active marks that are common across species (and across a wide range of tissues at different developmental stages) were found to extend beyond the previously defined regions. Thus, it appears that the active histone codes analyzed have a universality that has previously been underappreciated. Our results suggested that these universal codes, including those previously considered species-specific, could have an ancient origin, and are important in regulating animal gene expression abundance.

IntroductionCassava (Manihot esculenta) is an annual root crop which provides the major source of calories for over half a billion people around the world. Since its domestication ~10,000 years ago, cassava has been largely clonally propagated through stem cuttings. Minimal sexual recombination has led to an accumulation of deleterious mutations made evident by heavy inbreeding depression.MethodsTo locate and characterize these deleterious mutations, and to measure selection pressure across the cassava genome, we aligned 52 related Euphorbiaceae and other related species representing millions of years of evolution. With single base-pair resolution of genetic conservation, we used protein structure models, amino acid impact, and evolutionary conservation across the Euphorbiaceae to estimate evolutionary constraint. With known deleterious mutations, we aimed to improve genomic evaluations of plant performance through genomic prediction. We first tested this hypothesis through simulation utilizing multi-kernel GBLUP to predict simulated phenotypes across separate populations of cassava. ResultsSimulations showed a sizable increase of prediction accuracy when incorporating functional variants in the model when the trait was determined by

Transcription factors (TFs) bind combinatorially to cis-regulatory elements, orchestrating transcriptional programs. Although studies of chromatin state and chromosomal interactions have demonstrated dynamic neurodevelopmental cis-regulatory landscapes, parallel understanding of TF interactions lags. To elucidate combinatorial TF binding driving mouse basal ganglia development, we integrated chromatin immunoprecipitation sequencing (ChIP-seq) for twelve TFs, H3K4me3-associated enhancer-promoter interactions, chromatin and gene expression data, and functional enhancer assays. We identified sets of putative regulatory elements with shared TF binding (TF-pRE modules) that orchestrate distinct processes of GABAergic neurogenesis and suppress other cell fates. The majority of pREs were bound by one or two TFs; however, a small proportion were extensively bound. These sequences had exceptional evolutionary conservation and motif density, complex chromosomal interactions, and activity as in vivo enhancers. Our results provide insights into the combinatorial TF-pRE interactions that activate and repress expression programs during telencephalon neurogenesis and demonstrate the value of TF binding toward modeling developmental transcriptional wiring., Competing Interests: Declaration of interests J.L.R.R. is a co-founder, a stockholder, and currently on the scientific board of Neurona, a company studying the therapeutic use of interneuron transplantation., (Copyright © 2024 The Author(s). Published by Elsevier Inc. All rights reserved.)

The Mycobacterium tuberculosis genome harbours nine toxin-antitoxin (TA) systems of the mazEF family. These consist of two proteins, a toxin and an antitoxin, encoded in an operon. While the toxin has a conserved fold, the antitoxins are structurally diverse and the toxin binding region is typically intrinsically disordered before binding. We describe high throughput methodology for accurate mapping of interfacial residues and apply it to three MazEF complexes. The method involves screening one partner protein against a panel of chemically masked single cysteine mutants of its interacting partner, displayed on the surface of yeast cells. Such libraries have much lower diversity than those generated by saturation mutagenesis, simplifying library generation and data analysis. Further, because of the steric bulk of the masking reagent, labeling of virtually all exposed epitope residues should result in loss of binding, and buried residues are inaccessible to the labeling reagent. The binding residues are deciphered by probing the loss of binding to the labeled cognate partner by flow cytometry. Using this methodology, we have identified the interfacial residues for MazEF3, MazEF6 and MazEF9 TA systems of M. tuberculosis. In the case of MazEF9, where a crystal structure was available, there was excellent agreement between our predictions and the crystal structure, superior to those with AlphaFold2. We also report detailed biophysical characterization of the MazEF3 and MazEF9 TA systems and measured the relative affinities between cognate and non-cognate toxin–antitoxin partners in order to probe possible cross-talk between these systems.

Insects are able to detect a plethora of olfactory cues using a divergent family of odorant receptors (ORs). Despite the divergent nature of this family, related species frequently express several evolutionarily conserved OR orthologues. In the largest order of insects, Coleoptera, it remains unknown whether OR orthologues have conserved or divergent functions in different species. Using HEK293 cells, we addressed this question through functional characterization of two groups of OR orthologues in three species of the Curculionidae (weevil) family, the conifer‐feeding bark beetles Ips typographus L. ("Ityp") and Dendroctonus ponderosae Hopkins ("Dpon") (Scolytinae), and the pine weevil Hylobius abietis L. ("Habi"; Molytinae). The ORs of H. abietis were annotated from antennal transcriptomes. The results show highly conserved response specificities, with one group of orthologues (HabiOR3/DponOR8/ItypOR6) responding exclusively to 2‐phenylethanol (2‐PE), and the other group (HabiOR4/DponOR9/ItypOR5) responding to angiosperm green leaf volatiles (GLVs). Both groups of orthologues belong to the coleopteran OR subfamily 2B, and share a common ancestor with OR5 in the cerambycid Megacyllene caryae, also tuned to 2‐PE, suggesting a shared evolutionary history of 2‐PE receptors across two beetle superfamilies. The detected compounds are ecologically relevant for conifer‐feeding curculionids, and are probably linked to fitness, with GLVs being used to avoid angiosperm nonhost plants, and 2‐PE being important for intraspecific communication and/or playing a putative role in beetle–microbe symbioses. To our knowledge, this study is the first to reveal evolutionary conservation of OR functions across several beetle species and hence sheds new light on the functional evolution of insect ORs. [ABSTRACT FROM AUTHOR]

Proteins sample a multitude of different conformations by undergoing small‐ and large‐scale conformational changes that are often intrinsic to their functions. Information about these changes is often captured in the Protein Data Bank by the apparently redundant deposition of independent structural solutions of identical proteins. Here, we mine these data to examine the conservation of large‐scale conformational changes between homologous proteins. This is important for both practical reasons, such as predicting alternative conformations of a protein by comparative modeling, and conceptual reasons, such as understanding the extent of conservation of different features in evolution. To study this question, we introduce a novel approach to compare conformational changes between proteins by the comparison of their difference distance maps (DDMs). We found that proteins undergoing similar conformational changes have similar DDMs and that this similarity could be quantified by the correlation between the DDMs. By comparing the DDMs of homologous protein pairs, we found that large‐scale conformational changes show a high level of conservation across a broad range of sequence identities. This shows that conformational space is usually conserved between homologs, even relatively distant ones. [ABSTRACT FROM AUTHOR]

Background: Elevated temperature and drought stress have substantial impacts on fruit quality, especially in terms of sugar metabolism and content. β-Amylase (BAM) plays a critical role in regulating jujube fruit sugar levels and abiotic stress response. Nevertheless, little is known about the regulatory functions of the BAM genes in jujube fruit. Results: Nine jujube BAM genes were identified, clustered into four groups, and characterized to elucidate their structure, function, and distribution. Multiple sequence alignment and gene structure analysis showed that all ZjBAM genes contain Glu-186 and Glu-380 residues and are highly conserved. Phylogenetic and synteny analysis further indicated that the ZjBAM gene family is evolutionarily conserved and formed collinear pairs with the BAM genes of peach, apple, poplar, Arabidopsis thaliana, and cucumber. A single tandem gene pair was found within the ZjBAM gene family and is indicative of putative gene duplication events. We also explored the physicochemical properties, conserved motifs, and chromosomal and subcellular localization of ZjBAM genes as well as the interaction networks and 3D structures of ZjBAM proteins. A promoter cis-acting element analysis suggested that ZjBAM promoters comprise elements related to growth, development, phytohormones, and stress response. Furthermore, a metabolic pathways annotation analysis showed that ZjBAMs are significantly upregulated in the starch and sucrose metabolism, thereby controlling starch-maltose interconversion and hydrolyzing starch to maltose. Transcriptome and qRT-PCR analyses revealed that ZjBAMs respond positively to elevated temperature and drought stress. Specifically, ZjBAM1, ZjBAM2, ZjBAM5, and ZjBAM6 are significantly upregulated in response to severe drought. Bimolecular fluorescence complementation analysis demonstrated ZjBAM1-ZjAMY3, ZjBAM8-ZjDPE1, and ZjBAM7-ZjDPE1 protein interactions that were mainly present in the plasma membrane and nucleus. Conclusion: The jujube BAM gene family exhibits high evolutionary conservation. The various expression patterns of ZjBAM gene family members indicate that they play key roles in jujube growth, development, and abiotic stress response. Additionally, ZjBAMs interact with α-amylase and glucanotransferase. Collectively, the present study provides novel insights into the structure, evolution, and functions of the jujube BAM gene family, thus laying a foundation for further exploration of ZjBAM functional mechanisms in response to elevated temperature and drought stress, while opening up avenues for the development of economic forests in arid areas. [ABSTRACT FROM AUTHOR]

Abstract Background Transcription factors (TFs) bind specifically to TF binding sites (TFBSs) at cis-regulatory regions to control transcription. It is critical to locate these TF-DNA interactions to understand transcriptional regulation. Efforts to predict bona fide TFBSs benefit from the availability of experimental data mapping DNA binding regions of TFs (chromatin immunoprecipitation followed by sequencing - ChIP-seq). Results In this study, we processed ~ 10,000 public ChIP-seq datasets from nine species to provide high-quality TFBS predictions. After quality control, it culminated with the prediction of ~ 56 million TFBSs with experimental and computational support for direct TF-DNA interactions for 644 TFs in > 1000 cell lines and tissues. These TFBSs were used to predict > 197,000 cis-regulatory modules representing clusters of binding events in the corresponding genomes. The high-quality of the TFBSs was reinforced by their evolutionary conservation, enrichment at active cis-regulatory regions, and capacity to predict combinatorial binding of TFs. Further, we confirmed that the cell type and tissue specificity of enhancer activity was correlated with the number of TFs with binding sites predicted in these regions. All the data is provided to the community through the UniBind database that can be accessed through its web-interface ( https://unibind.uio.no/ ), a dedicated RESTful API, and as genomic tracks. Finally, we provide an enrichment tool, available as a web-service and an R package, for users to find TFs with enriched TFBSs in a set of provided genomic regions. Conclusions UniBind is the first resource of its kind, providing the largest collection of high-confidence direct TF-DNA interactions in nine species.

Germ layer formation is among the earliest differentiation events in metazoan embryos. In triploblasts, three germ layers are formed, among which the endoderm gives rise to the epithelial lining of the gut tube and associated organs including the liver, pancreas and lungs. In frogs (Xenopus), where early germ layer formation has been studied extensively, the process of endoderm specification involves the interplay of dozens of transcription factors. Here, we review the interactions between these factors, summarized in a transcriptional gene regulatory network (GRN). We highlight regulatory connections conserved between frog, fish, mouse, and human endodermal lineages. Especially prominent is the conserved role and regulatory targets of the Nodal signaling pathway and the T-box transcription factors, Vegt and Eomes. Additionally, we highlight network topologies and motifs, and speculate on their possible roles in development.