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

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.

Genome-Wide Identification and Characterization of Major RNAi Genes Highlighting Their Associated Factors in Cowpea (Vigna unguiculata (L.) Walp.)

In different regions of the world, cowpea (Vigna unguiculata (L.) Walp.) is an important vegetable and an excellent source of protein. It lessens the malnutrition of the underprivileged in developing nations and has some positive effects on health, such as a reduction in the prevalence of cancer and cardiovascular disease. However, occasionally, certain biotic and abiotic stresses caused a sharp fall in cowpea yield. Major RNA interference (RNAi) genes like Dicer-like (DCL), Argonaute (AGO), and RNA-dependent RNA polymerase (RDR) are essential for the synthesis of their associated factors like domain, small RNAs (sRNAs), transcription factors, micro-RNAs, and cis-acting factors that shield plants from biotic and abiotic stresses. In this study, applying BLASTP search and phylogenetic tree analysis with reference to the Arabidopsis RNAi (AtRNAi) genes, we discovered 28 VuRNAi genes, including 7 VuDCL, 14 VuAGO, and 7 VuRDR genes in cowpea. We looked at the domains, motifs, gene structures, chromosomal locations, subcellular locations, gene ontology (GO) terms, and regulatory factors (transcription factors, micro-RNAs, and cis-acting elements (CAEs)) to characterize the VuRNAi genes and proteins in cowpea in response to stresses. Predicted VuDCL1, VuDCL2(a, b), VuAGO7, VuAGO10, and VuRDR6 genes might have an impact on cowpea growth, development of the vegetative and flowering stages, and antiviral defense. The VuRNAi gene regulatory features miR395 and miR396 might contribute to grain quality improvement, immunity boosting, and pathogen infection resistance under salinity and drought conditions. Predicted CAEs from the VuRNAi genes might play a role in plant growth and development, improving grain quality and production and protecting plants from biotic and abiotic stresses. Therefore, our study provides crucial information about the functional roles of VuRNAi genes and their associated components, which would aid in the development of future cowpeas that are more resilient to biotic and abiotic stress. The manuscript is available as a preprint at this link: doi:10.1101/2023.02.15.528631v1.

Multi-omics for studying and understanding polar life

Polar ecosystems are experiencing amongst the most rapid rates of regional warming on Earth. Here, we discuss 'omics' approaches to investigate polar biodiversity, including the current state of the art, future perspectives and recommendations. We propose a community road map to generate and more fully exploit multi-omics data from polar organisms. These data are needed for the comprehensive evaluation of polar biodiversity and to reveal how life evolved and adapted to permanently cold environments with extreme seasonality. We argue that concerted action is required to mitigate the impact of warming on polar ecosystems via conservation efforts, to sustainably manage these unique habitats and their ecosystem services, and for the sustainable bioprospecting of novel genes and compounds for societal gain.

Species differences in hormonally mediated gene expression underlie the evolutionary loss of sexually dimorphic coloration in Sceloporus lizards

Phenotypic sexual dimorphism often involves the hormonal regulation of sex-biased expression for underlying genes. However, it is generally unknown whether the evolution of hormonally mediated sexual dimorphism occurs through upstream changes in tissue sensitivity to hormone signals, downstream changes in responsiveness of target genes, or both. Here, we use comparative transcriptomics to explore these possibilities in 2 species of Sceloporus lizards exhibiting different patterns of sexual dichromatism. Sexually dimorphic S. undulatus develops blue and black ventral coloration in response to testosterone, while sexually monomorphic S. virgatus does not, despite exhibiting similar sex differences in circulating testosterone levels. We administered testosterone implants to juveniles of each species and used RNAseq to quantify gene expression in ventral skin. Transcriptome-wide responses to testosterone were stronger in S. undulatus than in S. virgatus, suggesting species differences in tissue sensitivity to this hormone signal. Species differences in the expression of genes for androgen metabolism and sex hormone-binding globulin were consistent with this idea, but expression of the androgen receptor gene was higher in S. virgatus, complicating this interpretation. Downstream of androgen signaling, we found clear species differences in hormonal responsiveness of genes related to melanin synthesis, which were upregulated by testosterone in S. undulatus, but not in S. virgatus. Collectively, our results indicate that hormonal regulation of melanin synthesis pathways contributes to the development of sexual dimorphism in S. undulatus, and that changes in the hormonal responsiveness of these genes in S. virgatus contribute to the evolutionary loss of ventral coloration.

Genome-wide identification, classification, and characterization of lectin gene superfamily in sweet orange (Citrus sinensis L.)

