Exon size distribution and the origin of introns

Characterization of MKK family genes and their responses to temperature stress and Vibrio parahaemolyticus infection in noble scallop Chlamys nobilis

Mitogen-activated protein kinase kinase (MKK), the key element of the Mitogen-activated protein kinase (MAPK) signaling pathway, is crucial for the immune response to adverse environments in aquatic animals. Nevertheless, there is limited information regarding the role of the MKK gene family in mollusks. In our study, genome data and transcriptome were used to identify four MKK genes (CnMKK4, CnMKK5, CnMKK6, and CnMKK7) in the noble scallop. The result of the gene structure, motif analysis, and phylogenetic tree revealed that MKK genes are relatively conserved in bivalves. Moreover, four CnMKK genes were significantly highly expressed in immune-related tissues, suggesting that CnMKKs may related to bivalve immunity. Furthermore, CnMKK6 and CgMKK4 were significantly differentially expressed (P < 0.05) under 24 h of temperature stress, and all CnMKKs were significantly differentially expressed (P < 0.05) under 24 h of Vibrio parahaemolyticus infection. These results showed that the CnMKKs may have a significant impact under biotic and abiotic stresses. In conclusion, the result of the CnMKKs provides valuable insights into comprehending the function of MKK genes in mollusks.

Superstatistical and DNA sequence coding of the human genome

In this work, by considering superstatistics we investigate the short-range correlations (SRCs) and the fluctuations in the distribution of lengths of strings of nucleotides. To this end, a stochastic model provides the distributions of the size of the exons based on the q-Gamma and inverse q-Gamma distributions. Specifically, we define a time series for exon sizes to investigate the SRC and the fluctuations through the superstatistics distributions. To test the model's viability, we use the Project Ensembl database of genes to extract the time evolution of exon sizes, calculated in terms of the number of base pairs (bp) in these biological databases. Our findings show that, depending on the chromosome, both distributions are suitable for describing the length distribution of human DNA for lengths greater than 10 bp. In addition, we used Bayesian statistics to perform a selection model approach, which revealed weak evidence for the inverse q-Gamma distribution for a considerable number of chromosomes.

Coupled protein synthesis and ribosome-guided piRNA processing on mRNAs

PIWI-interacting small RNAs (piRNAs) protect the germline genome and are essential for fertility. piRNAs originate from transposable element (TE) RNAs, long non-coding RNAs, or 3´ untranslated regions (3´UTRs) of protein-coding messenger genes, with the last being the least characterized of the three piRNA classes. Here, we demonstrate that the precursors of 3´UTR piRNAs are full-length mRNAs and that post-termination 80S ribosomes guide piRNA production on 3´UTRs in mice and chickens. At the pachytene stage, when other co-translational RNA surveillance pathways are sequestered, piRNA biogenesis degrades mRNAs right after pioneer rounds of translation and fine-tunes protein production from mRNAs. Although 3´UTR piRNA precursor mRNAs code for distinct proteins in mice and chickens, they all harbor embedded TEs and produce piRNAs that cleave TEs. Altogether, we discover a function of the piRNA pathway in fine-tuning protein production and reveal a conserved piRNA biogenesis mechanism that recognizes translating RNAs in amniotes.

Shedding light on dark genes: enhanced targeted resequencing by optimizing the combination of enrichment technology and DNA fragment length

The exome contains many obscure regions difficult to explore with current short-read sequencing methods. Repetitious genomic regions prevent the unique alignment of reads, which is essential for the identification of clinically-relevant genetic variants. Long-read technologies attempt to resolve multiple-mapping regions, but they still produce many sequencing errors. Thus, a new approach is required to enlighten the obscure regions of the genome and rescue variants that would be otherwise neglected. This work aims to improve the alignment of multiple-mapping reads through the extension of the standard DNA fragment size. As Illumina can sequence fragments up to 550 bp, we tested different DNA fragment lengths using four major commercial WES platforms and found that longer DNA fragments achieved a higher genotypability. This metric, which indicates base calling calculated by combining depth of coverage with the confidence of read alignment, increased from hundreds to thousands of genes, including several associated with clinical phenotypes. While depth of coverage has been considered crucial for the assessment of WES performance, we demonstrated that genotypability has a greater impact in revealing obscure regions, with ~1% increase in variant calling in respect to shorter DNA fragments. Results confirmed that this approach enlightened many regions previously not explored.

