Arthur Chun-Chieh Shih

Uncovering Small RNA-Mediated Responses to Phosphate Deficiency in Arabidopsis by Deep Sequencing

Recent studies have demonstrated the important role of plant microRNAs (miRNAs) under nutrient deficiencies. In this study, deep sequencing of Arabidopsis (Arabidopsis thaliana) small RNAs was conducted to reveal miRNAs and other small RNAs that were differentially expressed in response to phosphate (Pi) deficiency. About 3.5 million sequence reads corresponding to 0.6 to 1.2 million unique sequence tags from each Pi-sufficient or Pi-deficient root or shoot sample were mapped to the Arabidopsis genome. We showed that upon Pi deprivation, the expression of miR156, miR399, miR778, miR827, and miR2111 was induced, whereas the expression of miR169, miR395, and miR398 was repressed. We found cross talk coordinated by these miRNAs under different nutrient deficiencies. In addition to miRNAs, we identified one Pi starvation-induced DICER-LIKE1-dependent small RNA derived from the long terminal repeat of a retrotransposon and a group of 19-nucleotide small RNAs corresponding to the 5′ end of tRNA and expressed at a high level in Pi-starved roots. Importantly, we observed an increased abundance of TAS4-derived trans-acting small interfering RNAs (ta-siRNAs) in Pi-deficient shoots and uncovered an autoregulatory mechanism of PAP1/MYB75 via miR828 and TAS4-siR81(−) that regulates the biosynthesis of anthocyanin. This finding sheds light on the regulatory network between miRNA/ta-siRNA and its target gene. Of note, a substantial amount of miR399* accumulated under Pi deficiency. Like miR399, miR399* can move across the graft junction, implying a potential biological role for miR399*. This study represents a comprehensive expression profiling of Pi-responsive small RNAs and advances our understanding of the regulation of Pi homeostasis mediated by small RNAs.

Simultaneous amino acid substitutions at antigenic sites drive influenza A hemagglutinin evolution.

The HA1 domain of HA, the major antigenic protein of influenza A viruses, contains all of the antigenic sites of HA and is under continual immune-driven selection. To resolve controversies on whether only a few or many residue sites of HA1 have undergone positive selection, whether positive selection at HA1 is continual or punctuated, and whether antigenic change is punctuated, we introduce an approach to analyze 2,248 HA1 sequences collected from 1968 to 2005. We identify 95 substitutions at 63 sites from 1968 to 2005 and show that each substitution occurred very rapidly. The rapid substitution and the fact that 57 of the 63 sites are antigenic sites indicate that hitchhiking plays a minor role and that most of these sites, many more than previously found, have undergone positive selection. Strikingly, 88 of the 95 substitutions occurred in groups, and multiple mutations at antigenic sites sped up the fixation process. Our results suggest that positive selection has been ongoing most of the time, not sporadic, and that multiple mutations at antigenic sites cumulatively enhance antigenic drift, indicating that antigenic change is less punctuated than recently proposed.

Characterizing Regulatory and Functional Differentiation between Maize Mesophyll and Bundle Sheath Cells by Transcriptomic Analysis

To study the regulatory and functional differentiation between the mesophyll (M) and bundle sheath (BS) cells of maize (Zea mays), we isolated large quantities of highly homogeneous M and BS cells from newly matured second leaves for transcriptome profiling by RNA sequencing. A total of 52,421 annotated genes with at least one read were found in the two transcriptomes. Defining a gene with more than one read per kilobase per million mapped reads as expressed, we identified 18,482 expressed genes; 14,972 were expressed in M cells, including 53 M-enriched transcription factor (TF) genes, whereas 17,269 were expressed in BS cells, including 214 BS-enriched TF genes. Interestingly, many TF gene families show a conspicuous BS preference in expression. Pathway analyses reveal differentiation between the two cell types in various functional categories, with the M cells playing more important roles in light reaction, protein synthesis and folding, tetrapyrrole synthesis, and RNA binding, while the BS cells specialize in transport, signaling, protein degradation and posttranslational modification, major carbon, hydrogen, and oxygen metabolism, cell division and organization, and development. Genes coding for several transporters involved in the shuttle of C4 metabolites and BS cell wall development have been identified, to our knowledge, for the first time. This comprehensive data set will be useful for studying M/BS differentiation in regulation and function.