Lectins are sugar-binding proteins found abundantly in plants. Lectin superfamily members have diverse roles, including plant growth, development, cellular processes, stress responses, and defense against microbes. However, the genome-wide identification and functional analysis of lectin genes in sweet orange (Citrus sinensis L.) remain unexplored. Therefore, we used integrated bioinformatics approaches (IBA) for in-depth genome-wide identification, characterization, and regulatory factor analysis of sweet orange lectin genes. Through genome-wide comparative analysis, we identified a total of 141 lectin genes distributed across 10 distinct gene families such as 68 CsB-Lectin, 13 CsLysin Motif (LysM), 4 CsChitin-Bind1, 1 CsLec-C, 3 CsGal-B, 1 CsCalreticulin, 3 CsJacalin, 13 CsPhloem, 11 CsGal-Lec, and 24 CsLectinlegB.This classification relied on characteristic domain and phylogenetic analysis, showing significant homology with Arabidopsis thaliana's lectin gene families. A thorough analysis unveiled common similarities within specific groups and notable variations across different protein groups. Gene Ontology (GO) enrichment analysis highlighted the predicted genes' roles in diverse cellular components, metabolic processes, and stress-related regulation. Additionally, network analysis of lectin genes with transcription factors (TFs) identified pivotal regulators like ERF, MYB, NAC, WRKY, bHLH, bZIP, and TCP. The cis-acting regulatory elements (CAREs) found in sweet orange lectin genes showed their roles in crucial pathways, including light-responsive (LR), stress-responsive (SR), hormone-responsive (HR), and more. These findings will aid in the in-depth molecular examination of these potential genes and their regulatory elements, contributing to targeted enhancements of sweet orange species in breeding programs.

Identification and characterization of Glycolate oxidase gene family in garden lettuce (Lactuca sativa cv. 'Salinas') and its response under various biotic, abiotic, and developmental stresses

Glycolate oxidase (GLO) is an FMN-containing enzyme localized in peroxisomes and performs in various molecular and biochemical mechanisms. It is a key player in plant glycolate and glyoxylate accumulation pathways. The role of GLO in disease and stress resistance is well-documented in various plant species. Although studies have been conducted regarding the role of GLO genes from spinach on a microbial level, the direct response of GLO genes to various stresses in short-season and leafy plants like lettuce has not been published yet. The genome of Lactuca sativa cultivar 'Salinas' (v8) was used to identify GLO gene members in lettuce by performing various computational analysis. Dual synteny, protein-protein interactions, and targeted miRNA analyses were conducted to understand the function of GLO genes. The identified GLO genes showed further clustering into two groups i.e., glycolate oxidase (GOX) and hydroxyacid oxidase (HAOX). Genes were observed to be distributed unevenly on three chromosomes, and syntenic analysis revealed that segmental duplication was prevalent. Thus, it might be the main reason for GLO gene diversity in lettuce. Almost all LsGLO genes showed syntenic blocks in respective plant genomes under study. Protein-protein interactions of LsGLO genes revealed various functional enrichments, mainly photorespiration, and lactate oxidation, and among biological processes oxidative photosynthetic carbon pathway was highly significant. Results of in-depth analyses disclosed the interaction of GLO genes with other members of the glycolate pathway and the activity of GLO genes in various organs and developmental stages in lettuce. The extensive genome evaluation of GLO gene family in garden lettuce is believed to be a reference for cloning and studying functional analyses of GLO genes and characterizing other members of glycolate/glyoxylate biosynthesis pathway in various plant species.

Improving yield-related traits by editing the promoter of the heading date gene Ehd1 in rice

KEY MESSAGE

We developed an efficient promoter editing method to create different weak Ehd1 alleles in elite japonica rice variety ZJ8 with slightly delayed heading and improved yield for use in breeding. Heading date is an important agronomic trait of rice (Oryza sativa) that determines the planting areas and cultivation seasons of different varieties, thus affecting final yield. Early heading date 1 (Ehd1) is a major rice integrator gene in the regulatory network of heading date whose expression level is negatively correlated with heading date and grain yield. Some elite japonica varieties such as Zhongjia 8 (ZJ8) show very early heading with poor agronomic traits when planted in South China. This problem can be addressed by downregulating the expression of Ehd1. In this study, we analyzed the cis-regulatory elements in the Ehd1 promoter region. We then used CRISPR/Cas9-mediated editing to modify the Ehd1 promoter at multiple target sites in ZJ8. We rapidly identified homozygous allelic mutations in the T2 generation via long-read sequencing. We obtained several Ehd1 promoter mutants with different degrees of lower Ehd1 expression, delayed heading date, and improved yield-related traits. We developed an efficient promoter editing method to create different weak Ehd1 alleles for breeding selection. Using this method, a series of heading date materials from elite varieties can be created to expand the planting area of rice and improve grain yields.

Idiosyncratic and dose-dependent epistasis drives variation in tomato fruit size

Epistasis between genes is traditionally studied with mutations that eliminate protein activity, but most natural genetic variation is in cis-regulatory DNA and influences gene expression and function quantitatively. In this study, we used natural and engineered cis-regulatory alleles in a plant stem-cell circuit to systematically evaluate epistatic relationships controlling tomato fruit size. Combining a promoter allelic series with two other loci, we collected over 30,000 phenotypic data points from 46 genotypes to quantify how allele strength transforms epistasis. We revealed a saturating dose-dependent relationship but also allele-specific idiosyncratic interactions, including between alleles driving a step change in fruit size during domestication. Our approach and findings expose an underexplored dimension of epistasis, in which cis-regulatory allelic diversity within gene regulatory networks elicits nonlinear, unpredictable interactions that shape phenotypes.

CRISPR-Cas-mediated unfolded protein response control for enhancing plant stress resistance

