Cis-regulatory elements: molecular mechanisms and evolutionary processes underlying divergence

DeepCORE: An interpretable multi-view deep neural network model to detect co-operative regulatory elements

Gene transcription is an essential process involved in all aspects of cellular functions with significant impact on biological traits and diseases. This process is tightly regulated by multiple elements that co-operate to jointly modulate the transcription levels of target genes. To decipher the complicated regulatory network, we present a novel multi-view attention-based deep neural network that models the relationship between genetic, epigenetic, and transcriptional patterns and identifies co-operative regulatory elements (COREs). We applied this new method, named DeepCORE, to predict transcriptomes in various tissues and cell lines, which outperformed the state-of-the-art algorithms. Furthermore, DeepCORE contains an interpreter that extracts the attention values embedded in the deep neural network, maps the attended regions to putative regulatory elements, and infers COREs based on correlated attentions. The identified COREs are significantly enriched with known promoters and enhancers. Novel regulatory elements discovered by DeepCORE showed epigenetic signatures consistent with the status of histone modification marks.

Genome-wide characterization of Related to ABI3/VP1 transcription factors among U's triangle Brassica species reveals a negative role for BnaA06.RAV3L in seed size

The transcription factors Related to ABI3/VP1 (RAV) are crucial for various plant processes and stress responses. Although the U's triangle Brassica species genomes have been released, the knowledge regarding the RAV family is still limited. In this study, we identified 123 putative RAV genes across the six U's triangle Brassica species (Brassica rapa, 14; Brassica oleracea, 14; Brassica nigra, 13; Brassica carinata, 27; Brassica juncea, 28; Brassica napus, 27). Phylogenetic analysis categorized them into three groups. The RAV genes exhibited diversity in both functional and structural aspects, particularly in gene structure and cis-acting elements within their promoters. The expression analysis revealed that BnaRAV genes in Group 1/2 exhibited diverse expression patterns across various tissues, while those in Group 3 did not show expression except for BnaRAV3L-2 and BnaRAV3L-6, which were exclusively expressed in seeds. Furthermore, the seed-specific expression of BnaA06. RAV3L (BnaRAV3L-2) was confirmed through promoter-GUS staining. Subcellular localization studies demonstrated that BnaA06.RAV3L is localized to the nucleus. The overexpression of BnaA06. RAV3L in Arabidopsis led to a remarkable inhibition of seed-specific traits such as seed width, seed length, seed area, and seed weight. This study provides insights into the functional evolution of the RAV gene family in U triangle Brassica species. It establishes a foundation for uncovering the molecular mechanisms underlying the negative role of RAV3L in seed development.

Genome-wide identification of MATE and ALMT genes and their expression profiling in mungbean (Vigna radiata L.) under aluminium stress

The Multidrug and toxic compound extrusion (MATE) and aluminium activated malate transporter (ALMT) gene families are involved in response to aluminium (Al) stress. In this study, we identified 48 MATE and 14 ALMT gene families in Vigna radiata genome and classified into 5 (MATE) and 3 (ALMT) clades by phylogenetic analysis. All the VrMATE and VrALMT genes were distributed across mungbean chromosomes. Tandem duplication was the main driving force for evolution and expansion of MATE gene family. Collinearity of mungbean with soybean indicated that MATE gene family is closely linked to Glycine max. Eight MATE transporters in clade 2 were found to be associated with previously characterized Al tolerance related MATEs in various plant species. Citrate exuding motif (CEM) was present in seven VrMATEs of clade 2. Promoter analysis revealed abundant plant hormone and stress responsive cis-elements. Results from quantitative real time-polymerase chain reaction (qRT-PCR) revealed that VrMATE19, VrMATE30 and VrALMT13 genes were markedly up-regulated at different time points under Al stress. Overall, this study offers a new direction for further molecular characterization of the MATE and ALMT genes in mungbean for Al tolerance.

Whole-genome analysis reveals distinct adaptation signatures to diverse environments in Chinese domestic pigs

BACKGROUND

Long-term natural and artificial selection has resulted in many genetic footprints within the genomes of pig breeds across distinct agroecological zones. Nevertheless, the mechanisms by which these signatures contribute to phenotypic diversity and facilitate environmental adaptation remain unclear.

RESULTS

Here, we leveraged whole-genome sequencing data from 82 individuals from 6 domestic pig breeds originating in tropical, high-altitude, and frigid regions. Population genetic analysis suggested that habitat isolation significantly shaped the genetic diversity and contributed to population stratification in local Chinese pig breeds. Analysis of selection signals revealed regions under selection for adaptation in tropical (55.5 Mb), high-altitude (43.6 Mb), and frigid (17.72 Mb) regions. The potential functions of the selective sweep regions were linked to certain complex traits that might play critical roles in different geographic environments, including fat coverage in frigid environments and blood indicators in tropical and high-altitude environments. Candidate genes under selection were significantly enriched in biological pathways involved in environmental adaptation. These pathways included blood circulation, protein degradation, and inflammation for adaptation to tropical environments; heart and lung development, hypoxia response, and DNA damage repair for high-altitude adaptation; and thermogenesis, cold-induced vasodilation (CIVD), and the cell cycle for adaptation to frigid environments. By examining the chromatin state of the selection signatures, we identified the lung and ileum as two candidate functional tissues for environmental adaptation. Finally, we identified a mutation (chr1: G246,175,129A) in the cis-regulatory region of ABCA1 as a plausible promising variant for adaptation to tropical environments.

CONCLUSIONS

In this study, we conducted a genome-wide exploration of the genetic mechanisms underlying the adaptability of local Chinese pig breeds to tropical, high-altitude, and frigid environments. Our findings shed light on the prominent role of cis-regulatory elements in environmental adaptation in pigs and may serve as a valuable biological model of human plateau-related disorders and cardiovascular diseases.

Integration of chromosome locations and functional aspects of enhancers and topologically associating domains in knowledge graphs enables versatile queries about gene regulation

Knowledge about transcription factor binding and regulation, target genes, cis-regulatory modules and topologically associating domains is not only defined by functional associations like biological processes or diseases but also has a determinative genome location aspect. Here, we exploit these location and functional aspects together to develop new strategies to enable advanced data querying. Many databases have been developed to provide information about enhancers, but a schema that allows the standardized representation of data, securing interoperability between resources, has been lacking. In this work, we use knowledge graphs for the standardized representation of enhancers and topologically associating domains, together with data about their target genes, transcription factors, location on the human genome, and functional data about diseases and gene ontology annotations. We used this schema to integrate twenty-five enhancer datasets and two domain datasets, creating the most powerful integrative resource in this field to date. The knowledge graphs have been implemented using the Resource Description Framework and integrated within the open-access BioGateway knowledge network, generating a resource that contains an interoperable set of knowledge graphs (enhancers, TADs, genes, proteins, diseases, GO terms, and interactions between domains). We show how advanced queries, which combine functional and location restrictions, can be used to develop new hypotheses about functional aspects of gene expression regulation.

Diverse Gene Regulatory Mechanisms Alter Rattlesnake Venom Gene Expression at Fine Evolutionary Scales

Understanding and predicting the relationships between genotype and phenotype is often challenging, largely due to the complex nature of eukaryotic gene regulation. A step towards this goal is to map how phenotypic diversity evolves through genomic changes that modify gene regulatory interactions. Using the Prairie Rattlesnake (Crotalus viridis) and related species, we integrate mRNA-seq, proteomic, ATAC-seq and whole-genome resequencing data to understand how specific evolutionary modifications to gene regulatory network components produce differences in venom gene expression. Through comparisons within and between species, we find a remarkably high degree of gene expression and regulatory network variation across even a shallow level of evolutionary divergence. We use these data to test hypotheses about the roles of specific trans-factors and cis-regulatory elements, how these roles may vary across venom genes and gene families, and how variation in regulatory systems drive diversity in venom phenotypes. Our results illustrate that differences in chromatin and genotype at regulatory elements play major roles in modulating expression. However, we also find that enhancer deletions, differences in transcription factor expression, and variation in activity of the insulator protein CTCF also likely impact venom phenotypes. Our findings provide insight into the diversity and gene-specificity of gene regulatory features and highlight the value of comparative studies to link gene regulatory network variation to phenotypic variation.

Identification of DNase I hypersensitive sites in the human genome by multiple sequence descriptors

DNase I hypersensitive sites (DHSs) are chromatin regions highly sensitive to DNase I enzymes. Studying DHSs is crucial for understanding complex transcriptional regulation mechanisms and localizing cis-regulatory elements (CREs). Numerous studies have indicated that disease-related loci are often enriched in DHSs regions, underscoring the importance of identifying DHSs. Although wet experiments exist for DHSs identification, they are often labor-intensive. Therefore, there is a strong need to develop computational methods for this purpose. In this study, we used experimental data to construct a benchmark dataset. Seven feature extraction methods were employed to capture information about human DHSs. The F-score was applied to filter the features. By comparing the prediction performance of various classification algorithms through five-fold cross-validation, random forest was proposed to perform the final model construction. The model could produce an overall prediction accuracy of 0.859 with an AUC value of 0.837. We hope that this model can assist scholars conducting DNase research in identifying these sites.

