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

OsTBP2.1, a TATA-Binding Protein, Alters the Ratio of OsNRT2.3b to OsNRT2.3a and Improves Rice Grain Yield

The OsNRT2.3a and OsNRT2.3b isoforms play important roles in the uptake and transport of nitrate during rice growth. However, it is unclear which cis-acting element controls the transcription of OsNRT2.3 into these specific isoforms. In this study, we used a yeast one-hybrid assay to obtain the TATA-box binding protein OsTBP2.1, which binds to the TATA-box of OsNRT2.3, and verified its important role through transient expression and RNA-seq. We found that the TATA-box of OsNRT2.3 mutants and binding protein OsTBP2.1 together increased the transcription ratio of OsNRT2.3b to OsNRT2.3a. The overexpression of OsTBP2.1 promoted nitrogen uptake and increased rice yield compared with the wild-type; however, the OsTBP2.1 T-DNA mutant lines exhibited the opposite trend. Detailed analyses demonstrated that the TATA-box was the key cis-regulatory element for OsNRT2.3 to be transcribed into OsNRT2.3a and OsNRT2.3b. Additionally, this key cis-regulatory element, together with the binding protein OsTBP2.1, promoted the development of rice and increased grain yield.

Gnocis: An integrated system for interactive and reproducible analysis and modelling of cis-regulatory elements in Python 3

Gene expression is regulated through cis-regulatory elements (CREs), among which are promoters, enhancers, Polycomb/Trithorax Response Elements (PREs), silencers and insulators. Computational prediction of CREs can be achieved using a variety of statistical and machine learning methods combined with different feature space formulations. Although Python packages for DNA sequence feature sets and for machine learning are available, no existing package facilitates the combination of DNA sequence feature sets with machine learning methods for the genome-wide prediction of candidate CREs. We here present Gnocis, a Python package that streamlines the analysis and the modelling of CRE sequences by providing extensible APIs and implementing the glue required for combining feature sets and models for genome-wide prediction. Gnocis implements a variety of base feature sets, including motif pair occurrence frequencies and the k-spectrum mismatch kernel. It integrates with Scikit-learn and TensorFlow for state-of-the-art machine learning. Gnocis additionally implements a broad suite of tools for the handling and preparation of sequence, region and curve data, which can be useful for general DNA bioinformatics in Python. We also present Deep-MOCCA, a neural network architecture inspired by SVM-MOCCA that achieves moderate to high generalization without prior motif knowledge. To demonstrate the use of Gnocis, we applied multiple machine learning methods to the modelling of D. melanogaster PREs, including a Convolutional Neural Network (CNN), making this the first study to model PREs with CNNs. The models are readily adapted to new CRE modelling problems and to other organisms. In order to produce a high-performance, compiled package for Python 3, we implemented Gnocis in Cython. Gnocis can be installed using the PyPI package manager by running ‘pip install gnocis’. The source code is available on GitHub, at https://github.com/bjornbredesen/gnocis.

The Genetic Basis of Gene Expression Divergence in Antennae of Two Closely Related Moth Species, Helicoverpa armigera and Helicoverpa assulta

The closely related species Helicoverpa armigera (H. armigera) and Helicoverpa assulta (H. assulta) have different host plant ranges and share two principal components of sex pheromones but with reversed ratios. The antennae are the main olfactory organ of insects and play a crucial role in host plant selection and mate seeking. However, the genetic basis for gene expression divergence in the antennae of the two species is unclear. We performed an allele-specific expression (ASE) analysis in the antennal transcriptomes of the two species and their F₁ hybrids, examining the connection between gene expression divergence and phenotypic differences. The results show that the proportion of genes classified as all cis was higher than that of all trans in males and reversed in females. The contribution of regulatory patterns to gene expression divergence in males was less than that in females, which explained the functional differentiation of male and female antennae. Among the five groups of F₁ hybrids, the fertile males from the cross of H. armigera female and H. assulta male had the lowest proportion of misexpressed genes, and the inferred regulatory patterns were more accurate. By using this group of F₁ hybrids, we discovered that cis-related regulations play a crucial role in gene expression divergence of sex pheromone perception-related proteins. These results are helpful for understanding how specific changes in the gene expression of olfactory-related genes can contribute to rapid evolutionary changes in important olfactory traits in closely related moths.

A regulatory network of Sox and Six transcription factors initiate a cell fate transformation during hearing regeneration in adult zebrafish

Using adult zebrafish inner ears as a model for sensorineural regeneration, we ablated the mechanosensory receptors and characterized the single-cell epigenome and transcriptome at consecutive time points during hair cell regeneration. We utilized deep learning on the regeneration-induced open chromatin sequences and identified cell-specific transcription factor (TF) motif patterns. Enhancer activity correlated with gene expression and identified potential gene regulatory networks. A pattern of overlapping Sox- and Six-family TF gene expression and binding motifs was detected, suggesting a combinatorial program of TFs driving regeneration and cell identity. Pseudotime analysis of single-cell transcriptomic data suggested that support cells within the sensory epithelium changed cell identity to a “progenitor” cell population that could differentiate into hair cells. We identified a 2.6 kb DNA enhancer upstream of the sox2 promoter that, when deleted, showed a dominant phenotype that resulted in a hair-cell-regeneration-specific deficit in both the lateral line and adult inner ear.

Jimenez et al. interrogate the epigenomic and transcriptomic landscape of regenerating adult zebrafish inner-ear sensory epithelia. They show that the support-cell population transitions to an intermediate “progenitor” cell state that becomes new hair cells, and they demonstrate that the cell fate decisions may be driven by the coordinate regulation and spatial co-binding of Sox and Six transcription factors. By functionally validating a predicted regeneration-responsive enhancer upstream of sox2, they show that precise timing of sox2 expression is critical for hearing regeneration in zebrafish.

Fine mapping and candidate gene analysis of gummy stem blight resistance in cucumber stem

Two candidate genes (Csa6G046210 and Csa6G046240) were identified by fine-mapping gsb-s6.2 for gummy stem blight resistance in cucumber stem. Gummy stem blight (GSB) is a serious fungal disease caused by Didymella bryoniae, that affects cucumber yield and quality worldwide. However, no GSB-resistant genes have been identified in cucumber cultivars. In this study, the wild cucumber accession 'PI 183967' was used as a source of resistance to GSB in adult stems. An F₂ population was mapped using resistant line 'LM189' and susceptible line 'LM6' derived from a cross between 'PI 183967' and '931'. By developing InDel and SNP markers, the gsb-s6.2 QTL on Chr. 6 was fine-mapped to a 34 kb interval harboring six genes. Gene Expression analysis after inoculation showed that two candidate genes (Csa6G046210 and Csa6G046240) were induced and differentially expressed between the resistant and susceptible parents, and may be involved in disease defense. Sequence alignment showed that Csa6G046210 encodes a multiple myeloma tumor-associated protein, and it harbored two nonsynonymous SNPs and one InDel in the third and the fourth exons, and two InDels in the TATA-box of the basal promoter region. Csa6G046240 encodes a MYB transcription factor with six variants in the AP2/ERF and MYB motifs in the promoter. These two candidate genes lay the foundation for revealing the mechanism of GSB resistance and may be useful for marker-assisted selection in cucumber disease-resistant breeding.

Characterization of sequence determinants of enhancer function using natural genetic variation