Anatomical and transcriptional dynamics of maize embryonic leaves during seed germination

Our anatomical analysis revealed that a dry maize seed contains four to five embryonic leaves at different developmental stages. Rudimentary kranz structure (KS) is apparent in the first leaf with a substantial density, but its density decreases toward younger leaves. Upon imbibition, leaf expansion occurs rapidly with new KSs initiated from the palisade-like ground meristem cells in the middle of the leaf. In parallel to the anatomical analysis, we obtained the time course transcriptomes for the embryonic leaves in dry and imbibed seeds every 6 h up to hour 72. Over this time course, the embryonic leaves exhibit transcripts of 30,255 genes at a level that can be regarded as “expressed.” In dry seeds, ∼25,500 genes are expressed, showing functional enrichment in transcription, RNA processing, protein synthesis, primary metabolic pathways, and calcium transport. During the 72-h time course, ∼13,900 genes, including 590 transcription factor genes, are differentially expressed. Indeed, by 30 h postimbibition, ∼2,200 genes expressed in dry seeds are already down-regulated, and ∼2,000 are up-regulated. Moreover, the top 1% expressed genes at 54 h or later are very different from those before 30 h, reflecting important developmental and physiological transitions. Interestingly, clusters of genes involved in hormone metabolism, signaling, and responses are differentially expressed at various time points and TF gene expression is also modular and stage specific. Our dataset provides an opportunity for hypothesizing the timing of regulatory actions, particularly in the context of KS development.

The prognostic significance of RUNX2 and miR-10a/10b and their inter-relationship in breast cancer

BackgroundThe major cancer related mortality is caused by metastasis and invasion. It is important to identify genes regulating metastasis and invasion in order to curtail metastatic spread of cancer cells.MethodsThis study investigated the association between RUNX2 and miR-10a/miR-10b and the risk of breast cancer relapse. Expression levels of RUNX2 and miR-10a/b in108 pairs of tumor and non-tumor tissue of breast cancer were assayed by quantitative PCR analysis and evaluated for their prognostic implications.ResultsThe median expression levels of RUNX2 and miR-10b in tumor tissue normalized using adjacent non-tumor tissue were significantly higher in relapsed patients than in relapse-free patients. Higher expression of these three genes were significantly correlated with the hazard ratio for breast cancer recurrence (RUNX2: 3.02, 95% CI = 1.50 ~ 6.07; miR-10a: 2.31, 95% CI = 1.00 ~ 5.32; miR-10b: 3.96, 95% CI = 1.21 ~ 12.98). The joint effect of higher expression of all three genes was associated with a hazard ratio of 12.37 (95% CI = 1.62 ~ 94.55) for relapse. In a breast cancer cell line, RUNX2 silencing reduced the expression of miR-10a/b and also impaired cell motility, while RUNX2 overexpression elicited opposite effects.ConclusionsThese findings indicate that higher expression of RUNX2 and miR-10a/b was associated with adverse outcome of breast cancer. Expression levels of RUNX2 and miR-10a/b individually or jointly are potential prognostic factors for predicting breast cancer recurrence. Data from in vitro studies support the notion that RUNX2 promoted cell motility by upregulating miR-10a/b.

Transcriptome dynamics of developing maize leaves and genomewide prediction of cis elements and their cognate transcription factors