Genome-wide identification, expression analysis of WRKY transcription factors in Citrus ichangensis and functional validation of CiWRKY31 in response to cold stress

BACKGROUND

Ichang papeda (Citrus ichangensis), a wild perennial plant of the Rutaceae family, is a cold-hardy plant. WRKY transcription factors are crucial regulators of plant growth and development as well as abiotic stress responses. However, the WRKY genes in C. ichangensis (CiWRKY) and their expression patterns under cold stress have not been thoroughly investigated, hindering our understanding of their role in cold tolerance.

RESULTS

In this study, a total of 52 CiWRKY genes identified in the genome of C. ichangensis were classified into three main groups and five subgroups based on phylogenetic analysis. Comprehensive analyses of motif features, conserved domains, and gene structures were performed. Segmental duplication plays a significant role in the CiWRKY gene family expansion. Cis-acting element analysis revealed the presence of various stress-responsive elements in the promoters of the majority of CiWRKYs. Gene ontology (GO) analysis and protein-protein interaction predictions indicate that the CiWRKYs exhibit crucial roles in regulation of both development and stress response. Expression profiling analysis demonstrates that 14 CiWRKYs were substantially induced under cold stress. Virus-induced gene silencing (VIGS) assay confirmed that CiWRKY31, one of the cold-induced WRKYs, functions positively in regulation of cold tolerance.

CONCLUSION

Sequence and protein properties of CiWRKYs were systematically analyzed. Among the 52 CiWRKY genes 14 members exhibited cold-responsive expression patterns, and CiWRKY31 was verified to be a positive regulator of cold tolerance. These findings pave way for future investigations to understand the molecular functions of CiWRKYs in cold tolerance and contribute to unravelling WRKYs that may be used for engineering cold tolerance in citrus.

Genomic determinants, architecture, and constraints in drought-related traits in Corymbia calophylla

BACKGROUND

Drought adaptation is critical to many tree species persisting under climate change, however our knowledge of the genetic basis for trees to adapt to drought is limited. This knowledge gap impedes our fundamental understanding of drought response and application to forest production and conservation. To improve our understanding of the genomic determinants, architecture, and trait constraints, we assembled a reference genome and detected ~ 6.5 M variants in 432 phenotyped individuals for the foundational tree Corymbia calophylla.

RESULTS

We found 273 genomic variants determining traits with moderate heritability (h2SNP = 0.26-0.64). Significant variants were predominantly in gene regulatory elements distributed among several haplotype blocks across all chromosomes. Furthermore, traits were constrained by frequent epistatic and pleiotropic interactions.

CONCLUSIONS

Our results on the genetic basis for drought traits in Corymbia calophylla have several implications for the ability to adapt to climate change: (1) drought related traits are controlled by complex genomic architectures with large haplotypes, epistatic, and pleiotropic interactions; (2) the most significant variants determining drought related traits occurred in regulatory regions; and (3) models incorporating epistatic interactions increase trait predictions. Our findings indicate that despite moderate heritability drought traits are likely constrained by complex genomic architecture potentially limiting trees response to climate change.

Defining cis-regulatory elements and transcription factors that control human cortical interneuron development

Although human cortical interneurons (cINs) are a minority population in the cerebral cortex, disruption of interneuron development is a frequent contributor to neurodevelopmental disorders. Here, we utilized a model for deriving cINs from human embryonic stem cells to profile chromatin state changes and generate an atlas of cis-regulatory elements (CREs) controlling human cIN development. We used these data to define candidate transcription factors (TFs) that may bind these CREs to regulate interneuron progenitor specification. Among these were RFX3 and RFX4, risk genes for autism spectrum disorder (ASD) with uncharacterized roles in human neuronal development. Using RFX3 and RFX4 knockdown models, we demonstrated new requirements for both genes in interneuron progenitor specification, with RFX3 deficiency causing precocious neuronal differentiation while RFX4 deficiency instead resulted in cessation of progenitor cell proliferation. Together, this work both defined central features of cis-regulatory control and identified new TF requirements for human interneuron development.

An oncogenic enhancer promotes melanoma progression via regulating ETV4 expression

BACKGROUND

Enhancers are important gene regulatory elements that promote the expression of critical genes in development and disease. Aberrant enhancer can modulate cancer risk and activate oncogenes that lead to the occurrence of various cancers. However, the underlying mechanism of most enhancers in cancer remains unclear. Here, we aim to explore the function and mechanism of a crucial enhancer in melanoma.

METHODS

Multi-omics data were applied to identify an enhancer (enh17) involved in melanoma progression. To evaluate the function of enh17, CRISPR/Cas9 technology were applied to knockout enh17 in melanoma cell line A375. RNA-seq, ChIP-seq and Hi-C data analysis integrated with luciferase reporter assay were performed to identify the potential target gene of enh17. Functional experiments were conducted to further validate the function of the target gene ETV4. Multi-omics data integrated with CUT&Tag sequencing were performed to validate the binding profile of the inferred transcription factor STAT3.

RESULTS

An enhancer, named enh17 here, was found to be aberrantly activated and involved in melanoma progression. CRISPR/Cas9-mediated deletion of enh17 inhibited cell proliferation, migration, and tumor growth of melanoma both in vitro and in vivo. Mechanistically, we identified ETV4 as a target gene regulated by enh17, and functional experiments further support ETV4 as a target gene that is involved in cancer-associated phenotypes. In addition, STAT3 acts as a transcription factor binding with enh17 to regulate the transcription of ETV4.

CONCLUSIONS

Our findings revealed that enh17 plays an oncogenic role and promotes tumor progression in melanoma, and its transcriptional regulatory mechanisms were fully elucidated, which may open a promising window for melanoma prevention and treatment.

Genome-wide identification and characterization of DTX family genes highlighting their locations, functions, and regulatory factors in banana (Musa acuminata)

The detoxification efflux carriers (DTX) are a significant group of multidrug efflux transporter family members that play diverse functions in all kingdoms of living organisms. However, genome-wide identification and characterization of DTX family transporters have not yet been performed in banana, despite its importance as an economic fruit plant. Therefore, a detailed genome-wide analysis of DTX family transporters in banana (Musa acuminata) was conducted using integrated bioinformatics and systems biology approaches. In this study, a total of 37 DTX transporters were identified in the banana genome and divided into four groups (I, II, III, and IV) based on phylogenetic analysis. The gene structures, as well as their proteins' domains and motifs, were found to be significantly conserved. Gene ontology (GO) annotation revealed that the predicted DTX genes might play a vital role in protecting cells and membrane-bound organelles through detoxification mechanisms and the removal of drug molecules from banana cells. Gene regulatory analyses identified key transcription factors (TFs), cis-acting elements, and post-transcriptional regulators (miRNAs) of DTX genes, suggesting their potential roles in banana. Furthermore, the changes in gene expression levels due to pathogenic infections and non-living factor indicate that banana DTX genes play a role in responses to both biotic and abiotic stresses. The results of this study could serve as valuable tools to improve banana quality by protecting them from a range of environmental stresses.

Osteoarthritis as an Enhanceropathy: Gene Regulation in Complex Musculoskeletal Disease

PURPOSE OF REVIEW

Osteoarthritis is a complex and highly polygenic disease. Over 100 reported osteoarthritis risk variants fall in non-coding regions of the genome, ostensibly conferring functional effects through the disruption of regulatory elements impacting target gene expression. In this review, we summarise the progress that has advanced our knowledge of gene enhancers both within the field of osteoarthritis and more broadly in complex diseases.

RECENT FINDINGS

Advances in technologies such as ATAC-seq have facilitated our understanding of chromatin states in specific cell types, bolstering the interpretation of GWAS and the identification of effector genes. Their application to osteoarthritis research has revealed enhancers as the principal regulatory element driving disease-associated changes in gene expression. However, tissue-specific effects in gene regulatory mechanisms can contribute added complexity to biological interpretation. Understanding gene enhancers and their altered activity in specific cell and tissue types is the key to unlocking the genetic complexity of osteoarthritis. The use of single-cell technologies in osteoarthritis research is still in its infancy. However, such tools offer great promise in improving our functional interpretation of osteoarthritis GWAS and the identification of druggable targets. Large-scale collaborative efforts will be imperative to understand tissue and cell-type specific molecular mechanisms underlying enhancer function in disease.

Genomic signatures associated with recurrent scale loss in cyprinid fish