Sequence variation in enhancers that control cell-type-specific gene transcription contributes significantly to phenotypic variation within human populations. However, it remains difficult to predict precisely the effect of any given sequence variant on enhancer function due to the complexity of DNA sequence motifs that determine transcription factor (TF) binding to enhancers in their native genomic context. Using F₁-hybrid cells derived from crosses between distantly related inbred strains of mice, we identified thousands of enhancers with allele-specific TF binding and/or activity. We find that genetic variants located within the central region of enhancers are most likely to alter TF binding and enhancer activity. We observe that the AP-1 family of TFs (Fos/Jun) are frequently required for binding of TEAD TFs and for enhancer function. However, many sequence variants outside of core motifs for AP-1 and TEAD also impact enhancer function, including sequences flanking core TF motifs and AP-1 half sites. Taken together, these data represent one of the most comprehensive assessments of allele-specific TF binding and enhancer function to date and reveal how sequence changes at enhancers alter their function across evolutionary timescales.

A temporal gradient of cytonuclear coordination of chaperonins and chaperones during RuBisCo biogenesis in allopolyploid plants

Ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCo), consisting of subunits encoded by nuclear and cytoplasmic genes, is a model for cytonuclear evolution in plant allopolyploids. To date, coordinated cytonuclear evolutionary responses of auxiliary cofactors involved in RuBisCo biogenesis remain unexplored. This study characterized and compared genomic and transcriptional cytonuclear coevolutionary responses of chaperonin/chaperones in RuBisCo folding and assembly processes across different allopolyploids. We discovered significant cytonuclear evolutionary responses in folding cofactors, with diminishing or attenuated responses later during assembly. Our results have general significance for understanding the unrecognized cytonuclear evolution of chaperonin/chaperone genes, structural and functional features of intermediate complexes, and the functioning stage of the Raf2 cofactor. Generally, the results reveal a hitherto unexplored dimension of allopolyploidy in plants.

Ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCo) has long been studied from many perspectives. As a multisubunit (large subunits [LSUs] and small subunits[SSUs]) protein encoded by genes residing in the chloroplast (rbcL) and nuclear (rbcS) genomes, RuBisCo also is a model for cytonuclear coevolution following allopolyploid speciation in plants. Here, we studied the genomic and transcriptional cytonuclear coordination of auxiliary chaperonin and chaperones that facilitate RuBisCo biogenesis across multiple natural and artificially synthesized plant allopolyploids. We found similar genomic and transcriptional cytonuclear responses, including respective paternal-to-maternal conversions and maternal homeologous biased expression, in chaperonin/chaperon-assisted folding and assembly of RuBisCo in different allopolyploids. One observation is about the temporally attenuated genomic and transcriptional cytonuclear evolutionary responses during early folding and later assembly process of RuBisCo biogenesis, which were established by long-term evolution and immediate onset of allopolyploidy, respectively. Our study not only points to the potential widespread and hitherto unrecognized features of cytonuclear evolution but also bears implications for the structural interaction interface between LSU and Cpn60 chaperonin and the functioning stage of the Raf2 chaperone.

Molecular basis of ocean acidification sensitivity and adaptation in Mytilus galloprovincialis

Predicting the potential for species adaption to climate change is challenged by the need to identify the physiological mechanisms that underpin species vulnerability. Here, we investigated the sensitivity to ocean acidification in marine mussels during early development, and specifically the trochophore stage. Using RNA and DNA sequencing and in situ RNA hybridization, we identified developmental processes associated with abnormal development and rapid adaptation to low pH. Trochophores exposed to low pH seawater exhibited 43 differentially expressed genes. Gene annotation and in situ hybridization of differentially expressed genes point to pH sensitivity of (1) shell field development and (2) cellular stress response. Five genes within these two processes exhibited shifts in allele frequencies indicative of a potential for rapid adaptation. This case study contributes direct evidence that protecting species’ existing genetic diversity is a critical management action to facilitate species resilience to climate change.

•

Marine mussel larval development and genetic adaptation in low pH seawater

•

RNA and DNA responses reveal impacts on shell field development and cell stress

•

Five genes exhibited both physiological sensitivity and long-term adaptive potential

•

Conserving standing genetic variation could bolster resilience to global change

Evolutionary changes in the chromatin landscape contribute to reorganization of a developmental gene network during rapid life history evolution in sea urchins

Chromatin configuration is highly dynamic during embryonic development in animals, exerting an important point of control in transcriptional regulation. Yet there exists remarkably little information about the role of evolutionary changes in chromatin configuration to the evolution of gene expression and organismal traits. Genome-wide assays of chromatin configuration, coupled with whole-genome alignments, can help address this gap in knowledge in several ways. In this study we present a comparative analysis of regulatory element sequences and accessibility throughout embryogenesis in three sea urchin species with divergent life histories: a lecithotroph Heliocidaris erythrogramma, a closely related planktotroph H. tuberculata, and a distantly related planktotroph Lytechinus variegatus. We identified distinct epigenetic and mutational signatures of evolutionary modifications to the function of putative cis-regulatory elements in H. erythrogramma that have accumulated nonuniformly throughout the genome, suggesting selection, rather than drift, underlies many modifications associated with the derived life history. Specifically, regulatory elements composing the sea urchin developmental gene regulatory network are enriched for signatures of positive selection and accessibility changes which may function to alter binding affinity and access of developmental transcription factors to these sites. Furthermore, regulatory element changes often correlate with divergent expression patterns of genes involved in cell type specification, morphogenesis, and development of other derived traits, suggesting these evolutionary modifications have been consequential for phenotypic evolution in H. erythrogramma. Collectively, our results demonstrate that selective pressures imposed by changes in developmental life history rapidly reshape the cis-regulatory landscape of core developmental genes to generate novel traits and embryonic programs.

Three distinct Atoh1 enhancers cooperate for sound receptor hair cell development

Cochlear hair cells (HCs) in the inner ear are responsible for sound detection. For HC fate specification, the master transcription factor Atoh1 is both necessary and sufficient. Atoh1 expression is dynamic and tightly regulated during development, but the cis-regulatory elements mediating this regulation remain unresolved. Unexpectedly, we found that deleting the only recognized Atoh1 enhancer, defined here as Eh1, failed to impair HC development. By using the assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq), we discovered two additional Atoh1 enhancers: Eh2 and Eh3. Notably, Eh2 deletion was sufficient for impairing HC development, and concurrent deletion of Eh1 and Eh2 or all three enhancers resulted in nearly complete absence of HCs. Lastly, we showed that Atoh1 binds to all three enhancers, consistent with its autoregulatory function. Our findings reveal that the cooperative action of three distinct enhancers underpins effective Atoh1 regulation during HC development, indicating potential therapeutic approaches for HC regeneration.

Genome-wide identification and expression analysis of MATE gene family in citrus fruit (Citrus clementina)

Multidrug and toxic compound extrusion (MATE) proteins are a class of secondary active multidrug transporters. In plants, this family has significantly expanded and is involved in numerous plant physiological processes. Although MATE proteins have been identified in an increasing number of species, the understanding about this family in citrus remains unclear. In this study, a total of 69 MATE transporters were identified in the citrus genome (Citrus clementina) and classified into four groups by phylogenetic analysis. Tandem and segmental duplication events were the main causes of the citrus MATE family expansion. RNA-seq and qRT-PCR analyses were performed during citrus fruit development. The results indicated that CitMATE genes showed specific expression profiles in citrus peels and flesh at different developmental stages. Combined with the variations of flavonoids and citrate levels in citrus fruit, we suggested that CitMATE43 and CitMATE66 may be involved in the transport process of flavonoids and citrate in citrus fruit, respectively. In addition, two flavonoids positive regulators, CitERF32 and CitERF33, both directly bind to and activated the CitMATE43 promoter. Our results provide comprehensive information on citrus MATE genes and valuable understanding for the flavonoids and citrate metabolism in citrus fruit.

Tissue Specificity of Gene Expression Evolves Across Mammal Species

