Intron Length Research Articles

Chloroplast genome sequencing and divergence analysis of 18 Pyrus species: insights into intron length polymorphisms and evolutionary processes

Pears constitute an essential temperate crop and are primarily produced through interspecific hybridization owing to self-incompatibility that complicates their breeding history. To address this, we sequenced the complete chloroplast (cp) genomes of 18 Pyrus and one Malus species using the Illumina HiSeq4000 platform. The cp genomes ranged from 159,885 bp to 160,153 bp and exhibited a conserved circular DNA structure with an average GC content of 36.5%. Each cp genome contained 127 genes, including 83 protein-coding, 36 tRNA, and 8 rRNA genes. Divergence analysis with mVISTA showed high conservation in the coding regions and notable variations in the non-coding regions. All species shared 17 intron-containing genes, with ycf3 and clpP each having two introns. Five intron-containing genes (ndhB, rpl2, rps12, trnA-UGC, and trnE-UUC) were located in the inverted repeat regions, while trnL-UAA was located in the large single-copy region, with conserved intron lengths across Pomoideae. We identified polymorphic intron sequences in the rpl22, petB, clpP, ndhA, and rps16 genes and designed primers for these regions. Notably, the two Pyrus ussuriensis accessions Doonggeullebae and Cheongdangrori showed intron-length polymorphisms despite being classified as the same species. Phylogenetic analysis of the cp genome sequences revealed two major clusters, indicating distinct maternal lineages and evolutionary origins. This study underscores the importance of cp gene polymorphisms in P. fauriei, P. calleryana, P. ussuriensis, and P. pyrifolia, providing valuable insights into Pyrus evolution as well as aiding in the conservation and breeding of pear germplasm.

Bioinformatics-Driven Refinement of the Commonly Used TPI Nonsense-Mediated Decay Reporter System.

The cellular nonsense-mediated decay (NMD) pathway recognizes and degrades mRNAs with unusual structural features, such as long 3' UTRs or overlapping reading frames, and therefore serves as a transcript quality control mechanism. A broad spectrum of today's knowledge about the nonsense-mediated mRNA decay pathway has been discovered using NMD reporter systems, mostly consisting of multiple exons, with a wild type (WT) and a premature termination codon (PTC) containing variant. In a preliminary NMD study, we used the seven-exon triose phosphate isomerase (TPI) reporter and observed that in this well-known NMD reporter, surprisingly, not all splice sites are used constitutively, but additional cryptic splice sites are used. As this is more frequently observed in the construction of minigenes, especially when unknown splicing regulatory elements are removed, e.g. by shortening introns, this may affect the reliability of such reporters. To demonstrate how such minigenes can be improved in general with respect to constitutive splice site recognition, we restored an intron length in the TPI reporter or made bioinformatic adjustments to splice regulatory elements (SREs) or intrinsic strength of the splice sites themselves. As a result, this NMD reporter could be made more robust and specific for the evaluation of NMD sensitivity within a single transcript. The modifications of the TPI reporter shown here as examples can generally be used for the transfer of cellular multiexon transcripts to minigenes.

Genome-Wide Identification and Functional Validation of Actin Depolymerizing Factor (ADF) Gene Family in Gossypium hirsutum L.

The Actin Depolymerizing Factor (ADF) protein, highly conserved among eukaryotes, is essential for plant growth, development, and stress responses. Cotton, a vital economic crop with applications spanning oilseed, textiles, and military sectors, has seen a limited exploration of its ADF gene family. This research has identified 118 unique ADF sequences across four principal cotton species: Gossypium hirsutum L., Gossypium barbadense Linn, Gossypium raimondii, and Asiatic cotton. The study found that the structural domains and physicochemical properties of these proteins are largely uniform across species. The ADF genes were classified into four subfamilies with a notable expansion in groups III and IV due to tandem and chromosomal duplication events. A thorough analysis revealed a high degree of conservation in gene structure, including exon counts and the lengths of introns and exons, with the majority of genes containing three exons, aligning with the characteristics of the ADF family. RNA-seq analysis uncovered a spectrum of responses by GhADFs to various abiotic stresses with GhADF19 showing the most significant reaction. Virus-induced gene silencing (VIGS) experiments were conducted to assess the role of GhADF19 in plant growth under abiotic stress. The results demonstrated that plants with silenced GhADF19 exhibited significantly slower growth rates and lower dry weights when subjected to cold, salt, and drought stress compared to the control group. This marked reduction in growth and dry weight under stress conditions highlights the potential importance of GhADF19 in stress tolerance mechanisms.

