Analysis of the Functional Relevance of Epigenetic Chromatin Marks in the First Intron Associated with Specific Gene Expression Patterns

Pervasive Conservation of Intron Number and Other Genetic Elements Revealed by a Chromosome-level Genome Assembly of the Hyper-polymorphic Nematode Caenorhabditis brenneri

With within-species genetic diversity estimates that span the gamut of that seen across the entirety of animals, the Caenorhabditis genus of nematodes holds unique potential to provide insights into how population size and reproductive strategies influence gene and genome organization and evolution. Our study focuses on Caenorhabditis brenneri, currently known as one of the most genetically diverse nematodes within its genus and, notably, across Metazoa. Here, we present a high-quality, gapless genome assembly and annotation for C. brenneri, revealing a common nematode chromosome arrangement characterized by gene-dense central regions and repeat-rich arms. A comparison of C. brenneri with other nematodes from the "Elegans" group revealed conserved macrosynteny but a lack of microsynteny, characterized by frequent rearrangements and low correlation of orthogroup size, indicative of high rates of gene turnover, consistent with previous studies. We also assessed genome organization within corresponding syntenic blocks in selfing and outcrossing species, affirming that selfing species predominantly experience loss of both genes and intergenic DNA. A comparison of gene structures revealed a strikingly small number of shared introns across species, yet consistent distributions of intron number and length, regardless of population size or reproductive mode, suggesting that their evolutionary dynamics are primarily reflective of functional constraints. Our study provides valuable insights into genome evolution and expands the nematode genome resources with the highly genetically diverse C. brenneri, facilitating research into various aspects of nematode biology and evolutionary processes.

Genes divided according to the relative position of the longest intron show increased representation in different KEGG pathways

Despite the fact that introns mean an energy and time burden for eukaryotic cells, they play an irreplaceable role in the diversification and regulation of protein production. As a common feature of eukaryotic genomes, it has been reported that in protein-coding genes, the longest intron is usually one of the first introns. The goal of our work was to find a possible difference in the biological function of genes that fulfill this common feature compared to genes that do not. Data on the lengths of all introns in genes were extracted from the genomes of six vertebrates (human, mouse, koala, chicken, zebrafish and fugu) and two other model organisms (nematode worm and arabidopsis). We showed that more than 40% of protein-coding genes have the relative position of the longest intron located in the second or third tertile of all introns. Genes divided according to the relative position of the longest intron were found to be significantly increased in different KEGG pathways. Genes with the longest intron in the first tertile predominate in a range of pathways for amino acid and lipid metabolism, various signaling, cell junctions or ABC transporters. Genes with the longest intron in the second or third tertile show increased representation in pathways associated with the formation and function of the spliceosome and ribosomes. In the two groups of genes defined in this way, we further demonstrated the difference in the length of the longest introns and the distribution of their absolute positions. We also pointed out other characteristics, namely the positive correlation between the length of the longest intron and the sum of the lengths of all other introns in the gene and the preservation of the exact same absolute and relative position of the longest intron between orthologous genes.

Pervasive conservation of intron number and other genetic elements revealed by a chromosome-level genomic assembly of the hyper-polymorphic nematode Caenorhabditis brenneri

With within-species genetic diversity estimates that span the gambit of that seen across the entirety of animals, the Caenorhabditis genus of nematodes holds unique potential to provide insights into how population size and reproductive strategies influence gene and genome organization and evolution. Our study focuses on Caenorhabditis brenneri, currently known as one of the most genetically diverse nematodes within its genus and metazoan phyla. Here, we present a high-quality gapless genome assembly and annotation for C. brenneri, revealing a common nematode chromosome arrangement characterized by gene-dense central regions and repeat rich peripheral parts. Comparison of C. brenneri with other nematodes from the 'Elegans' group revealed conserved macrosynteny but a lack of microsynteny, characterized by frequent rearrangements and low correlation iof orthogroup sizes, indicative of high rates of gene turnover. We also assessed genome organization within corresponding syntenic blocks in selfing and outcrossing species, affirming that selfing species predominantly experience loss of both genes and intergenic DNA. Comparison of gene structures revealed strikingly small number of shared introns across species, yet consistent distributions of intron number and length, regardless of population size or reproductive mode, suggesting that their evolutionary dynamics are primarily reflective of functional constraints. Our study provides valuable insights into genome evolution and expands the nematode genome resources with the highly genetically diverse C. brenneri, facilitating research into various aspects of nematode biology and evolutionary processes.