Plants consistently encounter environmental stresses that negatively affect their growth and development. To mitigate these challenges, plants have developed a range of adaptive strategies, including the unfolded protein response (UPR), which enables them to manage endoplasmic reticulum (ER) stress resulting from various adverse conditions. The CRISPR-Cas system has emerged as a powerful tool for plant biotechnology, with the potential to improve plant tolerance and resistance to biotic and abiotic stresses, as well as enhance crop productivity and quality by targeting specific genes, including those related to the UPR. This review highlights recent advancements in UPR signaling pathways and CRISPR-Cas technology, with a particular focus on the use of CRISPR-Cas in studying plant UPR. We also explore prospective applications of CRISPR-Cas in engineering UPR-related genes for crop improvement. The integration of CRISPR-Cas technology into plant biotechnology holds the promise to revolutionize agriculture by producing crops with enhanced resistance to environmental stresses, increased productivity, and improved quality traits.

Redundancy and the role of protein copy numbers in the cell polarization machinery of budding yeast

How can a self-organized cellular function evolve, adapt to perturbations, and acquire new sub-functions? To make progress in answering these basic questions of evolutionary cell biology, we analyze, as a concrete example, the cell polarity machinery of Saccharomyces cerevisiae. This cellular module exhibits an intriguing resilience: it remains operational under genetic perturbations and recovers quickly and reproducibly from the deletion of one of its key components. Using a combination of modeling, conceptual theory, and experiments, we propose that multiple, redundant self-organization mechanisms coexist within the protein network underlying cell polarization and are responsible for the module's resilience and adaptability. Based on our mechanistic understanding of polarity establishment, we hypothesize that scaffold proteins, by introducing new connections in the existing network, can increase the redundancy of mechanisms and thus increase the evolvability of other network components. Moreover, our work gives a perspective on how a complex, redundant cellular module might have evolved from a more rudimental ancestral form.

Gene editing of non-coding regulatory DNA and its application in crop improvement

The development of the clustered regularly interspaced short palindromic repeats/CRISPR-associated protein (CRISPR/Cas) system has provided precise and efficient strategies to edit target genes and generate transgene-free crops. Significant progress has been made in the editing of protein-coding genes; however, studies on the editing of non-coding DNA with regulatory roles lags far behind. Non-coding regulatory DNAs, including those which can be transcribed into long non-coding RNAs (lncRNAs), and miRNAs, together with cis-regulatory elements (CREs), play crucial roles in regulating plant growth and development. Therefore, the combination of CRISPR/Cas technology and non-coding regulatory DNA has great potential to generate novel alleles that affect various agronomic traits of crops, thus providing valuable genetic resources for crop breeding. Herein, we review recent advances in the roles of non-coding regulatory DNA, attempts to edit non-coding regulatory DNA for crop improvement, and potential application of novel editing tools in modulating non-coding regulatory DNA. Finally, the existing problems, possible solutions, and future applications of gene editing of non-coding regulatory DNA in modern crop breeding practice are also discussed.

Sexual conflict, heterochrony and tissue specificity as evolutionary problems of adaptive plasticity in development

Differential gene expression represents a fundamental cause and manifestation of phenotypic plasticity. Adaptive phenotypic plasticity in gene expression as a trait evolves when alleles that mediate gene regulation serve to increase organismal fitness by improving the alignment of variation in gene expression with variation in circumstances. Among the diverse circumstances that a gene encounters are distinct cell types, developmental stages and sexes, as well as an organism's extrinsic ecological environments. Consequently, adaptive phenotypic plasticity provides a common framework to consider diverse evolutionary problems by considering the shared implications of alleles that produce context-dependent gene expression. From this perspective, adaptive plasticity represents an evolutionary resolution to conflicts of interest that arise from any negatively pleiotropic effects of expression of a gene across ontogeny, among tissues, between the sexes, or across extrinsic environments. This view highlights shared properties within the general relation of fitness, trait expression and context that may nonetheless differ substantively in the grain of selection within and among generations to influence the likelihood of adaptive plasticity as an evolutionary response. Research programmes that historically have focused on these separate issues may use the insights from one another by recognizing their shared dependence on context-dependent gene regulatory evolution.

Deep flanking sequence engineering for efficient promoter design using DeepSEED

Designing promoters with desirable properties is essential in synthetic biology. Human experts are skilled at identifying strong explicit patterns in small samples, while deep learning models excel at detecting implicit weak patterns in large datasets. Biologists have described the sequence patterns of promoters via transcription factor binding sites (TFBSs). However, the flanking sequences of cis-regulatory elements, have long been overlooked and often arbitrarily decided in promoter design. To address this limitation, we introduce DeepSEED, an AI-aided framework that efficiently designs synthetic promoters by combining expert knowledge with deep learning techniques. DeepSEED has demonstrated success in improving the properties of Escherichia coli constitutive, IPTG-inducible, and mammalian cell doxycycline (Dox)-inducible promoters. Furthermore, our results show that DeepSEED captures the implicit features in flanking sequences, such as k-mer frequencies and DNA shape features, which are crucial for determining promoter properties.

Asymmetric Evolution of Protein Domains in the Leucine-Rich Repeat Receptor-Like Kinase Family of Plant Signaling Proteins