Climate-dependent distribution of Lemna minor, L. turionifera and L. × japonica (Lemnaceae) in temperate Eurasia and high variability of their genome size

In spite of a burst of duckweed biology studies, understanding their diversity is still in progress. We examined 174 reliably-identified (intron length polymorphism in TUBB2) samples of the Lemna minor complex from East Europe and North Asia (41 localities of L. minor, 58 localities of L. turionifera, and 75 localities of their hybrid, L. × japonica) and described clear patterns of their distribution. It is mainly determined by climatic preferences (i.e., the ability to survive unfavorable conditions). Namely, the ranges of L. minor and L. turionifera in the lowlands are separated by the Urals, with populations located west and east of these mountains respectively and coinciding with an eastward increase of climate continentality. However, historical (i.e. European origin of the species) and not ecological factors should explain the absence of L. minor in rather mild climate of the southern Russian Far East. The hybrid has an intermediate ability to tolerate continental climate, and co-occurs with one of its parents. In southern mountains (the Caucasus and Tian-Shan) the three taxa are separated by local climate conditions, determined by the terrain. We revealed high continuous variability of genome sizes within each taxon (two-fold in L. minor and L. turionifera; three-fold in L. × japonica), which is not explained by environmental factors. We propose an approach for distinguishing these three taxa in temperate Eurasia using morphology. We have not documented the invasive spread of L. aequinoctialis, L. minuta and L. turionifera in East Europe.

Design of potential Ilp markers for Hordeum vulgare and H. marinum Atg2, Atg7, Atg8 genes

Aim. Design of potential intron length polymorphism markers for ATG2, ATG7, ATG8 genes. Methods. Use of classical bioinformatics methods to search for loci encoding ATG2, ATG7, ATG8 in the Hordeum marinum genome based on homology with cultivated barley genes; analysis of exon-intron structure of genes, design of ILP markers. Results. The sequences of ATG2, ATG7, ATG8 genes in H. marinum were determined. Analyzed differences in the exon-intron structure of the ATG2, ATG7, ATG8 loci in cultivated (H. vulgare) and salt-resistant wild barley (H. marinum). 25 potential ILP markers were designed for introns that showed significant interspecies length polymorphism: 13 for ATG2, 6 for ATG7, 3 for ATG8A, 2 for ATG8B and 1 for ATG8D. Conclusions. Based on gene structure comparative analysis of the autophagosome formation proteins ATG2, ATG7, ATG8 of cultivated H. vulgare and salt-tolerasnt H. marinum, a number of introns are recommended for further use as a potential source of ILP markers for the study of interspecies genetic differentiation of cereals.

Cereal barcoding by assessing intron length polymorphism of γ-tubulin genes

Aim. The aim of the study was to obtain typical DNA profiles for common cereal plants using a method of assessing γ-tubulin intron length polymorphism and to evaluate the possibility of using this gene for plant species barcoding. Methods. Cereals, in particular wheat, aegilops, barley, and rice, have been studied. The method of γ-tubulin intron length polymorphism assessment was used. Amplified DNA fragments were separated via polyacrylamide gel electrophoresis under non-denaturating conditions and were visualized via silver nitrate staining. Results. The analyzed cereal species-specific DNA profiles of amplicons of the γ-tubulin gene introns were obtained, which made it possible to differentiate the species. Fingerprinting data was used for cluster analysis and dendrogram construction. Conclusions. The feasibility of using the proposed cereal barcoding method has been established. The studied plant DNA profiles can be effectively used in assessing the quality of food products such as flour.

Post-transfer adaptation of HGT-acquired genes and contribution to guanine metabolic diversification in land plants.

Horizontal gene transfer (HGT) is a major driving force in the evolution of prokaryotic and eukaryotic genomes. Despite recent advances in distribution and ecological importance, the extensive pattern, especially in seed plants, and post-transfer adaptation of HGT-acquired genes in land plants remain elusive. We systematically identified 1150 foreign genes in 522 land plant genomes that were likely acquired via at least 322 distinct transfers from nonplant donors and confirmed that recent HGT events were unevenly distributed between seedless and seed plants. HGT-acquired genes evolved to be more similar to native genes in terms of average intron length due to intron gains, and HGT-acquired genes containing introns exhibited higher expression levels than those lacking introns, suggesting that intron gains may be involved in the post-transfer adaptation of HGT in land plants. Functional validation of bacteria-derived gene GuaD in mosses and gymnosperms revealed that the invasion of foreign genes introduced a novel bypass of guanine degradation and resulted in the loss of native pathway genes in some gymnosperms, eventually shaping three major types of guanine metabolism in land plants. We conclude that HGT has played a critical role in land plant evolution.

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.