Stress-mediated convergence of splicing landscapes in male and female rock doves

BACKGROUND

The process of alternative splicing provides a unique mechanism by which eukaryotes are able to produce numerous protein products from the same gene. Heightened variability in the proteome has been thought to potentiate increased behavioral complexity and response flexibility to environmental stimuli, thus contributing to more refined traits on which natural and sexual selection can act. While it has been long known that various forms of environmental stress can negatively affect sexual behavior and reproduction, we know little of how stress can affect the alternative splicing associated with these events, and less still about how splicing may differ between sexes. Using the model of the rock dove (Columba livia), our team previously uncovered sexual dimorphism in the basal and stress-responsive gene transcription of a biological system necessary for facilitating sexual behavior and reproduction, the hypothalamic-pituitary-gonadal (HPG) axis. In this study, we delve further into understanding the mechanistic underpinnings of how changes in the environment can affect reproduction by testing the alternative splicing response of the HPG axis to an external stressor in both sexes.

RESULTS

This study reveals dramatic baseline differences in HPG alternative splicing between males and females. However, after subjecting subjects to a restraint stress paradigm, we found a significant reduction in these differences between the sexes. In both stress and control treatments, we identified a higher incidence of splicing activity in the pituitary in both sexes as compared to other tissues. Of these splicing events, the core exon event is the most abundant form of splicing and more frequently occurs in the coding regions of the gene. Overall, we observed less splicing activity in the 3'UTR (untranslated region) end of transcripts than the 5'UTR or coding regions.

CONCLUSIONS

Our results provide vital new insight into sex-specific aspects of the stress response on the HPG axis at an unprecedented proximate level. Males and females uniquely respond to stress, yet exhibit splicing patterns suggesting a convergent, optimal splicing landscape for stress response. This information has the potential to inform evolutionary theory as well as the development of highly-specific drug targets for stress-induced reproductive dysfunction.

Advances in Transcriptomics: Investigating Cardiovascular Disease at Unprecedented Resolution

Whole-genome transcriptional profiling has become a standard genomic approach to investigate biological processes. RNA sequencing (RNAseq) in particular has witnessed myriad applications in genetics and various biomedical fields. RNAseq involves a relatively simple experimental protocol of RNA extraction and cDNA library preparation and, because of decreasing next-generation sequencing cost and lower computational burden for data processing, has obtained a central role in the modern biology. The recent application of RNAseq methodology to single-cell transcriptional profiling has enabled the more precise characterization of cell lineage and cell state genetic profiles. The development of bioinformatic and statistical tools has provided for differential gene expression analysis, RNA isoform analysis, haplotype-specific analysis of gene expression (allele-specific expression), and analysis of expression quantitative trait loci. We give an overview of these and recent developments in RNAseq methodology with emphasis on quality control, read mapping, feature counting, differential gene expression, allele-specific expression and expression quantitative trait loci analysis, and fusion transcript detection. We describe utilization of RNAseq as a diagnostic tool in Mendelian diseases, complex phenotypes, and cancer and give an overview of long read RNAseq technology. Furthermore, we discuss in detail the recent revolution in single-cell transcriptomics that is reshaping modern biology.

P_RNA_scaffolder: a fast and accurate genome scaffolder using paired-end RNA-sequencing reads

BACKGROUND

Obtaining complete gene structures is one major goal of genome assembly. Some gene regions are fragmented in low quality and high-quality assemblies. Therefore, new approaches are needed to recover gene regions. Genomes are widely transcribed, generating messenger and non-coding RNAs. These widespread transcripts can be used to scaffold genomes and complete transcribed regions.

RESULTS

We present P_RNA_scaffolder, a fast and accurate tool using paired-end RNA-sequencing reads to scaffold genomes. This tool aims to improve the completeness of both protein-coding and non-coding genes. After this tool was applied to scaffolding human contigs, the structures of both protein-coding genes and circular RNAs were almost completely recovered and equivalent to those in a complete genome, especially for long proteins and long circular RNAs. Tested in various species, P_RNA_scaffolder exhibited higher speed and efficiency than the existing state-of-the-art scaffolders. This tool also improved the contiguity of genome assemblies generated by current mate-pair scaffolding and third-generation single-molecule sequencing assembly.