Scale morphology represents a fundamental feature of fish and a key evolutionary trait underlying fish diversification. Despite frequent and recurrent scale loss throughout fish diversification, comprehensive genome-wide analyses of the genomic signatures associated with scale loss in divergent fish lineages remain scarce. In the current study, we investigated genome-wide signatures, specifically convergent protein-coding gene loss, amino acid substitutions, and cis-regulatory sequence changes, associated with recurrent scale loss in two divergent Cypriniformes lineages based on large-scale genomic, transcriptomic, and epigenetic data. Results demonstrated convergent changes in many genes related to scale formation in divergent scaleless fish lineages, including loss of P/Q-rich scpp genes (e.g. scpp6 and scpp7), accelerated evolution of non-coding elements adjacent to the fgf and fgfr genes, and convergent amino acid changes in genes (e.g. snap29) under relaxed selection. Collectively, these findings highlight the existence of a shared genetic architecture underlying recurrent scale loss in divergent fish lineages, suggesting that evolutionary outcomes may be genetically repeatable and predictable in the convergence of scale loss in fish.

Rapid response of fly populations to gene dosage across development and generations

Although the effects of genetic and environmental perturbations on multicellular organisms are rarely restricted to single phenotypic layers, our current understanding of how developmental programs react to these challenges remains limited. Here, we have examined the phenotypic consequences of disturbing the bicoid regulatory network in early Drosophila embryos. We generated flies with two extra copies of bicoid, which causes a posterior shift of the network's regulatory outputs and a decrease in fitness. We subjected these flies to EMS mutagenesis, followed by experimental evolution. After only 8-15 generations, experimental populations have normalized patterns of gene expression and increased survival. Using a phenomics approach, we find that populations were normalized through rapid increases in embryo size driven by maternal changes in metabolism and ovariole development. We extend our results to additional populations of flies, demonstrating predictability. Together, our results necessitate a broader view of regulatory network evolution at the systems level.

Genome-Wide Identification and Expression Analysis of Bx Involved in Benzoxazinoids Biosynthesis Revealed the Roles of DIMBOA during Early Somatic Embryogenesis in Dimocarpus longan Lour

Benzoxazinoids (BXs) are tryptophan-derived indole metabolites and play a role in various physiological processes, such as auxin metabolism. Auxin is essential in the process of somatic embryogenesis (SE) in plants. In this study, we used bioinformatics, transcriptome data, exogenous treatment experiments, and qPCR analysis to study the evolutionary pattern of Bx genes in green plants, the regulatory mechanism of DlBx genes during early SE, and the effect of 2,4-dihydroxy-7-methoxy-1,4-benzoxazine-3-one (DIMBOA) on the early SE in Dimocarpus longan Lour. The results showed that 27 putative DlBxs were identified in the longan genome; the Bx genes evolved independently in monocots and dicots, and the main way of gene duplication for the DlBx was tandem duplication (TD) and the DlBx were strongly constrained by purification selection during evolution. The transcriptome data indicated varying expression levels of DlBx during longan early SE, and most DlBxs responded to light, temperature, drought stress, and 2,4-dichlorophenoxyacetic acid (2,4-D) treatment; qRT-PCR results showed DlBx1, DlBx6g and DlBx6h were responsive to auxin, and treatment with 0.1mg/L DIMBOA for 9 days significantly upregulated the expression levels of DlBx1, DlBx3g, DlBx6c, DlBx6f, DlB6h, DlBx7d, DlBx8, and DlBx9b. The correlation analysis showed a significantly negative correlation between the expression level of DlBx1 and the endogenous IAA contents; DIMBOA significantly promoted the early SE and significantly changed the endogenous IAA content, and the IAA content increased significantly at the 9th day and decreased significantly at the 13th day. Therefore, the results suggested that DIMBOA indirectly promote the early SE by changing the endogenous IAA content via affecting the expression level of DlBx1 and hydrogen peroxide (H2O2) content in longan.

Genome-wide identification of the walnut MYC gene family and functional characterization of Xinjiang wild walnut under low-temperature stress

Introduction: MYC transcription factors are the basic regulators of the jasmonic acid signaling pathway and play important roles in plant growth and development and the response to adverse stress. In recent years, severe winter freezing and late spring frost in the main planting area of walnut in Xinjiang have affected the growth and development of walnut, which has become a prominent problem restricting walnut production. Xinjiang wild walnut is the only remaining wild species of walnuts in China, which contains a lot of genes with excellent traits, and is important for the cultivation and breeding. Methods: In this paper, the physicochemical properties and bioinformatics of MYC transcription factor members in walnut were analyzed, and the nine MYC were screened from the transcriptome data under low temperature stress. At last, we study the subcellular localizations and the expression patterns of the nine MYC members in Xinjiang wild walnut. Results: The results revealed that 30 MYC members were identified from published walnut whole-genome data, and their evolutionary relationships with Arabidopsis and poplar were divided into six groups according to clustering analysis, among which JrMYC22 and JrMYC23 had high homology with PtrMYC2b, which is induced by jasmonic acid in response to low-temperature stress. Walnut MYC members are unevenly distributed on 12 chromosomes. The prediction of promoter cis-acting elements of walnut MYC transcription factor family members revealed that cis-acting elements related to jasmonic acid and lowtemperature stress were the ones with the greatest number of members, with 12. In addition, all nine MYC family members in Xinjiang wild walnut plants responding to low-temperature stress exhibited strong fluorescence responses in the nucleus. The expression levels of these members in response to low-temperature stress revealed that JrMYC28, JrMYC31, JrMYC33, JrMYC34, and JrMYC35 were highly expressed, and it was hypothesized that JrMYC28, JrMYC31, JrMYC33, JrMYC34, and JrMYC35 might play a key role in the response to lowtemperature stress. Discussion: The results of this study provide a theoretical basis for further research on the functional mechanisms of the MYC transcription factor family members in walnut.

Genes and gene networks underlying spatial cognition in food-caching chickadees

Substantial progress has been made in understanding the genetic architecture of phenotypes involved in a variety of evolutionary processes. Behavioral genetics remains, however, among the least understood. We explore the genetic architecture of spatial cognitive abilities in a wild passerine bird, the mountain chickadee (Poecile gambeli). Mountain chickadees cache thousands of seeds in the fall and require specialized spatial memory to recover these caches throughout the winter. We previously showed that variation in spatial cognition has a direct effect on fitness and has a genetic basis. It remains unknown which specific genes and developmental pathways are particularly important for shaping spatial cognition. To further dissect the genetic basis of spatial cognitive abilities, we combine experimental quantification of spatial cognition in wild chickadees with whole-genome sequencing of 162 individuals, a new chromosome-scale reference genome, and species-specific gene annotation. We have identified a set of genes and developmental pathways that play a key role in creating variation in spatial cognition and found that the mechanism shaping cognitive variation is consistent with selection against mildly deleterious non-coding mutations. Although some candidate genes were organized into connected gene networks, about half do not have shared regulation, highlighting that multiple independent developmental or physiological mechanisms contribute to variation in spatial cognitive abilities. A large proportion of the candidate genes we found are associated with synaptic plasticity, an intriguing result that leads to the hypothesis that certain genetic variants create antagonism between behavioral plasticity and long-term memory, each providing distinct benefits depending on ecological context.

A genome-wide genetic association study reveals SNPs significantly associated with environmental variables and specific leaf area in European beech

European beech is negatively affected by climate change and a further growth decline is predicted for large parts of its distribution range. Despite the importance of this species, little is known about its genetic adaptation and especially the genetic basis of its physiological traits. Here, we used genotyping by sequencing to identify SNPs in 43 German European beech populations growing under different environmental conditions. In total, 28 of these populations were located along a precipitation and temperature gradient in northern Germany, and single tree-based hydraulic and morphological traits were available. We obtained a set of 13,493 high-quality SNPs that were used for environmental and SNP-trait association analysis. In total, 22 SNPs were identified that were significantly associated with environmental variables or specific leaf area (SLA). Several SNPs were located in genes related to stress response. The majority of the significant SNPs were located in non-coding (intergenic and intronic) regions. These may be in linkage disequilibrium with the causative coding or regulatory regions. Our study gives insights into the genetic basis of abiotic adaptation in European beech, and provides genetic resources that can be used in future studies on this species. Besides clear patterns of local adaptation to environmental conditions of the investigated populations, the analyzed morphological and hydraulic traits explained most of the explainable genetic variation. Thus, they could successfully be altered in tree breeding programs, which may help to increase the adaptation of European beech to changing environmental conditions in the future.

Emx2 underlies the development and evolution of marsupial gliding membranes

