Kanako O. Koyanagi

Large-scale identification and characterization of alternative splicing variants of human gene transcripts using 56,419 completely sequenced and manually annotated full-length cDNAs.

We report the first genome-wide identification and characterization of alternative splicing in human gene transcripts based on analysis of the full-length cDNAs. Applying both manual and computational analyses for 56 419 completely sequenced and precisely annotated full-length cDNAs selected for the H-Invitational human transcriptome annotation meetings, we identified 6877 alternative splicing genes with 18 297 different alternative splicing variants. A total of 37 670 exons were involved in these alternative splicing events. The encoded protein sequences were affected in 6005 of the 6877 genes. Notably, alternative splicing affected protein motifs in 3015 genes, subcellular localizations in 2982 genes and transmembrane domains in 1348 genes. We also identified interesting patterns of alternative splicing, in which two distinct genes seemed to be bridged, nested or having overlapping protein coding sequences (CDSs) of different reading frames (multiple CDS). In these cases, completely unrelated proteins are encoded by a single locus. Genome-wide annotations of alternative splicing, relying on full-length cDNAs, should lay firm groundwork for exploring in detail the diversification of protein function, which is mediated by the fast expanding universe of alternative splicing variants.

Complete Genome Analysis of a Novel Intertypic Recombinant Human Adenovirus Causing Epidemic Keratoconjunctivitis in Japan

ABSTRACT For 4 months from September 2008, 102 conjunctival swab specimens were collected for surveillance purposes from patients across Japan suspected of having epidemic keratoconjunctivitis (EKC). Human adenovirus (HAdV) DNA was detected in 61 samples by PCR, though the HAdV type for 6 of the PCR-positive samples could not be determined by phylogenetic analysis using a partial hexon gene sequence. Moreover, for 2 months from January 2009, HAdV strains with identical sequences were isolated from five conjunctival swab samples obtained from EKC patients in five different regions of Japan. For the analyses of the 11 samples mentioned above, we determined the nucleotide sequences of the entire penton base, hexon, and fiber genes and early 3 (E3) region, which are variable regions among HAdV types, and compared them to those of other HAdV species D strains. The nucleotide sequences of loops 1 and 2 in the hexons of all 11 samples showed high degrees of identity with those of the HAdV type 15 (HAdV-15) and HAdV-29 prototype strains. However, the fiber gene and E3 region sequences showed high degrees of identity with those of HAdV-9, and the penton base gene sequence showed a high degree of identity with the penton base gene sequences of HAdV-9 and -26. Moreover, the complete genome sequence of the 2307-S strain, which was isolated by viral culture from 1 of the 11 samples, was determined. The 2307-S strain was a recombinant HAdV between HAdV-9, -15, -26, -29, and/or another HAdV type; however, the recombination sites in the genome were not obvious. We propose that this virus is a novel intertypic recombinant, HAdV-15/29/H9, and may be an etiological agent of EKC.

Evola : Ortholog database of all human genes in H-InvDB with manual curation of phylogenetic trees

