Kousuke Hanada

A comparison of the molecular clock of hepatitis C virus in the United States and Japan predicts that hepatocellular carcinoma incidence in the United States will increase over the next two decades

The prevalence of hepatitis C virus (HCV)-related hepatocellular carcinoma (HCC) is considerably lower in the U.S. than in Japan. To elucidate this difference, we determined the time origin of the HCV epidemic in each country by using molecularly clocked long-term serial samples obtained from HCV carriers of genotypes 1a and 1b. The molecular clock estimated that HCV genotype 1 first appeared in Japan in around 1882, whereas emergence in the U.S. was delayed until around 1910. In addition, by statistical analysis using coalescent theory, the major spread time for HCV infection in Japan occurred in the 1930s, whereas widespread dissemination of HCV in the U.S. occurred in the 1960s. These estimates of viral spread time are consistent with epidemiologic observations and predict that the burden of HCC in the U.S. will increase in the next two to three decades, possibly to equal that currently experienced in Japan.

Importance of Lineage-Specific Expansion of Plant Tandem Duplicates in the Adaptive Response to Environmental Stimuli

Plants have substantially higher gene duplication rates compared with most other eukaryotes. These plant gene duplicates are mostly derived from whole genome and/or tandem duplications. Earlier studies have shown that a large number of duplicate genes are retained over a long evolutionary time, and there is a clear functional bias in retention. However, the influence of duplication mechanism, particularly tandem duplication, on duplicate retention has not been thoroughly investigated. We have defined orthologous groups (OGs) between Arabidopsis (Arabidopsis thaliana) and three other land plants to examine the functional bias of retained duplicate genes during vascular plant evolution. Based on analysis of Gene Ontology categories, it is clear that genes in OGs that expanded via tandem duplication tend to be involved in responses to environmental stimuli, while those that expanded via nontandem mechanisms tend to have intracellular regulatory roles. Using Arabidopsis stress expression data, we further demonstrated that tandem duplicates in expanded OGs are significantly enriched in genes that are up-regulated by biotic stress conditions. In addition, tandem duplication of genes in an OG tends to be highly asymmetric. That is, expansion of OGs with tandem genes in one organismal lineage tends to be coupled with losses in the other. This is consistent with the notion that these tandem genes have experienced lineage-specific selection. In contrast, OGs with genes duplicated via nontandem mechanisms tend to experience convergent expansion, in which similar numbers of genes are gained in parallel. Our study demonstrates that the expansion of gene families and the retention of duplicates in plants exhibit substantial functional biases that are strongly influenced by the mechanism of duplication. In particular, genes involved in stress responses have an elevated probability of retention in a single-lineage fashion following tandem duplication, suggesting that these tandem duplicates are likely important for adaptive evolution to rapidly changing environments.

Evolutionary History and Stress Regulation of Plant Receptor-Like Kinase/Pelle Genes

Receptor-Like Kinase (RLK)/Pelle genes play roles ranging from growth regulation to defense response, and the dramatic expansion of this family has been postulated to be crucial for plant-specific adaptations. Despite this, little is known about the history of or the factors that contributed to the dramatic expansion of this gene family. In this study, we show that expansion coincided with the establishment of land plants and that RLK/Pelle subfamilies were established early in land plant evolution. The RLK/Pelle family expanded at a significantly higher rate than other kinases, due in large part to expansion of a few subfamilies by tandem duplication. Interestingly, these subfamilies tend to have members with known roles in defense response, suggesting that their rapid expansion was likely a consequence of adaptation to fast-evolving pathogens. Arabidopsis (Arabidopsis thaliana) expression data support the importance of RLK/Pelles in biotic stress response. We found that hundreds of RLK/Pelles are up-regulated by biotic stress. Furthermore, stress responsiveness is correlated with the degree of tandem duplication in RLK/Pelle subfamilies. Our findings suggest a link between stress response and tandem duplication and provide an explanation for why a large proportion of the RLK/Pelle gene family is found in tandem repeats. In addition, our findings provide a useful framework for potentially predicting RLK/Pelle stress functions based on knowledge of expansion pattern and duplication mechanism. Finally, we propose that the detection of highly variable molecular patterns associated with specific pathogens/parasites is the main reason for the up-regulation of hundreds of RLK/Pelles under biotic stress.