Phenotypic variation among species is a product of evolutionary changes to developmental programs1,2. However, how these changes generate novel morphological traits remains largely unclear. Here we studied the genomic and developmental basis of the mammalian gliding membrane, or patagium-an adaptative trait that has repeatedly evolved in different lineages, including in closely related marsupial species. Through comparative genomic analysis of 15 marsupial genomes, both from gliding and non-gliding species, we find that the Emx2 locus experienced lineage-specific patterns of accelerated cis-regulatory evolution in gliding species. By combining epigenomics, transcriptomics and in-pouch marsupial transgenics, we show that Emx2 is a critical upstream regulator of patagium development. Moreover, we identify different cis-regulatory elements that may be responsible for driving increased Emx2 expression levels in gliding species. Lastly, using mouse functional experiments, we find evidence that Emx2 expression patterns in gliders may have been modified from a pre-existing program found in all mammals. Together, our results suggest that patagia repeatedly originated through a process of convergent genomic evolution, whereby regulation of Emx2 was altered by distinct cis-regulatory elements in independently evolved species. Thus, different regulatory elements targeting the same key developmental gene may constitute an effective strategy by which natural selection has harnessed regulatory evolution in marsupial genomes to generate phenotypic novelty.

Genome-Wide Characterization and Functional Validation of the ACS Gene Family in the Chestnut Reveals Its Regulatory Role in Ovule Development

Ovule abortion significantly contributes to a reduction in chestnut yield. Therefore, an examination of the mechanisms underlying ovule abortion is crucial for increasing chestnut yield. In our previous study, we conducted a comprehensive multiomic analysis of fertile and abortive ovules and found that ACS genes in chestnuts (CmACS) play a crucial role in ovule development. Therefore, to further study the function of ACS genes, a total of seven CmACS members were identified, their gene structures, conserved structural domains, evolutionary trees, chromosomal localization, and promoter cis-acting elements were analyzed, and their subcellular localization was predicted and verified. The spatiotemporal specificity of the expression of the seven CmACS genes was confirmed via qRT-PCR analysis. Notably, CmACS7 was exclusively expressed in the floral organs, and its expression peaked during fertilization and decreased after fertilization. The ACC levels remained consistently greater in fertile ovules than in abortive ovules. The ACSase activity of CmACS7 was identified using the genetic transformation of chestnut healing tissue. Micro Solanum lycopersicum plants overexpressing CmACS7 had a significantly greater rate of seed failure than did wild-type plants. Our results suggest that ovule fertilization activates CmACS7 and increases ACC levels, whereas an overexpression of CmACS7 leads to an increase in ACC content in the ovule prior to fertilization, which can lead to abortion. In conclusion, the present study demonstrated that chestnut ovule abortion is caused by poor fertilization and not by nutritional competition. Optimization of the pollination and fertilization of female flowers is essential for increasing chestnut yield and reducing ovule abortion.

Plasticity and environment-specific relationships between gene expression and fitness in Saccharomyces cerevisiae

Phenotypic evolution is shaped by interactions between organisms and their environments. The environment influences how an organism's genotype determines its phenotype and how this phenotype affects its fitness. To better understand this dual role of the environment in the production and selection of phenotypic variation, we empirically determined and compared the genotype-phenotype-fitness relationship for mutant strains of the budding yeast Saccharomyces cerevisiae in four environments. Specifically, we measured how mutations in the promoter of the metabolic gene TDH3 modified its expression level and affected its growth on media with four different carbon sources. In each environment, we observed a clear relationship between TDH3 expression level and fitness, but this relationship differed among environments. Genetic variants with similar effects on TDH3 expression in different environments often had different effects on fitness and vice versa. Such environment-specific relationships between phenotype and fitness can shape the evolution of phenotypic plasticity. The set of mutants we examined also allowed us to compare the effects of mutations disrupting binding sites for key transcriptional regulators and the TATA box, which is part of the core promoter sequence. Mutations disrupting the binding sites for the transcription factors had more variable effects on expression among environments than mutations disrupting the TATA box, yet mutations with the most environmentally variable effects on fitness were located in the TATA box. This observation suggests that mutations affecting different molecular mechanisms are likely to contribute unequally to regulatory sequence evolution in changing environments.

Targeting the Epigenetic Reader ENL Inhibits Super-Enhancer-Driven Oncogenic Transcription and Synergizes with BET Inhibition to Suppress Tumor Progression

Epigenetic alterations at cis-regulatory elements (CRE) fine-tune transcriptional output. Epigenetic readers interact with CREs and can cooperate with other chromatin regulators to drive oncogene transcription. Here, we found that the YEATS domain-containing histone acetylation reader ENL (eleven-nineteen leukemia) acts as a key regulator of super-enhancers (SE), which are highly active distal CREs, across cancer types. ENL occupied the majority of SEs with substantially higher preference over typical enhancers, and the enrichment of ENL at SEs depended on its ability to bind acetylated histones. Rapid depletion of ENL by auxin-inducible degron tagging severely repressed the transcription of SE-controlled oncogenes, such as MYC, by inducing the decommissioning of their SEs, and restoring ENL protein expression largely reversed these effects. Additionally, ENL was indispensable for the rapid activation of SE-regulated immediate early genes in response to growth factor stimulation. Furthermore, ENL interacted with the histone chaperone FACT complex and was required for the deposition of FACT over CREs, which mediates nucleosome reorganization required for transcription initiation and elongation. Proper control of transcription by ENL and ENL-associated FACT was regulated by the histone reader BRD4. ENL was overexpressed in colorectal cancer and functionally contributed to colorectal cancer growth and metastasis. ENL degradation or inhibition synergized with BET inhibitors that target BRD4 in restraining colorectal cancer progression. These findings establish the essential role of epigenetic reader ENL in governing SE-driven oncogenic transcription and uncover the potential of ENL intervention to increase sensitivity to BET inhibition.

SIGNIFICANCE

ENL plays a key role in decoding epigenetic marks at highly active oncogenic super-enhancers and can be targeted in combination with BET inhibition as a promising synergistic strategy for optimizing cancer treatment.

Epigenomic insights into common human disease pathology

The epigenome-the chemical modifications and chromatin-related packaging of the genome-enables the same genetic template to be activated or repressed in different cellular settings. This multi-layered mechanism facilitates cell-type specific function by setting the local sequence and 3D interactive activity level. Gene transcription is further modulated through the interplay with transcription factors and co-regulators. The human body requires this epigenomic apparatus to be precisely installed throughout development and then adequately maintained during the lifespan. The causal role of the epigenome in human pathology, beyond imprinting disorders and specific tumour suppressor genes, was further brought into the spotlight by large-scale sequencing projects identifying that mutations in epigenomic machinery genes could be critical drivers in both cancer and developmental disorders. Abrogation of this cellular mechanism is providing new molecular insights into pathogenesis. However, deciphering the full breadth and implications of these epigenomic changes remains challenging. Knowledge is accruing regarding disease mechanisms and clinical biomarkers, through pathogenically relevant and surrogate tissue analyses, respectively. Advances include consortia generated cell-type specific reference epigenomes, high-throughput DNA methylome association studies, as well as insights into ageing-related diseases from biological 'clocks' constructed by machine learning algorithms. Also, 3rd-generation sequencing is beginning to disentangle the complexity of genetic and DNA modification haplotypes. Cell-free DNA methylation as a cancer biomarker has clear clinical utility and further potential to assess organ damage across many disorders. Finally, molecular understanding of disease aetiology brings with it the opportunity for exact therapeutic alteration of the epigenome through CRISPR-activation or inhibition.

The BTB-ZF gene Bm-mamo regulates pigmentation in silkworm caterpillars

The color pattern of insects is one of the most diverse adaptive evolutionary phenotypes. However, the molecular regulation of this color pattern is not fully understood. In this study, we found that the transcription factor Bm-mamo is responsible for black dilute (bd) allele mutations in the silkworm. Bm-mamo belongs to the BTB zinc finger family and is orthologous to mamo in Drosophila melanogaster. This gene has a conserved function in gamete production in Drosophila and silkworms and has evolved a pleiotropic function in the regulation of color patterns in caterpillars. Using RNAi and clustered regularly interspaced short palindromic repeats (CRISPR) technology, we showed that Bm-mamo is a repressor of dark melanin patterns in the larval epidermis. Using in vitro binding assays and gene expression profiling in wild-type and mutant larvae, we also showed that Bm-mamo likely regulates the expression of related pigment synthesis and cuticular protein genes in a coordinated manner to mediate its role in color pattern formation. This mechanism is consistent with the dual role of this transcription factor in regulating both the structure and shape of the cuticle and the pigments that are embedded within it. This study provides new insight into the regulation of color patterns as well as into the construction of more complex epidermal features in some insects.

Maize ZmLAZ1-3 gene negatively regulates drought tolerance in transgenic Arabidopsis

BACKGROUND

Molecular mechanisms in response to drought stress are important for the genetic improvement of maize. In our previous study, nine ZmLAZ1 members were identified in the maize genome, but the function of ZmLAZ1 was largely unknown.

RESULTS

