Phylogenetic analysis of endogenous viral elements in the rice genome reveals local chromosomal evolution in Oryza AA-genome species.

Introduction: Rice genomes contain endogenous viral elements homologous to rice tungro bacilliform virus (RTBV) from the pararetrovirus family Caulimoviridae. These viral elements, known as endogenous RTBV-like sequences (eRTBVLs), comprise five subfamilies, eRTBVL-A, -B, -C, -D, and -X. Four subfamilies (A, B, C, and X) are present to a limited degree in the genomes of the Asian cultivated rice Oryza sativa (spp. japonica and indica) and the closely related wild species Oryza rufipogon. Methods: The eRTBVL-D sequences are widely distributed within these and other Oryza AA-genome species. Fifteen eRTBVL-D segments identified in the japonica (Nipponbare) genome occur mostly at orthologous chromosomal positions in other AA-genome species. The eRTBVL-D sequences were inserted into the genomes just before speciation of the AA-genome species. Results and discussion: Ten eRTBVL-D segments are located at six loci, which were used for our evolutionary analyses during the speciation of the AA-genome species. The degree of genetic differentiation varied among the eRTBVL-D segments. Of the six loci, three showed phylogenetic trees consistent with the standard speciation pattern (SSP) of the AA-genome species (Type A), and the other three represented phylogenies different from the SSP (Type B). The atypical phylogenetic trees for the Type B loci revealed chromosome region-specific evolution among the AA-genome species that is associated with phylogenetic incongruences: complex genome rearrangements between eRTBVL-D segments, an introgression between the distant species, and low genetic diversity of a shared eRTBVL-D segment. Using eRTBVL-D as an indicator, this study revealed the phylogenetic incongruence of local chromosomal regions with different topologies that developed during speciation.

Inferring chromatin accessibility during murine hematopoiesis through phylogenetic analysis.

OBJECTIVE: Diversification of cell types and changes in epigenetic states during cell differentiation processes are important for understanding development. Recently, phylogenetic analysis using DNA methylation and histone modification information has been shown useful for inferring these processes. The purpose of this study was to examine whether chromatin accessibility data can help infer these processes in murine hematopoiesis. RESULTS: Chromatin accessibility data could partially infer the hematopoietic differentiation hierarchy. Furthermore, based on the ancestral state estimation of internal nodes, the open/closed chromatin states of differentiating progenitor cells could be predicted with a specificity of 0.86-0.99 and sensitivity of 0.29-0.72. These results suggest that the phylogenetic analysis of chromatin accessibility could offer important information on cell differentiation, particularly for organisms from which progenitor cells are difficult to obtain.

Translation of continuous artificial selection on phenotype into genotype during rice breeding programs.

Understanding genetic diversity among local populations is a primary goal of modern crop breeding programs. Here, we demonstrated the genetic relationships of rice varieties in Hokkaido, Japan, one of the northern limits of rice cultivation around the world. Furthermore, artificial selection during rice breeding programs has been characterized using genome sequences. We utilized 8,565 single nucleotide polymorphisms and insertion/deletion markers distributed across the genome in genotype-by-sequencing for genetic diversity analyses. Phylogenetics, genetic population structure, and principal component analysis showed that a total of 110 varieties were classified into four distinct clusters according to different populations geographically and historically. Furthermore, the genome sequences of 19 rice varieties along with historic representations in Hokkaido, nucleotide diversity and FST values in each cluster revealed that artificial selection of elite phenotypes focused on chromosomal regions. These results clearly demonstrated the history of the selections on agronomic traits as genome sequences among current rice varieties from Hokkaido.

Biosynthetic Short Neuropeptides: A Rational Theory Based on Experimental Results for the Missing Pain-Relief Opioid Endomorphin Precursor Gene.

Endomorphins are neuropeptides that bind strongly to µ-opioid receptors and are considered to play important roles in pain modulation and other biological functions. Two endomorphins have been identified, to date, endomorphine-1 and -2; both are tetrapeptides and differ by only a single amino acid in the third position. Both peptides were isolated from bovine brains; however, their precursor genes have not been identified. In this study, a nucleotide sequence corresponding to the endomorphin-1 peptide in an expressed sequence tag database has been found and a preproendomorphin-like precursor peptide from human brain complementary DNA (cDNA) has been cloned. The cDNA consists of nucleotide sequences of two already annotated predicted genes, and the putative peptide differs by one amino acid from the isolated endomorphin peptides. It is proposed herein that there is the possibility of unknown short proteins or peptide precursors being missed by automated gene prediction programs based on similarities of known protein sequences. A novel concept of how to produce endomorphins from a similar peptide is described. The oxidatively modified base might provide a clue for understanding discrepancies between nucleotide sequences on the genome and those on cDNAs.