Orthologs are genes in different species that evolved from a common ancestral gene by speciation. Currently, with the rapid growth of transcriptome data of various species, more reliable orthology information is prerequisite for further studies. However, detection of orthologs could be erroneous if pairwise distance-based methods, such as reciprocal BLAST searches, are utilized. Thus, as a sub-database of H-InvDB, an integrated database of annotated human genes (http://h-invitational.jp/), we constructed a fully curated database of evolutionary features of human genes, called ‘Evola’. In the process of the ortholog detection, computational analysis based on conserved genome synteny and transcript sequence similarity was followed by manual curation by researchers examining phylogenetic trees. In total, 18 968 human genes have orthologs among 11 vertebrates (chimpanzee, mouse, cow, chicken, zebrafish, etc.), either computationally detected or manually curated orthologs. Evola provides amino acid sequence alignments and phylogenetic trees of orthologs and homologs. In ‘dN/dS view’, natural selection on genes can be analyzed between human and other species. In ‘Locus maps’, all transcript variants and their exon/intron structures can be compared among orthologous gene loci. We expect the Evola to serve as a comprehensive and reliable database to be utilized in comparative analyses for obtaining new knowledge about human genes. Evola is available at http://www.h-invitational.jp/evola/.

Recombination analysis of intermediate human adenovirus type 53 in Japan by complete genome sequence.

Human adenovirus type 53 (HAdV-53) has commonly been detected in samples from epidemic keratoconjunctivitis (EKC) patients in Japan since 1996. HAdV-53 is an intermediate virus, containing hexon-chimeric, penton base and fiber structures similar to HAdV-22 and -37, HAdV-37 and HAdV-8, respectively. HAdV-53-like intermediate strains were first isolated from EKC samples in Japan in the 1980s. Here, the complete genome sequences of three such HAdV-53-like intermediate strains (870006C, 880249C and 890357C) and four HAdV-53 strains were determined, and their relationships were analysed. The seven HAdV strains were classified into three groups, 870006C/880249C, 890357C and the four HAdV-53 strains, on the basis of phylogenetic analyses of the partial and complete genome sequences. HAdV strains within the same group showed the highest nucleotide identities (99.87-100.00 %). Like HAdV-53, the hexon loop 1 and 2 regions of 870006C, 880249C and 890357C showed the highest identity with HAdV-22. However, these strains did not show a hexon-chimeric structure similar to HAdV-22 and -37, or a penton base similar to HAdV-37. The fiber genes of 870006C and 880249C were identical to that of HAdV-37, but not HAdV-8. Thus, the three intermediate HAdVs isolated in the 1980s were similar to each other but not to HAdV-53. The recombination breakpoints were inferred by the Recombination Detection Program (rdp) using whole-genome sequences of these seven HAdV and of 12 HAdV-D strains from GenBank. HAdV-53 may have evolved from intermediate HAdVs circulating in the 1980s, and from HAdV-8, -22 and -37, by recombination of sections cut at the putative breakpoints.

Highly Diversified Molecular Evolution of Downstream Transcription Start Sites in Rice and Arabidopsis

Alternative usage of transcription start sites (TSSs) is one of the key mechanisms to generate gene variation in eukaryotes. Here, we show diversified molecular evolution of TSSs in remotely related flowering plants, rice (Oryza sativa) and Arabidopsis (Arabidopsis thaliana), by comprehensive analyses of large collections of full-length cDNAs and genome sequences. We determined 45,917 representative TSSs within 23,445 loci of rice and 35,313 TSSs within 16,964 loci of Arabidopsis, about two TSSs per locus in either species. The nucleotide features around TSSs displayed distinct patterns when the most upstream TSSs were compared with downstream TSSs. We found that CG-skew and AT-skew were clearly different between upstream and downstream TSSs, and that this difference was commonly observed in rice and Arabidopsis. Relative entropy analysis revealed that the most upstream TSSs had retained canonical cis elements, whereas downstream TSSs showed atypical nucleotide features. Expression patterns were distinguishable between upstream and downstream TSSs. These results indicate that plant TSSs were generally diversified in downstream regions, resulting in the development of new gene expression patterns. Furthermore, our comparative analysis of TSS variation between the species showed a positive correlation between TSS number and gene conservation. Rice and Arabidopsis might have evolved novel TSSs in an independent manner, which led to diversification of these two species.

Evolutionary force of AT-rich repeats to trap genomic and episomal DNAs into the rice genome: lessons from endogenous pararetrovirus

In plant genomes, the incorporation of DNA segments is not a common method of artificial gene transfer. Nevertheless, various segments of pararetroviruses have been found in plant genomes in recent decades. The rice genome contains a number of segments of endogenous rice tungro bacilliform virus-like sequences (ERTBVs), many of which are present between AT dinucleotide repeats (ATrs). Comparison of genomic sequences between two closely related rice subspecies, japonica and indica, allowed us to verify the preferential insertion of ERTBVs into ATrs. In addition to ERTBVs, the comparative analyses showed that ATrs occasionally incorporate repeat sequences including transposable elements, and a wide range of other sequences. Besides the known genomic sequences, the insertion sequences also represented DNAs of unclear origins together with ERTBVs, suggesting that ATrs have integrated episomal DNAs that would have been suspended in the nucleus. Such insertion DNAs might be trapped by ATrs in the genome in a host-dependent manner. Conversely, other simple mono- and dinucleotide sequence repeats (SSR) were less frequently involved in insertion events relative to ATrs. Therefore, ATrs could be regarded as hot spots of double-strand breaks that induce non-homologous end joining. The insertions within ATrs occasionally generated new gene-related sequences or involved structural modifications of existing genes. Likewise, in a comparison between Arabidopsis thaliana and Arabidopsis lyrata, the insertions preferred ATrs to other SSRs. Therefore ATrs in plant genomes could be considered as genomic dumping sites that have trapped various DNA molecules and may have exerted a powerful evolutionary force.

Comparative genomic analysis of retrogene repertoire in two green algae Volvox carteri and Chlamydomonas reinhardtii

BackgroundRetroposition, one of the processes of copying the genetic material, is an important RNA-mediated mechanism leading to the emergence of new genes. Because the transcription controlling segments are usually not copied to the new location in this mechanism, the duplicated gene copies (retrocopies) become pseudogenized. However, few can still survive, e.g. by recruiting novel regulatory elements from the region of insertion. Subsequently, these duplicated genes can contribute to the formation of lineage-specific traits and phenotypic diversity. Despite the numerous studies of the functional retrocopies (retrogenes) in animals and plants, very little is known about their presence in green algae, including morphologically diverse species. The current availability of the genomes of both uni- and multicellular algae provides a good opportunity to conduct a genome-wide investigation in order to fill the knowledge gap in retroposition phenomenon in this lineage.ResultsHere we present a comparative genomic analysis of uni- and multicellular algae, Chlamydomonas reinhardtii and Volvox carteri, respectively, to explore their retrogene complements. By adopting a computational approach, we identified 141 retrogene candidates in total in both genomes, with their fraction being significantly higher in the multicellular Volvox. Majority of the retrogene candidates showed signatures of functional constraints, thus indicating their functionality. Detailed analyses of the identified retrogene candidates, their parental genes, and homologs of both, revealed that most of the retrogene candidates were derived from ancient retroposition events in the common ancestor of the two algae and that the parental genes were subsequently lost from the respective lineages, making many retrogenes ‘orphan’.ConclusionWe revealed that the genomes of the green algae have maintained many possibly functional retrogenes in spite of experiencing various molecular evolutionary events during a long evolutionary time after the retroposition events. Our first report about the retrogene set in the green algae provides a good foundation for any future investigation of the repertoire of retrogenes and facilitates the assessment of the evolutionary impact of retroposition on diverse morphological traits in this lineage.ReviewersThis article was reviewed by William Martin and Piotr Zielenkiewicz.

Genome-wide association mapping focusing on a rice population derived from rice breeding programs in a region.

Plant breeding programs in local regions may generate genetic variations that are desirable to local populations and shape adaptability during the establishment of local populations. To elucidate genetic bases for this process, we proposed a new approach for identifying the genetic bases for the traits improved during rice breeding programs; association mapping focusing on a local population. In the present study, we performed association mapping focusing on a local rice population, consisting of 63 varieties, in Hokkaido, the northernmost region of Japan and one of the northern limits of rice cultivation worldwide. Six and seventeen QTLs were identified for heading date and low temperature germinability, respectively. Of these, 13 were novel QTLs in this population and 10 corresponded to the QTLs previously reported based on QTL mapping. The identification of QTLs for traits in local populations including elite varieties may lead to a better understanding of the genetic bases of elite traits. This is of direct relevance for plant breeding programs in local regions.

Inferring Cell Differentiation Processes Based on Phylogenetic Analysis of Genome-Wide Epigenetic Information: Hematopoiesis as a Model Case

How cells divide and differentiate is a fundamental question in organismal development; however, the discovery of differentiation processes in various cell types is laborious and sometimes impossible. Phylogenetic analysis is typically used to reconstruct evolutionary processes based on inherent characters. It could also be used to reconstruct developmental processes based on the developmental changes that occur during cell proliferation and differentiation. In this study, DNA methylation information from differentiated hematopoietic cells was used to perform phylogenetic analyses. The results were assessed for their validity in inferring hierarchical differentiation processes of hematopoietic cells and DNA methylation processes of differentiating progenitor cells. Overall, phylogenetic analyses based on DNA methylation information facilitated inferences regarding hematopoiesis.

Rice genomes recorded ancient pararetrovirus activities: Virus genealogy and multiple origins of endogenization during rice speciation

Viral fossils in rice genomes are a best entity to understand ancient pararetrovirus activities through host plant history because of our advanced knowledge of the genomes and evolutionary history with rice and its related species. Here, we explored organization, geographic origins and genealogy of rice pararetroviruses, which were turned into endogenous rice tungro bacilliform virus-like (eRTBVL) sequences. About 300 eRTBVL sequences from three representative rice genomes were clearly classified into six families. Most of the endogenization events of the eRTBVLs were initiated before differentiation of the rice progenitor (> 160,000 years ago). We successfully followed the genealogy of old relic viruses during rice speciation, and inferred the geographical origins for these viruses. Possible virus genomic sequences were explained mostly by recombinations between different virus families. Interestingly, we discovered that only a few recombination events among the numerous occasions had determined the virus genealogy.