Organism: Brassica rapa ssp. pekinensis
We have obtained the sequences of BrGGP2 (Bra035554; ortholog to AT5G55120) and BrGPP (Bra040607; ortholog to AT3G02870) from the Brassica database (http://brassicadb.org/brad/). Upon multiple sequence alignment and further database search, we found that the sequences of Bra035554 and Bra040607 are experienced adjacent gene fusions. Bra035554 gained one domain (PF04885.8 – stigma specific domain; function unknown) coded by two exons at the N-terminus and Bra040607 gained one domain (PF12874.2 – zinc finger (C2H2 type) domain) coded by seven exons at the C-terminus. In order to validate these gene fusions by experiments, we performed PCR with different primer pairs to clone the full-length of Bra035554 and Bra040607 using cDNAs of different tissues. If the gene fusion is true, then we must get the amplification of full-length of Bra035554 and Bra040607 by using B. rapa cDNA samples. But, we did not get any amplification with any cDNAs. Intriguingly, when I ran PCR with N-terminus and C-terminus specific primers, we got amplification (see the attached figure). Meanwhile, I also searched NCBI-EST database and found no sequence to support gene fusion. Now, I can conclude either (i) the sequences of Bra035554 and Bra040607 provided in Brassica database may be wrong? (Brassica is a primary database for Brassica species, so the sequences must be genuine), or (ii) further research is required for the transcriptional regulation of Bra035554 and Bra040607. What else I should do to confirm whether the gene fusions in Bra035554 and Bra040607 are true?