A super minigene with a short promoter and truncated introns recapitulates essential features of transcription and splicing regulation of the SMN1 and SMN2 genes

Here we report a Survival Motor Neuron 2 (SMN2) super minigene, SMN2Sup, encompassing its own promoter, all exons, their flanking intronic sequences and the entire 3'-untranslated region. We confirm that the pre-mRNA generated from SMN2Sup undergoes splicing to produce a translation-competent mRNA. We demonstrate that mRNA generated from SMN2Sup produces more SMN than an identical mRNA generated from a cDNA clone. We uncover that overexpression of SMN triggers skipping of exon 3 of SMN1/SMN2. We define the minimal promoter and regulatory elements associated with the initiation and elongation of transcription of SMN2. The shortened introns within SMN2Sup preserved the ability of camptothecin, a transcription elongation inhibitor, to induce skipping of exons 3 and 7 of SMN2. We show that intron 1-retained transcripts undergo nonsense-mediated decay. We demonstrate that splicing factor SRSF3 and DNA/RNA helicase DHX9 regulate splicing of multiple exons in the context of both SMN2Sup and endogenous SMN1/SMN2. Prevention of SMN2 exon 7 skipping has implications for the treatment of spinal muscular atrophy (SMA). We validate the utility of the super minigene in monitoring SMN levels upon splicing correction. Finally, we demonstrate how the super minigene could be employed to capture the cell type-specific effects of a pathogenic SMN1 mutation.

The position of the longest intron is related to biological functions in some human genes

The evidence that introns can influence different levels of transfer of genetic information between DNA and the final product is increasing. Longer first introns were found to be a general property of eukaryotic gene structure and shown to contain a higher fraction of conserved sequence and different functional elements. Our work brings more precise information about the position of the longest introns in human protein-coding genes and possible connection with biological function and gene expression. According to our results, the position of the longest intron can be localized to the first third of introns in 64%, the second third in 19%, and the third in 17%, with notable peaks at the middle and last introns of approximately 5% and 6%, respectively. The median lengths of the longest introns decrease with increasing distance from the start of the gene from approximately 15,000 to 5,000 bp. We have shown that the position of the longest intron is in some cases linked to the biological function of the given gene. For example, DNA repair genes have the longest intron more often in the second or third. In the distribution of gene expression according to the position of the longest intron, tissue-specific profiles can be traced with the highest expression usually at the absolute positions of intron 1 and 2. In this work, we present arguments supporting the hypothesis that the position of the longest intron in a gene is another biological factor modulating the transmission of genetic information. The position of the longest intron is related to biological functions in some human genes.

Phosphorylation-dependent association of human chromatin protein PC4 to linker histone H1 regulates genome organization and transcription

Human Positive Coactivator 4 (PC4) is a multifaceted chromatin protein involved in diverse cellular processes including genome organization, transcription regulation, replication, DNA repair and autophagy. PC4 exists as a phospho-protein in cells which impinges on its acetylation by p300 and thereby affects its transcriptional co-activator functions via double-stranded DNA binding. Despite the inhibitory effects, the abundance of phosphorylated PC4 in cells intrigued us to investigate its role in chromatin functions in a basal state of the cell. We found that casein kinase-II (CKII)-mediated phosphorylation of PC4 is critical for its interaction with linker histone H1. By employing analytical ultracentrifugation and electron microscopy imaging of in vitro reconstituted nucleosomal array, we observed that phospho-mimic (PM) PC4 displays a superior chromatin condensation potential in conjunction with linker histone H1. ATAC-sequencing further unveiled the role of PC4 phosphorylation to be critical in inducing chromatin compaction of a wide array of coding and non-coding genes in vivo. Concordantly, phospho-PC4 mediated changes in chromatin accessibility led to gene repression and affected global histone modifications. We propose that the abundance of PC4 in its phosphorylated state contributes to genome compaction contrary to its co-activator function in driving several cellular processes like gene transcription and autophagy.

Functional Characterization and Molecular Marker Development of the Proenkephalin as Biomarker of Food Addiction in Food Habit Domestication of Mandarin Fish (Siniperca Chuatsi)