Significance Maize is a major crop and a model plant for studying C4 leaf development. However, its regulatory network of leaf development is poorly understood. We used transcriptomes of developing leaves to study gene-expression dynamics and coexpression to reveal functional transition during maize leaf development. More significantly, we developed methods to predict transcription factor-binding sites (TFBSs) and their cognate transcription factors (TFs) or to use the known Arabidopsis TF–TFBS pairs to predict the maize TF–TFBS pairs. In total, we predicted 1,340 novel TFBSs and 253 new TF–TFBS pairs in maize. Twelve predicted TF–TFBS interactions were validated by functional tests, suggesting that our methods perform well. Our study has significantly expanded our knowledge of the regulatory network of maize leaf development. Maize is a major crop and a model plant for studying C4 photosynthesis and leaf development. However, a genomewide regulatory network of leaf development is not yet available. This knowledge is useful for developing C3 crops to perform C4 photosynthesis for enhanced yields. Here, using 22 transcriptomes of developing maize leaves from dry seeds to 192 h post imbibition, we studied gene up- and down-regulation and functional transition during leaf development and inferred sets of strongly coexpressed genes. More significantly, we developed a method to predict transcription factor binding sites (TFBSs) and their cognate transcription factors (TFs) using genomic sequence and transcriptomic data. The method requires not only evolutionary conservation of candidate TFBSs and sets of strongly coexpressed genes but also that the genes in a gene set share the same Gene Ontology term so that they are involved in the same biological function. In addition, we developed another method to predict maize TF–TFBS pairs using known TF–TFBS pairs in Arabidopsis or rice. From these efforts, we predicted 1,340 novel TFBSs and 253 new TF–TFBS pairs in the maize genome, far exceeding the 30 TF–TFBS pairs currently known in maize. In most cases studied by both methods, the two methods gave similar predictions. In vitro tests of 12 predicted TF–TFBS interactions showed that our methods perform well. Our study has significantly expanded our knowledge on the regulatory network involved in maize leaf development.

Assembler for de novo assembly of large genomes

Significance Assembling a large genome faces three challenges: assembly quality, computer memory requirement, and execution time. Our developed assembler, JR-Assembler, uses (a) a strategy that selects good seeds for contig construction, (b) an extension strategy that uses whole sequencing reads to increase the chance to jump over repeats and to expedite extension, and (c) detecting misassemblies by remapping reads to assembled sequences. Compared with current assemblers, JR-Assembler achieves a better overall assembly quality, requires less execution time and requires, with one exception, less memory. The advantages of JR-Assembler in memory usage and execution time will increase slowly as the read length increases. Thus, contrary to the prevailing view, the extension approach seems superior to the de Bruijn graph approach. Assembling a large genome using next generation sequencing reads requires large computer memory and a long execution time. To reduce these requirements, we propose an extension-based assembler, called JR-Assembler, where J and R stand for “jumping” extension and read “remapping.” First, it uses the read count to select good quality reads as seeds. Second, it extends each seed by a whole-read extension process, which expedites the extension process and can jump over short repeats. Third, it uses a dynamic back trimming process to avoid extension termination due to sequencing errors. Fourth, it remaps reads to each assembled sequence, and if an assembly error occurs by the presence of a repeat, it breaks the contig at the repeat boundaries. Fifth, it applies a less stringent extension criterion to connect low-coverage regions. Finally, it merges contigs by unused reads. An extensive comparison of JR-Assembler with current assemblers using datasets from small, medium, and large genomes shows that JR-Assembler achieves a better or comparable overall assembly quality and requires lower memory use and less central processing unit time, especially for large genomes. Finally, a simulation study shows that JR-Assembler achieves a superior performance on memory use and central processing unit time than most current assemblers when the read length is 150 bp or longer, indicating that the advantages of JR-Assembler over current assemblers will increase as the read length increases with advances in next generation sequencing technology.

Phylo-mLogo: an interactive and hierarchical multiple-logo visualization tool for alignment of many sequences