The ZmLAZ1-3 gene was cloned from B73, and its drought-tolerant function was elucidated by expression analysis in transgenic Arabidopsis. The expression of ZmLAZ1-3 was upregulated by drought stress in different maize inbred lines. The driving activity of the ZmLAZ1-3 promoter was induced by drought stress and related to the abiotic stress-responsive elements such as MYB, MBS, and MYC. The results of subcellular localization indicated that the ZmLAZ1-3 protein localized on the plasma membrane and chloroplast. The ectopic expression of the ZmLAZ1-3 gene in Arabidopsis significantly reduced germination ratio and root length, decreased biomass, and relative water content, but increased relative electrical conductivity and malondialdehyde content under drought stress. Moreover, transcriptomics analysis showed that the differentially expressed genes between the transgenic lines and wild-type were mainly associated with response to abiotic stress and biotic stimulus, and related to pathways of hormone signal transduction, phenylpropanoid biosynthesis, mitogen-activated protein kinase signaling, and plant-pathogen interaction.

CONCLUSION

The study suggests that the ZmLAZ1-3 gene is a negative regulator in regulating drought tolerance and can be used to improve maize drought tolerance via its silencing or knockout.

From tradition to innovation: conventional and deep learning frameworks in genome annotation

Following the milestone success of the Human Genome Project, the 'Encyclopedia of DNA Elements (ENCODE)' initiative was launched in 2003 to unearth information about the numerous functional elements within the genome. This endeavor coincided with the emergence of numerous novel technologies, accompanied by the provision of vast amounts of whole-genome sequences, high-throughput data such as ChIP-Seq and RNA-Seq. Extracting biologically meaningful information from this massive dataset has become a critical aspect of many recent studies, particularly in annotating and predicting the functions of unknown genes. The core idea behind genome annotation is to identify genes and various functional elements within the genome sequence and infer their biological functions. Traditional wet-lab experimental methods still rely on extensive efforts for functional verification. However, early bioinformatics algorithms and software primarily employed shallow learning techniques; thus, the ability to characterize data and features learning was limited. With the widespread adoption of RNA-Seq technology, scientists from the biological community began to harness the potential of machine learning and deep learning approaches for gene structure prediction and functional annotation. In this context, we reviewed both conventional methods and contemporary deep learning frameworks, and highlighted novel perspectives on the challenges arising during annotation underscoring the dynamic nature of this evolving scientific landscape.

Assessing NGS-based computational methods for predicting transcriptional regulators with query gene sets

This article provides an in-depth review of computational methods for predicting transcriptional regulators with query gene sets. Identification of transcriptional regulators is of utmost importance in many biological applications, including but not limited to elucidating biological development mechanisms, identifying key disease genes, and predicting therapeutic targets. Various computational methods based on next-generation sequencing (NGS) data have been developed in the past decade, yet no systematic evaluation of NGS-based methods has been offered. We classified these methods into two categories based on shared characteristics, namely library-based and region-based methods. We further conducted benchmark studies to evaluate the accuracy, sensitivity, coverage, and usability of NGS-based methods with molecular experimental datasets. Results show that BART, ChIP-Atlas, and Lisa have relatively better performance. Besides, we point out the limitations of NGS-based methods and explore potential directions for further improvement.

Key points

An introduction to available computational methods for predicting functional TRs from a query gene set.A detailed walk-through along with practical concerns and limitations.A systematic benchmark of NGS-based methods in terms of accuracy, sensitivity, coverage, and usability, using 570 TR perturbation-derived gene sets.NGS-based methods outperform motif-based methods. Among NGS methods, those utilizing larger databases and adopting region-centric approaches demonstrate favorable performance. BART, ChIP-Atlas, and Lisa are recommended as these methods have overall better performance in evaluated scenarios.

Genomic dissection of the correlation between milk yield and various health traits using functional and evolutionary information about imputed sequence variants of 34,497 German Holstein cows

BACKGROUND

Over the last decades, it was subject of many studies to investigate the genomic connection of milk production and health traits in dairy cattle. Thereby, incorporating functional information in genomic analyses has been shown to improve the understanding of biological and molecular mechanisms shaping complex traits and the accuracies of genomic prediction, especially in small populations and across-breed settings. Still, little is known about the contribution of different functional and evolutionary genome partitioning subsets to milk production and dairy health. Thus, we performed a uni- and a bivariate analysis of milk yield (MY) and eight health traits using a set of ~34,497 German Holstein cows with 50K chip genotypes and ~17 million imputed sequence variants divided into 27 subsets depending on their functional and evolutionary annotation. In the bivariate analysis, eight trait-combinations were observed that contrasted MY with each health trait. Two genomic relationship matrices (GRM) were included, one consisting of the 50K chip variants and one consisting of each set of subset variants, to obtain subset heritabilities and genetic correlations. In addition, 50K chip heritabilities and genetic correlations were estimated applying merely the 50K GRM.

RESULTS

In general, 50K chip heritabilities were larger than the subset heritabilities. The largest heritabilities were found for MY, which was 0.4358 for the 50K and 0.2757 for the subset heritabilities. Whereas all 50K genetic correlations were negative, subset genetic correlations were both, positive and negative (ranging from -0.9324 between MY and mastitis to 0.6662 between MY and digital dermatitis). The subsets containing variants which were annotated as noncoding related, splice sites, untranslated regions, metabolic quantitative trait loci, and young variants ranked highest in terms of their contribution to the traits` genetic variance. We were able to show that linkage disequilibrium between subset variants and adjacent variants did not cause these subsets` high effect.

CONCLUSION

Our results confirm the connection of milk production and health traits in dairy cattle via the animals` metabolic state. In addition, they highlight the potential of including functional information in genomic analyses, which helps to dissect the extent and direction of the observed traits` connection in more detail.

Gene expression variation underlying tissue-specific responses to copper stress in Drosophila melanogaster

Copper is one of a handful of biologically necessary heavy metals that is also a common environmental pollutant. Under normal conditions, copper ions are required for many key physiological processes. However, in excess, copper results in cell and tissue damage ranging in severity from temporary injury to permanent neurological damage. Because of its biological relevance, and because many conserved copper-responsive genes respond to nonessential heavy metal pollutants, copper resistance in Drosophila melanogaster is a useful model system with which to investigate the genetic control of the heavy metal stress response. Because heavy metal toxicity has the potential to differently impact specific tissues, we genetically characterized the control of the gene expression response to copper stress in a tissue-specific manner in this study. We assessed the copper stress response in head and gut tissue of 96 inbred strains from the Drosophila Synthetic Population Resource using a combination of differential expression analysis and expression quantitative trait locus mapping. Differential expression analysis revealed clear patterns of tissue-specific expression. Tissue and treatment specific responses to copper stress were also detected using expression quantitative trait locus mapping. Expression quantitative trait locus associated with MtnA, Mdr49, Mdr50, and Sod3 exhibited both genotype-by-tissue and genotype-by-treatment effects on gene expression under copper stress, illuminating tissue- and treatment-specific patterns of gene expression control. Together, our data build a nuanced description of the roles and interactions between allelic and expression variation in copper-responsive genes, provide valuable insight into the genomic architecture of susceptibility to metal toxicity, and highlight candidate genes for future functional characterization.

Genome-wide identification of R2R3-MYB family genes and gene response to stress in ginger

Ginger (Zingiber officinale Roscoe) is an important plant used worldwide for medicine and food. The R2R3-MYB transcription factor (TF) family has essential roles in plant growth, development, and stresses resistance, and the number of genes in the family varies greatly among different types of plants. However, genome-wide discovery of ZoMYBs and gene responses to stresses have not been reported in ginger. Therefore, genome-wide analysis of R2R3-MYB genes in ginger was conducted in this study. Protein phylogenetic relations and conserved motifs and chromosome localization and duplication, structure, and cis-regulatory elements were analyzed. In addition, the expression patterns of selected genes were analyzed under two different stresses. A total of 299 candidate ZoMYB genes were discovered in ginger. Based on groupings of R2R3-MYB genes in the model plant Arabidopsis thaliana (L.) Heynh., ZoMYBs were divided into eight groups. Genes were distributed across 22 chromosomes at uneven densities. In gene duplication analysis, 120 segmental duplications were identified in the ginger genome. Gene expression patterns of 10 ZoMYBs in leaves of ginger under abscisic acid (ABA) and low-temperature stress treatments were different. The results will help to determine the exact roles of ZoMYBs in anti-stress responses in ginger.

Optimized barley phytase gene expression by focused FIND-IT screening for mutations in cis-acting regulatory elements

Introduction

Induced modification of plant gene expression is of both fundamental and applied importance. Cis-acting regulatory elements (CREs) are major determinants of the spatiotemporal strength of gene expression. Yet, there are few examples where induced genetic variation in predetermined CREs has been exploited to improve or investigate crop plants.

Methods

The digital PCR based FIND-IT technology was applied to discover barley mutants with CRE variants in the promoter of the nutritional important barley grain phytase (PAPhy_a) gene.

Results and discussion