Proenkephalin (PENK), as the precursor of endogenous opioid enkephalin (ENK), is widely present in the nervous system and plays an important role in animal food addiction and rewarding behavior. In our study, we intend to study the functional characterization and molecular marker development of the penk gene related to food habit domestication of mandarin fish. We found that the penk gene of mandarin fish had three types of endogenous opioid peptide sequences. Compared with other tissues, penk mRNA was highly expressed in the whole brain. Intracerebroventricular (ICV) injection of lysine or methionine significantly increased the expression of penk mRNA. The expression of penk mRNA in the brain of mandarin fish that could be easily domesticated from eating live prey fish to artificial diets was significantly higher than those that could not. After feeding with high-carbohydrate artificial diets, the expression of penk mRNA showed no significant difference between mandarin fish with hypophagia and those that still ate normally. A total of four single nucleotide polymorphisms (SNP) loci related to easy domestication toward eating artificial diets were screened from the mandarin fish population. Additionally, the TT genotype at one of the loci was significantly correlated with the food habit domestication of mandarin fish.

Introns control stochastic allele expression bias

Monoallelic expression (MAE) or extreme allele bias can account for incomplete penetrance, missing heritability and non-Mendelian diseases. In cancer, MAE is associated with shorter patient survival times and higher tumor grade. Prior studies showed that stochastic MAE is caused by stochastic epigenetic silencing, in a gene and tissue-specific manner. Here, we used C. elegans to study stochastic MAE in vivo. We found allele bias/MAE to be widespread within C. elegans tissues, presenting as a continuum from fully biallelic to MAE. We discovered that the presence of introns within alleles robustly decreases MAE. We determined that introns control MAE at distinct loci, in distinct cell types, with distinct promoters, and within distinct coding sequences, using a 5'-intron position-dependent mechanism. Bioinformatic analysis showed human intronless genes are significantly enriched for MAE. Our experimental evidence demonstrates a role for introns in regulating MAE, possibly explaining why some mutations within introns result in disease.

Noncoding sequence variants define a novel regulatory element in the first intron of the N-acetylglutamate synthase gene

N-acetylglutamate synthase deficiency is an autosomal recessive urea cycle disorder caused either by decreased expression of the NAGS gene or defective NAGS enzyme resulting in decreased production of N-acetylglutamate (NAG), an allosteric activator of carbamylphosphate synthetase 1 (CPS1). NAGSD is the only urea cycle disorder that can be effectively treated with a single drug, N-carbamylglutamate (NCG), a stable NAG analog, which activates CPS1 to restore ureagenesis. We describe three patients with NAGSD due to four novel noncoding sequence variants in the NAGS regulatory regions. All three patients had hyperammonemia that resolved upon treatment with NCG. Sequence variants NM_153006.2:c.427-222G>A and NM_153006.2:c.427-218A>C reside in the 547 bp-long first intron of NAGS and define a novel NAGS regulatory element that binds retinoic X receptor α. Sequence variants NC_000017.10:g.42078967A>T (NM_153006.2:c.-3065A>T) and NC_000017.10:g.42078934C>T (NM_153006.2:c.-3098C>T) reside in the NAGS enhancer, within known HNF1 and predicted glucocorticoid receptor binding sites, respectively. Reporter gene assays in HepG2 and HuH-7 cells demonstrated that all four substitutions could result in reduced expression of NAGS. These findings show that analyzing noncoding regions of NAGS and other urea cycle genes can reveal molecular causes of disease and identify novel regulators of ureagenesis.

A Global Vista of the Epigenomic State of the Mouse Submandibular Gland

The parotid, submandibular, and sublingual glands represent a trio of oral secretory glands whose primary function is to produce saliva, facilitate digestion of food, provide protection against microbes, and maintain oral health. While recent studies have begun to shed light on the global gene expression patterns and profiles of salivary glands, particularly those of mice, relatively little is known about the location and identity of transcriptional control elements. Here we have established the epigenomic landscape of the mouse submandibular salivary gland (SMG) by performing chromatin immunoprecipitation sequencing experiments for 4 key histone marks. Our analysis of the comprehensive SMG data sets and comparisons with those from other adult organs have identified critical enhancers and super-enhancers of the mouse SMG. By further integrating these findings with complementary RNA-sequencing based gene expression data, we have unearthed a number of molecular regulators such as members of the Fox family of transcription factors that are enriched and likely to be functionally relevant for SMG biology. Overall, our studies provide a powerful atlas of cis-regulatory elements that can be leveraged for better understanding the transcriptional control mechanisms of the mouse SMG, discovery of novel genetic switches, and modulating tissue-specific gene expression in a targeted fashion.

A novel genetic variant associated with benign paroxysmal positional vertigo within the LOXL1

Background