Tissue specificity of gene expression sheds light on the tissue-selective manifestation of hereditary disease despite the same DNA across all tissues. The evolutionary path of such tissue specificity provides essential information about the tissue-specific function of genes and the validity of disease animal models. With recent improvements of the sequencing technology, more and more large-scale transcriptomics studies have been conducted among different species across multiple tissues. In this study, we exploit existing transcriptomics resources of humans, cynomolgus macaques, rats, mice, and dogs across 13 tissues. We find that although tissue specificity of homologous gene expression is largely well conserved across species, a total of 380 genes shift or are in the process of shifting their tissue specificity. The tissue-specificity-shifting genes are less conserved than those preserving their tissue specificity or housekeeping genes. Interestingly, tissue-specificity-shifting genes tend to be less conserved at the third codon positions, likely due to their relaxed synonymous codon usage bias. Moreover, compared with genes, cassette exons are more likely to shift their tissue specificity of splicing across the five species.

Comparative Proteome and Cis-Regulatory Element Analysis Reveals Specific Molecular Pathways Conserved in Dog and Human Brains

Brain development and function are governed by precisely regulated protein expressions in different regions. To date, multiregional brain proteomes have been systematically analyzed only for adult human and mouse brains. To understand the underpinnings of brain development and function, we generated proteomes from six regions of the postnatal brain at three developmental stages of domestic dogs (Canis familiaris), which are special among animals in terms of their remarkable human-like social cognitive abilities. Quantitative analysis of the spatiotemporal proteomes identified region-enriched synapse types at different developmental stages and differential myelination progression in different brain regions. Through integrative analysis of inter-regional expression patterns of orthologous proteins and genome-wide cis-regulatory element frequencies, we found that proteins related with myelination and hippocampus were highly correlated between dog and human but not between mouse and human, although mouse is phylogenetically closer to human. Moreover, the global expression patterns of neurodegenerative disease and autism spectrum disorder-associated proteins in dog brain more resemble human brain than in mouse brain. The high similarity of myelination and hippocampus-related pathways in dog and human at both proteomic and genetic levels may contribute to their shared social cognitive abilities. The inter-regional expression patterns of disease-associated proteins in the brain of different species provide important information to guide mechanistic and translational study using appropriate animal models.

Genomic signatures of the evolution of a diurnal lifestyle in Strigiformes

Understanding the targets of selection associated with changes in behavioral traits represents an important challenge of current evolutionary research. Owls (Strigiformes) are a diverse group of birds, most of which are considered nocturnal raptors. However, a few owl species independently adopted a diurnal lifestyle in their recent evolutionary history. We searched for signals of accelerated rates of evolution associated with a diurnal lifestyle using a genome-wide comparative approach. We estimated substitution rates in coding and noncoding conserved regions of the genome of seven owl species, including three diurnal species. Substitution rates of the noncoding elements were more accelerated than those of protein-coding genes. We identified new, owl-specific conserved noncoding elements as candidates of parallel evolution during the emergence of diurnality in owls. Our results shed light on the molecular basis of adaptation to a new niche and highlight the importance of regulatory elements for evolutionary changes in behavior. These elements were often involved in the neuronal development of the brain.

Speciation and adaptation research meets genome editing

Understanding the genetic basis of reproductive isolation and adaptive traits in natural populations is one of the fundamental goals in evolutionary biology. Genome editing technologies based on CRISPR-Cas systems and site-specific recombinases have enabled us to modify a targeted genomic region as desired and thus to conduct functional analyses of target loci, genes and mutations even in non-conventional model organisms. Here, we review the technical properties of genome editing techniques by classifying them into the following applications: targeted gene knock-out for investigating causative gene functions, targeted gene knock-in of marker genes for visualizing expression patterns and protein functions, precise gene replacement for identifying causative alleles and mutations, and targeted chromosomal rearrangement for investigating the functional roles of chromosomal structural variations. We describe examples of their application to demonstrate functional analysis of naturally occurring genetic variations and discuss how these technologies can be applied to speciation and adaptation research. This article is part of the theme issue 'Genetic basis of adaptation and speciation: from loci to causative mutations'.

A reference single-cell regulomic and transcriptomic map of cynomolgus monkeys

Non-human primates are attractive laboratory animal models that accurately reflect both developmental and pathological features of humans. Here we present a compendium of cell types across multiple organs in cynomolgus monkeys (Macaca fascicularis) using both single-cell chromatin accessibility and RNA sequencing data. The integrated cell map enables in-depth dissection and comparison of molecular dynamics, cell-type compositions and cellular heterogeneity across multiple tissues and organs. Using single-cell transcriptomic data, we infer pseudotime cell trajectories and cell-cell communications to uncover key molecular signatures underlying their cellular processes. Furthermore, we identify various cell-specific cis-regulatory elements and construct organ-specific gene regulatory networks at the single-cell level. Finally, we perform comparative analyses of single-cell landscapes among mouse, monkey and human. We show that cynomolgus monkey has strikingly higher degree of similarities in terms of immune-associated gene expression patterns and cellular communications to human than mouse. Taken together, our study provides a valuable resource for non-human primate cell biology.

Non-human primates are attractive laboratory animal models that can accurately reflect some developmental and pathological features of humans. Here the authors chart a reference cell map of cynomolgus monkeys using both scATAC-seq and scRNA-seq data across multiple organs, providing insights into the molecular dynamics and cellular heterogeneity of this organism.

An enhancer of Agouti contributes to parallel evolution of cryptically colored beach mice

Identifying the genetic basis of repeatedly evolved traits provides a way to reconstruct their evolutionary history and ultimately investigate the predictability of evolution. Here, we focus on the oldfield mouse (Peromyscus polionotus), which occurs in the southeastern United States, where it exhibits considerable color variation. Dorsal coats range from dark brown in mainland mice to near white in mice inhabiting sandy beaches; this light pelage has evolved independently on Florida's Gulf and Atlantic coasts as camouflage from predators. To facilitate genomic analyses, we first generated a chromosome-level genome assembly of Peromyscus polionotus subgriseus. Next, in a uniquely variable mainland population (Peromyscus polionotus albifrons), we scored 23 pigment traits and performed targeted resequencing in 168 mice. We find that pigment variation is strongly associated with an ∼2-kb region ∼5 kb upstream of the Agouti signaling protein coding region. Using a reporter-gene assay, we demonstrate that this regulatory region contains an enhancer that drives expression in the dermis of mouse embryos during the establishment of pigment prepatterns. Moreover, extended tracts of homozygosity in this Agouti region indicate that the light allele experienced recent and strong positive selection. Notably, this same light allele appears fixed in both Gulf and Atlantic coast beach mice, despite these populations being separated by >1,000 km. Together, our results suggest that this identified Agouti enhancer allele has been maintained in mainland populations as standing genetic variation and from there, has spread to and been selected in two independent beach mouse lineages, thereby facilitating their rapid and parallel evolution.

Exploring high-resolution chromatin interaction changes and functional enhancers of myogenic marker genes during myogenic differentiation

Skeletal muscle differentiation (myogenesis) is a complex and highly coordinated biological process regulated by a series of myogenic marker genes. Chromatin interactions between gene’s promoters and their enhancers have an important role in transcriptional control. However, the high-resolution chromatin interactions of myogenic genes and their functional enhancers during myogenesis remain largely unclear. Here, we used circularized chromosome conformation capture coupled with next generation sequencing (4C-seq) to investigate eight myogenic marker genes in C2C12 myoblasts (C2C12-MBs) and C2C12 myotubes (C2C12-MTs). We revealed dynamic chromatin interactions of these marker genes during differentiation and identified 163 and 314 significant interaction sites (SISs) in C2C12-MBs and C2C12-MTs, respectively. The interacting genes of SISs in C2C12-MTs were mainly involved in muscle development, and histone modifications of the SISs changed during differentiation. Through functional genomic screening, we also identified 25 and 41 putative active enhancers in C2C12-MBs and C2C12-MTs, respectively. Using luciferase reporter assays for putative enhancers of Myog and Myh3, we identified eight activating enhancers. Furthermore, dCas9-KRAB epigenome editing and RNA-Seq revealed a role for Myog enhancers in the regulation of Myog expression and myogenic differentiation in the native genomic context. Taken together, this study lays the groundwork for understanding 3D chromatin interaction changes of myogenic genes during myogenesis and provides insights that contribute to our understanding of the role of enhancers in regulating myogenesis.