Mutants with higher or lower gene expression and ultimately higher or lower mature grain phytase activity (MGPA), respectively, were discovered. Field trials and inositol phosphate profiling during germination showed that PAPhy_a does not influence agronomic performance under the trial conditions but it does shorten the lag time of phosphate mobilization during germination. Higher endogenous MGPA is an improvement of grain quality for feed use as it improves the phosphate bioavailability for monogastric animals. Moreover, as the targeted CRE motifs of the PAPhy_a promoter are shared with a range of seed expressed genes like key cereal and legume storage genes, the current results demonstrates a concept for modulating individual gene expression levels of a range of seed genes.

Extreme restructuring of cis-regulatory regions controlling a deeply conserved plant stem cell regulator

A striking paradox is that genes with conserved protein sequence, function and expression pattern over deep time often exhibit extremely divergent cis-regulatory sequences. It remains unclear how such drastic cis-regulatory evolution across species allows preservation of gene function, and to what extent these differences influence how cis-regulatory variation arising within species impacts phenotypic change. Here, we investigated these questions using a plant stem cell regulator conserved in expression pattern and function over ~125 million years. Using in-vivo genome editing in two distantly related models, Arabidopsis thaliana (Arabidopsis) and Solanum lycopersicum (tomato), we generated over 70 deletion alleles in the upstream and downstream regions of the stem cell repressor gene CLAVATA3 (CLV3) and compared their individual and combined effects on a shared phenotype, the number of carpels that make fruits. We found that sequences upstream of tomato CLV3 are highly sensitive to even small perturbations compared to its downstream region. In contrast, Arabidopsis CLV3 function is tolerant to severe disruptions both upstream and downstream of the coding sequence. Combining upstream and downstream deletions also revealed a different regulatory outcome. Whereas phenotypic enhancement from adding downstream mutations was predominantly weak and additive in tomato, mutating both regions of Arabidopsis CLV3 caused substantial and synergistic effects, demonstrating distinct distribution and redundancy of functional cis-regulatory sequences. Our results demonstrate remarkable malleability in cis-regulatory structural organization of a deeply conserved plant stem cell regulator and suggest that major reconfiguration of cis-regulatory sequence space is a common yet cryptic evolutionary force altering genotype-to-phenotype relationships from regulatory variation in conserved genes. Finally, our findings underscore the need for lineage-specific dissection of the spatial architecture of cis-regulation to effectively engineer trait variation from conserved productivity genes in crops.

Genome-Wide Comprehensive Identification and In Silico Characterization of Lectin Receptor-Like Kinase Gene Family in Barley (Hordeum vulgare L.)

Lectin receptor-like kinases (LecRLKs) are a significant subgroup of the receptor-like kinases (RLKs) protein family. They play crucial roles in plant growth, development, immune responses, signal transduction, and stress tolerance. However, the genome-wide identification and characterization of LecRLK genes and their regulatory elements have not been explored in a major cereal crop, barley (Hordeum vulgare L.). Therefore, in this study, integrated bioinformatics tools were used to identify and characterize the LecRLK gene family in barley. Based on the phylogenetic tree and domain organization, a total of 113 LecRLK genes were identified in the barley genome (referred to as HvlecRLK) corresponding to the LecRLK genes of Arabidopsis thaliana. These putative HvlecRLK genes were classified into three groups: 62 G-type LecRLKs, 1 C-type LecRLK, and 50 L-type LecRLKs. They were unevenly distributed across eight chromosomes, including one unknown chromosome, and were predominantly located in the plasma membrane (G-type HvlecRLK (96.8%), C-type HvlecRLK (100%), and L-type HvlecRLK (98%)). An analysis of motif composition and exon-intron configuration revealed remarkable homogeneity with the members of AtlecRLK. Notably, most of the HvlecRLKs (27 G-type, 43 L-type) have no intron, suggesting their rapid functionality. The Ka/Ks and syntenic analysis demonstrated that HvlecRLK gene pairs evolved through purifying selection and gene duplication was the major factor for the expansion of the HvlecRLK gene family. Exploration of gene ontology (GO) enrichment indicated that the identified HvlecRLK genes are associated with various cellular processes, metabolic pathways, defense mechanisms, kinase activity, catalytic activity, ion binding, and other essential pathways. The regulatory network analysis identified 29 transcription factor families (TFFs), with seven major TFFs including bZIP, C2H2, ERF, MIKC_MADS, MYB, NAC, and WRKY participating in the regulation of HvlecRLK gene functions. Most notably, eight TFFs were found to be linked to the promoter region of both L-type HvleckRLK64 and HvleckRLK86. The promoter cis-acting regulatory element (CARE) analysis of barley identified a total of 75 CARE motifs responsive to light responsiveness (LR), tissue-specific (TS), hormone responsiveness (HR), and stress responsiveness (SR). The maximum number of CAREs was identified in HvleckRLK11 (25 for LR), HvleckRLK69 (17 for TS), and HvleckRLK80 (12 for HR). Additionally, HvleckRLK14, HvleckRLK16, HvleckRLK33, HvleckRLK50, HvleckRLK52, HvleckRLK56, and HvleckRLK110 were predicted to exhibit higher responses in stress conditions. In addition, 46 putative miRNAs were predicted to target 81 HvlecRLK genes and HvlecRLK13 was the most targeted gene by 8 different miRNAs. Protein-protein interaction analysis demonstrated higher functional similarities of 63 HvlecRLKs with 7 Arabidopsis STRING proteins. Our overall findings provide valuable information on the LecRLK gene family which might pave the way to advanced research on the functional mechanism of the candidate genes as well as to develop new barley cultivars in breeding programs.

Experimental evidence for cancer resistance in a bat species

Mammals exhibit different rates of cancer, with long-lived species generally showing greater resistance. Although bats have been suggested to be resistant to cancer due to their longevity, this has yet to be systematically examined. Here, we investigate cancer resistance across seven bat species by activating oncogenic genes in their primary cells. Both in vitro and in vivo experiments suggest that Myotis pilosus (MPI) is particularly resistant to cancer. The transcriptomic and functional analyses reveal that the downregulation of three genes (HIF1A, COPS5, and RPS3) largely contributes to cancer resistance in MPI. Further, we identify the loss of a potential enhancer containing the HIF1A binding site upstream of COPS5 in MPI, resulting in the downregulation of COPS5. These findings not only provide direct experimental evidence for cancer resistance in a bat species but also offer insights into the natural mechanisms of cancer resistance in mammals.

Cell-type-specific cis-regulatory divergence in gene expression and chromatin accessibility revealed by human-chimpanzee hybrid cells

Although gene expression divergence has long been postulated to be the primary driver of human evolution, identifying the genes and genetic variants underlying uniquely human traits has proven to be quite challenging. Theory suggests that cell-type-specific cis-regulatory variants may fuel evolutionary adaptation due to the specificity of their effects. These variants can precisely tune the expression of a single gene in a single cell-type, avoiding the potentially deleterious consequences of trans-acting changes and non-cell type-specific changes that can impact many genes and cell types, respectively. It has recently become possible to quantify human-specific cis-acting regulatory divergence by measuring allele-specific expression in human-chimpanzee hybrid cells-the product of fusing induced pluripotent stem (iPS) cells of each species in vitro. However, these cis-regulatory changes have only been explored in a limited number of cell types. Here, we quantify human-chimpanzee cis-regulatory divergence in gene expression and chromatin accessibility across six cell types, enabling the identification of highly cell-type-specific cis-regulatory changes. We find that cell-type-specific genes and regulatory elements evolve faster than those shared across cell types, suggesting an important role for genes with cell-type-specific expression in human evolution. Furthermore, we identify several instances of lineage-specific natural selection that may have played key roles in specific cell types, such as coordinated changes in the cis-regulation of dozens of genes involved in neuronal firing in motor neurons. Finally, using novel metrics and a machine learning model, we identify genetic variants that likely alter chromatin accessibility and transcription factor binding, leading to neuron-specific changes in the expression of the neurodevelopmentally important genes FABP7 and GAD1. Overall, our results demonstrate that integrative analysis of cis-regulatory divergence in chromatin accessibility and gene expression across cell types is a promising approach to identify the specific genes and genetic variants that make us human.

Identification and expression analysis of ATP-binding cassette (ABC) transporters revealed its role in regulating stress response in pear (Pyrus bretchneideri)

BACKGROUND

ATP-binding cassette (ABC) transporter proteins constitute a plant gene superfamily crucial for growth, development, and responses to environmental stresses. Despite their identification in various plants like maize, rice, and Arabidopsis, little is known about the information on ABC transporters in pear. To investigate the functions of ABC transporters in pear development and abiotic stress response, we conducted an extensive analysis of ABC gene family in the pear genome.

RESULTS