The coding sequences of developmental genes are expected to be deeply conserved, with cis-regulatory change driving the modulation of gene function. In contrast, proteins with roles in defense are expected to evolve rapidly, in molecular arms races with pathogens. However, some gene families include both developmental and defense genes. In these families, does the tempo and mode of evolution differ between genes with divergent functions, despite shared ancestry and structure? The leucine-rich repeat receptor-like kinase (LRR-RLKs) protein family includes members with roles in plant development and defense, thus providing an ideal system for answering this question. LRR-RLKs are receptors that traverse plasma membranes. LRR domains bind extracellular ligands; RLK domains initiate intracellular signaling cascades in response to ligand binding. In LRR-RLKs with roles in defense, LRR domains evolve faster than RLK domains. To determine whether this asymmetry extends to LRR-RLKs that function primarily in development, we assessed evolutionary rates and tested for selection acting on 11 subfamilies of LRR-RLKs, using deeply sampled protein trees. To assess functional evolution, we performed heterologous complementation assays in Arabidopsis thaliana (Arabidopsis). We found that the LRR domains of all tested LRR-RLK proteins evolved faster than their cognate RLK domains. All tested subfamilies of LRR-RLKs had strikingly similar patterns of molecular evolution, despite divergent functions. Heterologous transformation experiments revealed that multiple mechanisms likely contribute to the evolution of LRR-RLK function, including escape from adaptive conflict. Our results indicate specific and distinct evolutionary pressures acting on LRR versus RLK domains, despite diverse organismal roles for LRR-RLK proteins.

Developmental regulation of conserved non-coding element evolution provides insights into limb loss in squamates

Limb loss shows recurrent phenotypic evolution across squamate lineages. Here, based on three de novo-assembled genomes of limbless lizards from different lineages, we showed that divergence of conserved non-coding elements (CNEs) played an important role in limb development. These CNEs were associated with genes required for limb initiation and outgrowth, and with regulatory signals in the early stage of limb development. Importantly, we identified the extensive existence of insertions and deletions (InDels) in the CNEs, with the numbers ranging from 111 to 756. Most of these CNEs with InDels were lineage-specific in the limbless squamates. Nearby genes of these InDel CNEs were important to early limb formation, such as Tbx4, Fgf10, and Gli3. Based on functional experiments, we found that nucleotide mutations and InDels both affected the regulatory function of the CNEs. Our study provides molecular evidence underlying limb loss in squamate reptiles from a developmental perspective and sheds light on the importance of regulatory element InDels in phenotypic evolution.

Genome-wide identification and expression analysis of the trihelix transcription factor family in sesame (Sesamum indicum L.) under abiotic stress

BACKGROUND

The plant trihelix gene family is among the earliest discovered transcription factor families, and it is vital in modulating light, plant growth, and stress responses.

METHODS

The identification and characterization of trihelix family members in the sesame genome were analyzed by bioinformatics methods, and the expression patterns of sesame trihelix genes were assessed by quantitative real-time PCR.

RESULTS

There were 34 trihelix genes discovered in the genome of sesame, which were irregularly distributed among 10 linkage groups. Also, the genome contained 5 duplicate gene pairs. The 34 trihelix genes were divided into six sub-families through a phylogenetic study. A tissue-specific expression revealed that SiTH genes exhibited spatial expression patterns distinct from other trihelix genes in the same subfamily. The cis-element showed that the SiTHs gene promoter contained various elements associated with responses to hormones and multiple abiotic stresses. Additionally, the expression patterns of 8 SiTH genes in leaves under abiotic stresses demonstrated that all selected genes were significantly upregulated or downregulated at least once in the stress period. Furthermore, the SiTH4 gene was significantly induced in response to drought and salt stress, showing that SiTH genes may be engaged in the stress response mechanisms of sesame.

CONCLUSION

These findings establish a foundation for further investigation of the trihelix gene-mediated response to abiotic stress in sesame.

Advances in cis-element- and natural variation-mediated transcriptional regulation and applications in gene editing of major crops

Transcriptional regulation is crucial to control of gene expression. Both spatio-temporal expression patterns and expression levels of genes are determined by the interaction between cis-acting elements and trans-acting factors. Numerous studies have focused on the trans-acting factors that mediate transcriptional regulatory networks. However, cis-acting elements, such as enhancers, silencers, transposons, and natural variations in the genome, are also vital for gene expression regulation and could be utilized by clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9)-mediated gene editing to improve crop quality and yield. In this review, we discuss current understanding of cis-element-mediated transcriptional regulation in major crops, including rice (Oryza sativa), wheat (Triticum aestivum), and maize (Zea mays), as well as the latest advancements in gene editing techniques and their applications in crops to highlight prospective strategies for crop breeding.

Genome-wide analysis of the Pleurotus eryngii laccase gene (PeLac) family and functional identification of PeLac5

