BACKGROUND: Winter wheat undergoes vernalization, a process activated by prolonged exposure to low temperatures. During this phase, flowering signals are generated and transported to the apical meristems, stimulating the transition to the inflorescence meristem while inhibiting tiller bud elongation. Although some vernalization genes have been identified, the key cis-regulatory elements and precise mechanisms governing this process in wheat remain largely unknown.
RESULTS: In this study, we construct extensive epigenomic and transcriptomic profiling across multiple tissues-leaf, axillary bud, and shoot apex-during the vernalization of winter wheat. Epigenetic modifications play a crucial role in eliciting tissue-specific responses and sub-genome-divergent expressions during vernalization. Notably, we observe that H3K27me3 primarily regulates vernalization-induced genes and has limited influence on vernalization-repressed genes. The integration of these datasets enables the identification of 10,600 putative vernalization-related regulatory elements including distal accessible chromatin regions (ACRs) situated 30Kb upstream of VRN3, contributing to the construction of a comprehensive regulatory network. Furthermore, we discover that TaSPL7/15, integral components of the aging-related flowering pathway, interact with the VRN1 promoter and VRN3 distal regulatory elements. These interactions finely regulate their expressions, consequently impacting the vernalization process and flowering.
CONCLUSIONS: Our study offers critical insights into wheat vernalization\'s epigenomic dynamics and identifies the putative regulatory elements crucial for developing wheat germplasm with varied vernalization characteristics. It also establishes a vernalization-related transcriptional network, and uncovers that TaSPL7/15 from the aging pathway participates in vernalization by directly binding to the VRN1 promoter and VRN3 distal regulatory elements.

Malaria transmission and endemicity in Africa remains hugely disproportionate compared to the rest of the world. The complex life cycle of P. falciparum (Pf) between the vertebrate human host and the anopheline vector results in differential expression of genes within and between hosts. An in-depth understanding of Pf interaction with various human genes through regulatory elements will pave way for identification of newer tools in the arsenal for malaria control. Therefore, the regulatory elements (REs) involved in the over- or under-expression of various host immune genes hold the key to elucidating alternative control measures that can be applied for disease surveillance, prompt diagnosis and treatment. We carried out an RNAseq analysis to identify differentially expressed genes and network elucidation of non-coding RNAs and target genes associated with immune response in individuals with different clinical outcomes. Raw RNAseq datasets, retrieved for analyses include individuals with severe (Gambia-20), symptomatic (Burkina Faso-15), asymptomatic (Mali-16) malaria as well as uninfected controls (Tanzania-20; Mali-36). Of the total 107 datasets retrieved, we identified 5534 differentially expressed genes (DEGs) among disease and control groups. A peculiar pattern of DEGs was observed, with individuals presenting with severe/symptomatic malaria having the highest and most diverse upregulated genes, while a reverse phenomenon was recorded among asymptomatic and uninfected individuals. In addition, we identified 141 differentially expressed micro RNA (miRNA), of which 78 and 63 were upregulated and downregulated respectively. Interactome analysis revealed a moderate interaction between DEGs and miRNAs. Of all identified miRNA, five were unique (hsa-mir-32, hsa-mir-25, hsa-mir-221, hsa-mir-29 and hsa-mir-148) because of their connectivity to several genes, including hsa-mir-221 connected to 16 genes. Six-hundred and eight differentially expressed long non coding RNA (lncRNA) were also identified, including SLC7A11, LINC01524 among the upregulated ones. Our study provides important insight into host immune genes undergoing differential expression under different malaria conditions. It also identified unique miRNAs and lncRNAs that modify and/or regulate the expression of various immune genes. These regulatory elements we surmise, have the potential to serve a diagnostic purpose in discriminating between individuals with severe/symptomatic malaria and those with asymptomatic infection or uninfected, following further clinical validation from field isolates.