In this study, 177 ABC transporter genes were successfully identified in the pear genome, classified into seven subfamilies: 8 ABCAs, 40 ABCBs, 24 ABCCs, 8 ABCDs, 9 ABCEs, 8 ABCFs, and 80 ABCGs. Ten motifs were common among all ABC transporter proteins, while distinct motif structures were observed for each subfamily. Distribution analysis revealed 85 PbrABC transporter genes across 17 chromosomes, driven primarily by WGD and dispersed duplication. Cis-regulatory element analysis of PbrABC promoters indicated associations with phytohormones and stress responses. Tissue-specific expression profiles demonstrated varied expression levels across tissues, suggesting diverse functions in development. Furthermore, several PbrABC genes responded to abiotic stresses, with 82 genes sensitive to salt stress, including 40 upregulated and 23 downregulated genes. Additionally, 91 genes were responsive to drought stress, with 22 upregulated and 36 downregulated genes. These findings highlight the pivotal role of PbrABC genes in abiotic stress responses.

CONCLUSION

This study provides evolutionary insights into PbrABC transporter genes, establishing a foundation for future research on their functions in pear. The identified motifs, distribution patterns, and stress-responsive expressions contribute to understanding the regulatory mechanisms of ABC transporters in pear. The observed tissue-specific expression profiles suggest diverse roles in developmental processes. Notably, the significant responses to salt and drought stress emphasize the importance of PbrABC genes in mediating adaptive responses. Overall, our study advances the understanding of PbrABC transporter genes in pear, opening avenues for further investigations in plant molecular biology and stress physiology.

Evolution of chemosensory tissues and cells across ecologically diverse Drosophilids

Chemosensory tissues exhibit significant between-species variability, yet the evolution of gene expression and cell types underlying this diversity remain poorly understood. To address these questions, we conducted transcriptomic analyses of five chemosensory tissues from six Drosophila species and integrated the findings with single-cell datasets. While stabilizing selection predominantly shapes chemosensory transcriptomes, thousands of genes in each tissue have evolved expression differences. Genes that have changed expression in one tissue have often changed in multiple other tissues but at different past epochs and are more likely to be cell type-specific than unchanged genes. Notably, chemosensory-related genes have undergone widespread expression changes, with numerous species-specific gains/losses including novel chemoreceptors expression patterns. Sex differences are also pervasive, including a D. melanogaster-specific excess of male-biased expression in sensory and muscle cells in its forelegs. Together, our analyses provide new insights for understanding evolutionary changes in chemosensory tissues at both global and individual gene levels.

Natural variation in a K+ -preferring HKT transporter contributes to wheat shoot K+ accumulation and salt tolerance

Soil salinity can adversely affect crop growth and yield, and an improved understanding of the genetic factors that confer salt tolerance could inform breeding strategies to engineer salt-tolerant crops and improve productivity. Here, a group of K+ -preferring HKT transporters, TaHKT8, TaHKT9 and TaHKT10, were identified and negatively regulate the wheat shoot K+ accumulation and salt tolerance. A genome-wide association study (GWAS) and candidate gene association analysis further revealed that TaHKT9-B substantially underlies the natural variation of wheat shoot K+ accumulation under saline soil conditions. Specifically, an auxin responsive element (ARE) within an 8-bp insertion in the promoter of TaHKT9-B is strongly associated with shoot K+ content among wheat accessions. This ARE can be directly bound by TaARF4 for transcriptional activation of TaHKT9-B, which subsequently attenuates shoot K+ accumulation and salt tolerance. Moreover, the tae-miR390/TaTAS3/TaARF4 pathway was identified to regulate the salt-induced root development and salt tolerance in wheat. Taken together, our study describes the genetic basis and accompanying mechanism driving phenotypic variation in wheat shoot K+ accumulation and salt tolerance. The identified tae-miR390/TaTAS3/TaARF4/TaHKT9-B module is an important regulator in wheat subjected to salt stress, which provides the potentially important genetic resources for breeders to improve wheat salt tolerance.

Comprehensive analysis of annexin gene family and its expression in response to branching architecture and salt stress in crape myrtle

BACKGROUND

Annexin (ANN) is calcium (Ca2+)-dependent and phospholipid binding protein family, which is involved in plant growth and development and response to various stresses. However, little known about ANN genes were identified from crape myrtle, an ornamental horticultural plant widely cultivated in the world.

RESULTS

Here, 9 LiANN genes were identified from Lagerstroemia indica, and their characterizations and functions were investigated in L. indica for the first time. The LiANN genes were divided into 2 subfamilies. The gene structure, chromosomal location, and collinearity relationship were also explored. In addition, the GO annotation analysis of these LiANNs indicated that they are enriched in molecular functions, cellular components, and biological processes. Moreover, transcription factors (TFs) prediction analysis revealed that bHLH, MYB, NAC, and other TFs can interact with the LiANN promoters. Interestingly, the LiANN2/4/6-9 were demonstrated to play critical roles in the branching architecture of crape myrtle. Furthermore, the LiANN2/6/8/9 were differentially expressed under salt treatment, and a series of TFs regulating LiANN2/6/8/9 expression were predicted to play essential roles in salt resistance.

CONCLUSIONS

These results shed light on profile and function of the LiANN gene family, and lay a foundation for further studies of the LiANN genes.

Expression analysis of genes encoding extracellular leucine-rich repeat proteins in Arabidopsis thaliana

Leucine-rich repeat (LRR)-containing proteins have been identified in diverse species, including plants. The diverse intracellular and extracellular LRR variants are responsible for numerous biological processes. We analyzed the expression patterns of Arabidopsis thaliana extracellular LRR (AtExLRR) genes, 10 receptor-like proteins, and 4 additional genes expressing the LRR-containing protein by a promoter: β-glucuronidase (GUS) study. According to in silico expression studies, several AtExLRR genes were expressed in a tissue- or stage-specific and abiotic/hormone stress-responsive manner, indicating their potential participation in specific biological processes. Based on the promoter: GUS assay, AtExLRRs were expressed in different cells and organs. A quantitative real-time PCR investigation revealed that the expressions of AtExLRR3 and AtExLRR9 were distinct under various abiotic stress conditions. This study investigated the potential roles of extracellular LRR proteins in plant growth, development, and response to various abiotic stresses.

Quantification and modeling of turnover dynamics of de novo transcripts in Drosophila melanogaster

Most of the transcribed eukaryotic genomes are composed of non-coding transcripts. Among these transcripts, some are newly transcribed when compared to outgroups and are referred to as de novo transcripts. De novo transcripts have been shown to play a major role in genomic innovations. However, little is known about the rates at which de novo transcripts are gained and lost in individuals of the same species. Here, we address this gap and estimate the de novo transcript turnover rate with an evolutionary model. We use DNA long reads and RNA short reads from seven geographically remote samples of inbred individuals of Drosophila melanogaster to detect de novo transcripts that are gained on a short evolutionary time scale. Overall, each sampled individual contains around 2500 unspliced de novo transcripts, with most of them being sample specific. We estimate that around 0.15 transcripts are gained per year, and that each gained transcript is lost at a rate around 5× 10-5 per year. This high turnover of transcripts suggests frequent exploration of new genomic sequences within species. These rate estimates are essential to comprehend the process and timescale of de novo gene birth.

Deciphering the Role of Rapidly Evolving Conserved Elements in Primate Brain Development and Exploring Their Potential Involvement in Alzheimer's Disease

Although previous studies have identified human-specific accelerated regions as playing a key role in the recent evolution of the human brain, the characteristics and cellular functions of rapidly evolving conserved elements (RECEs) in ancestral primate lineages remain largely unexplored. Here, based on large-scale primate genome assemblies, we identify 888 RECEs that have been highly conserved in primates that exhibit significantly accelerated substitution rates in the ancestor of the Simiiformes. This primate lineage exhibits remarkable morphological innovations, including an expanded brain mass. Integrative multiomic analyses reveal that RECEs harbor sequences with potential cis-regulatory functions that are activated in the adult human brain. Importantly, genes linked to RECEs exhibit pronounced expression trajectories in the adult brain relative to the fetal stage. Furthermore, we observed an increase in the chromatin accessibility of RECEs in oligodendrocytes from individuals with Alzheimer's disease (AD) compared to that of a control group, indicating that these RECEs may contribute to brain aging and AD. Our findings serve to expand our knowledge of the genetic underpinnings of brain function during primate evolution.

Genome-wide analysis of the ABC gene family in almond and functional predictions during flower development, freezing stress, and salt stress

ABC (ATP-binding cassette) transporter proteins are one of the most extensive protein families known to date and are ubiquitously found in animals, plants, and microorganisms. ABCs have a variety of functions, such as plant tissue development regulation, hormone transport, and biotic and abiotic stress resistance. However, the gene characterization and function of the ABC gene family in almond (Prunus dulcis) have not been thoroughly studied. In this study, we identified 117 PdABC genes using the whole genome of 'Wanfeng' almond obtained by sequencing and explored their protein characterization. The PdABC family members were classified into eight subfamilies. The members of the same subfamily had conserved motifs but poorly conserved numbers of exons and introns and were unevenly distributed among the eight subfamilies and on the eight chromosomes. Expression patterns showed that PdABC family members were significantly differentially expressed during almond development, dormant freezing stress, and salt stress. We found that PdABC59 and PdABC77 had extremely high expression levels in pollen. PdABC63 and PdABC64 had high expression levels during almond petal development and multiple stages of flower development. PdABC98 was highly expressed in annual dormant branches after six temperature-freezing stress treatments. PdABC29, PdABC69, and PdABC98 were highly expressed under different concentrations of salt stress. This study preliminarily investigated the expression characteristics of ABC genes in different tissues of almond during flower development, freezing stress and salt stress, and the results will provide a reference for further in-depth research and breeding of almond in the future.