Inferring changes in histone modification during cell differentiation by ancestral state estimation based on phylogenetic trees of cell types: Human hematopoiesis as a model case.

Revealing the landscape of epigenetic changes in cells during differentiation is important for understanding the development of organisms. In this study, to infer such epigenetic changes during human hematopoiesis, ancestral state estimation based on a phylogenetic tree was applied to map the epigenomic changes in six kinds of histone modifications onto the hierarchical cell differentiation process of hematopoiesis using epigenomes of eight types of differentiated hematopoietic cells. The histone modification changes inferred during hematopoiesis showed that changes that occurred on the branches separating different cell types reflected the characteristics of hematopoiesis in terms of genomic position and gene function. These results suggested that ancestral state estimation based on phylogenetic analysis of histone modifications in differentiated hematopoietic cells could reconstruct an appropriate landscape of histone modification changes during hematopoiesis. Since integration of the inferred changes of different histone modifications could reveal genes with specific histone marks such as active histone marks and bivalent histone marks on each internal branch of cell-type trees, this approach could provide valuable information for understanding the cell differentiation steps of each cell lineage.

Inferring changes in histone modification during cell differentiation by ancestral state estimation based on phylogenetic trees of cell types: Human hematopoiesis as a model case.

Revealing the landscape of epigenetic changes in cells during differentiation is important for understanding the development of organisms. In this study, to infer such epigenetic changes during human hematopoiesis, ancestral state estimation based on a phylogenetic tree was applied to map the epigenomic changes in six kinds of histone modifications onto the hierarchical cell differentiation process of hematopoiesis using epigenomes of eight types of differentiated hematopoietic cells. The histone modification changes inferred during hematopoiesis showed that changes that occurred on the branches separating different cell types reflected the characteristics of hematopoiesis in terms of genomic position and gene function. These results suggested that ancestral state estimation based on phylogenetic analysis of histone modifications in differentiated hematopoietic cells could reconstruct an appropriate landscape of histone modification changes during hematopoiesis. Since integration of the inferred changes of different histone modifications could reveal genes with specific histone marks such as active histone marks and bivalent histone marks on each internal branch of cell-type trees, this approach could provide valuable information for understanding the cell differentiation steps of each cell lineage.

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

BACKGROUND: Retroposition, 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. RESULTS: Here 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'. CONCLUSION: We 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. REVIEWERS: This article was reviewed by William Martin and Piotr Zielenkiewicz.

Interregional Coevolution Analysis Revealing Functional and Structural Interrelatedness between Different Genomic Regions in Human Mastadenovirus D.

UNLABELLED: Human mastadenovirus D (HAdV-D) is exceptionally rich in type among the seven human adenovirus species. This feature is attributed to frequent intertypic recombination events that have reshuffled orthologous genomic regions between different HAdV-D types. However, this trend appears to be paradoxical, as it has been demonstrated that the replacement of some of the interacting proteins for a specific function with other orthologues causes malfunction, indicating that intertypic recombination events may be deleterious. In order to understand why the paradoxical trend has been possible in HAdV-D evolution, we conducted an interregional coevolution analysis between different genomic regions of 45 different HAdV-D types and found that ca. 70% of the genome has coevolved, even though these are fragmented into several pieces via short intertypic recombination hot spot regions. Since it is statistically and biologically unlikely that all of the coevolving fragments have synchronously recombined between different genomes, it is probable that these regions have stayed in their original genomes during evolution as a platform for frequent intertypic recombination events in limited regions. It is also unlikely that the same genomic regions have remained almost untouched during frequent recombination events, independently, in all different types, by chance. In addition, the coevolving regions contain the coding regions of physically interacting proteins for important functions. Therefore, the coevolution of these regions should be attributed at least in part to natural selection due to common biological constraints operating on all types, including protein-protein interactions for essential functions. Our results predict additional unknown protein interactions. IMPORTANCE: Human mastadenovirus D, an exceptionally type-rich human adenovirus species and causative agent of different diseases in a wide variety of tissues, including that of ocular region and digestive tract, as well as an opportunistic infection in immunocompromised patients, is known to have highly diverged through frequent intertypic recombination events; however, it has also been demonstrated that the replacement of a component protein of a multiprotein system with a homologous protein causes malfunction. The present study solved this apparent paradox by looking at which genomic parts have coevolved using a newly developed method. The results revealed that intertypic recombination events have occurred in limited genomic regions and been avoided in the genomic regions encoding proteins that physically interact for a given function. This approach detects purifying selection against recombination events causing the replacement of partial components of multiprotein systems and therefore predicts physical and functional interactions between different proteins and/or genomic elements.

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.

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.

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.