CONCLUSIONS

The P_RNA_scaffolder can improve the contiguity of genome assembly and benefit gene prediction. This tool is available at http://www.fishbrowser.org/software/P_RNA_scaffolder .

Modeling the evolution dynamics of exon-intron structure with a general random fragmentation process

Background

Most eukaryotic genes are interrupted by spliceosomal introns. The evolution of exon-intron structure remains mysterious despite rapid advance in genome sequencing technique. In this work, a novel approach is taken based on the assumptions that the evolution of exon-intron structure is a stochastic process, and that the characteristics of this process can be understood by examining its historical outcome, the present-day size distribution of internal translated exons (exon). Through the combination of simulation and modeling the size distribution of exons in different species, we propose a general random fragmentation process (GRFP) to characterize the evolution dynamics of exon-intron structure. This model accurately predicts the probability that an exon will be split by a new intron and the distribution of novel insertions along the length of the exon.

Results

As the first observation from this model, we show that the chance for an exon to obtain an intron is proportional to its size to the 3rd power. We also show that such size dependence is nearly constant across gene, with the exception of the exons adjacent to the 5^′ UTR. As the second conclusion from the model, we show that intron insertion loci follow a normal distribution with a mean of 0.5 (center of the exon) and a standard deviation of 0.11. Finally, we show that intron insertions within a gene are independent of each other for vertebrates, but are more negatively correlated for non-vertebrate. We use simulation to demonstrate that the negative correlation might result from significant intron loss during evolution, which could be explained by selection against multi-intron genes in these organisms.

Conclusions

The GRFP model suggests that intron gain is dynamic with a higher chance for longer exons; introns are inserted into exons randomly with the highest probability at the center of the exon. GRFP estimates that there are 78 introns in every 10 kb coding sequences for vertebrate genomes, agreeing with empirical observations. GRFP also estimates that there are significant intron losses in the evolution of non-vertebrate genomes, with extreme cases of around 57% intron loss in Drosophila melanogaster, 28% in Caenorhabditis elegans, and 24% in Oryza sativa.

Detection of Alternatively Spliced or Processed RNAs in Cancer Using Oligonucleotide Microarray

Deregulation of gene expression plays a pivotal role in tumorigenesis, so the ability to detect RNA alterations is of great value in cancer diagnosis and management. DNA microarrays have been used to measure changes in mRNA or microRNA level, but less often the change of RNA isoforms. Here we appraise the utilization of microarray in detecting alternatively processed RNAs, which have alternative splice forms, retained introns, or altered 3' untranslated regions. We cover the methodology and focus on cancer studies. Recent development in parallel or deep sequencing used in transcriptome analysis is also discussed.

Evolution of cyclin B3 shows an abrupt three-fold size increase, due to the extension of a single exon in placental mammals, allowing for new protein-protein interactions

Cyclin B3 evolution has the unique peculiarity of an abrupt 3-fold increase of the protein size in the mammalian lineage due to the extension of a single exon. We have analyzed the evolution of the gene to define the modalities of this event and the possible consequences on the function of the protein. Database searches can trace the appearance of the gene to the origin of metazoans. Most introns were already present in early metazoans, and the intron-exon structure as well as the protein size were fairly conserved in invertebrates and nonmammalian vertebrates. Although intron gains are considered as rare events, we identified two cases, one at the prochordate-chordate transition and one in murids, resulting from different mechanisms. At the emergence of mammals, the gene was relocated from chromosome 6 of platypus to the X chromosome in marsupials, but the exon extension occurred only in placental mammals. A repetitive structure of 18 amino acids, of uncertain origin, is detectable in the 3,000-nt mammalian exon-encoded sequence, suggesting an extension by multiple internal duplications, some of which are still detectable in the primate lineage. Structure prediction programs suggest that the repetitive structure has no associated three-dimensional structure but rather a tendency for disorder. Splice variant isoforms were detected in several mammalian species but without conserved pattern, notably excluding the constant coexistence of premammalian-like transcripts, without the extension. The yeast two-hybrid method revealed that, in human, the extension allowed new interactions with ten unrelated proteins, most of them with specific three-dimensional structures involved in protein-protein interactions, and some highly expressed in testis, as is cyclin B3. The interactions with activator of cAMP-responsive element modulator in testis (ACT), germ cell-less homolog 1, and chromosome 1 open reading frame 14 remain to be verified in vivo since they may not be expressed at the same stages of spermatogenesis as cyclin B3.

