Skip to main content

Gene fusion as an important mechanism to generate new genes in the genus Oryza

Abstract

Background

Events of gene fusion have been reported in several organisms. However, the general role of gene fusion as part of new gene origination remains unknown.

Results

We conduct genome-wide interrogations of four Oryza genomes by designing and implementing novel pipelines to detect fusion genes. Based on the phylogeny of ten plant species, we detect 310 fusion genes across four Oryza species. The estimated rate of origination of fusion genes in the Oryza genus is as high as 63 fusion genes per species per million years, which is fixed at 16 fusion genes per species per million years and much higher than that in flies. By RNA sequencing analysis, we find more than 44% of the fusion genes are expressed and 90% of gene pairs show strong signals of purifying selection. Further analysis of CRISPR/Cas9 knockout lines indicates that newly formed fusion genes regulate phenotype traits including seed germination, shoot length and root length, suggesting the functional significance of these genes.

Conclusions

We detect new fusion genes that may drive phenotype evolution in Oryza. This study provides novel insights into the genome evolution of Oryza.

Background

New genes, which have been shown to be critical in the understanding of phenotypic evolution [1, 2], can be generated by a number of molecular mechanisms [3, 4], including gene duplication, de novo evolution, and gene fusion. Thus, investigation of the magnitude, origination, and corresponding evolutionary processes underlying the formation and fixation of new genes is critical to our understanding of the evolution of genome complexity that contributes to protein diversity. New genes formed by gene fusion can occur between 2 ancestrally neighboring genes, and sometimes lead to the evolution of novel complex domain structures [4,5,6,7]. In terms of the parental materials (DNA or RNA sequence) that produced to process, the fusion event results into fusion gene and retro-fusion gene, respectively.

Gene fusion was found to be a potent process for evolutionary novelties in bacteria. For example, in 18 bacteria species, including Escherichia coli, it was found that 0.2–1% of the genes in each genome were fusion genes [8]. However, it has been a challenge to understand the role of gene fusion in multicellular eukaryotic organisms. Early medical genetic studies revealed that some cancer-genesis mutations have arisen from the fusion of 2 adjacent genes in the genomes of cancer patients [9], implicating the role of gene fusion as a mechanism for new gene evolution, even though such fused genes may have deleterious consequences upon their birth. Analyses of the origination of fusion genes in humans and other hominoids have shown that they can form at both the DNA (i.e., alternative splicing site skip or mutation, transposon element related movement or recombination) and RNA levels (i.e., retroposition with subsequent flanking sequence recruitment [5, 9, 10]), thereby revealing the diverse molecular processes that can lead to gene fusions.

Whether or not gene fusion is a general mechanism for new gene evolution in plants remains unknown. However, recent sequencing efforts in several model organisms have predicted that gene fusion is likely an active molecular process in eukaryotes. For example, in human it has been shown that at least 4–5% of tandemly duplicated genes are transcribed into single putative fusion transcripts [11, 12]. With alternative splicing and deletion mechanisms, a fusion gene that encodes a chimeric protein can be generated during evolution [5]. This prediction was supported by the high proportion of chimeric proteins identified in several organisms [7, 10, 13], and although it was unknown how often the gene fusion mechanisms were invoked, ~30% of new genes recruited various genomic sequences and formed chimeric gene structures in Drosophila [14]. In rice, ~50% of new genes on the short arm of chromosome 3 may have been formed by chimeric mechanisms [15]. These chimeric new genes, which could have been formed by any number of the recombination-based molecular processes that result in gene fusion, inspired a search for the global distribution of gene fusion events across the genomes of a set of closely related species.

The major technical hurdles that have prevented the detection and study of gene fusions in the past have been (1) the lack of transcriptome data derived from multiple and recently diverged species; (2) the inclusion of low expressed transcripts; and (3) the lack of a suitable out-group to infer fusion gene origination processes. The low expression associated with recently evolved genes often makes the identification of new genes difficult, including those created by gene fusion mechanisms [16]. With the advent of high coverage RNA sequencing (RNA-seq), however, it recently became feasible to identify and investigate fusion genes at the transcriptome level in multiple species. The DNA and RNA sequencing of multiple Oryza species which have diverged over a recent evolutionary time frame provide an excellent opportunity to study new gene evolution through gene fusion.