The laccase gene family encodes multiple isozymes that are crucial for the degradation of substrates and the regulation of developmental processes in fungi. Pleurotus eryngii is an important edible and medicinal fungus belonging to the Basidiomycota phylum and can grow on a variety of natural substrates. In the present study, genome-wide profiling of P. eryngii identified 10 genes encoding its laccase isoenzymes. Conservative sequence analysis demonstrated that all PeLacs possess classical laccase structural domains. Phylogenetic analysis yielded four major subgroups, the members of which are similar with respect to conserved gene organization, protein domain architecture, and consensus motifs. The 10 PeLacs formed three groups together with 12 PoLacs in Pleurotus ostreatus, indicating that they share a high level of evolutionary homology. Cis-responsive element analysis implied that PeLacs genes play a role in growth and development and lignocellulose degradation. Targeted overexpression of PeLac5 reduced the time to primordia formation and their development to fruiting bodies. Gene expression patterns in the presence of different lignocellulosic substrates indicate that three PeLacs genes (2, 4, and 9) are key to lignocellulose degradation. This work presents the first inventory of laccase genes in P. eryngii and preliminarily explores their functions, which may help to uncover the manner by which these proteins utilize substrates.

Environmentally robust cis-regulatory changes underlie rapid climatic adaptation

Changes in gene expression are thought to play a major role in adaptive evolution. While it is known that gene expression is highly sensitive to the environment, very few studies have determined the influence of genetic and environmental effects on adaptive gene expression differences in natural populations. Here, we utilize allele-specific expression to characterize cis and trans gene regulatory divergence in temperate and tropical house mice in two metabolic tissues under two thermal conditions. First, we show that gene expression divergence is pervasive between populations and across thermal conditions, with roughly 5 to 10% of genes exhibiting genotype-by-environment interactions. Second, we found that most expression divergence was due to cis-regulatory changes that were stable across temperatures. In contrast, patterns of expression plasticity were largely attributable to trans-effects, which showed greater sensitivity to temperature. Nonetheless, we found a small subset of temperature-dependent cis-regulatory changes, thereby identifying loci underlying expression plasticity. Finally, we performed scans for selection in wild house mice to identify genomic signatures of rapid adaptation. Genomic outliers were enriched in genes with evidence for cis-regulatory divergence. Notably, these genes were associated with phenotypes that affected body weight and metabolism, suggesting that cis-regulatory changes are a possible mechanism for adaptive body size evolution between populations. Our results show that gene expression plasticity, largely controlled in trans, may facilitate the colonization of new environments, but that evolved changes in gene expression are largely controlled in cis, illustrating the genetic and nongenetic mechanisms underlying the establishment of populations in new environments.

Genome-wide analysis of the GRAS gene family in Liriodendron chinense reveals the putative function in abiotic stress and plant development

Introduction

GRAS genes encode plant-specific transcription factors that play essential roles in plant growth and development. However, the members and the function of the GRAS gene family have not been reported in Liriodendron chinense. L. chinense, a tree species in the Magnolia family that produces excellent timber for daily life and industry. In addition, it is a good relict species for plant evolution research.

Methods

Therefore, we conducted a genome-wide study of the LcGRAS gene family and identified 49 LcGRAS genes in L. chinense.

Results

We found that LcGRAS could be divided into 13 sub-groups, among which there is a unique branch named HAM-t. We carried out RNA sequencing analysis of the somatic embryos from L. chinense and found that LcGRAS genes are mainly expressed after heart-stage embryo development, suggesting that LcGRAS may have a function during somatic embryogenesis. We also investigated whether GRAS genes are responsive to stress by carrying out RNA sequencing (RNA-seq) analysis, and we found that the genes in the PAT subfamily were activated upon stress treatment, suggesting that these genes may help plants survive stressful environments. We found that PIF was downregulated and COR was upregulated after the transient overexpression of PATs, suggesting that PAT may be upstream regulators of cold stress.

Discussion

Collectively, LcGRAS genes are conserved and play essential roles in plant development and adaptation to abiotic stress.

Sequencing 4.3 million mutations in wheat promoters to understand and modify gene expression

Wheat is an important contributor to global food security, and further improvements are required to feed a growing human population. Functional genetics and genomics tools can help us to understand the function of different genes and to engineer beneficial changes. In this study, we used a promoter capture assay to sequence 2-kb regions upstream of all high-confidence annotated genes from 1,513 mutagenized plants from the tetraploid wheat variety Kronos. We identified 4.3 million induced mutations with an accuracy of 99.8%, resulting in a mutation density of 41.9 mutations per kb. We also remapped Kronos exome capture reads to Chinese Spring RefSeq v1.1, identified 4.7 million mutations, and predicted their effects on annotated genes. Using these predictions, we identified 59% more nonsynonymous substitutions and 49% more truncation mutations than in the original study. To show the biological value of the promoter dataset, we selected two mutations within the promoter of the VRN-A1 vernalization gene. Both mutations, located within transcription factor binding sites, significantly altered VRN-A1 expression, and one reduced the number of spikelets per spike. These publicly available sequenced mutant datasets provide rapid and inexpensive access to induced variation in the promoters and coding regions of most wheat genes. These mutations can be used to understand and modulate gene expression and phenotypes for both basic and commercial applications, where limited governmental regulations can facilitate deployment. These mutant collections, together with gene editing, provide valuable tools to accelerate functional genetic studies in this economically important crop.

A Quantitative, Genome-Wide Analysis in Drosophila Reveals Transposable Elements' Influence on Gene Expression Is Species-Specific