Integrative analyses on the ciliates Colpoda illuminate the life history evolution of soil microorganisms.

Microorganisms play a central role in sustaining soil ecosystems and agriculture, and these functions are usually associated with their complex life history. Yet, the regulation and evolution of life history have remained enigmatic and poorly understood, especially in protozoa, the third most abundant group of organisms in the soil. Here, we explore the life history of a cosmopolitan species-Colpoda steinii. Our analysis has yielded a high-quality macronuclear genome for C. steinii, with size of 155 Mbp and 37,123 protein-coding genes, as well as mean intron length of ~93 bp, longer than most other studied ciliates. Notably, we identify two possible whole-genome duplication events in C. steinii, which may account for its genome being about twice the size of C. inflata's, another co-existing species. We further resolve the gene expression profiles in diverse life stages of C. steinii, which are also corroborated in C. inflata. During the resting cyst stage, genes associated with cell death and vacuole formation are upregulated, and translation-related genes are downregulated. While the translation-related genes are upregulated during the excystment of resting cysts. Reproductive cysts exhibit a significant reduction in cell adhesion. We also demonstrate that most genes expressed in specific life stages are under strong purifying selection. This study offers a deeper understanding of the life history evolution that underpins the extraordinary success and ecological functions of microorganisms in soil ecosystems.IMPORTANCEColpoda species, as a prominent group among the most widely distributed and abundant soil microorganisms, play a crucial role in sustaining soil ecosystems and promoting plant growth. This investigation reveals their exceptional macronuclear genomic features, including significantly large genome size, long introns, and numerous gene duplications. The gene expression profiles and the specific biological functions associated with the transitions between various life stages are also elucidated. The vast majority of genes linked to life stage transitions are subject to strong purifying selection, as inferred from multiple natural strains newly isolated and deeply sequenced. This substantiates the enduring and conservative nature of Colpoda's life history, which has persisted throughout the extensive evolutionary history of these highly successful protozoa in soil. These findings shed light on the evolutionary dynamics of microbial eukaryotes in the ever-fluctuating soil environments. This integrative research represents a significant advancement in understanding the life histories of these understudied single-celled eukaryotes.

HnRNPM protects against the dsRNA-mediated interferon response by repressing LINE-associated cryptic splicing

RNA splicing is pivotal in post-transcriptional gene regulation, yet the exponential expansion of intron length in humans poses a challenge for accurate splicing. Here, we identify hnRNPM as an essential RNA-binding protein that suppresses cryptic splicing through binding to deep introns, maintaining human transcriptome integrity. Long interspersed nuclear elements (LINEs) in introns harbor numerous pseudo splice sites. hnRNPM preferentially binds at intronic LINEs to repress pseudo splice site usage for cryptic splicing. Remarkably, cryptic exons can generate long dsRNAs through base-pairing of inverted ALU transposable elements interspersed among LINEs and consequently trigger an interferon response, a well-known antiviral defense mechanism. Significantly, hnRNPM-deficient tumors show upregulated interferon-associated pathways and elevated immune cell infiltration. These findings unveil hnRNPM as a guardian of transcriptome integrity by repressing cryptic splicing and suggest that targeting hnRNPM in tumors may be used to trigger an inflammatory immune response, thereby boosting cancer surveillance.

Genome-Wide Identification and Expression Pattern Analysis of TIFY Family Genes Reveal Their Potential Roles in Phalaenopsis aphrodite Flower Opening.

The TIFY gene family (formerly known as the zinc finger proteins expressed in inflorescence meristem (ZIM) family) not only functions in plant defense responses but also are widely involved in regulating plant growth and development. However, the identification and functional analysis of TIFY proteins remain unexplored in Orchidaceae. Here, we identified 19 putative TIFY genes in the Phalaenopsis aphrodite genome. The phylogenetic tree classified them into four subfamilies: 14 members from JAZ, 3 members from ZML, and 1 each from PPD and TIFY. Sequence analysis revealed that all Phalaenopsis TIFY proteins contained a TIFY domain. Exon-intron analysis showed that the intron number and length of Phalaenopsis TIFY genes varied, whereas the same subfamily and subgroup genes had similar exon or intron numbers and distributions. The most abundant cis-elements in the promoter regions of the 19 TIFY genes were associated with light responsiveness, followed by MeJA and ABA, indicating their potential regulation by light and phytohormones. The 13 candidate TIFY genes screened from the transcriptome data exhibited two types of expression trends, suggesting their different roles in cell proliferation and cell expansion of floral organ growth during Phalaenopsis flower opening. Overall, this study serves as a background for investigating the underlying roles of TIFY genes in floral organ growth in Phalaenopsis.