Here we report a computational search method, assembled into a pipeline named GriffinDetector (scripts and documentation are available online: http://longlab.uchicago.edu/?q=GriffinDetector), which was designed to identify fusion genes in several recently diverged species. We applied this pipeline to interrogate 8 Oryza species/varietal groups and 2 additional out-group species, to detect the presence of fusion genes in four focus genomes Oryza sativa L. v.g. japonica, O. sativa L. v.g. indica, O. barthii, and O. glaberrima. O. glaberrima is a cultivated African rice species which originated from O. barthii that diverged from each other within 0.40 million years [17,18,19,20]. Using these comprehensive genomic resources and GriffinDetector, we detected a surprisingly high number of new gene fusion events that led to a high rate of fused protein formation in the recent history of Oryza, suggesting frequent adaptive events in these species are triggered by new fusion genes.

Results

Predicting new gene fusions

We developed a set of computational methods to process the genome sequence data of four recently diverged Oryza species, O. sativa L v.g. japonica, O. sativa L v.g. indica, O. glaberrima, and O. barthii. Generally, a fusion gene has at least 2 distinct parts originated from 2 independent parental genes. The algorithm used to identify fusion genes needs to resolve the major question of whether 2 different gene parts fuse to become one single gene. In single genome studies, analyses are limited to the identification of fusion genes through comparisons of the candidate gene to its paralogous genes [7] or to different stages of transcriptome data [21]. Without out-group support, it can be hard to distinguish gene fusion events from gene fission events—which occur when a longer homologous parental gene is split into 2 short homologous daughter genes. Benefiting from next-generation sequencing technology, it is possible to detect fusion genes based on multiple genome comparisons [4, 15]. Our pipeline, GriffinDetector, detects fusion genes using protein sequences annotated in several genomes and considers the phylogenetic relationships of homologous gene structures.

The basic design of GriffinDetector is to classify BLAST alignment hits of a query gene from a focus species (the species where we want to detect fusion genes) to several species into 2 groups, the long homologous gene group, i.e., potential fusion genes (e.g., two exons structure in Fig. 1 below the in-group), and the short homologous gene group, i.e., potential parent genes (e.g., one exon per gene in Fig. 1). If the query gene (e.g., from focus species A) has at least 2 short homologous gene hits across all species (all groups) while it has the long homologous gene only in the closely related species (species B and C in in-group in Fig. 1), then, most likely the short homologous genes have been integrated into one fusion gene along the lineage (gray dot in Fig. 1). Thus, we assess the query gene as a fusion gene in focus species A.

Fig. 1
figure 1

The basic design for detecting fusion genes. The three species in gray box have both long homologous gene hits (fusion gene candidates) and short homologous gene hits (parental genes), while the other six species in the mid-group and out-group only have the parental genes. Thus, we can deduce a fusion event (gray dot) occurred in the lineage leading to the in-group

GriffinDetector consists of 3 main parts. First, the species used in analysis (input data including out-group species, phylogeny, and focus species—which is defined in the control file) are categorized into three groups: “in-group,” “mid-group,” and “out-group” (Fig. 2 A–C). Next, all query genes in the focus species are aligned to the other species with BLASTP [22], followed by an alignment assortment into long homologous copies and short homologous copies for each species (Fig. 2 D, E) and group (Fig. 2 A, F). Finally, a query gene is classified as a candidate fusion gene if it meets additional threshold criteria as outlined in more detail in the “Methods” section.

Fig. 2
figure 2

GriffinDetector pipeline flowchart. The light blue ovals indicate the data input or final output data, the green diamonds indicate the thresholds used, the purple cylinders indicate the data generated, and the red boxes indicate data processing. A, F Whether the gene is present in the species /group or not. B, C All the species are categorized into three groups according to their phylogeny: species in the out-group are fixed based on the control file (red dashed box) while the other two groups are dynamic (green and yellow dashed box); the focus species is the initial species in the in-group, then, the closest node or clade will be added into the in-group gradually; the remaining species belong to the mid-group; the process will be stop when only one species remains in the mid-group. D, E When the species belonging to the three groups are clear, BLASTP hits from the focus species query gene are categorized into long hit copies or short hit copies for each species: hits having more than 80% sequence coverage of the query gene are recorded as long copies (line d); the remaining longest hits which have no overlap with each other indicate short copies (lines b and c); other hits (lines ai) are ignored

Abundance of fusion genes in the Oryza genomes

Using the OGE/IOMAP annotation dataset for O. sativa v.g. japonica, O. sativa v.g. indica, and Arabidopsis thaliana as the most distant out-group species within the phylogeny (Fig. 3, Zea mays and Sorghum bicolor are also used as out-groups in the pipeline’s robust analysis, more detail in the “Methods” section), we detected 157, 107, 122, and 100 candidate fusion genes (isoforms) in japonica, indica, O. barthii, and O. glaberrima, respectively (Additional file 1), which correspond to 114, 100, 109, and 82 genomic loci, respectively. Since some of the genes among these species may be orthologous to one another, we removed all redundant copies by comparing the candidate fusion genes via a BLASTP alignment [22], where hits having more than 80% identity over 100 bp or greater were classified as orthologous gene pairs. This resulted in a final candidate gene fusion list of 310 across the 4 focused genomes (Table 1, more details in Additional file 2, notes that the orthologous genes in different species are considered as one fusion gene), most of which were species-specific genes.

Fig. 3
figure 3

The phylogeny and relative annotation versions used in the analysis. Species in the light gray box are focus species, fusion genes are detected in all these four species for all six combinations (Additional file 2). The dark box represents the most distantly related species

Table 1 Fusion gene number detected in four focus species

A high rate of fusion gene origination

Of the 310 candidate fusion genes detected in this study, 58 were shared between 2 and 4 species by simple cross intersections (IDs in Additional file 2). However, when phylogeny was take into consideration, only 31 shared fusion genes could be properly matched onto the species topology tree, suggesting the episodic origination of these fusion genes (Fig. 3). The other 252 fusion genes were categorized as being specie-specific and must have originated very recently in evolutionary time (<1 million years, MY) [17,18,19]. Thus, the simple estimation of the fusion gene birth rate can be as high as 63 (252/4) per MY per species (Table 1 and Fig. 3). However, for individual species, they showed a great variation of fixation rate (Additional file 3). Analysis of resequencing data from 110 O. glaberrima accessions showed a mean fixation frequency of 93.3%, with 35 of the 43 O. glaberrima-specific fusion genes being 100% fixed (Additional file 3), while 94 O. barthii accessions held a mean fixation frequency of 21.5%, with only 2 genes being 100% fixed. On the other hand, by retrieving CDS from vcf files in the 3000 Rice Genomes Project (http://iric.irri.org/resources/3000-genomes-project), O. sativa v.g. japonica and O. sativa v.g. indica were estimated to have mean fixation frequencies of 31.8% and 15.4%, respectively (Additional file 3). If we take a frequency of 80% as our gene fixation threshold, then we estimate that the fusion gene fixation rate in rice can be as high as 16 (i.e., 64 species-specific fusion genes/4 species) new fusion genes per MY per species, which is nearly 1~2 orders of magnitude greater than that documented for flies and humans [23].

One third of fusion genes show specific expression patterns, while another one third of fusion genes exhibit similar expression patterns to their parental genes

To determine if some of the candidate fusion genes identified in this study are potentially functional, we analyzed 157 japonica fusion gene isoforms for evidence of gene expression and signatures of selection.

Using RNA-seq data from root, leaf, and mixed stage panicle tissues of O. sativa v.g. japonica, O. barthii, O. glaberrima, O. punctata, O. brachyantha, and L. perrieri [24], fragments per kilobase of exon per million fragments (FPKM) [25] values were retrieved for the 157 O. sativa v.g. japonica candidate fusion gene isoforms and their homologous copies (see “Methods” section, Additional file 4). First, we observed that the long homologous copies of 37 fusion gene isoforms (23.6%), identified in all six species, did not produce detectable FPKM signals. However, of the short homologous copies, all but one group of genes (related to fusion gene isoform Osjap05g14570.1) produced detectable FPKM signals (Additional file 4, Additional file 5: Figure S1). Next, among the 69 fusion genes (loci) that had detectable expression signals (FPKM value >0) in O. sativa v.g. japonica, we noticed 2 fusion genes (Osjap05g04120 and Osjap01g42770) that shared the same expression patterns as their short homologous copies in all three tissues (Additional file 6, Pearson’s r range from 0.864 to 0.996, p value range from 0.071 to 7.55E−05), 16 fusion genes that shared the same expression patterns as one of their short homologous copies in all three tissues with significant p values (Additional file 6, p<0.05), as well as an additional 6 with weak p values (Additional file 6, 0.05<p<0.1). Lastly, we selected the maximum FPKM value from three tissues for each species to represent the expression levels of both the long homologous and short homologous copies. The expression level from the former was significantly less than the latter (Additional file 7, Additional file 5: Figure S2, p=2.57 × 10−8); this narrowed expression level of new fusion gene is in accord with the phenomenon in de novo genes [26].

Based on FPKM expression data for the 69 fusion gene pairs, 24 gene pairs shared the same/similar expression patterns as their parental genes, 26 gene pairs had novel expression patterns, and the remaining 19 lacked data for their parental copies (Additional file 6). These results indicate that over one third of the candidate fusion genes may have undergone neofunctionalization or subfunctionalization, which has also been suggested by previous studies in human and Drosophila [7, 27]. Previous studies support that the new genes in Drosophila and human tend to move off their sex chromosomes and tend to be expressed in testis [23, 28,29,30]. Since the rice tribe does not have sex chromosomes, but does have reproductive tissues (i.e., panicles), we compared FPKM expression data of the fusion genes across three tissues and detected signals for 48, 27, and 52 genes in root, leaf, and panicle tissues, respectively (Additional file 6). There was no significant difference between panicle and root tissue (Fisher’s exact test, p = 0.4948), but the number of expressed fusion genes in leaf tissue were significantly lower than in root and panicle (p = 1.22 × 10−4 and 2.00 × 10−3). These expression patterns suggest that fusion genes may play specific roles in root and panicle biology as compared with leaf tissue.

Fusion genes may drive phenotype evolution

To further investigate whether these expressed fusion genes are associated with adaptive phenotypes, we chose a new gene that represented “neofunctionalization” (i.e., Osjap07g28390 have a novel expression pattern, whose expression level did not correlate with two parental genes), and another that represented “subfunctionalization” (i.e., Osjap09g15430, which has a highly similar expression pattern with that of one of its parental genes), for gene editing experiments using CRISPR/Cas9.

A total of 30 lines were obtained for the Osjap07g28390 gene editing experiment, of which 14 were shown to be edited, and included 3 homozygous mutants. One Osjap07g28390#1 edited line had a single-nucleotide deletion (A) for the sgRNA1 and a two nucleotide (AA) deletion at sgRNA2; transformant Osjap07g28390#6 contained a three nucleotide (AGC) deletion at the sgRNA1; while transformant Osjap07g28390#23 contained a large fragment deletion between the two sgRNAs (Fig. 4 A).

Fig 4
figure 4

Location of sgRNA on the chimeric diagram and editing informations of three homozygous lines for these two fusion genes. Chimeric structure was drawn based on pairwise blast results, two sgRNAs of each fusion gene were marked by green triangle and their edit types for three selected lines were aligned chromatographically

A total of 29 putative transgenic lines were generated for the Osjap09g15430 fusion gene, of which 13 were shown to be edited. Selected homozygous lines are as follows: Osjap09g15430#13 contained a 5-nucleotide deletion at sgRNA1 and a 15-nucleotide deletion at sgRNA2; Osjap09g15430#19 contained a 2-nucleotide (TT) insertion at sgRNA1, and single-nucleotide (A) deletion at sgRNA2; while Osjap09g15430#26 contained a one-nucleotide (G) insertion at sgRNA1 and one-nucleotide (T) insertion at sgRNA2 (Fig. 4 B).

We found that the germination rate and root length of the three Osjap07g28390 lines was much lower that of the same Nipponbare line and that root and shoot length of the three Osjap09g15430 edited lines were shorter (Fig. 5) than the same Nipponbare line grown under identical growing conditions. These results are consistent with the expression preponderance in panicle and root. When we searched these two genes against rice stress database BAR (http://bar.utoronto.ca/efprice/cgi-bin/efpWeb.cgi), they both mainly response to drought stress (Additional file 5: Figure S3). As reported in other species, drought often reduced germination rate, root length, and shoot length [31, 32]; thus, the shorter root and lower germination in CRISPR mutant possibly reflect a result of drought adaptation.

Fig. 5
figure 5

Phenotypic comparison between CRISPR/Cas9 knockout mutants and wild type plants. A Morphological characteristics donated by knockout of two fusion genes, the photo were taken at the fifth day after sowing. B Germination rate of different genotype. C Root length of germinated seedlings. D Shoot length of germinated seedlings. All these traits were observed at five days after sowing, with 3 homozygous lines of each fusion gene showed similar phenotype. The same Nipponbare line was used as control. The statistical analysis between transgenic lines and Nipponbare was done by Wilcox test. Asterisks indicate significant at p< 0.05

Duplication of parental gene subsequent with fusion event is the main origination structural pattern

To better understand how gene fusions that are formed in rice, we investigated the structures of the 69 japonica expressed fusion genes we identified. First, the 69 expressed fusion genes were categorized into 3 major types based on the parental gene presence and locus in the genome of the focus species [33]. Type I fusion genes have both parental genes that are detectable (Type I) and are formed under the following 3 scenarios where the parental genes are located in the same region of the chromosome, and a segmental duplication generated a novel copy which was then fused together by alternative splicing site mutations or skip (Fig. 6 A); the parental genes are located in different regions of the genome, but duplicate copies are generated, inserted into an adjacent region, and then fused (Fig. 6B); or, one of the parental genes generated a tandemly duplicated copy while the other parental gene (unlinked) generated a duplicate copy that inserted nearby and then fused (Fig. 6C). Type II fusion genes are ones when only one parental gene is detectable. Here, one of the parental genes generates a novel duplicated copy that is inserted adjacent to another parental gene, which has not been duplicated and they fuse together (Fig. 6 D); or one of the parental genes generates a tandemly duplicated copy, and another parental gene, elsewhere in the genome, is moved to an adjacent region and then they fuse together (Fig. 6E). (3) Type III fusion genes are ones where no original parental gene is detectable. Here, one parental gene is transposed to an adjacent region of another parental gene, and then becomes a fusion gene (Fig. 6G), or the parental genes are located in the same region, and an alternative splicing site mutation or skip leads to the fusion of the 2 parental genes (Fig. 6F).

Fig 6
figure 6

Origination patterns of fusion genes. A–C Parental copies are detectable in focus species (Type I); D, E only one parental copy is detectable (Type II); F, G none of the parental genes are detectable (Type III)

Following the guidelines above, we categorized the 69 expressed fusion genes into 35 Type I, 19 Type II, and 15 Type III (Additional file 8). In Type I, 2 of 35 fusion genes had parental genes in the same region in the focus species (Fig. 6A), 4 were located near one of the parental genes (Fig. 6C), and the other 29 were located in different regions of the genome (Fig. 6B). Of the 19 Type II fusion genes, only one was located in the adjacent region to one parental gene (Fig. 6E) while the others followed the hypotheses outlined in Fig. 6D. Lastly, of the 15 Type III fusion genes, and considering the location of parental genes in other species, 11 of them had the parental homologous copies in the same syntenic regions (Fig. 6G), while parental genes in the remaining 4 cases were located in different chromosomal regions (Fig. 6F).

Our data suggests there are 13 fusion genes that have only one short homologous copy in L. perrieri and A. thaliana while they have two parental copies in other species (Additional file 8). This result reveals that some of these parental genes are also newly evolved. Since the lack of homologous copies that have more than 20% coverage of the fusion gene in L. perrieri and A. thaliana, these young parental genes may have originated de novo. These results are consistent with our previous works that showed that the genus Oryza has a high rate of de novo gene [26] origination and the neofunctionalization in gene fission and gene fusion is often dependent upon a duplication of parental gene [33, 34].

Discussion

Here, we presented a new pipeline, GriffinDetector, to scan for candidate fusion genes in focus species utilizing several closely related species, a given phylogeny, and three group species. We then applied GriffinDetector to the OGE/IOMAP rice tribe genome data set, which showed to be on the top of published genome data [24] especially for the analysis of new gene evolution [26]. We detected 310 fusion genes in 4 focus species, O. sativa L. v.g. japonica, O. sativa L. v.g. indica, O. barthii, and O. glaberrima (Additional file 2). The abundance of fusion genes detected in these AA genome species reveals a high birth rate of functional fusion genes, of which around 1/3 share the same expression pattern with at least one parental gene, and another 1/3 developed totally new expression patterns. Further, almost 1/2 of the fusion genes detected were formed first by gene duplication followed by gene fusion. CRISPR/Cas9 knockout mutants of 2 randomly selected fusion genes demonstrated that fusion genes contribute to phenotypes such as root length and seed germination, irrespective of expression novelty, which implies that fusion genes may play important roles in rice genome evolution even if they share redundancy with their parental genes. As Ohno (1970) pointed out, gene duplication serve as the most important mechanism for origination of new gene copies. Compared with other patterns of new gene like orphan gene, fusion gene derived from gene duplication or exon shuffling could make use of existing nucleotide sequences and was estimated to affect a wider range of phenotypic outcomes for the potential greater genetic diversity [13]. Works in humans and Drosophila also suggested that ~80% of new genes are formed by DNA-based duplication, 5–10% by de novo duplication, and ~10% by retroposition [4].

As a pipeline, GriffinDetector has several novel attributes compared with previous analytical tools used to detect gene fusions [35,36,37]. First, GriffinDetector is the first pipeline that is able to integrate data from several genomes to scan for fusion genes automatically. Second, GriffinDetector allows the use of a maximum amount of genomic information to make judgments about whether a query gene is absent from one group or not. This strategy largely reduces the number of false positives that can arise from imperfections of an assembly or annotation. Third, since GriffinDetector works on several closely related genomes, possible biased data present in one genome will not have a significant effect on the final results, unless the bias is in the focus species. Lastly, our results show that GriffinDetector has great potential to become a gateway tool that will inspire innovative projects for those working on integrated analyses of multiple, closely related, genomes to detect and interrogate fusion genes.

GriffinDetector is an efficient pipeline. Our analyses indicates that changing the out-group species, or annotation of a sister species, will not significantly affect the results, thereby demonstrating the robustness of the detection methodology (see “Methods”). We observed that the vast majority (i.e., 76.7%) of the fusion genes detected by GriffinDetector were shared among different out-group and non-focus annotation combinations (Table 2, Fig. 7, “Methods”). However, we also observed that utilization of alternative annotations of focus genomes had a significant impact on the detection of fusion genes (Table 2). For example, when we focused on the detection of fusion genes in O. sativa L. v.g. japonica rice, and then switched the OGE/IOMAP-MAKER annotation to the MSU or RAP-DB annotation, the fusion gene number decreased from 157 to 24 or 7, respectively. Such large variations are not unexpected. For example, when Wang et al. [38] and Sakai et al. [39]  scanned the indica and japonica genome assemblies for retrocopies, they detected 1235 and 150, respectively. One major reason is that the public annotations, like the MSU- and RAP-DB annotations, tend to use very strict parameters to maintain a conservative gene pool, which was shown to be unsuitable for new gene analysis [16]. In our analysis, we chose the OGE/IOMAP annotations for all subsequent analysis, since most of the other species used in the analysis were annotated in a similar manner [24, 26] which used the newest MAKER-P [40, 41] annotation pipeline to generate a reliable and consistent annotation data set.

Table 2. Number of candidate fusion gene isoforms
Fig. 7
figure 7

Comparisons of candidate fusion gene isoforms among six combinations. Three-letter vectors denote the various annotations and most distantly related species. Letters in the first position: O represents OGE/IOMAP, U for MSU, R for RAP-DB annotations. In the second position, O stands for OGE/IOMAP, and G for GLEAN; Letters in third position: T stands for TAIR, S stands for Sorghum bicolor, and Z stands for Zea mays. Venn diagrams were used to display the number of fusion genes that were shared or unique across different species or annotations. A, B Most distantly related species shift, A sativa v.g. japonica IDs, B glaberrima IDs, 155 and 99 candidate fusion gene isoforms are shared by three combinations respectively; C, E sativa v.g. indica annotation shift, C sativa v.g. japonica IDs, E glaberrima IDs, 150 and 94 candidate fusion gene isoforms are shared by two combinations respectively; D sativa v.g. japonica annotation shift, D glaberrima IDs (sativa v.g. japonica IDs are not comparable), 92 candidate fusion gene isoforms are shared by the three combinations

High rate of fusion gene origination in Oryza

Many groups have demonstrated that fusion genes play an important role in eukaryotic new gene origination which has led to major contributions of adaptive evolutionary novelties [42]. Here, we demonstrated that several species in the Oryza genus exhibit high rates of fusion gene origination.

The total number of fusion genes detected by GriffinDetector was 310 among 4 Oryza species (Additional file 2) which diverged very recently (<1 MY) [17,18,19]. Excluding fusion genes that may have originated earlier, and fusion genes that may not be fixed in the population, we estimated that the fixation rate of fusion genes in Oryza was as high as 16 fusion gene per MY per species. This rate is much higher than that found in other eukaryotic systems such as Drosophila, where it has been estimated to be fixed at around 0.16 fusion genes per MY per species [7]. Among them, O. glaberrima showed an especially high fixation frequency (Additional file 3), which is consistent with the fact that O. glaberrima underwent a genetic bottleneck effect leading to its low level of genetic diversity [43,44,45,46].

Among the 69 fusion genes in O. sativa L. v.g. japonica that have an expression signal in the tissues tested, only 2 cases were confirmed to have an intron-loss structure (Fig. 8), Osjap01g42770 was shared with indica, while Osjap11g25180 was japonica specific. Thus, we estimate the retro-fusion gene birth rate in rice to be around 2 per MY (1.77 per MY if weighted, 1 + 1 × 33/43 = 1.77), which is very close to 1.7 (197/(1.5/2/6.5 × 103) = 1.7) per MY estimated in a previous study [38]. This rate is more than 10 times the rate found in human and Drosophila, which was estimated to be 0.14 [23, 30] and 0.16 [47] per MY, respectively, and two times the rate in zebrafish, estimated at 1 per MY [48]. It is interesting to note that the estimates of birth rates of new genes seem to increase over time. Our rate (16 fusion genes per MY and 2 retro-fusion genes per MY) is far higher than previous studies [7, 23, 30], but close to the results reported for rice [38] and zebrafish [48]. This phenomenon can be explained in two ways: (1) due to a fast decay rate, many old fusion genes will be eliminated from the genome, while very young fusion genes can still be observed. Both our data and previous studies support this hypothesis. Our data suggests that among the 4 focus species which diverged within 1 MY, at least 11 fusion gene loss events have occurred (Table 1). Previous studies indicate that the estimated half-life of fusion genes was 0.44 MY [7], which means that older fusion genes have a much higher chance of being deleted. The distributions of fusion genes according to Ks values (or age) clearly show that more fusion genes are enriched at a young age [7, 48]. Wang et al. [38] estimated the birth rate of retro-fusion genes using young genes which were less than 1 MY (Ks ≤ 0.013), with a rate as high as 7 genes per MY, but when all the 197 genes are used (the Ks ≤ 1.5, which indicate the oldest age is 115.4 MY), this rate decreases to 1.7 gene per MY. (2) The mechanisms or strategies for genome evolution of plants, primates, fish, and insects may have large differences, so fusion genes might play a more important role in rice and zebrafish compared with their roles in human or Drosophila. Our recent work suggests that plants do have large differences in gene structure evolution as compared to human and fly systems [49], which therefore may contribute to the differences we observed in gene fusion birth rates.

Fig. 8
figure 8

Two retro-fusion gene structures. Osjap0312940 was retroposed to be the second exon of Osjap11g25180; Osjap01g08100 was retroposed to be the first exon of Osjap01g42770, while Osjap03g27560 was retroposed to become the second and third exon. The second intron of Osjap01g42770 may be retained from the parental gene or is newly evolved. Tandem repeats were detected for both gene pairs (dash lines)

The patterns of fusion gene formation

The progenitor of cereal crops was shown to be an ancient aneuploid that experienced a whole genome duplication ~70 MY ago [18]. It has been suggested that genome content and gene expression quickly change after whole genome duplication events (WGDs) [50]. Since all the fusion genes we detected were very young, most of them being less than 1 MY old, these fusion genes have likely originated on a case-by-case basis and are not related to this large-scale duplication event(s).

In our analysis, we classified the formation of fusion genes into 3 categories: (type I) both parental copies are detectable in the focus species; (type II) only one parental copy is detectable in the focus species; and (type III) none of parental short copies are detectable. Overall, these 3 categories correspond with 3 types of molecular evolution events: (1) sequences with in the fusion gene are copies of parental genes; (2) one part of a fusion gene sequence was from a parental gene duplication, and the other part of the fusion gene was a parental gene, itself; and finally, (3) the fusion gene sequence was a fused unit of 2 parental genes. We identified 35 fusion genes belong to the type I category, 19 to the type II category, and 15 to the type III category. Although many previous studies have pointed out that the rate of duplicated gene loss is lower than expected, and several theories have been proposed to explain these lower rates of gene loss [51], we cannot guarantee that all fusion genes detected in this study are indexed in the correct group, possibly because of visible recent gene loss in the focus species. For example, Additional file 2 and Fig. 6 show that fusion gene loss events have occurred 11 times. However, this obviously does not change our estimate for the high origination rate of new fusion genes, given the high number of young fusion genes that were detected herein.

It is well known that the two primary mechanisms responsible for the majority of gene fusion events in eukaryotes are DNA-level recombination [52] and retroposition [23, 30, 38]. Our Oryza fusion gene data did not reveal the presence of a large amount of TEs (23 in total, Additional file 9) like LINEs, LTRs, Pack-Mules, and Helitrons, which previous studies have suggested may contribute to fusion gene origination [53,54,55]. Here, we only detected 2 cases of retro-fusion genes (Fig. 8) out of a total of 69 fusion genes which have transcription evidence. The low retro-fusion gene birth rate observed in our study coincides with previous studies [56] and may be due to the finding that retrogenes favor the recruitment of promoter regions, and not additional protein coding genes/exons, en route to fusion gene formation [4]. We note here that GriffinDetector was designed specifically to focus on protein coding gene fusions.

Williford and Betrán recently classified the formation of a fusion gene into two phases [57]. First, sequences from two sources were placed in adjacent regions, and second, sequence insertion, deletion, mutation, or recruitment of noncoding regions occurs, which provides suitable splicing sites for the two sequence sources to become a fusion gene. This simple but important statement may help us to explain the data we observed, such as the case of Osjap04g02830 (Additional file 8), where the parental genes in O. sativa L. v.g. japonica were located near to one another, while in the other species they were not, which means the parental gene may have been re-located recently and, in a very short time, allowed for a new fusion gene to form.

In an analysis of the formation of fusion genes in Drosophila, fusion genes were classified into 3 major types: retro-fusion genes, ectopic recombination fusion genes, and tandem duplicated fusion genes [13]. These types are supported by previous studies: Rogers and colleagues detected 14 fusion genes, 8 of them related to the tandem duplication mechanisms; Shuang Yang and Arguello et al. detected 17 duplicated genes which were generated through ectopic recombination in a 12 MY time frame, most of which are functional and evolved diverse functions and chimerical structural [58]. In contrast, our data showed very little evidence of tandem duplication fusion genes (Additional file 8), which raises interesting new questions of origination to explore in future studies.

Phenotypic effects of new genes and biological diversity

Phenotypic evolution in morphology and development leads to biological diversity. Whether or not the gene evolution participates in the evolution of genetic basis underlying phenotypic evolution is a new direction to explore, besides the better known role of regulatory systems for genes [3]. Huang et al. [59] show that a species-specific gene in Arabidopsis thaliana evolved important phenotypic effects in development and morphology. In this study, we used the CRISPR knockout technique to create mutant lines for the fusion genes Osjap07g28390 and Osjap09g15430. We found that both genes evolved significant phenotypic effects in three phenotypic traits, germination rates, root length, and shoot length. Regardless of possible neofunctionalization (Osjap07g28390) or subfunctionalization (Osjap09g15430) in molecular functions, they both are involved in drought resistance, an important trait adaptive to frequent environmental changes where rice often encountered. These data reveal that the genetic basis of the morphology and organs in Oryza evolved rapidly by acquiring new genes recently.

Conclusions

We detect new fusion genes that may drive phenotype evolution in Oryza. This study provides novel insights into the genome evolution of Oryza.

Methods

Development of GriffinDetector

GriffinDetector is a pipeline script written in Perl and requires two assistant programs (makeblastdb and blastp) from the BLAST+ package [22]. GriffinDetector is available at http://longlab.uchicago.edu/?q=GriffinDetector with detailed documentation, which also can be download at GitHub repository (https://github.com/zhangcj2022/GriffinDetector).

The GriffinDetector pipeline uses a control file to set the main parameters including information regarding the necessary input data (you can get the example in package). It divides all of the species listed in the control file into 3 groups: in-group, mid-group, and out-group. The out-group is fixed (e.g., Fig. 2C, the three species in red dashed box), while the in-group and mid-group are dynamic (e.g., Fig. 2C, the species in green and yellow dashed boxes). In Fig. 2C, from left to the right, the number of species belonging to the in-group gradually increases, while the number of species in mid-group decreases. In the initial round, GriffinDetector will take one species from the focus species listed in the control file, as the initial in-group (Fig. 2C 1st) and in the following rounds, GriffinDetector will add one species (Fig. 2C 2nd, 3th, 5th) or one clade (Fig. 2C 4th) which is the nearest one to the current in-group clade into the in-group to generate a new in-group. This dynamic progress will not stop until only one species remains in the mid-group. Thus, GriffinDetector sort alternative combinations of species with associated with each focus species.

Based on the in-group combinations, protein sequences of each focus species are aligned using BLASTP [22] to protein sequences of other species (including itself). GriffinDetector will inspect the BLASTP results and categorize the hits into short homologous copies and long homologous copies (Fig. 2D, E). Hits that have more than 80% coverage of the query gene will be assigned as “long homologous copies” immediately (Fig. 2E line d), then GriffinDetector will compare the remaining hits to each other and filter the shorter hits (Fig. 2E line ai) which overlap with other longer hits (Fig. 2E line c). Finally, hits that overlap more than 20% with the query gene (Fig. 2E line b & c) will be assigned as short homologous copies. If the final remaining short homologous copies in one species are more than two (including two) and the total coverage of the query gene is larger than 80%, GriffinDetector identifies that species as having the short copies. Thus, GriffinDetector is able to determine whether each species has the long homologous copy or not, as well as short copies.

If one homologous copy is present in one group, it will be denoted as 1; otherwise, it will be denoted as 0. GriffinDetector defines the homologous copy as present in one group if the homologous copy is present in all species of the group, but if the group has two or more species, the homologous copy is allowed to be absent in one species at most (Fig. 2A). For example, if the query gene has a long homologous copy present in the in-group and mid-group, but is absent from the out-group, we use a vector to denote it as L-110. If the short homologous copies are present in all three groups, we record it as S-111 (e.g., Fig. 1). Generally, we have 23 combinations for each long homologous copy and 23 combinations for each short homologous copy, where the total number of possible combinations is 64. GriffinDetector only recognizes three combinations as fusion genes: L-110 and S-111, L-100 and S-111, L-100 and S-110.

Data and resources

In our analysis, we used 12 species and 15 annotation versions. All related data, including protein sequences and GFF files of Oryza species, can be accessed through OGE/IOMAP website (http://oge.gramene.org/, more detail in http://bioinfor.kib.ac.cn/?q=node/10), which are required by the pipeline. The fragments per kilobase of transcript per million mapped reads (FPKM) values were downloaded from the OGE/IOMAP website, as well as the TE annotation of O. japonica.

RNA-seq data from leaf root and mixed stage panicle tissues were as described in Stein et al. [24]. The RNA-seq reads were aligned to the respected genomes using Tophat software [60]. The alignment files were inputted into Cufflink software [61] to generate the FPKM value.

The OGE/IOMAP annotation was generated using the MAKER annotation pipeline v2.2.8a [40, 41] as describe in Stein et al. [24]. The de novo transcript assembly, reference guide transcript assembly, the cDNA, and protein annotations from O. sativa L v.g. japonica, O. glaberrima, Brachypodium distachyon [62], and Oryza cDNA database (NCBI taxonomy ID 4527) were used as evidence for annotation.

Efficiency of GriffinDetector

Since GriffinDetector is the first automated pipeline designed to detect fusion genes across multiple genomes, it is important to determine its efficiency in detecting such genes. Here, we applied the GriffinDetector to the rice tribe and its closely related out-group species to estimate its robustness. First, we set O. sativa L. v.g. japonica, O. sativa L. v.g. indica, O. barthii, and O. glaberrima as focus species respectively; O. brachyantha, Leersia perrieri, and Arabidopsis thaliana as out group species; and O. rufipogon, O. longistaminata, and O. punctate as other species, to search for fusion genes in the 4 focus species respectively based on their phylogeny (Fig. 3). In addition, we utilized different versions of annotations for the japonica and indica genome, and adjusted the most distantly related out-group from A. thaliana to Sorghum bicolor and Zea mays to estimate the influences of annotation and out-group distance on gene fusion detection, respectively. For the other 7 species, we only used the OGE/IOMAP MAKER annotations [24]. We use three-letter vectors to denote the 6 different combinations of annotations and out-groups utilized for these analyses: O-O-T, O-O-S, O-O-M, U-O-T, R-O-T, O-G-T (Table 2). The letter in the first position indicates the annotation version for japonica: O = OGE/IOMAP, U = MSU [63, 64], and R = RAP-DB [64, 65]. The second letter position indicates the annotation version for indica: O = OGE/IOMAP and G = GLEAN [66]. The third letter position indicates the most distantly related out-group species: T = A. thaliana TAIR annotation [67], S = S. bicolor [68], and Z = Z. mays [69].

The fusion gene data sets identified using these 6 combinations were compared to each other (Fig. 7), and we reached the following conclusions: (1) changing the most distantly related out-group species did not significantly change the results (i.e., >95% of the fusion genes detected were shared); (2) changing the annotation of which is close to focus species due to the phylogeny will slightly affect the final results (i.e., >76.7% of the fusion genes were shared); and (3) changing the annotation version of the focus species will greatly affect the results (i.e., the number of fusion genes detected decreased from 157 to 7). In Table 2, columns O-O-T, O-O-S, and O-O-M, which display the results from the analyses in which the most distantly related species was changed to either A. thaliana, S. bicolor, or Z. mays, the final number of candidate fusion gene isoforms barely changed in japonica, indica, O. barthii, and O. glaberrima. A more detailed comparison between candidate fusion gene isoforms of japonica (Fig. 7A) and O. glaberrima (Fig. 7B) indicates that most (>95%) of the detected fusion genes are shared among these three combinations, which means changing the most distantly related out-group species hardly has any effect on the fusion genes identified in each focus species. The O-O-T, U-O-T, and R-O-T columns, which differed in the japonica annotation used, showed a slight change in the number of fusion gene candidates detected for indica (N = 6), O. barthii (N = 15), and O. glaberrima (N = 13). Similar results were seen by varying the indica genome annotations (i.e., columns O-O-T and O-G-T), where the number of fusion gene candidates in japonica, O. barthii, and O. glaberrima changed by 15, 8, and 4. A more detailed comparison between candidate fusion gene isoforms of japonica (Fig. 7C) and O. glaberrima (Fig. 7D, E) indicate that over three fourths (>76.7%) of the results are shared among these combinations; this rate (>76.7%) is lower than the previous one (>95%) and indicates that changing closely related species (e.g., while treating japonica as focus species which means we focus on detecting fusion genes in japonica, but changing the indica annotation from Version A to B, the final results of detected fusion genes will be different) annotations will somewhat affect the final results. On the other hand, when the japonica annotation changed from OGE/IOMAP to MSU to RAP-DB, the candidate fusion gene isoforms in japonica decreased from 157 to 24 to 7. Further, when the indica annotation changed from OGE/IOMAP to GLEAN, the candidate fusion gene isoforms in indica increased from 107 to 178. Both results show that the annotation used plays an important role in fusion gene detection.

Since the annotations of most species used were generated by OGE/ IOMAP (Fig. 3), we decide to use the O-O-T combination for all subsequent analyses which yield representative results as described above.

The origination rate and fixation frequency of fusion genes

Origination of Oryza species were based on species-specific fusion gene. For the calculation of fixation frequency on species-specific fusion genes, we collected O. glaberrima and O. barthii population SRA data from Meyer et al. [20] and Wang et al. [45], and the SNPs were called by Genome Analysis Toolkit (GATK 4.1.9.0) [70]. The variation of coding sequence was retrieved by vcf2fasta (https://github.com/santiagosnchez/vcf2fasta). On the other hand, we directly download vcf files for 100 O. japonica individuals and 101 O. indica individuals from Rice SNP-seek database (https://snpseek.irri.org/). Based on the re-annotated gff of species-specific fusion genes in four focus species (Additional file 10), we extracted coding sequence (CDS) of each fusion gene.

We mapped the SNPs on the reference and checked the exact sites in the CDS region, if the SNPs changed the frame of the CDS or it become a stop codon, then we defined the fusion gene was not fixed in the individual. Thus, we counted the total individual numbers of fusion genes not fixed among each species and defined the fixation frequency as the percentage of fixed individual to total individual number in each species. Finally, we define the fixation rate as the total fixed fusion gene number against species number.

RNA-seq data analysis

RNA-seq data derived from three tissues (i.e., leaf, root, mixed state panicles) from 6 species: japonica, O. barthii, O. glaberrima, O. punctata, O. brachyantha, and L. perrieri, were generated by the OGE/IOMAP project [24]. Using the same reference genome coordinate, if the region for an FPKM value had more than 50% sequence overlap to a fusion gene or homologous short copies, the FPKM value was used as the expression signal measurement for the relative genes. The FPKM values for the 157 fusion genes (isoforms) in japonica, as well as for all the homologous short copies, are presented in Additional file 4. The reference gene LOC_Os08g03290, which was recommended by previous work [71], was also used as a control to present the expression levels. Among these, 69 fusion genes (loci) have expression signals and were used in the following analysis.

Using the natural logarithm of the FPKM values and two R libraries, limma and ggplot2, we generated the heat maps (Additional file 5: Figure S4) to visualize the gene expression data employed in the analysis. This reference gene was the only gene studied previously in two separate analyses [71, 72] and has expression signals in all 3 tissues in japonica. For cases where there was more than one homologous copy in one species, we took the highest FPKM value as the representative for the expression level of the short homologous copies.

Knockout mutant generated by CRISPR/Cas9

Vector system of CRISPR/Cas9 was kindly obtained from Professor yaoguang Liu’s lab. Two candidate sgRNA were designed for each selected fusion gene. The two sgRNA were tandemly fused by paired adapters Uctcg-B1’&gRctga-B2 and Uctga-B2’&gRcggt-BL, and taken effect under the control of rice promoter OsU3 and OsU6a, respectively. All these constructs were introduced into E. coli and confirmed by sequencing. Positive single clone of each fusion gene construct were rejuvenated and used for plasmid extraction. Qualified plasmids were introduced into Agrobacterium tumefaciens EHA105 with electroporation method, and then transformed into Japonica rice Nipponbare for genome editing.

For the Agrobacterium-mediated transformation in Japonica rice variety Nipponbare efficiently, tissue culture-based method was adopted. Embryonic callus cells induced from the mature seeds scutella serve as starting material (Additional file 5: Figure S5A), then high potential ones were inoculated onto subculture medium for callus propagation for 20 days (Additional file 5: Figure S5B). After shaking one single colony Agrobacterium harboring CRISPR/Cas9 vector in suspension medium with antibiotics, glucose, and acetosyringone, target callus source co-cultivated with the fresh activated agrobacterium for transformation and then dried callus particles were transferred to co-culture medium for another 3 days (Additional file 5: Figure S5C). The transformed callus was rinsed with ddH2O and subjected to resistance screen on selective culture medium containing carbenicillin and hygromycin B (Additional file 5: Figure S5D). The potential positive transformants were further selected for proliferation on differentiation medium (Additional file 5: Figure S5E). Finally, the regenerated seedlings were processed for rooting before transplantation (Additional file 5: Figure S5F).

The mutants potentially edited by Cas9 protein were first screened by hygromycin gene-specific primer. Editing information of T0 transformants were obtained by PCR sequencing with gene-specific primers flanking sgRNA, and the homozygous lines were obtained by analyzing each chromatograms against reference sequence in DSDecodeM (http://skl.scau.edu.cn/dsdecode/). Seeds of homozygous transformants were propagated to T2 generation for phenotype observation.

Availability of data and materials

All the data feed into pipeline was available at website http://bioinfor.kib.ac.cn/?q=node/10. The raw RNA-seq data are available in NCBI under the accession number PRJNA13141, PRJNA13765, PRJNA13770, PRJNA163065, PRJNA30379, and PRJNA70533. The generated vcf files for O. barthii and O. glaberrima population data were deposited at https://data.cyverse.org/dav-anon/iplant/home/zhouyanli11/ricedata. Other datasets generated during the current study are included in this article and its supplementary information files. Additional information and resources in the manuscript are available from the corresponding authors on reasonable request. The codes of GrrifinDetector are available under a GPL-3.0 License on GitHub https://github.com/zhangcj2022/GriffinDetector.git [73] and Zenodo [74].

References

  1. Burns DM, Horn V, Paluh J, Yanofsky C. Evolution of the tryptophan synthetase of fungi. Analysis of experimentally fused Escherichia coli tryptophan synthetase alpha and beta chains. J Biol Chem. 1990;265(4):2060–9. https://doi.org/10.1016/S0021-9258(19)39940-5.

    CAS  Article  PubMed  Google Scholar 

  2. Lang D, Thoma R, Henn-Sax M, Sterner R, Wilmanns M. Structural evidence for evolution of the beta/alpha barrel scaffold by gene duplication and fusion. Science. 2000;289(5484):1546–50. https://doi.org/10.1126/science.289.5484.1546.

    CAS  Article  PubMed  Google Scholar 

  3. Chen SD, Krinsky BH, Long MY: New genes as drivers of phenotypic evolution(vol 14, pg 645, 2013). Nat Rev Genet 2013, 14:745-745.

  4. Long MY, VanKuren NW, Chen SD, Vibranovski MD. New gene evolution: little did we know. Annu Rev Genet. 2013;47:307–33.

    CAS  Article  Google Scholar 

  5. Long MY: A new function evolved from gene fusion. Genome Res 2001, 11:308-308, 11, 1657, DOI: https://doi.org/10.1101/gr.165700.

  6. Akiva P, Toporik A, Edelheit S, Peretz Y, Diber A, Shemesh R, et al. Transcription-mediated gene fusion in the human genome. Genome Res. 2006;16(1):30–6. https://doi.org/10.1101/gr.4137606.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  7. Rogers RL, Bedford T, Hardl DL. Formation and longevity of chimeric and duplicate genes in Drosophila melanogaster. Genetics. 2009;181(1):313–22. https://doi.org/10.1534/genetics.108.091538.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  8. Snel B, Bork P, Huynen M. Genome evolution - gene fusion versus gene fission. Trends Genet. 2000;16(1):9–11. https://doi.org/10.1016/S0168-9525(99)01924-1.

    CAS  Article  PubMed  Google Scholar 

  9. Thomson TM, Lozano JJ, Loukili N, Carrio R, Serras F, Cormand B, et al. Fusion of the human gene for the polyubiquitination coeffector UEV1 with Kua, a newly identified gene. Genome Res. 2000;10(11):1743–56. https://doi.org/10.1101/gr.GR-1405R.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  10. Zhang Y, Lu SJ, Zhao SQ, Zheng XF, Long MY, Wei LP. Positive selection for the male functionality of a co-retroposed gene in the hominoids. BMC Evol Biol. 2009;9(1). https://doi.org/10.1186/1471-2148-9-252.

  11. Parra G, Reymond A, Dabbouseh N, Dermitzakis ET, Castelo R, Thomson TM, et al. Tandem chimerism as a means to increase protein complexity in the human genome. Genome Res. 2006;16(1):37–44. https://doi.org/10.1101/gr.4145906.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  12. Goodstadt L, Ponting CP. Phylogenetic reconstruction of orthology, paralogy, and conserved synteny for dog and human. PLoS Comput Biol. 2006;2(9):e133. https://doi.org/10.1371/journal.pcbi.0020133.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  13. Rogers RL, Hartl DL. Chimeric genes as a source of rapid evolution in Drosophila melanogaster. Mol Biol Evol. 2012;29(2):517–29. https://doi.org/10.1093/molbev/msr184.

    CAS  Article  PubMed  Google Scholar 

  14. Zhou Q, Zhang G, Zhang Y, Xu S, Zhao R, Zhan Z, et al. On the origin of new genes in Drosophila. Genome Res. 2008;18(9):1446–55. https://doi.org/10.1101/gr.076588.108.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  15. Zhang CJ, Wang J, Marowsky NC, Long MY, Wing RA, Fan CZ. High occurrence of functional new chimeric genes in survey of rice chromosome 3 short arm genome sequences. Genome Biol Evol. 2013;5(5):1038–48. https://doi.org/10.1093/gbe/evt071.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  16. Zhang YE, Landback P, Vibranovski M, Long MY. New genes expressed in human brains: implications for annotating evolving genomes. Bioessays. 2012;34(11):982–91. https://doi.org/10.1002/bies.201200008.

    CAS  Article  PubMed  Google Scholar 

  17. Tang L, Zou XH, Achoundong G, Potgieter C, Second G, Zhang DY, et al. Phylogeny and biogeography of the rice tribe (Oryzeae): evidence from combined analysis of 20 chloroplast fragments. Mol Phylogenet Evol. 2010;54(1):266–77. https://doi.org/10.1016/j.ympev.2009.08.007.

    CAS  Article  PubMed  Google Scholar 

  18. Wang XY, Shi XL, Hao BL, Ge S, Luo JC. Duplication and DNA segmental loss in the rice genome: implications for diploidization. New Phytol. 2005;165(3):937–46. https://doi.org/10.1111/j.1469-8137.2004.01293.x.

    CAS  Article  PubMed  Google Scholar 

  19. Ma JX, Bennetzen JL. Rapid recent growth and divergence of rice nuclear genomes. Proc Natl Acad Sci U S A. 2004;101(34):12404–10. https://doi.org/10.1073/pnas.0403715101.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  20. Meyer RS, Choi JY, Sanches M, Plessis A, Flowers JM, Amas J, et al. Domestication history and geographical adaptation inferred from a SNP map of African rice. Nat Genet. 2016;48(9):1083–8. https://doi.org/10.1038/ng.3633.

    CAS  Article  PubMed  Google Scholar 

  21. Kaessmann H, Vinckenbosch N, Long MY. RNA-based gene duplication: mechanistic and evolutionary insights. Nat Rev Genet. 2009;10(1):19–31. https://doi.org/10.1038/nrg2487.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  22. Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, et al. BLAST+: architecture and applications. BMC Bioinformatics. 2009;10(1):421. https://doi.org/10.1186/1471-2105-10-421.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  23. Vinckenbosch N, Dupanloup I, Kaessmann H. Evolutionary fate of retroposed gene copies in the human genome. Proc Natl Acad Sci U S A. 2006;103(9):3220–5. https://doi.org/10.1073/pnas.0511307103.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  24. Stein JC, Yu Y, Copetti D, Zwickl DJ, Zhang L, Zhang C, et al. Genomes of 13 domesticated and wild rice relatives highlight genetic conservation, turnover and innovation across the genus Oryza. Nat Genet. 2018;50(2):285–96. https://doi.org/10.1038/s41588-018-0040-0.

    CAS  Article  PubMed  Google Scholar 

  25. Mortazavi A, Williams BA, Mccue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5(7):621–8. https://doi.org/10.1038/nmeth.1226.

    CAS  Article  PubMed  Google Scholar 

  26. Zhang L, Ren Y, Yang T, Li G, Chen J, Gschwend AR, et al. Rapid evolution of protein diversity by de novo origination in Oryza. Nat Ecol Evol. 2019;3(4):679–90. https://doi.org/10.1038/s41559-019-0822-5.

    Article  PubMed  Google Scholar 

  27. Prince VE, Pickett FB. Splitting pairs: The diverging fates of duplicated genes. Nat Rev Genet. 2002;3(11):827–37. https://doi.org/10.1038/nrg928.

    CAS  Article  PubMed  Google Scholar 

  28. Betran E, Thornton K, Long M. Retroposed new genes out of the X in Drosophila. Genome Res. 2002;12(12):1854–9. https://doi.org/10.1101/gr.604902.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  29. Vibranovski MD, Lopes HF, Karr TL, Long MY. Stage-specific expression profiling of Drosophila spermatogenesis suggests that meiotic sex chromosome inactivation drives genomic relocation of testis-expressed genes. Plos Genet. 2009;5(11):e1000731. https://doi.org/10.1371/journal.pgen.1000731.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  30. Marques AC, Dupanloup I, Vinckenbosch N, Reymond A, Kaessmann H. Emergence of young human genes after a burst of retroposition in primates. Plos Biol. 2005;3(11):1970–9. https://doi.org/10.1371/journal.pbio.0030357.

    CAS  Article  Google Scholar 

  31. Khodarahmpour Z. Effect of drought stress induced by polyethylene glycol (PEG) on germination indices in corn (Zea mays L.) hybrids. African J Biotechnol. 2011;10(79):18222–7. https://doi.org/10.5897/AJB11.2639.

    CAS  Article  Google Scholar 

  32. Queiroz MS, Oliveira CE, Steiner F, Zuffo AM, Zoz T, Vendruscolo EP, et al. Drought stresses on seed germination and early growth of maize and sorghum. J Agricult Sci. 2019;11(2):310–8. https://doi.org/10.5539/jas.v11n2p310.

    Article  Google Scholar 

  33. Wang W, Yu H, Long M. Duplication-degeneration as a mechanism of gene fission and the origin of new genes in Drosophila species. Nat Genet. 2004;36(5):523–7. https://doi.org/10.1038/ng1338.

    CAS  Article  PubMed  Google Scholar 

  34. Wang W, Lan H. Rapid and parallel chromosomal number reductions in Muntjac deer inferred from mitochondrial DNA phylogeny. Mol Biol Evol. 2000;17(9):1326–33. https://doi.org/10.1093/oxfordjournals.molbev.a026416.

    CAS  Article  PubMed  Google Scholar 

  35. Sboner A, Habegger L, Pflueger D, Terry S, Chen DZ, Rozowsky JS, et al. FusionSeq: a modular framework for finding gene fusions by analyzing paired-end RNA-sequencing data. Genome Biol. 2010;11(10):R104. https://doi.org/10.1186/gb-2010-11-10-r104.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  36. McPherson A, Hormozdiari F, Zayed A, Giuliany R, Ha G, Sun MG, et al. deFuse: an algorithm for gene fusion discovery in tumor RNA-Seq data. PLoS Comput Biol. 2011;7(5):e1001138. https://doi.org/10.1371/journal.pcbi.1001138.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  37. Asmann YW, Hossain A, Necela BM, Middha S, Kalari KR, Sun Z, et al. A novel bioinformatics pipeline for identification and characterization of fusion transcripts in breast cancer and normal cell lines. Nucleic Acids Res. 2011;39(15):e100. https://doi.org/10.1093/nar/gkr362.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  38. Wang W, Zheng HK, Fan CZ, Li J, Shi JJ, Cai ZQ, et al. High rate of chimeric gene origination by retroposition in plant genomes. Plant Cell. 2006;18(8):1791–802. https://doi.org/10.1105/tpc.106.041905.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  39. Sakai H, Mizuno H, Kawahara Y, Wakimoto H, Ikawa H, Kawahigashi H, Kanamori H, Matsumoto T, Itoh T, Gaut BS. Retrogenes in Rice (Oryza sativa L. ssp. japonica) Exhibit Correlated Expression with Their Source Genes. Genome Biology and Evolution. 2011;3:1357–68.

  40. Holt C, Yandell M. MAKER2: an annotation pipeline and genome-database management tool for second-generation genome projects. BMC Bioinformatics. 2011;12(1). https://doi.org/10.1186/1471-2105-12-491.

  41. Campbell MS, Law M, Holt C, Stein JC, Moghe GD, Hufnagel DE, et al. MAKER-P: a tool kit for the rapid creation, management, and quality control of plant genome annotations. Plant Physiol. 2014;164(2):513–24. https://doi.org/10.1104/pp.113.230144.

    CAS  Article  PubMed  Google Scholar 

  42. Kaessmann H. Origins, evolution, and phenotypic impact of new genes. Genome Res. 2010;20(10):1313–26. https://doi.org/10.1101/gr.101386.109.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  43. Veltman MA, Flowers JM, van Andel TR, Schranz ME. Origins and geographic diversification of African rice (Oryza glaberrima). PLOS ONE. 2019;14(3):e0203508. https://doi.org/10.1371/journal.pone.0203508.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  44. Semon M, Nielsen R, Jones MP, SR MC. The population structure of African cultivated rice Oryza glaberrima (Steud.): evidence for elevated levels of linkage disequilibrium caused by admixture with O. sativa and ecological adaptation. Genetics. 2005;169:1639–47. https://doi.org/10.1534/genetics.104.033175.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  45. Wang M, Yu Y, Haberer G, Marri PR, Fan C, Goicoechea JL, et al. The genome sequence of African rice (Oryza glaberrima) and evidence for independent domestication. Nat Genet. 2014;46(9):982–8. https://doi.org/10.1038/ng.3044.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  46. Nabholz B, Sarah G, Sabot F, Ruiz M, Adam H, Nidelet S, et al. Transcriptome population genomics reveals severe bottleneck and domestication cost in the African rice (Oryza glaberrima). Mol Ecol. 2014;23(9):2210–27. https://doi.org/10.1111/mec.12738.

    CAS  Article  PubMed  Google Scholar 

  47. Bai Y, Casola C, Feschotte C, Betran E. Comparative genomics reveals a constant rate of origination and convergent acquisition of functional retrogenes in Drosophila. Genome Biol. 2007;8(1):R11. https://doi.org/10.1186/gb-2007-8-1-r11.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  48. Fu B, Chen M, Zou M, Long M, He S. The rapid generation of chimerical genes expanding protein diversity in zebrafish. BMC Genomics. 2010;11(1):657. https://doi.org/10.1186/1471-2164-11-657.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  49. Zhang CJ, Gschwend AR, Ouyang YD, Long MY. Evolution of gene structural complexity: an alternative-splicing-based model accounts for intron-containing retrogenes. Plant Physiol. 2014;165(1):412–23. https://doi.org/10.1104/pp.113.231696.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  50. Semon M, Wolfe KH. Consequences of genome duplication. Current Opinion in Genetics & Development. 2007;17(6):505–12. https://doi.org/10.1016/j.gde.2007.09.007.

    CAS  Article  Google Scholar 

  51. Innan H, Kondrashov F. The evolution of gene duplications: classifying and distinguishing between models. Nature Reviews Genetics. 2010;11(2):97–108. https://doi.org/10.1038/nrg2689.

    CAS  Article  PubMed  Google Scholar 

  52. Arguello JR, Chen Y, Yang S, Wang W, Long M. Origination of an X-linked testes chimeric gene by illegitimate recombination in Drosophila. PLoS Genet. 2006;2(5):e77. https://doi.org/10.1371/journal.pgen.0020077.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  53. Bennetzen JL. Transposable elements, gene creation and genome rearrangement in flowering plants. Curr Opin Genet Dev. 2005;15(6):621–7. https://doi.org/10.1016/j.gde.2005.09.010.

    CAS  Article  PubMed  Google Scholar 

  54. Jiang N, Bao ZR, Zhang XY, Eddy SR, Wessler SR. Pack-MULE transposable elements mediate gene evolution in plants. Nature. 2004;431(7008):569–73. https://doi.org/10.1038/nature02953.

    CAS  Article  PubMed  Google Scholar 

  55. Morgante M, Brunner S, Pea G, Fengler K, Zuccolo A, Rafalski A. Gene duplication and exon shuffling by helitron-like transposons generate intraspecies diversity in maize. Nat Genet. 2005;37(9):997–1002. https://doi.org/10.1038/ng1615.

    CAS  Article  PubMed  Google Scholar 

  56. Buljan M, Frankish A, Bateman A. Quantifying the mechanisms of domain gain in animal proteins. Genome Biol. 2010;11(7):R74. https://doi.org/10.1186/gb-2010-11-7-r74.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  57. Williford A, Betrán E: Gene fusion. In eLS. John Wiley & Sons, Ltd; 2013, Gene Fusion

  58. Yang SA, Arguello JR, Li X, Ding Y, Zhou Q, Chen Y, et al. Repetitive element-mediated recombination as a mechanism for new gene origination in Drosophila. Plos Genet. 2008;4(1):e3. https://doi.org/10.1371/journal.pgen.0040003.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  59. Huang Y, Chen J, Dong C, Sosa D, Xia S, Ouyang Y, et al. Species-specific partial gene duplication in Arabidopsis thaliana evolved novel phenotypic effects on morphological traits under strong positive selection. Plant Cell. 2021;34(2):802–17. https://doi.org/10.1093/plcell/koab291.

    Article  PubMed Central  Google Scholar 

  60. Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25(9):1105–11. https://doi.org/10.1093/bioinformatics/btp120.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  61. Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28(5):511–U174. https://doi.org/10.1038/nbt.1621.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  62. Vogel JP, Garvin DF, Mockler TC, Schmutz J, Rokhsar D, Bevan MW, et al. Genome sequencing and analysis of the model grass Brachypodium distachyon. Nature. 2010;463(7282):763–8. https://doi.org/10.1038/nature08747.

    CAS  Article  Google Scholar 

  63. Ouyang S, Zhu W, Hamilton J, Lin H, Campbell M, Childs K, Thibaud-Nissen F, Malek RL, Lee Y, Zheng L, Orvis J, Haas B, Wortman J, Buell CR: The TIGR rice genome annotation resource: improvements and new features. Nucleic Acids Res 2007, 35:D883-D887, Database, DOI: https://doi.org/10.1093/nar/gkl976.

  64. Kawahara Y, de la Bastide M, Hamilton JP, Kanamori H, McCombie WR, Ouyang S, et al. Improvement of the Oryza sativa Nipponbare reference genome using next generation sequence and optical map data. Rice. 2013;6(1). https://doi.org/10.1186/1939-8433-6-4.

  65. Sakai H, Lee SS, Tanaka T, Numa H, Kim J, Kawahara Y, Wakimoto H, Yang C, Iwamoto M, Abe T, et al: Rice Annotation Project Database (RAP-DB): an integrative and interactive database for rice genomics. Plant Cell Physiol 2013, 54:E6-+.

  66. Gao ZY, Zhao SC, He WM, Guo LB, Peng YL, Wang JJ, et al. Dissecting yield-associated loci in super hybrid rice by resequencing recombinant inbred lines and improving parental genome sequences. Proc Natl Acad Sci U S A. 2013;110(35):14492–7. https://doi.org/10.1073/pnas.1306579110.

    Article  PubMed  PubMed Central  Google Scholar 

  67. Lamesch P, Berardini TZ, Li DH, Swarbreck D, Wilks C, Sasidharan R, et al. The Arabidopsis Information Resource (TAIR): improved gene annotation and new tools. Nucleic Acids Res. 2012;40(D1):D1202–10. https://doi.org/10.1093/nar/gkr1090.

    CAS  Article  PubMed  Google Scholar 

  68. Paterson AH, Bowers JE, Bruggmann R, Dubchak I, Grimwood J, Gundlach H, et al. The Sorghum bicolor genome and the diversification of grasses. Nature. 2009;457(7229):551–6. https://doi.org/10.1038/nature07723.

    CAS  Article  PubMed  Google Scholar 

  69. Schnable PS. The B73 maize genome: complexity, diversity, and dynamics (November, pg 1112, 2009). Science. 2012;337:1040.

    CAS  Google Scholar 

  70. McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20(9):1297–303. https://doi.org/10.1101/gr.107524.110.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  71. Wang L, Xie WB, Chen Y, Tang WJ, Yang JY, Ye RJ, et al. A dynamic gene expression atlas covering the entire life cycle of rice. Plant J. 2010;61(5):752–66. https://doi.org/10.1111/j.1365-313X.2009.04100.x.

    CAS  Article  PubMed  Google Scholar 

  72. Jain M. Genome-wide identification of novel internal control genes for normalization of gene expression during various stages of development in rice. Plant Sci. 2009;176(5):702–6. https://doi.org/10.1016/j.plantsci.2009.02.001.

    CAS  Article  Google Scholar 

  73. Zhang C. Gene fusion as an important mechanism to generate new genes in the genus Oryza (GriffinDetector). Github. 2022; https://github.com/zhangcj2022/GriffinDetector.git.

  74. Zhang C. Gene fusion as an important mechanism to generate new genes in the genus Oryza (GriffinDetector). Zenodo. 2022. https://doi.org/10.5281/zenodo.6480086.

Download references

Acknowledgements

We are grateful to Yidan Oyang for helpful discussions.

Review history

The review history is available as Additional file 11.

Peer review information

Wenjing She was the primary editor of this article and managed its editorial process and peer review in collaboration with the rest of the editorial team.

Funding

This research was financially funded by National Natural Science Foundation of China (31571311, 32000180), Youth Innovation Promotion Association CAS (2021394), Yunnan Fundamental Research Projects (202001AU070075), Yunnan Young & Elite Talents Project (YNWR-QNBJ-2019-268), and CAS Scholarship. ML was supported by NSF1026200 and NSF1051826; RW by NSF1026200.

Author information

Authors and Affiliations

Authors

Contributions

CZ, ML, and RW designed the idea and participated in revision and discussion; CZ, YZ, LZ, QY, and NL executed data analysis; JZ and MW helped in data collection and method; YZ, WF, and GL conducted experiment; YL and CZ drafted the manuscript with significant participation of ML. All author(s) read and approved the final manuscript.

Corresponding authors

Correspondence to Chengjun Zhang, Manyuan Long or Rod A. Wing.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Chimeric gene output of four focus species.

Additional file 2.

Candidate fusion gene presented in four focus species.

Additional file 3.

Fixation frequency of species-specific fusion gene included in four focus species.

Additional file 4.

FPKM of fusion gene and short copy homologs in O. sativa v.g. japonica, O. barthii, O. glaberrima, O. punctate, O. branchyantha and Leersia.

Additional file 5: Supplemental figures.

Figure S1. Expression profiles of fusion genes and their parent genes in O. sativa v.g. japonica and other sibling species like O. barthii, O. glaberrima, O. punctata, O. brchyantha and O. Leersia; Figure S2. Statistics of expression levels between long homolog and short homolog; Figure S3. Expression heat map of two selected fusion genes in response to abiotic stress; Figure S4. Spatial expression pattern of 69 expressed fusion genes; Figure S5. Generation of transgenic rice plants. (A) Embryogenic callus derived from the scutella of mature seeds. (B) Subculture elite callus for proliferation. (C) Inoculation of callus with Agrobacterium carrying CRISPR/Cas9 vector by co-culture. (D) Screen hygromycin-resistant callus by transferring to selective medium. (E) Regeneration of hygromycin-resistant callus on MS medium for 1 month. (F) Rooting and growth of transgenic rice plants.

Additional file 6.

Correlation of O. japonica fusion gene and its short copies FPKM value in different tissues.

Additional file 7.

Correlation between the maximum FPKM value of long-homologous copies and the short-homologous copies in four Oryza species.

Additional file 8.

The generation patterns of fusion genes.

Additional file 9.

Transposable element retrieved from O. sativa v.g. japonica fusion genes.

Additional file 10.

Gene structure re-annotation of species-specific fusion gene in four focus species.

Additional file 11.

Review history.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Zhou, Y., Zhang, C., Zhang, L. et al. Gene fusion as an important mechanism to generate new genes in the genus Oryza. Genome Biol 23, 130 (2022). https://doi.org/10.1186/s13059-022-02696-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13059-022-02696-w

Keywords

  • Fusion gene
  • Oryza
  • Evolutionary patterns
  • Phenotype