A Large Variation in the Rates of Synonymous Substitution for RNA Viruses and Its Relationship to a Diversity of Viral Infection and Transmission Modes

Abstract RNA viruses successfully adapt to various environments by repeatedly producing new mutants, often through generating a number of nucleotide substitutions. To estimate the degree of variation in mutation rates of RNA viruses and to understand the source of such variation, we studied the synonymous substitution rate because synonymous substitution is exempt from functional constraints at the protein level, and its rate reflects the mutation rate to a great extent. We estimated the synonymous substitution rates for a total of 49 different species of RNA viruses, and we found that the rates had tremendous variation by 5 orders of magnitude (from 1.3 × 10−7 to 6.2 × 10−2 /synonymous site/year). Comparing the synonymous substitution rates with the replication frequencies and replication error rates for the RNA viruses, we found that the main source of the rate variation was differences in the replication frequency because the rates of replication error were roughly constant over different RNA viruses. Moreover, we examined a relationship between viral life strategies and synonymous substitution rates to understand which viral life strategies affect replication frequencies. The results show that the variation of synonymous substitution rates has been influenced most by either the difference in the infection modes or the differences in the transmission modes. In conclusion, the variation of mutation rates for RNA viruses is caused by different replication frequencies, which are affected strongly by the infection and transmission modes.

An evolutionary view of functional diversity in family 1 glycosyltransferases.

Glycosyltransferases (GTs) (EC 2.4.x.y) catalyze the transfer of sugar moieties to a wide range of acceptor molecules, such as sugars, lipids, proteins, nucleic acids, antibiotics and other small molecules, including plant secondary metabolites. These enzymes can be classified into at least 92 families, of which family 1 glycosyltransferases (GT1), often referred to as UDP glycosyltransferases (UGTs), is the largest in the plant kingdom. To understand how UGTs expanded in both number and function during evolution of land plants, we screened genome sequences from six plants (Physcomitrella patens, Selaginella moellendorffii, Populus trichocarpa, Oryza sativa, Arabidopsis thaliana and Arabidopsis lyrata) for the presence of a conserved UGT protein domain. Phylogenetic analyses of the UGT genes revealed a significant expansion of UGTs, with lineage specificity and a higher duplication rate in vascular plants after the divergence of Physcomitrella. The UGTs from the six species fell into 24 orthologous groups that contained genes derived from the common ancestor of these six species. Some orthologous groups contained multiple UGT families with known functions, suggesting that UGTs discriminate compounds as substrates in a lineage-specific manner. Orthologous groups containing only a single UGT family tend to play a crucial role in plants, suggesting that such UGT families may have not expanded because of evolutionary constraints.

Genome-wide suppression of aberrant mRNA-like noncoding RNAs by NMD in Arabidopsis

The nonsense-mediated mRNA decay (NMD) pathway is a well-known eukaryotic surveillance mechanism that eliminates aberrant mRNAs that contain a premature termination codon (PTC). The UP-Frameshift (UPF) proteins, UPF1, UPF2, and UPF3, are essential for normal NMD function. Several NMD substrates have been identified, but detailed information on NMD substrates is lacking. Here, we noticed that, in Arabidopsis, most of the mRNA-like nonprotein-coding RNAs (ncRNAs) have the features of an NMD substrate. We examined the expression profiles of 2 Arabidopsis mutants, upf1-1 and upf3-1, using a whole-genome tiling array. The results showed that expression of not only protein-coding transcripts but also many mRNA-like ncRNAs (mlncRNAs), including natural antisense transcript RNAs (nat-RNAs) transcribed from the opposite strands of the coding strands, were up-regulated in both mutants. The percentage of the up-regulated mlncRNAs to all expressed mlncRNAs was much higher than that of the up-regulated protein-coding transcripts to all expressed protein- coding transcripts. This finding demonstrates that one of the most important roles of NMD is the genome-wide suppression of the aberrant mlncRNAs including nat-RNAs.