Evolution of a cis-Acting SNP That Controls Type VI Secretion in Vibrio cholerae

Mutations in regulatory mechanisms that control gene expression contribute to phenotypic diversity and thus facilitate the adaptation of microbes and other organisms to new niches. Comparative genomics can be used to infer rewiring of regulatory architecture based on large effect mutations like loss or acquisition of transcription factors but may be insufficient to identify small changes in noncoding, intergenic DNA sequence of regulatory elements that drive phenotypic divergence. In human-derived Vibrio cholerae, the response to distinct chemical cues triggers production of multiple transcription factors that can regulate the type VI secretion system (T6), a broadly distributed weapon for interbacterial competition. However, to date, the signaling network remains poorly understood because no regulatory element has been identified for the major T6 locus. Here we identify a conserved cis-acting single nucleotide polymorphism (SNP) controlling T6 transcription and activity. Sequence alignment of the T6 regulatory region from diverse V. cholerae strains revealed conservation of the SNP that we rewired to interconvert V. cholerae T6 activity between chitin-inducible and constitutive states. This study supports a model of pathogen evolution through a noncoding cis-regulatory mutation and preexisting, active transcription factors that confers a different fitness advantage to tightly regulated strains inside a human host and unfettered strains adapted to environmental niches. IMPORTANCE Organisms sense external cues with regulatory circuits that trigger the production of transcription factors, which bind specific DNA sequences at promoters ("cis" regulatory elements) to activate target genes. Mutations of transcription factors or their regulatory elements create phenotypic diversity, allowing exploitation of new niches. Waterborne pathogen Vibrio cholerae encodes the type VI secretion system "nanoweapon" to kill competitor cells when activated. Despite identification of several transcription factors, no regulatory element has been identified in the promoter of the major type VI locus, to date. Combining phenotypic, genetic, and genomic analysis of diverse V. cholerae strains, we discovered a single nucleotide polymorphism in the type VI promoter that switches its killing activity between a constitutive state beneficial outside hosts and an inducible state for constraint in a host. Our results support a role for noncoding DNA in adaptation of this pathogen.

Functional Variants Associated With CMPK2 and in ASB16 Influence Bovine Digital Dermatitis

Bovine digital dermatitis (BDD) is an infectious disease of the hoof in cattle with multifactorial etiology and a polygenic influence on susceptibility. With our study, we identified genomic regions with the impact on occurrence and development of BDD. We used 5,040 genotyped animals with phenotype information based on the M-stage system for genome-wide association. Significant associations for single-nucleotide polymorphisms were found near genes CMPK2 (chromosome 11) and ASB16 (chromosome 19) both being implicated in immunological processes. A sequence analysis of the chromosomal regions revealed rs208894039 and rs109521151 polymorphisms as having significant influence on susceptibility to the disease. Specific genotypes were significantly more likely to be affected by BDD and developed chronic lesions. Our study provides an insight into the genomic background for a genetic predisposition related to the pathogenesis of BDD. Results might be implemented in cattle-breeding programs and could pave the way for the establishment of a BDD prescreening test.

Role of CNSs Conserved Distal Cis-Regulatory Elements in CD4 + T Cell Development and Differentiation

Naïve CD4⁺ T cells differentiate into diverse subsets of effector cells and perform various homeostatic and immune functions. The differentiation and maintenance of these different subsets are controlled through the upregulation and silencing of master genes. Mechanistic studies of the regulation of these master genes identified conserved and distal intronic regulatory elements, which are accessible subsets of conserved non-coding sequences (CNSs), acting as cis-regulatory elements in a lineage-specific manner that controls the function of CD4⁺ T cells. Abnormal CNS activity is associated with incorrect expression of master genes and development of autoimmune diseases or immune suppression. Here, we describe the function of several conserved, distal cis-regulatory elements at the Foxp3, Rorc, Il-4, Il-10 and Il-17 gene locus were shown to play important roles in CD4⁺ T cells differentiation. Together, this review briefly outlines currently known CNSs, with a focus on their regulations and functions in complexes modulating the differentiation and maintenance of various CD4⁺ T cells subsets, in health and disease contexts, as well as during the conversion of T regulatory cells to T helper 17 cells. This article will provide a comprehensive view of CNSs conserved distal cis-regulatory elements at a few loci that control aspects of CD4⁺ T cells function.

Evolutionary footprint of plant immunity

There are pieces of evidence from genomic footprints and fossil records indicating that plants have co-evolved with microbes after terrestrialization for more than 407 million years. Therefore, to truly comprehend plant evolution, we need to understand the co-evolutionary process and history between plants and microbes. Recent developments in genomes and transcriptomes of a vast number of plant species as well as microbes have greatly expanded our knowledge of the evolution of the plant immune system. In this review, we summarize recent advances in the co-evolution between plants and microbes with emphasis on the plant side and point out future research needed for understanding plant-microbial co-evolution. Knowledge of the evolution and variation of the plant immune system will better equip us on designing crops with boosted performance in agricultural fields.

The relationship between sexual dimorphism and androgen response element proliferation in primate genomes

In the males of many vertebrate species, sexual selection has led to the evolution of sexually dimorphic traits, which often are developmentally controlled by androgen signaling involving androgen response elements (AREs). Evolutionary changes in the number and genomic locations of AREs can modify patterns of receptor regulation and potentially alter gene expression. Here, we use recently sequenced primate genomes to evaluate the hypothesis that the strength of sexual selection is related to the genome-wide number of AREs in a diversifying lineage. In humans, we find a higher incidence of AREs near male-biased genes and androgen-responsive genes when compared to randomly selected genes from the genome. In a set of primates, we find that gains or losses of AREs proximal to genes are correlated with changes in male expression levels and the degree of sex-biased expression of those genes. In a larger set of primates, we find that increases in indicators of sexual selection are correlated with genome-wide ARE counts. Our results suggest that the responsiveness of the genome to androgens in humans and their close relatives has been shaped by sexual selection that arises from competition among males for mating access to females.

A synthetic synthesis to explore animal evolution and development

Identifying the general principles by which genotypes are converted into phenotypes remains a challenge in the post-genomic era. We still lack a predictive understanding of how genes shape interactions among cells and tissues in response to signalling and environmental cues, and hence how regulatory networks generate the phenotypic variation required for adaptive evolution. Here, we discuss how techniques borrowed from synthetic biology may facilitate a systematic exploration of evolvability across biological scales. Synthetic approaches permit controlled manipulation of both endogenous and fully engineered systems, providing a flexible platform for investigating causal mechanisms in vivo. Combining synthetic approaches with multi-level phenotyping (phenomics) will supply a detailed, quantitative characterization of how internal and external stimuli shape the morphology and behaviour of living organisms. We advocate integrating high-throughput experimental data with mathematical and computational techniques from a variety of disciplines in order to pursue a comprehensive theory of evolution.

This article is part of the theme issue ‘Genetic basis of adaptation and speciation: from loci to causative mutations’.

Genome-wide analysis of cis-regulatory changes underlying metabolic adaptation of cavefish