Hold out the genome: a roadmap to solving the cis-regulatory code

Gene expression is regulated by transcription factors that work together to read cis-regulatory DNA sequences. The 'cis-regulatory code' - how cells interpret DNA sequences to determine when, where and how much genes should be expressed - has proven to be exceedingly complex. Recently, advances in the scale and resolution of functional genomics assays and machine learning have enabled substantial progress towards deciphering this code. However, the cis-regulatory code will probably never be solved if models are trained only on genomic sequences; regions of homology can easily lead to overestimation of predictive performance, and our genome is too short and has insufficient sequence diversity to learn all relevant parameters. Fortunately, randomly synthesized DNA sequences enable testing a far larger sequence space than exists in our genomes, and designed DNA sequences enable targeted queries to maximally improve the models. As the same biochemical principles are used to interpret DNA regardless of its source, models trained on these synthetic data can predict genomic activity, often better than genome-trained models. Here we provide an outlook on the field, and propose a roadmap towards solving the cis-regulatory code by a combination of machine learning and massively parallel assays using synthetic DNA.

Exploring the role of FBXL fbxl gene family in Soybean: Implications for plant height and seed size regulation

F-box proteins constitute a significant family in eukaryotes and, as a component of the Skp1p-cullin-F-box complex, are considered critical for cellular protein degradation and other biological processes in plants. Despite their importance, the functions of F-box proteins, particularly those with C-terminal leucine-rich repeat (LRR) domains, remain largely unknown in plants. Therefore, the present study conducted genome-wide identification and in silico characterization of F-BOX proteins with C-terminal LRR domains in soybean (Glycine max L.) (GmFBXLs). A total of 45 GmFBXLs were identified. The phylogenetic analysis showed that GmFBXLs could be subdivided into ten subgroups and exhibited a close relationship with those from Arabidopsis thaliana, Cicer aretineum, and Medicago trunculata. It was observed that most cis-regulatory elements in the promoter regions of GmFBXLs are involved in hormone signalling, stress responses, and developmental stages. In silico transcriptome data illustrated diverse expression patterns of the identified GmFBXLs across various tissues, such as shoot apical meristem, flower, green pods, leaves, nodules, and roots. Overexpressing (OE) GmFBXL12 in Tianlong No.1 cultivar resulted in a significant difference in seed size, number of pods, and number of seeds per plant, indicated a potential increase in yield compared to wild type. This study offers valuable perspectives into the role of FBXLs in soybean, serving as a foundation for future research. Additionally, the identified OE lines represent valuable genetic resources for enhancing seed-related traits in soybean.

The C2H2 zinc finger transcription factor CF2-II regulates multi-insecticide resistance-related gut-predominant ABC transporters in Aphis gossypii Glover

Clarifying the molecular mechanisms of cotton aphid resistance to various insecticides is crucial for the long-term safe application of insecticides in chemical control. ATP-binding cassette (ABC) transporters mediate the membrane transport of various substrates (including exogenous substances). Experiments confirmed that ABCB5, ABCF2, and MRP12 contributed to high levels of resistance to spirotetramat, cyantraniliprole, thiamethoxam or imidacloprid. Binding sites of the C2H2 zinc finger transcription factor CF2-II was predicted to be located in the promoters of ABCB5, ABCF2, and MRP12. The expression levels of ABCB5, ABCF2, and MRP12 were significantly upregulated after silencing CF2-II. The results of dual-luciferase reporter assays demonstrated a negative regulatory relationship between CF2-II and ABC transporter promoters. Furthermore, yeast one-hybrid (Y1H) and electrophoresis mobility shift assays (EMSAs) revealed that CF2-II inhibited the expression of ABC transporter genes through interaction with binding sites [ABCF2.p (-1149/-1140) or MRP12.p (-1189/-1181)]. The above results indicated that ABCB5, ABCF2, and MRP12 were negatively regulated by the transcription factor CF2-II, which will help us further understand the mechanism of transcriptional adaption of multi-insecticides resistant related ABC transporters in response to xenobiotics.

Selected humanization of yeast U1 snRNP leads to global suppression of pre-mRNA splicing and mitochondrial dysfunction in the budding yeast

The recognition of 5' splice site (5' ss) is one of the earliest steps of pre-mRNA splicing. To better understand the mechanism and regulation of 5' ss recognition, we selectively humanized components of the yeast U1 snRNP to reveal the function of these components in 5' ss recognition and splicing. We targeted U1C and Luc7, two proteins that interact with and stabilize the yeast U1 (yU1) snRNA and the 5' ss RNA duplex. We replaced the Zinc-Finger (ZnF) domain of yU1C with its human counterpart, which resulted in cold-sensitive growth phenotype and moderate splicing defects. Next, we added an auxin-inducible degron to yLuc7 protein and found that Luc7-depleted yU1 snRNP resulted in the concomitant loss of PRP40 and Snu71 (two other essential yeast U1 snRNP proteins), and further biochemical analyses suggest a model of how these three proteins interact with each other in the U1 snRNP. The loss of these proteins resulted in a significant growth retardation accompanied by a global suppression of pre-mRNA splicing. The splicing suppression led to mitochondrial dysfunction as revealed by a release of Fe 2+ into the growth medium and an induction of mitochondrial reactive oxygen species. Together, these observations indicate that the human U1C ZnF can substitute that of yeast, Luc7 is essential for the incorporation of the Luc7-Prp40-Snu71 trimer into yeast U1 snRNP, and splicing plays a major role in the regulation of mitochondria function in yeast.

Evolutionary analysis of conserved non-coding elements subsequent to whole-genome duplication in opium poppy

Whole-genome duplication (WGD) leads to the duplication of both coding and non-coding sequences within an organism's genome, providing an abundant supply of genetic material that can drive evolution, ultimately contributing to plant adaptation and speciation. Although non-coding sequences contain numerous regulatory elements, they have been understudied compared to coding sequences. In order to address this gap, we explored the evolutionary patterns of regulatory sequences, coding sequences and transcriptomes using conserved non-coding elements (CNEs) as regulatory element proxies following the recent WGD event in opium poppy (Papaver somniferum). Our results showed similar evolutionary patterns in subgenomes of regulatory and coding sequences. Specifically, the biased or unbiased retention of coding sequences reflected the same pattern as retention levels in regulatory sequences. Further, the divergence of gene expression patterns mediated by regulatory element variations occurred at a more rapid pace than that of gene coding sequences. However, gene losses were purportedly dependent on relaxed selection pressure in coding sequences. Specifically, the rapid evolution of tissue-specific benzylisoquinoline alkaloid production in P. somniferum was associated with regulatory element changes. The origin of a novel stem-specific ACR, which utilized ancestral cis-elements as templates, is likely to be linked to the evolutionary trajectory behind the transition of the PSMT1-CYP719A21 cluster from high levels of expression solely in P. rhoeas root tissue to its elevated expression in P. somniferum stem tissue. Our findings demonstrate that rapid regulatory element evolution can contribute to the emergence of new phenotypes and provide valuable insights into the high evolvability of regulatory elements.

DNA hydroxymethylation differences underlie phenotypic divergence of somatic growth in Nile tilapia reared in common garden

Phenotypic plasticity of metabolism and growth are essential for adaptation to new environmental conditions, such as those experienced during domestication. Epigenetic regulation plays a key role in this process but the underlying mechanisms are poorly understood, especially in the case of hydroxymethylation. Using reduced representation 5-hydroxymethylcytosine profiling, we compared the liver hydroxymethylomes in full-sib Nile tilapia with distinct growth rates (3.8-fold difference) and demonstrated that DNA hydroxymethylation is strongly associated with phenotypic divergence of somatic growth during the early stages of domestication. The 2677 differentially hydroxymethylated cytosines between fast- and slow-growing fish were enriched within gene bodies (79%), indicating a pertinent role in transcriptional regulation. Moreover, they were found in genes involved in biological processes related to skeletal system and muscle structure development, and there was a positive association between somatic growth and 5hmC levels in genes coding for growth factors, kinases and receptors linked to myogenesis. Single nucleotide polymorphism analysis revealed no genetic differentiation between fast- and slow-growing fish. In addition to unveiling a new link between DNA hydroxymethylation and epigenetic regulation of growth in fish during the initial stages of domestication, this study suggests that epimarkers may be applied in selective breeding programmes for superior phenotypes.