The evolution of gene expression responses are a critical component of adaptation to variable environments. Predicting how DNA sequence influences expression is challenging because the genotype to phenotype map is not well resolved for cis regulatory elements, transcription factor binding, regulatory interactions, and epigenetic features, not to mention how these factors respond to environment. We tested if flexible machine learning models could learn some of the underlying cis-regulatory genotype to phenotype map. We tested this approach using cold-responsive transcriptome profiles in 5 diverse Arabidopsis thaliana accessions. We first tested for evidence that cis regulation plays a role in environmental response, finding 14 and 15 motifs that were significantly enriched within the up- and down-stream regions of cold-responsive differentially regulated genes (DEGs). We next applied convolutional neural networks (CNNs), which learn de novo cis-regulatory motifs in DNA sequences to predict expression response to environment. We found that CNNs predicted differential expression with moderate accuracy, with evidence that predictions were hindered by biological complexity of regulation and the large potential regulatory code. Overall, DEGs between specific environments can be predicted based on variation in cis-regulatory sequences, although more information needs to be incorporated and better models may be required.

The CRISPR/Cas9 technique applied to modify the cattle genome has value in increasing animal health and welfare. Here, we established a simple, fast, and efficient cloning-free CRISPR/Cas9 protocol for large deletions of genomic loci in the frequently used model bovine MDBK cell line. The main advantages of our protocol are as follows: (i) pre-screening of the sgRNA efficiency with a fast and simple cleavage assay, (ii) reliable detection of genomic edits primarily by PCR and confirmed by DNA sequencing, and (iii) single cell sorting with FACS providing specific genetic information from modified cells of interest. Therefore, our method could be successfully applied in different studies, including functional validation of any genetic or regulatory elements.

Genetic predisposition to cardiac arrhythmias has been a field of intense investigation. Research initially focused on rare hereditary arrhythmias, but over the last two decades, the role of genetic variation (single nucleotide polymorphisms) in heart rate, rhythm, and arrhythmias has been taken into consideration as well. In particular, genome-wide association studies have identified hundreds of genomic loci associated with quantitative electrocardiographic traits, atrial fibrillation, and less common arrhythmias such as Brugada syndrome. A significant number of associated variants have been found to systematically localize in non-coding regulatory elements that control the tissue-specific and temporal transcription of genes encoding transcription factors, ion channels, and other proteins. However, the identification of causal variants and the mechanism underlying their impact on phenotype has proven difficult due to the complex tissue-specific, time-resolved, condition-dependent, and combinatorial function of regulatory elements, as well as their modest conservation across different model species. In this review, we discuss research efforts aimed at identifying and characterizing-trait-associated variant regulatory elements and the molecular mechanisms underlying their impact on heart rate or rhythm.

Recombination is responsible for breaking up haplotypes, influencing genetic variability, and the efficacy of selection. Bird genomes lack the protein PR domain-containing protein 9, a key determinant of recombination dynamics in most metazoans. Historical recombination maps in birds show an apparent stasis in positioning recombination events. This highly conserved recombination pattern over long timescales may constrain the evolution of recombination in birds. At the same time, extensive variation in recombination rate is observed across the genome and between different species of birds. Here, we characterize the fine-scale historical recombination map of an iconic migratory songbird, the Eurasian blackcap (Sylvia atricapilla), using a linkage disequilibrium-based approach that accounts for population demography. Our results reveal variable recombination rates among and within chromosomes, which associate positively with nucleotide diversity and GC content and negatively with chromosome size. Recombination rates increased significantly at regulatory regions but not necessarily at gene bodies. CpG islands are associated strongly with recombination rates, though their specific position and local DNA methylation patterns likely influence this relationship. The association with retrotransposons varied according to specific family and location. Our results also provide evidence of heterogeneous intrachromosomal conservation of recombination maps between the blackcap and its closest sister taxon, the garden warbler. These findings highlight the considerable variability of recombination rates at different scales and the role of specific genomic features in shaping this variation. This study opens the possibility of further investigating the impact of recombination on specific population-genomic features.