The Origin and Evolution of Porcine Reproductive and Respiratory Syndrome Viruses

Abstract Porcine reproductive and respiratory syndrome viruses (PRRSV) are divided into North American and European types, which show about 40% difference in their amino acid sequences. The divergence time of these two types has been estimated to be about 1980 from epidemiological data. This suggested that PRRSV have evolved at a higher evolutionary rate (order of 10−2/site/year) compared with other RNA viruses of 10−3 to 10−5/site/year. Here, to test the evolutionary history of PRRSV speculated by the epidemiological background, we estimated the divergence time and evolutionary rate of PRRSV with molecular evolutionary analysis. Estimated divergence time (1972–1988) corresponded well to that estimated by the epidemiological data, and the evolutionary rate (4.71–9.8) × 10−2 of PRRSV was indeed the highest among RNA viruses so far reported. Furthermore, we inferred important sites for the adaptation in order to examine how PRRSV have adapted to swine since they emerged. The adaptive sites were located not only in the epitopes related to immunity but also in the transmembrane regions including a signal peptide. In particular, the adaptive sites in the transmembrane regions were considered to affect compatibility to the host cell membrane. We conclude that PRRSV were transmitted from another host species to swine in about 1980 and have adapted to swine by altering the transmembrane regions.

Small open reading frames associated with morphogenesis are hidden in plant genomes

It is likely that many small ORFs (sORFs; 30–100 amino acids) are missed when genomes are annotated. To overcome this limitation, we identified ∼8,000 sORFs with high coding potential in intergenic regions of the Arabidopsis thaliana genome. However, the question remains as to whether these coding sORFs play functional roles. Using a designed array, we generated an expression atlas for 16 organs and 17 environmental conditions among 7,901 identified coding sORFs. A total of 2,099 coding sORFs were highly expressed under at least one experimental condition, and 571 were significantly conserved in other land plants. A total of 473 coding sORFs were overexpressed; ∼10% (49/473) induced visible phenotypic effects, a proportion that is approximately seven times higher than that of randomly chosen known genes. These results indicate that many coding sORFs hidden in plant genomes are associated with morphogenesis. We believe that the expression atlas will contribute to further study of the roles of sORFs in plants.

Two glycosyltransferases involved in anthocyanin modification delineated by transcriptome independent component analysis in Arabidopsis thaliana

To identify candidate genes involved in Arabidopsis flavonoid biosynthesis, we applied transcriptome coexpression analysis and independent component analyses with 1388 microarray data from publicly available databases. Two glycosyltransferases, UGT79B1 and UGT84A2 were found to cluster with anthocyanin biosynthetic genes. Anthocyanin was drastically reduced in ugt79b1 knockout mutants. Recombinant UGT79B1 protein converted cyanidin 3-O-glucoside to cyanidin 3-O-xylosyl(1→2)glucoside. UGT79B1 recognized 3-O-glucosylated anthocyanidins/flavonols and uridine diphosphate (UDP)-xylose, but not 3,5-O-diglucosylated anthocyanidins, indicating that UGT79B1 encodes anthocyanin 3-O-glucoside: 2′′-O-xylosyltransferase. UGT84A2 is known to encode sinapic acid: UDP-glucosyltransferase. In ugt84a2 knockout mutants, a major sinapoylated anthocyanin was drastically reduced. A comparison of anthocyanin profiles in ugt84a knockout mutants indicated that UGT84A2 plays a major role in sinapoylation of anthocyanin, and that other UGT84As contribute the production of 1-O-sinapoylglucose to a lesser extent. These data suggest major routes from cyanidin 3-O-glucoside to the most highly modified cyanidin in the potential intricate anthocyanin modification pathways in Arabidopsis.

sORF finder

SUMMARY sORF finder is a program package for identifying small open reading frames (sORFs) with high-coding potential. This application allows the identification of coding sORFs according to the nucleotide composition bias among coding sequences and the potential functional constraint at the amino acid level through evaluation of synonymous and non-synonymous substitution rates. AVAILABILITY Online tools and source codes are freely available at http://evolver.psc.riken.jp/. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.