Transposable elements (TEs) are parasite DNA sequences that are able to move and multiply along the chromosomes of all genomes. They can be controlled by the host through the targeting of silencing epigenetic marks, which may affect the chromatin structure of neighboring sequences, including genes. In this study, we used transcriptomic and epigenomic high-throughput data produced from ovarian samples of several Drosophila melanogaster and Drosophila simulans wild-type strains, in order to finely quantify the influence of TE insertions on gene RNA levels and histone marks (H3K9me3 and H3K4me3). Our results reveal a stronger epigenetic effect of TEs on ortholog genes in D. simulans compared with D. melanogaster. At the same time, we uncover a larger contribution of TEs to gene H3K9me3 variance within genomes in D. melanogaster, which is evidenced by a stronger correlation of TE numbers around genes with the levels of this chromatin mark in D. melanogaster. Overall, this work contributes to the understanding of species-specific influence of TEs within genomes. It provides a new light on the considerable natural variability provided by TEs, which may be associated with contrasted adaptive and evolutionary potentials.

A genome‑wide approach to the systematic and comprehensive analysis of LIM gene family in sorghum (Sorghum bicolor L.)

The LIM domain-containing proteins are dominantly found in plants and play a significant role in various biological processes such as gene transcription as well as actin cytoskeletal organization. Nevertheless, genome-wide identification as well as functional analysis of the LIM gene family have not yet been reported in the economically important plant sorghum (Sorghum bicolor L.). Therefore, we conducted an in silico identification and characterization of LIM genes in S. bicolor genome using integrated bioinformatics approaches. Based on phylogenetic tree analysis and conserved domain, we identified five LIM genes in S. bicolor (SbLIM) genome corresponding to Arabidopsis LIM (AtLIM) genes. The conserved domain, motif as well as gene structure analyses of the SbLIM gene family showed the similarity within the SbLIM and AtLIM members. The gene ontology (GO) enrichment study revealed that the candidate LIM genes are directly involved in cytoskeletal organization and various other important biological as well as molecular pathways. Some important families of regulating transcription factors such as ERF, MYB, WRKY, NAC, bZIP, C2H2, Dof, and G2-like were detected by analyzing their interaction network with identified SbLIM genes. The cis-acting regulatory elements related to predicted SbLIM genes were identified as responsive to light, hormones, stress, and other functions. The present study will provide valuable useful information about LIM genes in sorghum which would pave the way for the future study of functional pathways of candidate SbLIM genes as well as their regulatory factors in wet-lab experiments.

A digital twin reproducing gene regulatory network dynamics of early Ciona embryos indicates robust buffers in the network

How gene regulatory networks (GRNs) encode gene expression dynamics and how GRNs evolve are not well understood, although these problems have been studied extensively. We created a digital twin that accurately reproduces expression dynamics of 13 genes that initiate expression in 32-cell ascidian embryos. We first showed that gene expression patterns can be manipulated according to predictions by this digital model. Next, to simulate GRN rewiring, we changed regulatory functions that represented their regulatory mechanisms in the digital twin, and found that in 55 of 100 cases, removal of a single regulator from a conjunctive clause of Boolean functions did not theoretically alter qualitative expression patterns of these genes. In other words, we found that more than half the regulators gave theoretically redundant temporal or spatial information to target genes. We experimentally substantiated that the expression pattern of Nodal was maintained without one of these factors, Zfpm, by changing the upstream regulatory sequence of Nodal. Such robust buffers of regulatory mechanisms may provide a basis of enabling developmental system drift, or rewiring of GRNs without changing expression patterns of downstream genes, during evolution.

The Ancestral Modulation Hypothesis: Predicting Mechanistic Control of Sexually Heteromorphic Traits Using Evolutionary History

AbstractAcross the animal kingdom there are myriad forms within a sex across, and even within, species, rendering concepts of universal sex traits moot. The mechanisms that regulate the development of these trait differences are varied, although in vertebrates, common pathways involve gonadal steroid hormones. Gonadal steroids are often associated with heteromorphic trait development, where the steroid found at higher circulating levels is the one involved in trait development for that sex. Occasionally, there are situations in which a gonadal steroid associated with heteromorphic trait development in one sex is involved in heteromorphic or monomorphic trait development in another sex. We propose a verbal hypothesis, the ancestral modulation hypothesis (AMH), that uses the evolutionary history of the trait-particularly which sex ancestrally possessed higher trait values-to predict the regulatory pathway that governs trait expression. The AMH predicts that the genomic architecture appears first to resolve sexual conflict in an initially monomorphic trait. This architecture takes advantage of existing sex-biased signals, the gonadal steroid pathway, to generate trait heteromorphism. In cases where the other sex experiences evolutionary pressure for the new phenotype, that sex will co-opt the existing architecture by altering its signal to match that of the original high-trait-value sex. We describe the integrated levels needed to produce this pattern and what the expected outcomes will be given the evolutionary history of the trait. We present this framework as a testable hypothesis for the scientific community to investigate and to create further engagement and analysis of both ultimate and proximate approaches to sexual heteromorphism.

Genome-wide identification and molecular evolution of NAC gene family in Dendrobium nobile