BackgroundWhen aligning several hundreds or thousands of sequences, such as epidemic virus sequences or homologous/orthologous sequences of some big gene families, to reconstruct the epidemiological history or their phylogenies, how to analyze and visualize the alignment results of many sequences has become a new challenge for computational biologists. Although there are several tools available for visualization of very long sequence alignments, few of them are applicable to the alignments of many sequences.ResultsA multiple-logo alignment visualization tool, called Phylo-mLogo, is presented in this paper. Phylo-mLogo calculates the variabilities and homogeneities of alignment sequences by base frequencies or entropies. Different from the traditional representations of sequence logos, Phylo-mLogo not only displays the global logo patterns of the whole alignment of multiple sequences, but also demonstrates their local homologous logos for each clade hierarchically. In addition, Phylo-mLogo also allows the user to focus only on the analysis of some important, structurally or functionally constrained sites in the alignment selected by the user or by built-in automatic calculation.ConclusionWith Phylo-mLogo, the user can symbolically and hierarchically visualize hundreds of aligned sequences simultaneously and easily check the changes of their amino acid sites when analyzing many homologous/orthologous or influenza virus sequences. More information of Phylo-mLogo can be found at URL http://biocomp.iis.sinica.edu.tw/phylomlogo.

Editing site analysis in a gymnosperm mitochondrial genome reveals similarities with angiosperm mitochondrial genomes.

Sequence analysis of organelle genomes and comprehensive analysis of C-to-U editing sites from flowering and non-flowering plants have provided extensive sequence information from diverse taxa. This study includes the first comprehensive analysis of RNA editing sites from a gymnosperm mitochondrial genome, and utilizes informatics analyses to determine conserved features in the RNA sequence context around editing sites. We have identified 565 editing sites in 21 full-length and 4 partial cDNAs of the 39 protein-coding genes identified from the mitochondrial genome of Cycas taitungensis. The information profiles and RNA sequence context of C-to-U editing sites in the Cycas genome exhibit similarity in the immediate flanking nucleotides. Relative entropy analyses indicate that similar regions in the 5′ flanking 20 nucleotides have information content compared to angiosperm mitochondrial genomes. These results suggest that evolutionary constraints exist on the nucleotide sequences immediately adjacent to C-to-U editing sites, and similar regions are utilized in editing site recognition.

SinicView: A visualization environment for comparisons of multiple nucleotide sequence alignment tools

BackgroundDeluged by the rate and complexity of completed genomic sequences, the need to align longer sequences becomes more urgent, and many more tools have thus been developed. In the initial stage of genomic sequence analysis, a biologist is usually faced with the questions of how to choose the best tool to align sequences of interest and how to analyze and visualize the alignment results, and then with the question of whether poorly aligned regions produced by the tool are indeed not homologous or are just results due to inappropriate alignment tools or scoring systems used. Although several systematic evaluations of multiple sequence alignment (MSA) programs have been proposed, they may not provide a standard-bearer for most biologists because those poorly aligned regions in these evaluations are never discussed. Thus, a tool that allows cross comparison of the alignment results obtained by different tools simultaneously could help a biologist evaluate their correctness and accuracy.ResultsIn this paper, we present a versatile alignment visualization system, called SinicView, (for Sequence-aligning INnovative and Interactive Comparison VIEWer), which allows the user to efficiently compare and evaluate assorted nucleotide alignment results obtained by different tools. SinicView calculates similarity of the alignment outputs under a fixed window using the sum-of-pairs method and provides scoring profiles of each set of aligned sequences. The user can visually compare alignment results either in graphic scoring profiles or in plain text format of the aligned nucleotides along with the annotations information. We illustrate the capabilities of our visualization system by comparing alignment results obtained by MLAGAN, MAVID, and MULTIZ, respectively.ConclusionWith SinicView, users can use their own data sequences to compare various alignment tools or scoring systems and select the most suitable one to perform alignment in the initial stage of sequence analysis.