The majority of disease-associated variants identified through genome-wide association studies are located outside of protein-coding regions. Prioritizing candidate regulatory variants and gene targets to identify potential biological mechanisms for further functional experiments can be challenging. To address this challenge, we developed FORGEdb ( https://forgedb.cancer.gov/ ; https://forge2.altiusinstitute.org/files/forgedb.html ; and https://doi.org/10.5281/zenodo.10067458 ), a standalone and web-based tool that integrates multiple datasets, delivering information on associated regulatory elements, transcription factor binding sites, and target genes for over 37 million variants. FORGEdb scores provide researchers with a quantitative assessment of the relative importance of each variant for targeted functional experiments.

DNA regulatory elements intricately control when, where, and how genes are activated. Therefore, understanding the function of these elements could unveil the complexity of the genetic regulation network. Genome-wide significant variants are predominantly found in non-coding regions of DNA, so comprehending the predicted functional regulatory elements is crucial for understanding the biological context of these genomic markers, which can be incorporated into breeding programs. The emergence of CRISPR technology has provided a powerful tool for studying non-coding regulatory elements in genomes. In this study, we leveraged epigenetic data from the Functional Annotation of Animal Genomes project to identify promoter and putative enhancer regions associated with three genes (HBBA, IRF7, and PPARG) in the chicken genome. To identify the enhancer regions, we designed guide RNAs targeting the promoter and candidate enhancer regions and utilized CRISPR activation (CRISPRa) with dCas9-p300 and dCas9-VPR as transcriptional activators in chicken DF-1 cells. By comparing the expression levels of target genes between the promoter activation and the co-activation of the promoter and putative enhancers, we were able to identify functional enhancers that exhibited augmented upregulation. In conclusion, our findings demonstrate the remarkable efficiency of CRISPRa in precisely manipulating the expression of endogenous genes by targeting regulatory elements in the chicken genome, highlighting its potential for functional validation of non-coding regions.

Cellular identity is determined partly by cell type-specific epigenomic profiles that regulate gene expression. In neuroscience, there is a pressing need to isolate and characterize the epigenomes of specific CNS cell types in health and disease. In this study, we developed an in vivo tagging mouse model (Camk2a-NuTRAP) for paired isolation of neuronal DNA and RNA without cell sorting and then used this model to assess epigenomic regulation, DNA modifications in particular, of gene expression between neurons and glia.
After validating the cell-specificity of the Camk2a-NuTRAP model, we performed TRAP-RNA-Seq and INTACT-whole genome oxidative bisulfite sequencing (WGoxBS) to assess the neuronal translatome and epigenome in the hippocampus of young mice (4 months old). WGoxBS findings were validated with enzymatic methyl-Seq (EM-Seq) and nanopore sequencing. Comparing neuronal data to microglial and astrocytic data from NuTRAP models, microglia had the highest global mCG levels followed by astrocytes and then neurons, with the opposite pattern observed for hmCG and mCH. Differentially modified regions between cell types were predominantly found within gene bodies and distal intergenic regions, rather than proximal promoters. Across cell types there was a negative correlation between DNA modifications (mCG, mCH, hmCG) and gene expression at proximal promoters. In contrast, a negative correlation of gene body mCG and a positive relationship between distal promoter and gene body hmCG with gene expression was observed. Furthermore, we identified a neuron-specific inverse relationship between mCH and gene expression across promoter and gene body regions.
Neurons, astrocytes, and microglia demonstrate different genome-wide levels of mCG, hmCG, and mCH that are reproducible across analytical methods. However, modification-gene expression relationships are conserved across cell types. Enrichment of differential modifications across cell types in gene bodies and distal regulatory elements, but not proximal promoters, highlights epigenomic patterning in these regions as potentially greater determinants of cell identity. These findings also demonstrate the importance of differentiating between mC and hmC in neuroepigenomic analyses, as up to 30% of what is conventionally interpreted as mCG can be hmCG, which often has a different relationship to gene expression than mCG.

BACKGROUND: To identify the disease-causing gene in a Chinese family affected with congenital aniridia.
METHODS: Patients underwent systematic ophthalmic examinations such as anterior segment photography, fundus photography, optical coherence tomography, and fundus fluorescein angiography. The proband was screened for pathogenic variants by whole exome sequencing (WES) and copy number variant (CNV) analysis. Real-time quantitative PCR (RT-qPCR) was applied to confirm the CNV results. Breakpoints were identified by long-range PCR followed by Sanger sequencing.
RESULTS: All seven members of this Chinese family, including four patients and three normal individuals, were recruited for this study. All patients showed bilateral congenital aniridia with nystagmus, except the son of the proband, who presented with bilateral partial coloboma of the iris. A novel heterozygous deletion (chr11:31,139,019-31,655,997) containing the 3\' regulatory enhancers of the PAX6 gene was detected in this family. We also reviewed the reported microdeletions downstream of PAX6 in patients with aniridia.
CONCLUSIONS: We identified a novel microdeletion, 517 kb in size located about 133 kb downstream of the PAX6 gene, responsible for congenital aniridia in this Chinese family, which expands the spectrum of aniridia-associated mutations in PAX6.