NAC transcription factors are an important genes that regulate plant growth and development, and can regulate functions such as fruit ripening in plants. Based on genome data of Dendrobium nobile, the NAC gene family was identified and analyzed by bioinformatics methods. In this study, we identified 85 NAC genes in Dendrobium nobile genome, and systematically analyzed the NAC gene family. We found that they were distributed unevenly in the nineteen chromosomes. The amino acid length of D. nobile NAC gene family (DnoNACs) ranged from 80 to 1065, molecular weight ranged from 22.17 to 119.02 kD, and isoelectric point ranged from 4.61~9.26. Its promoter region contains multiple stress responsive elements, including light responsive, gibberellin-responsive, abscisic acid responsiveness, MeJA-responsiveness and drought-inducibility elements. Phylogenetic analysis indicates that the D. nobile NAC gene family is most closely related to Dendrobium catenatum and Dendrobium chrysotoxum. Analysis of SSR loci indicates that the fraction of mononucleotide repeats was the largest, as was the frequency of A/T. Non-coding RNA analysis showed that these 85 NAC genes contain 397 miRNAs. The collinearity analysis shows that 9 collinear locis were found on the chromosomes of D. nobile with Arabidopsis thaliana, and 75 collinear locis with D.chrysotoxum. QRT-PCR experiment under different salt concentration and temperature conditions verified the response mechanism of DnoNAC gene family under stress conditions. Most DnoNAC genes are sensitive to salt stress and temperature stress. The results of this study provide a reference for further understanding the function of NAC gene in D. nobile.

A genome-wide genetic screen uncovers determinants of human pigmentation

Skin color, one of the most diverse human traits, is determined by the quantity, type, and distribution of melanin. In this study, we leveraged the light-scattering properties of melanin to conduct a genome-wide screen for regulators of melanogenesis. We identified 169 functionally diverse genes that converge on melanosome biogenesis, endosomal transport, and gene regulation, of which 135 represented previously unknown associations with pigmentation. In agreement with their melanin-promoting function, the majority of screen hits were up-regulated in melanocytes from darkly pigmented individuals. We further unraveled functions of KLF6 as a transcription factor that regulates melanosome maturation and pigmentation in vivo, and of the endosomal trafficking protein COMMD3 in modulating melanosomal pH. Our study reveals a plethora of melanin-promoting genes, with broad implications for human variation, cell biology, and medicine.

Machine-guided design of synthetic cell type-specific cis-regulatory elements

Cis-regulatory elements (CREs) control gene expression, orchestrating tissue identity, developmental timing, and stimulus responses, which collectively define the thousands of unique cell types in the body. While there is great potential for strategically incorporating CREs in therapeutic or biotechnology applications that require tissue specificity, there is no guarantee that an optimal CRE for an intended purpose has arisen naturally through evolution. Here, we present a platform to engineer and validate synthetic CREs capable of driving gene expression with programmed cell type specificity. We leverage innovations in deep neural network modeling of CRE activity across three cell types, efficient in silico optimization, and massively parallel reporter assays (MPRAs) to design and empirically test thousands of CREs. Through in vitro and in vivo validation, we show that synthetic sequences outperform natural sequences from the human genome in driving cell type-specific expression. Synthetic sequences leverage unique sequence syntax to promote activity in the on-target cell type and simultaneously reduce activity in off-target cells. Together, we provide a generalizable framework to prospectively engineer CREs and demonstrate the required literacy to write regulatory code that is fit-for-purpose in vivo across vertebrates.

Divergent combinations of cis-regulatory elements control the evolution of phenotypic plasticity

The widespread occurrence of phenotypic plasticity across all domains of life demonstrates its evolutionary significance. However, how plasticity itself evolves and how it contributes to evolution is poorly understood. Here, we investigate the predatory nematode Pristionchus pacificus with its feeding structure plasticity using recombinant-inbred-line and quantitative-trait-locus (QTL) analyses between natural isolates. We show that a single QTL at a core developmental gene controls the expression of the cannibalistic morph. This QTL is composed of several cis-regulatory elements. Through CRISPR/Cas-9 engineering, we identify copy number variation of potential transcription factor binding sites that interacts with a single intronic nucleotide polymorphism. Another intronic element eliminates gene expression altogether, mimicking knockouts of the locus. Comparisons of additional isolates further support the rapid evolution of these cis-regulatory elements. Finally, an independent QTL study reveals evidence for parallel evolution at the same locus. Thus, combinations of cis-regulatory elements shape plastic trait expression and control nematode cannibalism.

Responsiveness to endurance training can be partly explained by the number of favorable single nucleotide polymorphisms an individual possesses