Intertypic modular exchanges of genomic segments by homologous recombination at universally conserved segments in human adenovirus species D.

Human adenovirus species D (HAdV-D), which is composed of clinically and epidemiologically important pathogens worldwide, contains more taxonomic "types" than any other species of the genus Mastadenovirus, although the mechanisms accounting for the high level of diversity remain to be disclosed. Recent studies of known and new types of HAdV-D have indicated that intertypic recombination between distant types contributes to the increasing diversity of the species. However, such findings raise the question as to how homologous recombination events occur between diversified types since homologous recombination is suppressed as nucleotide sequences diverge. In order to address this question, we investigated the distribution of the recombination boundaries in comparison with the landscape of intergenomic sequence conservation assessed according to the synonymous substitution rate (dS). The results revealed that specific genomic segments are conserved between even the most distantly related genomes; we call these segments "universally conserved segments" (UCSs). These findings suggest that UCSs facilitate homologous recombination, resulting in intergenomic segmental exchanges of UCS-flanking genomic regions as recombination modules. With the aid of such a mechanism, the haploid genomes of HAdV-Ds may have been reshuffled, resulting in chimeric genomes out of diversified repertoires in the HAdV-D population analogous to the MHC region reshuffled via crossing over in vertebrates. In addition, some HAdVs with chimeric genomes may have had the opportunity to avoid host immune responses thereby causing epidemics.

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.

In silico and in vivo identification of the intermediate filament vimentin that is downregulated downstream of Brachyury during Xenopus embryogenesis.

Brachyury, a member of the T-box transcription family, has been suggested to be essential for morphogenetic movements in various processes of animal development. However, little is known about its critical transcriptional targets. In order to identify targets of Brachyury and understand the molecular mechanisms underlying morphogenetic movements, we first searched the genome sequence of Xenopus tropicalis, the only amphibian genomic sequence available, for Brachyury-binding sequences known as T-half sites, and then screened for the ones conserved between vertebrate genomes. We found three genes that have evolutionarily conserved T-half sites in the promoter regions and examined these genes experimentally to determine whether their expressions were regulated by Brachyury, using the animal cap system of Xenopus laevis embryos. Eventually, we obtained evidence that vimentin, encoding an intermediate filament protein, was a potential target of Brachyury. This is the first report to demonstrate that Brachyury might affect the cytoskeletal structure through regulating the expression of an intermediate filament protein, vimentin.

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.

Complete genome analysis of a novel intertypic recombinant human adenovirus causing epidemic keratoconjunctivitis in Japan.

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.

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.

The Rice Annotation Project Database (RAP-DB): 2008 update.

The Rice Annotation Project Database (RAP-DB) was created to provide the genome sequence assembly of the International Rice Genome Sequencing Project (IRGSP), manually curated annotation of the sequence, and other genomics information that could be useful for comprehensive understanding of the rice biology. Since the last publication of the RAP-DB, the IRGSP genome has been revised and reassembled. In addition, a large number of rice-expressed sequence tags have been released, and functional genomics resources have been produced worldwide. Thus, we have thoroughly updated our genome annotation by manual curation of all the functional descriptions of rice genes. The latest version of the RAP-DB contains a variety of annotation data as follows: clone positions, structures and functions of 31 439 genes validated by cDNAs, RNA genes detected by massively parallel signature sequencing (MPSS) technology and sequence similarity, flanking sequences of mutant lines, transposable elements, etc. Other annotation data such as Gnomon can be displayed along with those of RAP for comparison. We have also developed a new keyword search system to allow the user to access useful information. The RAP-DB is available at: http://rapdb.dna.affrc.go.jp/ and http://rapdb.lab.nig.ac.jp/.