Benign paroxysmal positional vertigo (BPPV) is a common, self‐limited, and favorable prognostic peripheral vestibular disorder. BPPV is transmitted in an autosomal dominant fashion, but most cases occur sporadically. Little research has been reported regarding the mutation spectrum of sporadic BPPV in a large cohort. This study attempted to identify the causative candidate variants associated with BPPV in VDR, LOXL1, and LOXL1‐AS1.

Methods

An amplicon‐targeted next‐generation sequencing (NGS) method for VDR, LOXL1, and LOXL1‐AS1, was completed in 726 BPPV patients and 502 normal controls. A total of 30 variants (20 variants from VDR, nine variants from LOXL1, seven variants from LOXL1‐AS1) were identified in these two groups.

Results

Three of 30 variants were nonsynonymous mutations, but no significant difference was found between the BPPV group and the control group via association analysis. A single nucleotide variant (SNV), rs1078967, was identified that is located in intron 1 of LOXL1. The allelic frequency distribution differed significantly between the BPPV group and the control group (p = 0.002). Genotypic frequency was also significantly different (p = 0.006), as determined by gene‐based analyses.

Conclusion

This report is the first to analyze the variant spectrum of BPPV in a large Chinese population.

First Come, First Served: Sui Generis Features of the First Intron

Most of the transcribed genes in eukaryotic cells are interrupted by intervening sequences called introns that are co-transcriptionally removed from nascent messenger RNA through the process of splicing. In Arabidopsis, 79% of genes contain introns and more than 60% of intron-containing genes undergo alternative splicing (AS), which ostensibly is considered to increase protein diversity as one of the intrinsic mechanisms for fitness to the varying environment or the internal developmental program. In addition, recent findings have prevailed in terms of overlooked intron functions. Here, we review recent progress in the underlying mechanisms of intron function, in particular by focusing on unique features of the first intron that is located in close proximity to the transcription start site. The distinct deposition of epigenetic marks and nucleosome density on the first intronic DNA sequence, the impact of the first intron on determining the transcription start site and elongation of its own expression (called intron-mediated enhancement, IME), translation control in 5′-UTR, and the new mechanism of the trans-acting function of the first intron in regulating gene expression at the post-transcriptional level are summarized.

The Histone Marks Signature in Exonic and Intronic Regions Is Relevant in Early Response of Tomato Genes to Botrytis cinerea and in miRNA Regulation

Research into the relationship between epigenetic regulation and resistance to biotic stresses provides alternatives for plant protection and crop improvement. To unravel the mechanisms underlying tomato responses to Botrytis cinerea, we performed a chromatin immunoprecipitation (ChIP) analysis showing the increase in H3K9ac mark along the early induced genes SlyDES, SlyDOX1, and SlyLoxD encoding oxylipin-pathway enzymes, and SlyWRKY75 coding for a transcriptional regulator of hormonal signaling. This histone mark showed a more distinct distribution than the previously studied H3K4me3. The RNAPol-ChIP analysis reflected the actual gene transcription associated with increased histone modifications. A different pattern of marks in the oxylipin-related genes against P. syringae supported a pathogen-specific profile, while no significant differences occurred in SlyWRKY75. The epigenetic regulation of SlyWRKY75 by the intron-binding miR1127-3p was supported by the presence of SlyWRKY75 pre-mRNA in control plants. Interestingly, mRNA was found to be accumulated in response to B. cinerea and P. syringae, while reduction in miRNA only occurred against B. cinerea. The intronic region presented a similar pattern of marks than the rest of the gene in both pathosystems, except for H3K4me3 in the miRNA binding site upon B. cinerea. We located the gene encoding Sly-miR1127-3p, which presented reduced H3K4me3 on its promoter against B. cinerea.

Enrichment of rare alleles within epigenetic chromatin marks in the first intron

In previous studies, we demonstrated that some sites in the first intron likely regulate gene expression. In the present work, we sought to further confirm the functional relevance of first intron sites by estimating the quantity of rare alleles in the first intron. A basic hypothesis posited herein is that genomic regions carrying more functionally important sites will have a higher proportion of rare alleles. We estimated the proportions of rare single nucleotide polymorphisms with a minor allele frequency < 0.01 located in several histone marks in the first introns of various genes, and compared them with those in other introns and those in 2-kb upstream regions. As expected, rare alleles were found to be significantly enriched in most of the regulatory sites located in the first introns. Meanwhile, transcription factor binding sites were significantly more enriched in the 2-kb upstream regions (i.e., the regions of putative promoters of genes) than in the first introns. These results strongly support our proposal that the first intron sites of genes may have important regulatory functions in gene expression independent of promoters.