A versatile gene trap to visualize and interrogate the function of the vertebrate proteome

We report a multifunctional gene-trapping approach, which generates full-length Citrine fusions with endogenous proteins and conditional mutants from a single integration event of the FlipTrap vector. We identified 170 FlipTrap zebrafish lines with diverse tissue-specific expression patterns and distinct subcellular localizations of fusion proteins generated by the integration of an internal citrine exon. Cre-mediated conditional mutagenesis is enabled by heterotypic lox sites that delete Citrine and "flip" in its place mCherry with a polyadenylation signal, resulting in a truncated fusion protein. Inducing recombination with Cerulean-Cre results in fusion proteins that often mislocalize, exhibit mutant phenotypes, and dramatically knock down wild-type transcript levels. FRT sites in the vector enable targeted genetic manipulation of the trapped loci in the presence of Flp recombinase. Thus, the FlipTrap captures the functional proteome, enabling the visualization of full-length fluorescent fusion proteins and interrogation of function by conditional mutagenesis and targeted genetic manipulation.

Constructing Physical and Genomic Maps for Puccinia striiformis f. sp. tritici, the Wheat Stripe Rust Pathogen, by Comparing Its EST Sequences to the Genomic Sequence of P. graminis f. sp. tritici, the Wheat Stem Rust Pathogen

The wheat stripe rust fungus, Puccinia striiformis f. sp. tritici (Pst), does not have a known alternate host for sexual reproduction, which makes it impossible to study gene linkages through classic genetic and molecular mapping approaches. In this study, we compared 4,219 Pst expression sequence tags (ESTs) to the genomic sequence of P. graminis f. sp. tritici (Pgt), the wheat stem rust fungus, using BLAST searches. The percentages of homologous genes varied greatly among different Pst libraries with 54.51%, 51.21%, and 13.61% for the urediniospore, germinated urediniospore, and haustorial libraries, respectively, with an average of 33.92%. The 1,432 Pst genes with significant homology with Pgt sequences were grouped into physical groups corresponding to 237 Pgt supercontigs. The physical relationship was demonstrated by 12 pairs (57%), out of 21 selected Pst gene pairs, through PCR screening of a Pst BAC library. The results indicate that the Pgt genome sequence is useful in constructing Pst physical maps.

A chemokine targets the nucleus: Cxcl12-gamma isoform localizes to the nucleolus in adult mouse heart

Chemokines are extracellular mediators of complex regulatory circuits involved principally in cell-to-cell communication. Most studies to date of the essential chemokine Cxcl12 (Sdf-1) have focused on the ubiquitously expressed secreted isoforms α and β. Here we show that, unlike these isoforms and all other known chemokines, the alternatively transcribed γ isoform is an intracellular protein that localizes to the nucleolus in differentiated mouse Cardiac tissue. Our results demonstrate that nucleolar transportation is encoded by a nucleolar-localization signal in the unique carboxy-terminal region of Sdf-1γ, and is competent both in vivo and in vitro. The molecular mechanism underlying these unusual chemokine properties involves cardiac-specific transcription of an mRNA containing a unique short-leader sequence lacking the signal peptide and translation from a non-canonical CUG codon. Our results provide an example of genome economy even for essential and highly conserved genes such as Cxcl12, and suggest that chemokines can exert tissue specific functions unrelated to cell-to-cell communication.

Spontaneous symmetry breaking in genome evolution

The quest for evolutionary mechanisms providing separation between the coding (exons) and noncoding (introns) parts of genomic DNA remains an important focus of genetics. This work combines an analysis of the most recent achievements of genomics and fundamental concepts of random processes to provide a novel point of view on genome evolution. Exon sizes in sequenced genomes show a lognormal distribution typical of a random Kolmogoroff fractioning process. This implies that the process of intron incretion may be independent of exon size, and therefore could be dependent on intron–exon boundaries. All genomes examined have two distinctive classes of exons, each with different evolutionary histories. In the framework proposed in this article, these two classes of exons can be derived from a hypothetical ancestral genome by (spontaneous) symmetry breaking. We note that one of these exon classes comprises mostly alternatively spliced exons.