Cardiorespiratory fitness is a key component of health-related fitness. It is a necessary focus of improvement, especially for those that have poor fitness and are classed as untrained. However, much research has shown individuals respond differentially to identical training programs, suggesting the involvement of a genetic component in individual exercise responses. Previous research has focused predominantly on a relatively low number of candidate genes and their overall influence on exercise responsiveness. However, examination of gene-specific alleles may provide a greater level of understanding. Accordingly, this study aimed to investigate the associations between cardiorespiratory fitness and an individual's genotype following a field-based endurance program within a previously untrained population. Participants (age: 29 ± 7 years, height: 175 ± 9 cm, mass: 79 ± 21 kg, body mass index: 26 ± 7 kg/m2) were randomly assigned to either a training (n = 21) or control group (n = 24). The training group completed a periodized running program for 8-weeks (duration: 20-30-minutes per session, intensity: 6-7 Borg Category-Ratio-10 scale rating, frequency: 3 sessions per week). Both groups completed a Cooper 12-minute run test to estimate cardiorespiratory fitness at baseline, mid-study, and post-study. One thousand single nucleotide polymorphisms (SNPs) were assessed via saliva sample collections. Cooper run distance showed a significant improvement (0.23 ± 0.17 km [11.51 ± 9.09%], p < 0.001, ES = 0.48 [95%CI: 0.16-0.32]), following the 8-week program, whilst controls displayed no significant changes (0.03 ± 0.15 km [1.55 ± 6.98%], p = 0.346, ES = 0.08, [95%CI: -0.35-0.95]). A significant portion of the inter-individual variation in Cooper scores could be explained by the number of positive alleles a participant possessed (r = 0.92, R2 = 0.85, p < 0.001). These findings demonstrate the relative influence of key allele variants on an individual's responsiveness to endurance training.

Challenges and advances towards the rational design of microalgal synthetic promoters in Chlamydomonas reinhardtii

Microalgae hold enormous potential to provide a safe and sustainable source of high-value compounds, acting as carbon-fixing biofactories that could help to mitigate rapidly progressing climate change. Bioengineering microalgal strains will be key to optimizing and modifying their metabolic outputs, and to render them competitive with established industrial biotechnology hosts, such as bacteria or yeast. To achieve this, precise and tuneable control over transgene expression will be essential, which would require the development and rational design of synthetic promoters as a key strategy. Among green microalgae, Chlamydomonas reinhardtii represents the reference species for bioengineering and synthetic biology; however, the repertoire of functional synthetic promoters for this species, and for microalgae generally, is limited in comparison to other commercial chassis, emphasizing the need to expand the current microalgal gene expression toolbox. Here, we discuss state-of-the-art promoter analyses, and highlight areas of research required to advance synthetic promoter development in C. reinhardtii. In particular, we exemplify high-throughput studies performed in other model systems that could be applicable to microalgae, and propose novel approaches to interrogating algal promoters. We lastly outline the major limitations hindering microalgal promoter development, while providing novel suggestions and perspectives for how to overcome them.

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 quickly results in cell and tissue damage that can range in severity from temporary injury to permanent neurological damage. Because of its biological relevance, and because many conserved copper-responsive genes also respond to other non-essential heavy metal pollutants, copper resistance in Drosophila melanogaster is a useful model system with which to investigate the genetic control of the response to heavy metal stress. 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 (DSPR) using a combination of differential expression analysis and expression quantitative trait locus (eQTL) mapping. Differential expression analysis revealed clear patterns of tissue-specific expression, primarily driven by a more pronounced gene expression response in gut tissue. eQTL mapping of gene expression under control and copper conditions as well as for the change in gene expression following copper exposure (copper response eQTL) revealed hundreds of genes with tissue-specific local cis-eQTL and many distant trans-eQTL. eQTL associated with MtnA, Mdr49, Mdr50, and Sod3 exhibited genotype by environment effects on gene expression under copper stress, illuminating several 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 many candidate genes for future functional characterization.

Cis-regulatory landscapes in the evolution and development of the mammalian skull

Extensive morphological variation found in mammals reflects the wide spectrum of their ecological adaptations. The highest morphological diversity is present in the craniofacial region, where geometry is mainly dictated by the bony skull. Mammalian craniofacial development represents complex multistep processes governed by numerous conserved genes that require precise spatio-temporal control. A central question in contemporary evolutionary biology is how a defined set of conserved genes can orchestrate formation of fundamentally different structures, and therefore how morphological variability arises. In principle, differential gene expression patterns during development are the source of morphological variation. With the emergence of multicellular organisms, precise regulation of gene expression in time and space is attributed to cis-regulatory elements. These elements contribute to higher-order chromatin structure and together with trans-acting factors control transcriptional landscapes that underlie intricate morphogenetic processes. Consequently, divergence in cis-regulation is believed to rewire existing gene regulatory networks and form the core of morphological evolution. This review outlines the fundamental principles of the genetic code and genomic regulation interplay during development. Recent work that deepened our comprehension of cis-regulatory element origin, divergence and function is presented here to illustrate the state-of-the-art research that uncovered the principles of morphological novelty. This article is part of the theme issue 'The mammalian skull: development, structure and function'.

Speciation and development

Understanding general principles about the origin of species remains one of the foundational challenges in evolutionary biology. The genomic divergence between groups of individuals can spawn hybrid inviability and hybrid sterility, which presents a tantalizing developmental problem. Divergent developmental programs may yield either conserved or divergent phenotypes relative to ancestral traits, both of which can be responsible for reproductive isolation during the speciation process. The genetic mechanisms of developmental evolution involve cis- and trans-acting gene regulatory change, protein-protein interactions, genetic network structures, dosage, and epigenetic regulation, all of which also have roots in population genetic and molecular evolutionary processes. Toward the goal of demystifying Darwin's "mystery of mysteries," this review integrates microevolutionary concepts of genetic change with principles of organismal development, establishing explicit links between population genetic process and developmental mechanisms in the production of macroevolutionary pattern. This integration aims to establish a more unified view of speciation that binds process and mechanism.