An interpretable model of pre-mRNA splicing for animal and plant genes.

Pre-mRNA splicing is a fundamental step in gene expression, conserved across eukaryotes, in which the spliceosome recognizes motifs at the 3' and 5' splice sites (SSs), excises introns, and ligates exons. SS recognition and pairing is often influenced by protein splicing factors (SFs) that bind to splicing regulatory elements (SREs). Here, we describe SMsplice, a fully interpretable model of pre-mRNA splicing that combines models of core SS motifs, SREs, and exonic and intronic length preferences. We learn models that predict SS locations with 83 to 86% accuracy in fish, insects, and plants and about 70% in mammals. Learned SRE motifs include both known SF binding motifs and unfamiliar motifs, and both motif classes are supported by genetic analyses. Our comparisons across species highlight similarities between non-mammals, increased reliance on intronic SREs in plant splicing, and a greater reliance on SREs in mammalian splicing.

Chromosome-level genome assembly and characterization of the Calophaca sinica genome.

Calophaca sinica is a rare plant endemic to northern China which belongs to the Fabaceae family and possesses rich nutritional value. To support the preservation of the genetic resources of this plant, we have successfully generated a high-quality genome of C. sinica (1.06 Gb). Notably, transposable elements (TEs) constituted ~73% of the genome, with long terminal repeat retrotransposons (LTR-RTs) dominating this group of elements (~54% of the genome). The average intron length of the C. sinica genome was noticeably longer than what has been observed for closely related species. The expansion of LTR-RTs and elongated introns emerged had the largest influence on the enlarged genome size of C. sinica in comparison to other Fabaceae species. The proliferation of TEs could be explained by certain modes of gene duplication, namely, whole genome duplication (WGD) and dispersed duplication (DSD). Gene family expansion, which was found to enhance genes associated with metabolism, genetic maintenance, and environmental stress resistance, was a result of transposed duplicated genes (TRD) and WGD. The presented genomic analysis sheds light on the genetic architecture of C. sinica, as well as provides a starting point for future evolutionary biology, ecology, and functional genomics studies centred around C. sinica and closely related species.

Positioning of pyrimidine motifs around cassette exons defines their PTB-dependent splicing in Arabidopsis.

Alternative splicing (AS) is a complex process that generates transcript variants from a single pre-mRNA and is involved in numerous biological functions. Many RNA-binding proteins are known to regulate AS; however, little is known about the underlying mechanisms, especially outside the mammalian clade. Here, we show that polypyrimidine tract binding proteins (PTBs) from Arabidopsis thaliana regulate AS of cassette exons via pyrimidine (Py)-rich motifs close to the alternative splice sites. Mutational studies on three PTB-dependent cassette exon events revealed that only some of the Py motifs in this region are critical for AS. Moreover, in vitro binding of PTBs did not reflect a motif's impact on AS in vivo. Our mutational studies and bioinformatic investigation of all known PTB-regulated cassette exons from A. thaliana and human suggested that the binding position of PTBs relative to a cassette exon defines whether its inclusion or skipping is induced. Accordingly, exon skipping is associated with a higher frequency of Py stretches within the cassette exon, and in human also upstream of it, whereas exon inclusion is characterized by increased Py motif occurrence downstream of said exon. Enrichment of Py motifs downstream of PTB-activated 5' splice sites is also seen for PTB-dependent intron removal and alternative 5' splice site events from A. thaliana, suggesting this is a common step of exon definition. In conclusion, the position-dependent AS regulatory mechanism by PTB homologs has been conserved during the separate evolution of plants and mammals, while other critical features, in particular intron length, have considerably changed.

Assessment of genetic diversity and population structure of carrot accessions with DcSto Intron Length Polymorphism markers

Assessment of genetic diversity and population structure of carrot accessions with <i>DcSto</i> Intron Length Polymorphism markers

Reexamination of Aspergillus cristatus phylogeny in dark tea: Characteristics of the mitochondrial genome.