Cis-regulatory changes are key drivers of adaptative evolution. However, their contribution to the metabolic adaptation of organisms is not well understood. Here, we used a unique vertebrate model, Astyanax mexicanus-different morphotypes of which survive in nutrient-rich surface and nutrient-deprived cave waters-to uncover gene regulatory networks underlying metabolic adaptation. We performed genome-wide epigenetic profiling in the liver tissues of Astyanax and found that many of the identified cis-regulatory elements (CREs) have genetically diverged and have differential chromatin features between surface and cave morphotypes, while retaining remarkably similar regulatory signatures between independently derived cave populations. One such CRE in the hpdb gene harbors a genomic deletion in cavefish that abolishes IRF2 repressor binding and derepresses enhancer activity in reporter assays. Selection of this mutation in multiple independent cave populations supports its importance in cave adaptation, and provides novel molecular insights into the evolutionary trade-off between loss of pigmentation and adaptation to food-deprived caves.

PCRMS: a database of predicted cis-regulatory modules and constituent transcription factor binding sites in genomes

More accurate and more complete predictions of cis-regulatory modules (CRMs) and constituent transcription factor (TF) binding sites (TFBSs) in genomes can facilitate characterizing functions of regulatory sequences. Here, we developed a database predicted cis-regulatory modules (PCRMS) (https://cci-bioinfo.uncc.edu) that stores highly accurate and unprecedentedly complete maps of predicted CRMs and TFBSs in the human and mouse genomes. The web interface allows the user to browse CRMs and TFBSs in an organism, find the closest CRMs to a gene, search CRMs around a gene and find all TFBSs of a TF. PCRMS can be a useful resource for the research community to characterize regulatory genomes. Database URL: https://cci-bioinfo.uncc.edu/.

Targeting a gene regulatory element enhances rice grain yield by decoupling panicle number and size

Crop genetic improvement requires balancing complex tradeoffs caused by gene pleiotropy and linkage drags, as exemplified by IPA1 (Ideal Plant Architecture 1), a typical pleiotropic gene in rice that increases grains per panicle but reduces tillers. In this study, we identified a 54-base pair cis-regulatory region in IPA1 via a tiling-deletion-based CRISPR-Cas9 screen that, when deleted, resolves the tradeoff between grains per panicle and tiller number, leading to substantially enhanced grain yield per plant. Mechanistic studies revealed that the deleted fragment is a target site for the transcription factor An-1 to repress IPA1 expression in panicles and roots. Targeting gene regulatory regions should help dissect tradeoff effects and provide a rich source of targets for breeding complementary beneficial traits.

Transcriptome shock in interspecific F1 allotriploid hybrids between Brassica species

Interspecific hybridization drives the evolution of angiosperms and can be used to introduce novel alleles for important traits or to activate heterosis in crop breeding. Hybridization brings together gene expression networks from two different species, potentially causing global alterations of gene expression in the F1 plants which is called 'transcriptome shock'. Here, we explored such a transcriptome shock in allotriploid Brassica hybrids. We generated interspecific F1 allotriploid hybrids between the allotetraploid species Brassica napus and three accessions of the diploid species Brassica rapa. RNA-seq of the F1 hybrids and the parental plants revealed that 26.34-30.89% of genes were differentially expressed between the parents. We also analyzed expression level dominance and homoeolog expression bias between the parents and the F1 hybrids. The expression-level dominance biases of the Ar, An, and Cn subgenomes was genotype and stage dependent, whereas significant homoeolog expression bias was observed among three subgenomes from different parents. Furthermore, more genes were involved in trans regulation than in cis regulation in allotriploid F1 hybrids. Our findings provide new insights into the transcriptomic responses of cross-species hybrids and hybrids showing heterosis, as well as a new method for promoting the breeding of desirable traits in polyploid Brassica species.

Neuroepigenetic Mechanisms of Action of Ultrashort Peptides in Alzheimer's Disease

Epigenetic regulation of gene expression is necessary for maintaining higher-order cognitive functions (learning and memory). The current understanding of the role of epigenetics in the mechanism of Alzheimer’s disease (AD) is focused on DNA methylation, chromatin remodeling, histone modifications, and regulation of non-coding RNAs. The pathogenetic links of this disease are the misfolding and aggregation of tau protein and amyloid peptides, mitochondrial dysfunction, oxidative stress, impaired energy metabolism, destruction of the blood–brain barrier, and neuroinflammation, all of which lead to impaired synaptic plasticity and memory loss. Ultrashort peptides are promising neuroprotective compounds with a broad spectrum of activity and without reported side effects. The main aim of this review is to analyze the possible epigenetic mechanisms of the neuroprotective action of ultrashort peptides in AD. The review highlights the role of short peptides in the AD pathophysiology. We formulate the hypothesis that peptide regulation of gene expression can be mediated by the interaction of short peptides with histone proteins, cis- and transregulatory DNA elements and effector molecules (DNA/RNA-binding proteins and non-coding RNA). The development of therapeutic agents based on ultrashort peptides may offer a promising addition to the multifunctional treatment of AD.

Multidimensional chromatin profiling of zebrafish pancreas to uncover and investigate disease-relevant enhancers

The pancreas is a central organ for human diseases. Most alleles uncovered by genome-wide association studies of pancreatic dysfunction traits overlap with non-coding sequences of DNA. Many contain epigenetic marks of cis-regulatory elements active in pancreatic cells, suggesting that alterations in these sequences contribute to pancreatic diseases. Animal models greatly help to understand the role of non-coding alterations in disease. However, interspecies identification of equivalent cis-regulatory elements faces fundamental challenges, including lack of sequence conservation. Here we combine epigenetic assays with reporter assays in zebrafish and human pancreatic cells to identify interspecies functionally equivalent cis-regulatory elements, regardless of sequence conservation. Among other potential disease-relevant enhancers, we identify a zebrafish ptf1a distal-enhancer whose deletion causes pancreatic agenesis, a phenotype previously found to be induced by mutations in a distal-enhancer of PTF1A in humans, further supporting the causality of this condition in vivo. This approach helps to uncover interspecies functionally equivalent cis-regulatory elements and their potential role in human disease.

Alterations in cis-regulatory elements (CREs) can contribute to pancreatic diseases. Here the authors combine chromatin profiling and interaction points with in vivo reporter assays in zebrafish to uncover functionally equivalent human CREs, helping to predict disease-relevant enhancers.

Genome-Wide Identification and Expression Analysis of XTH Gene Family during Flower-Opening Stages in Osmanthus fragrans

Osmanthus fragrans is an aromatic plant which is widely used in landscaping and garden greening in China. However, the process of flower opening is significantly affected by ambient temperature changes. Cell expansion in petals is the primary factor responsible for flower opening. Xyloglucan endoglycolase/hydrolase (XTH) is a cell-wall-loosening protein involved in cell expansion or cell-wall weakening. Through whole-genome analysis, 38 OfXTH genes were identified in O. fragrans which belong to the four main phylogenetic groups. The gene structure, chromosomal location, synteny relationship, and cis-acting elements prediction and expression patterns were analyzed on a genome-wide scale. The expression patterns showed that most OfXTHs were closely associated with the flower-opening period of O. fragrans. At the early flower-opening stage (S1 and S2), transcriptome and qRT-PCR analysis revealed the expression of OfXTH24, 27, 32, 35, and 36 significantly increased under low ambient temperature (19 °C). It is speculated that the five genes might be involved in the regulation of flower opening by responding to ambient temperature changes. Our results provide solid foundation for the functional analysis of OfXTH genes and help to explore the mechanism of flower opening responding to ambient temperature in O. fragrans.

Transcriptional states and chromatin accessibility during bovine myoblasts proliferation and myogenic differentiation

Objectives

Although major advances have been made in bovine epigenome study, the epigenetic basis for fetal skeletal muscle development still remains poorly understood. The aim is to recapitulated the time course of fetal skeletal muscle development in vitro, and explore the dynamic changes of chromatin accessibility and gene expression during bovine myoblasts proliferation and differentiation.

Methods

PDGFR‐ cells were isolated from bovine fetal skeletal muscle, then cultured and induced myogenic differentiation in vitro in a time‐course study (P, D0, D2,and D4). The assay for transposase‐accessible chromatin sequencing (ATAC‐seq) and RNA sequencing (RNA‐seq) were performed.

Results

Among the enriched transcriptional factors with high variability, we determined the effects of MAFF, ZNF384, and KLF6 in myogenesis using RNA interference (RNAi). In addition, we identified both stage‐specific genes and chromatin accessibility regions to reveal the sequential order of gene expression, transcriptional regulatory, and signal pathways involved in bovine skeletal muscle development. Further investigation integrating chromatin accessibility and transcriptome data was conducted to explore cis‐regulatory regions in line with gene expression. Moreover, we combined bovine GWAS results of growth traits with regulatory regions defined by chromatin accessibility, providing a suggestive means for a more precise annotation of genetic variants of bovine growth traits.

Conclusion

Overall, these findings provide valuable information for understanding the stepwise regulatory mechanisms in skeletal muscle development and conducting beef cattle genetic improvement programs.

Genome Editing and Designer Crops for the Future

Domestication spanning over thousands of years led to the evolution of crops that are being cultivated in recent times. Later, selective breeding methods were practiced by human to produce improved cultivars/germplasm. Classical breeding was further transformed into molecular- and genomics-assisted breeding strategies, however, these approaches are labor-intensive and time-consuming. The advent of omics technologies has facilitated the identification of genes and genetic determinants that regulate particular traits allowing the direct manipulation of target genes and genomic regions to achieve desirable phenotype. Recently, genome editing technologies such as meganucleases (MN), zinc-finger nucleases (ZFNs), transcription activator-like effector nucleases (TALENs), and CRISPR (clustered regularly interspaced short palindromic repeats)/CRISPR-Associated protein 9 (Cas9) have gained popularity for precise editing of genes to develop crop varieties with superior agronomic, physiological, climate-resilient, and nutritional traits. Owing to the efficiency and precision, genome editing approaches have been widely used to design the crops that can survive the challenges posed by changing climate, and also cater the food and nutritional requirements for ever-growing population. Here, we briefly review different genome editing technologies deployed for crop improvement, and the fundamental differences between GE technology and transgene-based approach. We also summarize the recent advances in genome editing and how this radical expansion can complement the previously established technologies along with breeding for creating designer crops.

Genomic Survey and Cold-Induced Expression Patterns of bHLH Transcription Factors in Liriodendron chinense (Hemsl) Sarg

bHLH transcription factors play an animated role in the plant kingdom during growth and development, and responses to various abiotic stress. In this current study, we conducted, the genome-wide survey of bHLH transcription factors in Liriodendron chinense (Hemsl) Sarg., 91 LcbHLH family members were identified. Identified LcbHLH gene family members were grouped into 19 different subfamilies based on the conserved motifs and phylogenetic analysis. Our results showed that LcbHLH genes clustered in the same subfamily exhibited a similar conservative exon-intron pattern. Hydrophilicity value analysis showed that all LcbHLH proteins were hydrophilic. The Molecular weight (Mw) of LcbHLH proteins ranged from 10.19 kD (LcbHLH15) to 88.40 kD (LcbHLH50). A greater proportion, ~63%, of LcbHLH proteins had a theoretical isoelectric point (pI) less than seven. Additional analysis on the collinear relationships within species and among dissimilar species illustrated that tandem and fragment duplication are the foremost factors of amplification of this family in the evolution process, and they are all purified and selected. RNA-seq and real-time quantitative PCR analysis of LcbHLH members showed that the expression of LcbHLH35, 55, and 86 are up-regulated, and the expression of LcbHLH9, 20, 39, 54, 56, and 69 is down-regulated during cold stress treatments while the expression of LcbHLH24 was up-regulated in the short term and then later down-regulated. From our results, we concluded that LcbHLH genes might participate in cold-responsive processes of L. chinense. These findings provide the basic information of bHLH gene in L. chinense and their regulatory roles in plant development and cold stress response.

Rapid Identification of Pollen- and Anther-Specific Genes in Response to High-Temperature Stress Based on Transcriptome Profiling Analysis in Cotton

Anther indehiscence and pollen sterility caused by high temperature (HT) stress have become a major problem that decreases the yield of cotton. Pollen- and anther-specific genes play a critical role in the process of male reproduction and the response to HT stress. In order to identify pollen-specific genes that respond to HT stress, a comparative transcriptome profiling analysis was performed in the pollen and anthers of Gossypium hirsutum HT-sensitive Line H05 against other tissue types under normal temperature (NT) conditions, and the analysis of a differentially expressed gene was conducted in the pollen of H05 under NT and HT conditions. In total, we identified 1111 pollen-specific genes (PSGs), 1066 anther-specific genes (ASGs), and 833 pollen differentially expressed genes (DEGs). Moreover, we found that the late stage of anther included more anther- and pollen-specific genes (APSGs). Stress-related cis-regulatory elements (CREs) and hormone-responsive CREs are enriched in the promoters of APSGs, suggesting that APSGs may respond to HT stress. However, 833 pollen DEGs had only 10 common genes with 1111 PSGs, indicating that PSGs are mainly involved in the processes of pollen development and do not respond to HT stress. Promoters of these 10 common genes are enriched for stress-related CREs and MeJA-responsive CREs, suggesting that these 10 common genes are involved in the process of pollen development while responding to HT stress. This study provides a pathway for rapidly identifying cotton pollen-specific genes that respond to HT stress.

Nonadditive gene expression is correlated with nonadditive phenotypic expression in interspecific triploid hybrids of willow (Salix spp.)

Many studies have highlighted the complex and diverse basis for heterosis in inbred crops. Despite the lack of a consensus model, it is vital that we turn our attention to understanding heterosis in undomesticated, heterozygous, and polyploid species, such as willow (Salix spp.). Shrub willow is a dedicated energy crop bred to be fast-growing and high yielding on marginal land without competing with food crops. A trend in willow breeding is the consistent pattern of heterosis in triploids produced from crosses between diploid and tetraploid species. Here, we test whether differentially expressed genes are associated with heterosis in triploid families derived from diploid Salix purpurea, diploid Salix viminalis, and tetraploid Salix miyabeana parents. Three biological replicates of shoot tips from all family progeny and parents were collected after 12 weeks in the greenhouse and RNA extracted for RNA-Seq analysis. This study provides evidence that nonadditive patterns of gene expression are correlated with nonadditive phenotypic expression in interspecific triploid hybrids of willow. Expression-level dominance was most correlated with heterosis for biomass yield traits and was highly enriched for processes involved in starch and sucrose metabolism. In addition, there was a global dosage effect of parent alleles in triploid hybrids, with expression proportional to copy number variation. Importantly, differentially expressed genes between family parents were most predictive of heterosis for both field and greenhouse collected traits. Altogether, these data will be used to progress models of heterosis to complement the growing genomic resources available for the improvement of heterozygous perennial bioenergy crops.

Exploration of glutathione reductase for abiotic stress response in bread wheat (Triticum aestivum L.)

A total of seven glutathione reductase (GR) genes were identified in Triticum aestivum, which were used for comparative structural characterization, phylogenetic analysis and expression profiling with the GR genes of other cereal plants. The modulated gene expression and enzyme activity revealed the role of GRs in abiotic stress response in T. aestivum. Glutathione reductase (GR) is an enzymatic antioxidant that converts oxidized glutathione (GSSG) into reduced glutathione (GSH) through the ascorbate-glutathione cycle. In this study, a total of seven GR genes forming two homeologous groups were identified in the allohexaploid genome of bread wheat (Triticum aestivum). Besides, we identified three GR genes in each Aegilops tauschii, Brachypodium distachyon, Triticum urartu and Sorghum bicolor, which were used for comparative characterization. Phylogenetic analysis revealed the clustering of GR proteins into two groups; class I and class II, which were predicted to be localized in cytoplasm and chloroplast, respectively. The exon-intron and conserved motif patterns were almost conserved in each group, in which a maximum of 10 and 17 exons were present in chloroplastic and cytoplasmic GRs, respectively. The protein structure analysis confirmed the occurrence of conserved pyridine nucleotide disulfide oxidoreductase (Pyr_redox) and pyridine nucleotide disulfide oxidoreductase dimerization (Pyr_redox_dim) domains in each GR. The active site of GR proteins consisted of two conserved cysteine residues separated by four amino acid residues. Promoter analysis revealed the occurrence of growth and stress-related cis-active elements. Tissue-specific expression profiling suggested the involvement of GRs in both vegetative and reproductive tissue development in various plants. The differential expression of TaGR genes and enhanced GR enzyme activity suggested their roles under drought, heat, salt and arsenic stress. Interaction of GRs with other proteins and chemical compounds of the ascorbate-glutathione cycle revealed their coordinated functioning. The current study will provide a foundation for the validation of the precise role of each GR gene in future studies.

In vivo CRISPR-Cas9-mediated DNA chop identifies a cochlear outer hair cell-specific enhancer

Cochlear outer hair cells (OHCs) are essential for hearing. A short, OHC-specific enhancer is necessary but not yet available for gene therapeutic applications in OHC damage. Such damage is a major cause of deafness. Prestin is a motor protein exclusively expressed in OHCs. We hypothesized that the cis-regulatory DNA fragment deletion of Slc26a5 would affect its expression. We tested this hypothesis by conducting CRISPR/Cas9-mediated large DNA fragment deletion of mouse Slc26a5 intron regions. First, starting from a ~13 kbp fragment, step-by-step, we narrowed down the sequence to a 1.4 kbp segment. By deleting either a 13 kbp or 1.4 kbp fragment, we observed delayed Prestin expression. Second, we showed that 1.4 kbp was an OHC-specific enhancer because enhanced green fluorescent protein (EGFP) was highly and specifically expressed in OHCs in a transgenic mouse where EGFP was driven by the 1.4 kbp segment. More importantly, specific EGFP was also driven by its homologous 398 bp fragment in human Slc26a5. This suggests that the enhancer is likely to be evolutionarily conserved across different species.

The evolution, evolvability and engineering of gene regulatory DNA

Mutations in non-coding regulatory DNA sequences can alter gene expression, organismal phenotype and fitness^1-3. Constructing complete fitness landscapes, in which DNA sequences are mapped to fitness, is a long-standing goal in biology, but has remained elusive because it is challenging to generalize reliably to vast sequence spaces^4-6. Here we build sequence-to-expression models that capture fitness landscapes and use them to decipher principles of regulatory evolution. Using millions of randomly sampled promoter DNA sequences and their measured expression levels in the yeast Saccharomyces cerevisiae, we learn deep neural network models that generalize with excellent prediction performance, and enable sequence design for expression engineering. Using our models, we study expression divergence under genetic drift and strong-selection weak-mutation regimes to find that regulatory evolution is rapid and subject to diminishing returns epistasis; that conflicting expression objectives in different environments constrain expression adaptation; and that stabilizing selection on gene expression leads to the moderation of regulatory complexity. We present an approach for using such models to detect signatures of selection on expression from natural variation in regulatory sequences and use it to discover an instance of convergent regulatory evolution. We assess mutational robustness, finding that regulatory mutation effect sizes follow a power law, characterize regulatory evolvability, visualize promoter fitness landscapes, discover evolvability archetypes and illustrate the mutational robustness of natural regulatory sequence populations. Our work provides a general framework for designing regulatory sequences and addressing fundamental questions in regulatory evolution.

MAP4K4 controlled transcription factor POUM1 regulates PxABCG1 expression influencing Cry1Ac resistance in Plutella xylostella (L.)

Deciphering the molecular mechanisms of insect resistance to Bacillus thuringiensis (Bt) based biotechnology products including Bt sprays and Bt crops is critical for the long-term application of Bt technology. Previously, we established that down-regulation of the ABC transporter gene PxABCG1, trans-regulated by the MAPK signaling pathway, contributed to high-level resistance to Bt Cry1Ac toxin in diamondback moth, Plutella xylostella (L.). However, the underlying transcriptional regulatory mechanism was unknown. Herein, we identified putative binding sites (PBSs) of the transcription factor (TF) POUM1 in the PxABCG1 promoter and used a dual-luciferase reporter assay (DLRA) and yeast one-hybrid (Y1H) assay to reveal that POUM1 activates PxABCG1 via interaction with one of these sites. The expression of POUM1 was significantly decreased in the midgut tissue of Cry1Ac-resistant P. xylostella strains compared to a Cry1Ac-susceptible P. xylostella strain. Silencing of POUM1 expression resulted in reduced expression of the PxABCG1 gene and an increase in larval tolerance to Bt Cry1Ac toxin in the Cry1Ac-susceptible P. xylostella strain. Furthermore, silencing of PxMAP4K4 expression increased the expression of both POUM1 and PxABCG1 genes in the Cry1Ac-resistant P. xylostella strain. These results indicate that the POUM1 induces PxABCG1 expression, while the activated MAPK cascade represses PxABCG1 expression thus reducing Cry1Ac susceptibility in P. xylostella. This result deepens our understanding of the transcriptional regulatory mechanism of midgut Cry receptor genes and the molecular basis of the evolution of Bt resistance in insects.

Pangenomics in crop improvement-from coding structural variations to finding regulatory variants with pangenome graphs

Since the first reported crop pangenome in 2014, advances in high-throughput and cost-effective DNA sequencing technologies facilitated multiple such studies including the pangenomes of oilseed rape (Brassica napus L.), soybean [Glycine max (L.) Merr.], rice (Oryza sativa L.), wheat (Triticum aestivum L.), and barley (Hordeum vulgare L.). Compared with single-reference genomes, pangenomes provide a more accurate representation of the genetic variation present in a species. By combining the genomic data of multiple accessions, pangenomes allow for the detection and annotation of complex DNA polymorphisms such as structural variations (SVs), one of the major determinants of genetic diversity within a species. In this review we summarize the current literature on crop pangenomics, focusing on their application to find candidate SVs involved in traits of agronomic interest. We then highlight the potential of pangenomes in the discovery and functional characterization of noncoding regulatory sequences and their variations. We conclude with a summary and outlook on innovative data structures representing the complete content of plant pangenomes including annotations of coding and noncoding elements and outcomes of transcriptomic and epigenomic experiments.

The power of the (imperfect) palindrome: Sequence-specific roles of palindromic motifs in gene regulation

In human languages, a palindrome reads the same forward as backward (e.g., 'madam'). In regulatory DNA, a palindrome is an inverted sequence repeat that allows a transcription factor to bind as a homodimer or as a heterodimer with another type of transcription factor. Regulatory palindromes are typically imperfect, that is, the repeated sequences differ in at least one base pair, but the functional significance of this asymmetry remains poorly understood. Here, we review the use of imperfect palindromes in Drosophila photoreceptor differentiation and mammalian steroid receptor signaling. Moreover, we discuss mechanistic explanations for the predominance of imperfect palindromes over perfect palindromes in these two gene regulatory contexts. Lastly, we propose to elucidate whether specific imperfectly palindromic variants have specific regulatory functions in steroid receptor signaling and whether such variants can help predict transcriptional outcomes as well as the response of individual patients to drug treatments.

Gene losses may contribute to subterranean adaptations in naked mole-rat and blind mole-rat

Background

Naked mole-rats (Heterocephalus glaber, NMRs) and blind mole-rats (Spalax galili, BMRs) are representative subterranean rodents that have evolved many extraordinary traits, including hypoxia tolerance, longevity, and cancer resistance. Although multiple candidate loci responsible for these traits have been uncovered by genomic studies, many of them are limited to functional changes to amino acid sequence and little is known about the contributions of other genetic events. To address this issue, we focused on gene losses (unitary pseudogenes) and systematically analyzed gene losses in NMRs and BMRs, aiming to elucidate the potential roles of pseudogenes in their adaptation to subterranean lifestyle.

Results

We obtained the pseudogene repertoires in NMRs and BMRs, as well as their respective aboveground relatives, guinea pigs and rats, on a genome-wide scale. As a result, 167, 139, 341, and 112 pseudogenes were identified in NMRs, BMRs, guinea pigs, and rats, respectively. Functional enrichment analysis identified 4 shared and 2 species-specific enriched functional groups (EFGs) in subterranean lineages. Notably, the pseudogenes in these EFGs might be associated with either regressive (e.g., visual system) or adaptive (e.g., altered DNA damage response) traits. In addition, several pseudogenes including TNNI3K and PDE5A might be associated with specific cardiac features observed in subterranean lineages. Interestingly, we observed 20 convergent gene losses in NMRs and BMRs. Given that the functional investigations of these genes are generally scarce, we provided functional evidence that independent loss of TRIM17 in NMRs and BMRs might be beneficial for neuronal survival under hypoxia, supporting the positive role of eliminating TRIM17 function in hypoxia adaptation. Our results also suggested that pseudogenes, together with positively selected genes, reinforced subterranean adaptations cooperatively.

Conclusions

Our study provides new insights into the molecular underpinnings of subterranean adaptations and highlights the importance of gene losses in mammalian evolution.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01243-0.

Transcriptomic Variations and Network Hubs Controlling Seed Size and Weight During Maize Seed Development

To elucidate the mechanisms underlying seed development in maize, comprehensive RNA-seq analyses were conducted on Zhengdan1002 (ZD1002), Zhengdan958 (ZD958), and their parental lines during seven seed developmental stages. We found that gene expression levels were largely nonadditive in hybrids and that cis-only or trans × cis pattern played a large role in hybrid gene regulation during seed developmental stage. Weighted gene co-expression network (WGCNA) analysis showed that 36 modules were highly correlated (r = -0.90-0.92, p < 0.05) with kernel weight, length, and width during seed development. Forty-five transcription factors and 38 ribosomal protein genes were identified as major hub genes determining seed size/weight. We also described a network hub, Auxin Response Factor 12 of maize (ZmARF12), a member of a family of transcription factor that mediate gene expression in response to auxin, potentially links auxin signal pathways, cell division, and the size of the seeds. The ZmARF12 mutant exhibited larger seed size and higher grain weight. ZmARF12 transcription was negatively associated with cell division during seed development, which was confirmed by evaluating the yield of protoplasts that isolated from the kernels of the mutant and other inbred lines. Transient knock-down of ZmARF12 in maize plants facilitated cell expansion and division, whereas transient silencing of its potential interactor ZmIAA8 impaired cell division. ZmIAA8 expression was repressed in the ZmARF12 over-expressed protoplasts. The mutant phenotype and the genetics studies presented here illustrated evidence that ZmARF12 is a cell division repressor, and potentially determines the final seed size.

Variation in cis-regulation of a NAC transcription factor contributes to drought tolerance in wheat

Drought is a major environmental factor limiting wheat production worldwide, and developing drought-tolerant cultivars is a central challenge for wheat breeders globally. Therefore, it is important to identify genetic components determining drought tolerance in wheat. In this study, we identified a wheat NAC gene (TaNAC071-A) that is tightly associated with drought tolerance by a genome-wide association study. Knockdown of TaNAC071-A in wheat attenuated plant drought tolerance, whereas its overexpression significantly enhanced drought tolerance through improved water-use efficiency and increased expression of stress-responsive genes. This heightened water-saving mechanism mitigated the yield loss caused by water deficit. Further candidate gene association analysis showed that a 108-bp insertion in the promoter of TaNAC071-A alters its expression level and contributes to variation in drought tolerance among wheat accessions. This insertion contains two MYB cis-regulatory elements (CREs) that can be directly bound by the MYB transcription activator, TaMYBL1, thereby leading to increased TaNAC071-A expression and plant drought tolerance. Importantly, introgression of this 108-bp insertion allele, TaNAC071-A^In-693, into drought-sensitive cultivars could improve their drought tolerance, demonstrating that it is a valuable genetic resource for wheat breeding. Taken together, our findings highlight a major breakthrough in determining the genetic basis underlying phenotypic variation in wheat drought tolerance and showcase the potential of exploiting CRE-containing indels for improving important agronomical traits.

Mouse genomic and cellular annotations

Mice have emerged as one of the most popular and valuable model organisms in the research of human biology. This is due to their genetic and physiological similarity to humans, short generation times, availability of genetically homologous inbred strains, and relatively easy laboratory maintenance. Therefore, following the release of the initial human reference genome, the generation of the mouse reference genome was prioritised and represented an important scientific resource for the mouse genetics community. In 2002, the Mouse Genome Sequencing Consortium published an initial draft of the mouse reference genome which contained ~ 96% of the euchromatic genome of female C57BL/6 J mice. Almost two decades on from the publication of the initial draft, sequencing efforts have continued to increase the completeness and accuracy of the C57BL/6 J reference genome alongside advances in genome annotation. Additionally new sequencing technologies have provided a wealth of data that has added to the repertoire of annotations associated with traditional genomic annotations. Including but not limited to advances in regulatory elements, the 3D genome and individual cellular states. In this review we focus on the reference genome C57BL/6 J and summarise the different aspects of genomic and cellular annotations, as well as their relevance to mouse genetic research. We denote a genomic annotation as a functional unit of the genome. Cellular annotations are annotations of cell type or state, defined by the transcriptomic expression profile of a cell. Due to the wide-ranging number and diversity of annotations describing the mouse genome, we focus on gene, repeat and regulatory element annotation as well as two relatively new technologies; 3D genome architecture and single-cell sequencing outlining their utility in genetic research and their current challenges.

Modeling temporal and hormonal regulation of plant transcriptional response to wounding

Plants respond to wounding stress by changing gene expression patterns and inducing the production of hormones including jasmonic acid. This wounding transcriptional response activates specialized metabolism pathways such as the glucosinolate pathways in Arabidopsis thaliana. While the regulatory factors and sequences controlling a subset of wound-response genes are known, it remains unclear how wound response is regulated globally. Here, we how these responses are regulated by incorporating putative cis-regulatory elements, known transcription factor binding sites, in vitro DNA affinity purification sequencing, and DNase I hypersensitive sites to predict genes with different wound-response patterns using machine learning. We observed that regulatory sites and regions of open chromatin differed between genes upregulated at early and late wounding time-points as well as between genes induced by jasmonic acid and those not induced. Expanding on what we currently know, we identified cis-elements that improved model predictions of expression clusters over known binding sites. Using a combination of genome editing, in vitro DNA-binding assays, and transient expression assays using native and mutated cis-regulatory elements, we experimentally validated four of the predicted elements, three of which were not previously known to function in wound-response regulation. Our study provides a global model predictive of wound response and identifies new regulatory sequences important for wounding without requiring prior knowledge of the transcriptional regulators.