To enhance our understanding of Aspergillus cristatus, an important functional microorganism, the characteristics of its mitochondrial genome were analyzed and compared with related species. The mitochondrial genome of A. cristatus was determined to be 77,649 bp in length, with 15 protein-coding regions. Notably, its length surpassed that of the other species, primarily attributable to the intron length. Gene order exhibited significant variations, with greater conservation observed in the genus Penicillium compared to Aspergillus. Phylogenetic tree analyses indicated that the genera Aspergillus and Penicillium are closely related but monophyletic. Furthermore, the phylogenetic tree constructed based on protein-coding genes effectively distinguished all strains with high branching confidence. This approach provides a robust reflection of the evolutionary relationship between A. cristatus and its related species, offering potential for the development of molecular markers suitable for Aspergillus and Penicillium.

Patterns of structural variants within TP53 introns and relocation of the TP53 promoter: a commentary†.

Gene disruption from double-strand DNA breaks within introns is a mechanism of inactivating the tumor suppressor TP53. This occurs more frequently in osteosarcoma and biliary adenocarcinoma compared with other cancer types. The patterns of intron breakpoints within TP53 do not correlate with prevalence, intron length, or overall genome-wide levels of rearrangements. Therefore, these breakpoints appear to be selected for reasons other than to disrupt TP53. A recent article published by Saba et al in The Journal of Pathology illustrates a benefit to having breakpoints within intron 1 using high-quality matched genomic and transcriptomic osteosarcoma sequencing data as well as in vitro validation. The authors describe how the rearrangement results in relocation of the TP53 promoter region to regions upstream of genes that encode members of cartilage, growth plate development, osteoclast formation, and other TP53-related pathways. The upregulation of these genes by the TP53 promoter are gain-of-function events that are likely to promote tumor development and growth. Therefore, this article presents a potential new paradigm in which a single mutation would result in both the loss of a tumor suppressor and the gain of an oncogenic program. © 2024 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Development of Genome-Wide Intron Length Polymorphism (ILP) Markers in Tea Plant (Camellia sinensis) and Related Applications for Genetics Research.

The market value of tea is largely dependent on the tea species and cultivar. Therefore, it is important to develop efficient molecular markers covering the entire tea genome that can be used for the identification of tea varieties, marker-assisted breeding, and mapping important quantitative trait loci for beneficial traits. In this study, genome-wide molecular markers based on intron length polymorphism (ILP) were developed for tea trees. A total of 479, 1393, and 1342 tea ILP markers were identified using the PCR method in silico from the 'Shuchazao' scaffold genome, the chromosome-level genome of 'Longjing 43', and the ancient tea DASZ chromosome-level genome, respectively. A total of 230 tea ILP markers were used to amplify six tea tree species. Among these, 213 pairs of primers successfully characterize products in all six species, with 112 primer pairs exhibiting polymorphism. The polymorphism rate of primer pairs increased with the improvement in reference genome assembly quality level. The cross-species transferability analysis of 35 primer pairs of tea ILP markers showed an average amplification rate of 85.17% through 11 species in 6 families, with high transferability in Camellia reticulata and tobacco. We also used 40 pairs of tea ILP primers to evaluate the genetic diversity and population structure of C. tetracocca with 176 plants from Puan County, Guizhou Province, China. These genome-wide markers will be a valuable resource for genetic diversity analysis, marker-assisted breeding, and variety identification in tea, providing important information for the tea industry.

Transcription and splicing variations of SR genes accompany with genome-wide accumulation of long-introns in pine

Most of mRNAs in Eukaryote were matured after the removal of introns in their pre-mRNA transcripts. Serine/arginine-rich (SR) proteins are a group of splicing regulators regulating the splicing processes globally. Expressions of SR proteins themselves were extensively regulated, at both transcription and splicing levels, under different environmental conditions, specially heat stress conditions. The pine genome is characterized by super-long and easily methylated introns in a large number of genes that derived from the extensive accumulation of transposons (TEs). Here, we identified and analyzed the phylogenetic characteristics of 24 SR proteins and their encoding genes from the pine genome. Then we explored transcription and pre-mRNA splicing expression patterns of SR genes in P. massoniana seedlings under normal and heat stress temperature conditions. Our results showed that the transcription patterns of SR genes in pine exhibited significant changes compared to other plant species, and these changes were not strictly correlated with the intron length and DNA methylation intensity of the SR genes. Interestingly, none of the long introns of SR genes underwent alternative splicing (AS) in our experiment. Furthermore, the intensity of AS regulation may be related to the potential DNA methylation intensity of SR genes. Taken together, this study explores for the first time the characteristics of significant variations in the transcription and splicing patterns of SR proteins in a plant species with an over-accumulation of super-long introns.