NF-κB-driven suppression of FOXO3a contributes to EGFR mutation-independent gefitinib resistance.

Therapy with epidermal growth factor receptor (EGFR) tyrosine kinase inhibitors (EGFR-TKIs, such as gefitinib or erlotinib) significantly prolongs survival time for patients with tumors harboring an activated mutation on EGFR; however, up to 40% of lung cancer patients exhibit acquired resistance to EGFR-TKIs with an unknown mechanism. FOXO3a, a transcription factor of the forkhead family, triggers apoptosis, but the mechanistic details involved in EGFR-TKI resistance and cancer stemness remain largely unclear. Here, we observed that a high level of FOXO3a was correlated with EGFR mutation-independent EGFR-TKI sensitivity, the suppression of cancer stemness, and better progression-free survival in lung cancer patients. The suppression of FOXO3a obviously increased gefitinib resistance and enhanced the stem-like properties of lung cancer cells; consistent overexpression of FOXO3a in gefitinib-resistant lung cancer cells reduced these effects. Moreover, we identified that miR-155 targeted the 3'UTR of FOXO3a and was transcriptionally regulated by NF-κB, leading to repressed FOXO3a expression and increased gefitinib resistance, as well as enhanced cancer stemness of lung cancer in vitro and in vivo. Our findings indicate that FOXO3a is a significant factor in EGFR mutation-independent gefitinib resistance and the stemness of lung cancer, and suggest that targeting the NF-κB/miR-155/FOXO3a pathway has potential therapeutic value in lung cancer with the acquisition of resistance to EGFR-TKIs.

Are all of the human exons alternatively spliced?

Alternative mRNA splicing (AS) is a major mechanism for increasing regulatory complexity. A key concept in AS is the distinction between alternatively and constitutively spliced exons (ASEs and CSEs, respectively). ASEs and CSEs have been reported to be differentially regulated, and to have distinct biological properties. However, the recent flood of RNA-sequencing data has obscured the boundary between ASEs and CSEs. Researchers are beginning to question whether 'authentic CSEs' do exist, and whether the ASE/CSE distinction is biologically invalid. Here, I examine the influences of increasing transcriptome data on the human ASE/CSE classification and our past understanding of the properties of these two types of exons. Interestingly, although the percentage of human ASEs has increased dramatically in recent years, the overall distinction between ASEs and CSEs remain valid. For example, CSEs are longer, evolve more slowly, and less frequently correspond to intrinsically disordered protein regions than ASEs. In addition, only a relatively small number of human genes have their transcripts composed entirely of ASEs despite the large amount of high-throughput transcriptome information. Therefore, the 'backbone' concept of AS, in which CSEs constitute the invariant part and ASEs the flexible part of the transcript, appears to be generally true despite the increasing percentage of ASEs in the human exome.

DNA methylation is associated with an increased level of conservation at nondegenerate nucleotides in mammals.

DNA methylation at CpG dinucleotides can significantly increase the rate of cytosine-to-thymine mutations and the level of sequence divergence. Although the correlations between DNA methylation and genomic sequence evolution have been widely studied, an unaddressed yet fundamental question is how DNA methylation is associated with the conservation of individual nucleotides in different sequence contexts. Here, we demonstrate that in mammalian exons, the correlations between DNA methylation and the conservation of individual nucleotides are dependent on the type of exonic sequence (coding or untranslated), the degeneracy of coding nucleotides, background selection pressure, and the relative position (first or nonfirst exon in the transcript) where the nucleotides are located. For untranslated and nonzero-fold degenerate nucleotides, methylated sites are less conserved than unmethylated sites regardless of background selection pressure and the relative position of the exon. For zero-fold degenerate (or nondegenerate) nucleotides, however, the reverse trend is observed in nonfirst coding exons and first coding exons that are under stringent background selection pressure. Furthermore, cytosine-to-thymine mutations at methylated zero-fold degenerate nucleotides are predicted to be more detrimental than those that occur at unmethylated nucleotides. As zero-fold and nonzero-fold degenerate nucleotides are very close to each other, our results suggest that the "functional resolution" of DNA methylation may be finer than previously recognized. In addition, the positive correlation between CpG methylation and the level of conservation at zero-fold degenerate nucleotides implies that CpG methylation may serve as an "indicator" of functional importance of these nucleotides.

Position-dependent correlations between DNA methylation and the evolutionary rates of mammalian coding exons.

DNA cytosine methylation is a central epigenetic marker that is usually mutagenic and may increase the level of sequence divergence. However, methylated genes have been reported to evolve more slowly than unmethylated genes. Hence, there is a controversy on whether DNA methylation is correlated with increased or decreased protein evolutionary rates. We hypothesize that this controversy has resulted from the differential correlations between DNA methylation and the evolutionary rates of coding exons in different genic positions. To test this hypothesis, we compare human-mouse and human-macaque exonic evolutionary rates against experimentally determined single-base resolution DNA methylation data derived from multiple human cell types. We show that DNA methylation is significantly related to within-gene variations in evolutionary rates. First, DNA methylation level is more strongly correlated with C-to-T mutations at CpG dinucleotides in the first coding exons than in the internal and last exons, although it is positively correlated with the synonymous substitution rate in all exon positions. Second, for the first exons, DNA methylation level is negatively correlated with exonic expression level, but positively correlated with both nonsynonymous substitution rate and the sample specificity of DNA methylation level. For the internal and last exons, however, we observe the opposite correlations. Our results imply that DNA methylation level is differentially correlated with the biological (and evolutionary) features of coding exons in different genic positions. The first exons appear more prone to the mutagenic effects, whereas the other exons are more influenced by the regulatory effects of DNA methylation.

The evolution of the coding exome of the Arabidopsis species--the influences of DNA methylation, relative exon position, and exon length.

BACKGROUND: The evolution of the coding exome is a major driving force of functional divergence both between species and between protein isoforms. Exons at different positions in the transcript or in different transcript isoforms may (1) mutate at different rates due to variations in DNA methylation level; and (2) serve distinct biological roles, and thus be differentially targeted by natural selection. Furthermore, intrinsic exonic features, such as exon length, may also affect the evolution of individual exons. Importantly, the evolutionary effects of these intrinsic/extrinsic features may differ significantly between animals and plants. Such inter-lineage differences, however, have not been systematically examined. RESULTS: Here we examine how DNA methylation at CpG dinucleotides (CpG methylation), in the context of intrinsic exonic features (exon length and relative exon position in the transcript), influences the evolution of coding exons of Arabidopsis thaliana. We observed fairly different evolutionary patterns in A. thaliana as compared with those reported for animals. Firstly, the mutagenic effect of CpG methylation is the strongest for internal exons and the weakest for first exons despite the stringent selective constraints on the former group. Secondly, the mutagenic effect of CpG methylation increases significantly with length in first exons but not in the other two exon groups. Thirdly, CpG methylation level is correlated with evolutionary rates (dS, dN, and the dN/dS ratio) with markedly different patterns among the three exon groups. The correlations are generally positive, negative, and mixed for first, last, and internal exons, respectively. Fourthly, exon length is a CpG methylation-independent indicator of evolutionary rates, particularly for dN and the dN/dS ratio in last and internal exons. Finally, the evolutionary patterns of coding exons with regard to CpG methylation differ significantly between Arabidopsis species and mammals. CONCLUSIONS: Our results suggest that intrinsic features, including relative exonic position in the transcript and exon length, play an important role in the evolution of A. thaliana coding exons. Furthermore, CpG methylation is correlated with exonic evolutionary rates differentially between A. thaliana and animals, and may have served different biological roles in the two lineages.

Alternative RNA structure-coupled gene regulations in tumorigenesis.

Alternative RNA structures (ARSs), or alternative transcript isoforms, are critical for regulating cellular phenotypes in humans. In addition to generating functionally diverse protein isoforms from a single gene, ARS can alter the sequence contents of 5'/3' untranslated regions (UTRs) and intronic regions, thus also affecting the regulatory effects of these regions. ARS may introduce premature stop codon(s) into a transcript, and render the transcript susceptible to nonsense-mediated decay, which in turn can influence the overall gene expression level. Meanwhile, ARS can regulate the presence/absence of upstream open reading frames and microRNA targeting sites in 5'UTRs and 3'UTRs, respectively, thus affecting translational efficiencies and protein expression levels. Furthermore, since ARS may alter exon-intron structures, it can influence the biogenesis of intronic microRNAs and indirectly affect the expression of the target genes of these microRNAs. The connections between ARS and multiple regulatory mechanisms underline the importance of ARS in determining cell fate. Accumulating evidence indicates that ARS-coupled regulations play important roles in tumorigenesis. Here I will review our current knowledge in this field, and discuss potential future directions.

Independent effects of alternative splicing and structural constraint on the evolution of mammalian coding exons.

Alternative splicing (AS) is known to significantly affect exon-level protein evolutionary rates in mammals. Particularly, alternatively spliced exons (ASEs) have a higher nonsynonymous-to-synonymous substitution rate (dN/dS) ratio than constitutively spliced exons (CSEs), possibly because the former are required only occasionally for normal biological functions. Meanwhile, intrinsically disordered regions (IDRs), the protein regions lacking fixed 3D structures, are also reported to have an increased evolutionary rate due to lack of structural constraint. Interestingly, IDRs tend to be located in alternative protein regions. Yet which of these two factors is the major determinant of the increased dN/dS in mammalian ASEs remains unclear. By comparing human-macaque and human-mouse one-to-one orthologous genes, we demonstrate that AS and protein structural disorder have independent effects on mammalian exon evolution. We performed analyses of covariance to demonstrate that the slopes of the (dN/dS-percentage of IDR) regression lines differ significantly between CSEs and ASEs. In other words, the dN/dS ratios of both ASEs and CSEs increase with the proportion of IDR (PIDR), whereas ASEs have higher dN/dS ratios than CSEs when they have similar PIDRs. Since ASEs and IDRs may less frequently overlap with protein domains (which also affect dN/dS), we also examined the correlations between dN/dS ratio and exon type/PIDR by controlling for the density of protein domain. We found that the effects of exon type and PIDR on dN/dS are both independent of domain density. Our results imply that nature can select for different biological features with regard to ASEs and IDRs, even though the two biological features tend to be localized in the same protein regions.

Assessing determinants of exonic evolutionary rates in mammals.

From studies investigating the differences in evolutionary rates between genes, gene compactness and gene expression level have been identified as important determinants of gene-level protein evolutionary rate, as represented by nonsynonymous to synonymous substitution rate (d(N)/d(S)) ratio. However, the causes of exon-level variances in d(N)/d(S) are less understood. Here, we use principal component regression to examine to what extent 13 exon features explain the variance in d(N), d(S), and the d(N)/d(S) ratio of human-rhesus macaque or human-mouse orthologous exons. The exon features were grouped into six functional categories: expression features, mRNA splicing features, structural-functional features, compactness features, exon duplicability, and other features, including G + C content and exon length. Although expression features are important for determining d(N) and d(N)/d(S) between exons of different genes, structural-functional features and splicing features explained more of the variance for exons of the same genes. Furthermore, we show that compactness features can explain only a relatively small percentage of variance in exon-level d(N) or d(N)/d(S) in either between-gene or within-gene comparison. By contrast, d(S) yielded inconsistent results in the human-mouse comparison and the human-rhesus macaque comparison. This inconsistency may suggest rapid evolutionary changes of the mutation landscape in mammals. Our results suggest that between-gene and within-gene variation in d(N)/d(S) (and d(N)) are driven by different evolutionary forces and that the role of mRNA splicing in causing the variation in evolutionary rates of coding sequences may be underappreciated.

GEMSiRV: a software platform for GEnome-scale metabolic model simulation, reconstruction and visualization.

MOTIVATION: Genome-scale metabolic network models have become an indispensable part of the increasingly important field of systems biology. Metabolic systems biology studies usually include three major components-network model construction, objective- and experiment-guided model editing and visualization, and simulation studies based mainly on flux balance analyses. Bioinformatics tools are required to facilitate these complicated analyses. Although some of the required functions have been served separately by existing tools, a free software resource that simultaneously serves the needs of the three major components is not yet available. RESULTS: Here we present a software platform, GEMSiRV (GEnome-scale Metabolic model Simulation, Reconstruction and Visualization), to provide functionalities of easy metabolic network drafting and editing, amenable network visualization for experimental data integration and flux balance analysis tools for simulation studies. GEMSiRV comes with downloadable, ready-to-use public-domain metabolic models, reference metabolite/reaction databases and metabolic network maps, all of which can be input into GEMSiRV as the starting materials for network construction or simulation analyses. Furthermore, all of the GEMSiRV-generated metabolic models and analysis results, including projects in progress, can be easily exchanged in the research community. GEMSiRV is a powerful integrative resource that may facilitate the development of systems biology studies. AVAILABILITY: The software is freely available on the web at http://sb.nhri.org.tw/GEMSiRV.

An experimentally validated genome-scale metabolic reconstruction of Klebsiella pneumoniae MGH 78578, iYL1228.

Klebsiella pneumoniae is a Gram-negative bacterium of the family Enterobacteriaceae that possesses diverse metabolic capabilities: many strains are leading causes of hospital-acquired infections that are often refractory to multiple antibiotics, yet other strains are metabolically engineered and used for production of commercially valuable chemicals. To study its metabolism, we constructed a genome-scale metabolic model (iYL1228) for strain MGH 78578, experimentally determined its biomass composition, experimentally determined its ability to grow on a broad range of carbon, nitrogen, phosphorus and sulfur sources, and assessed the ability of the model to accurately simulate growth versus no growth on these substrates. The model contains 1,228 genes encoding 1,188 enzymes that catalyze 1,970 reactions and accurately simulates growth on 84% of the substrates tested. Furthermore, quantitative comparison of growth rates between the model and experimental data for nine of the substrates also showed good agreement. The genome-scale metabolic reconstruction for K. pneumoniae presented here thus provides an experimentally validated in silico platform for further studies of this important industrial and biomedical organism.

Changes in transcriptional orientation are associated with increases in evolutionary rates of enterobacterial genes.

BACKGROUND: Changes in transcriptional orientation ("CTOs") occur frequently in prokaryotic genomes. Such changes usually result from genomic inversions, which may cause a conflict between the directions of replication and transcription and an increase in mutation rate. However, CTOs do not always lead to the replication-transcription confrontation. Furthermore, CTOs may cause deleterious disruptions of operon structure and/or gene regulations. The currently existing CTOs may indicate relaxation of selection pressure. Therefore, it is of interest to investigate whether CTOs have an independent effect on the evolutionary rates of the affected genes, and whether these genes are subject to any type of selection pressure in prokaryotes. METHODS: Three closely related enterbacteria, Escherichia coli, Klebsiella pneumoniae and Salmonella enterica serovar Typhimurium, were selected for comparisons of synonymous (dS) and nonsynonymous (dN) substitution rate between the genes that have experienced changes in transcriptional orientation (changed-orientation genes, "COGs") and those that do not (same-orientation genes, "SOGs"). The dN/dS ratio was also derived to evaluate the selection pressure on the analyzed genes. Confounding factors in the estimation of evolutionary rates, such as gene essentiality, gene expression level, replication-transcription confrontation, and decreased dS at gene terminals were controlled in the COG-SOG comparisons. RESULTS: We demonstrate that COGs have significantly higher dN and dS than SOGs when a series of confounding factors are controlled. However, the dN/dS ratios are similar between the two gene groups, suggesting that the increase in dS can sufficiently explain the increase in dN in COGs. Therefore, the increases in evolutionary rates in COGs may be mainly mutation-driven. CONCLUSIONS: Here we show that CTOs can increase the evolutionary rates of the affected genes. This effect is independent of the replication-transcription confrontation, which is suggested to be the major cause of inversion-associated evolutionary rate increases. The real cause of such evolutionary rate increases remains unclear but is worth further explorations.

The genomic features that affect the lengths of 5' untranslated regions in multicellular eukaryotes.

BACKGROUND: The lengths of 5'UTRs of multicellular eukaryotes have been suggested to be subject to stochastic changes, with upstream start codons (uAUGs) as the major constraint to suppress 5'UTR elongation. However, this stochastic model cannot fully explain the variations in 5'UTR length. We hypothesize that the selection pressure on a combination of genomic features is also important for 5'UTR evolution. The ignorance of these features may have limited the explanatory power of the stochastic model. Furthermore, different selective constraints between vertebrates and invertebrates may lead to differences in the determinants of 5'UTR length, which have not been systematically analyzed. METHODS: Here we use a multiple linear regression model to delineate the correlation between 5'UTR length and the combination of a series of genomic features (G+C content, observed-to-expected (OE) ratios of uAUGs, upstream stop codons (uSTOPs), methylation-related CG/UG dinucleotides, and mRNA-destabilizing UU/UA dinucleotides) in six vertebrates (human, mouse, rat, chicken, African clawed frog, and zebrafish) and four invertebrates (fruit fly, mosquito, sea squirt, and nematode). The relative contributions of each feature to the variation of 5'UTR length were also evaluated. RESULTS: We found that 14%~33% of the 5'UTR length variations can be explained by a linear combination of the analyzed genomic features. The most important genomic features are the OE ratios of uSTOPs and G+C content. The surprisingly large weightings of uSTOPs highlight the importance of selection on upstream open reading frames (which include both uAUGs and uSTOPs), rather than on uAUGs per se. Furthermore, G+C content is the most important determinants for most invertebrates, but for vertebrates its effect is second to uSTOPs. We also found that shorter 5'UTRs are affected more by the stochastic process, whereas longer 5'UTRs are affected more by selection pressure on genomic features. CONCLUSIONS: Our results suggest that upstream open reading frames may be the real target of selection, rather than uAUGs. We also show that the selective constraints on genomic features of 5'UTRs differ between vertebrates and invertebrates, and between longer and shorter 5'UTRs. A more comprehensive model that takes these findings into consideration is needed to better explain 5'UTR length evolution.

Exploring the selective constraint on the sizes of insertions and deletions in 5' untranslated regions in mammals.

BACKGROUND: Small insertions and deletions ("indels" with size >/= 100 bp) whose lengths are not multiples of three (non-3n) are strongly constrained and depleted in protein-coding sequences. Such a constraint has never been reported in noncoding genomic regions. In 5'untranslated regions (5'UTRs) in mammalian genomes, upstream start codons (uAUGs) and upstream open reading frames (uORFs) can regulate protein translation. The presence of non-3n indels in uORFs can potentially disrupt the functions of these regulatory elements. We thus hypothesize that natural selection disfavors non-3n indels in 5'UTRs when these regulatory elements are present. RESULTS: We design the Indel Selection Index to measure the selective constraint on non-3n indels in 5'UTRs. The index controls for the genomic compositions of the analyzed 5'UTRs and measures the probability of non-3n indel depletion downstream of uAUGs. By comparing the experimentally supported transcripts of human-mouse orthologous genes, we demonstrate that non-3n indels downstream of two types of uAUGs (alternative translation initiation sites and the uAUGs of coding sequence-overlapping uORFs) are underrepresented. The results hold well regardless of differences in alignment tool, gene structures between human and mouse, or the criteria in selecting alternatively spliced isoforms used for the analysis. CONCLUSIONS: To our knowledge, this is the first study to demonstrate selective constraints on non-3n indels in 5'UTRs. Such constraints may be associated with the regulatory functions of uAUGs/uORFs in translational regulation or the generation of protein isoforms. Our study thus brings a new perspective to the evolution of 5'UTRs in mammals.

Phosphorylated and nonphosphorylated serine and threonine residues evolve at different rates in mammals.

Protein phosphorylation plays an important role in the regulation of protein function. Phosphorylated residues are generally assumed to be subject to functional constraint, but it has recently been suggested from a comparison of distantly related vertebrate species that most phosphorylated residues evolve at the rates consistent with the surrounding regions. To resolve the controversy, we infer the ancestral phosphoproteome of human and mouse to compare the evolutionary rates of phosphorylated and nonphosphorylated serine (S), threonine (T), and tyrosine (Y) residues. This approach enables accurate estimation of evolutionary rates as it does not assume deep conservation of phosphorylated residues. We show that phosphorylated S/T residues tend to evolve more slowly than nonphosphorylated S/T residues not only in disordered but also in ordered protein regions, indicating evolutionary conservation of phosphorylated S/T residues in mammals. Thus, phosphorylated S/T residues tend to be subject to stronger functional constraint than nonphosphorylated residues regardless of the protein regions in which they reside. In contrast, phosphorylated Y residues evolve at similar rates as nonphosphorylated ones. We also find that the human lineage has gained more phosphorylated T residues and lost fewer phosphorylated Y residues than the mouse lineage. The cause of the gain/loss imbalance remains a mystery but should be worth exploring.

Gene family size conservation is a good indicator of evolutionary rates.

The evolution of duplicate genes has been a topic of broad interest. Here, we propose that the conservation of gene family size is a good indicator of the rate of sequence evolution and some other biological properties. By comparing the human-chimpanzee-macaque orthologous gene families with and without family size conservation, we demonstrate that genes with family size conservation evolve more slowly than those without family size conservation. Our results further demonstrate that both family expansion and contraction events may accelerate gene evolution, resulting in elevated evolutionary rates in the genes without family size conservation. In addition, we show that the duplicate genes with family size conservation evolve significantly more slowly than those without family size conservation. Interestingly, the median evolutionary rate of singletons falls in between those of the above two types of duplicate gene families. Our results thus suggest that the controversy on whether duplicate genes evolve more slowly than singletons can be resolved when family size conservation is taken into consideration. Furthermore, we also observe that duplicate genes with family size conservation have the highest level of gene expression/expression breadth, the highest proportion of essential genes, and the lowest gene compactness, followed by singletons and then by duplicate genes without family size conservation. Such a trend accords well with our observations of evolutionary rates. Our results thus point to the importance of family size conservation in the evolution of duplicate genes.

CNVVdb: a database of copy number variations across vertebrate genomes.

SUMMARY: CNVVdb is a web interface for identification of putative copy number variations (CNVs) among 16 vertebrate species using the-same-species self-alignments and cross-species pairwise alignments. By querying genomic coordinates in the target species, all the potential paralogous/orthologous regions that overlap > or = 80-100% (adjustable) of the query sequences with user-specified sequence identity (> or = 60% ~ > or = 90%) are returned. Additional information is also given for the genes that are included in the returned regions, including gene description, alternatively spliced transcripts, gene ontology descriptions and other biologically important information. CNVVdb also provides information of pseudogenes and single nucleotide polymorphisms (SNPs) for the CNV-related genomic regions. Moreover, multiple sequence alignments of shared CNVs across species are also provided. With the combination of CNV, SNP, pseudogene and functional information, CNVVdb can be very useful for comparative and functional studies in vertebrates. AVAILABILITY: CNVVdb is freely accessible at (http://CNVVdb.genomics.sinica.edu.tw).

CAPIH: a Web interface for comparative analyses and visualization of host-HIV protein-protein interactions.

BACKGROUND: The Human Immunodeficiency Virus type one (HIV-1) is the major causing pathogen of the Acquired Immune Deficiency Syndrome (AIDS). A large number of HIV-1-related studies are based on three non-human model animals: chimpanzee, rhesus macaque, and mouse. However, the differences in host-HIV-1 interactions between human and these model organisms have remained unexplored. DESCRIPTION: Here we present CAPIH (Comparative Analysis of Protein Interactions for HIV-1), the first web-based interface to provide comparative information between human and the three model organisms in the context of host-HIV-1 protein interactions. CAPIH identifies genetic changes that occur in HIV-1-interacting host proteins. In a total of 1,370 orthologous protein sets, CAPIH identifies approximately 86,000 amino acid substitutions, approximately 21,000 insertions/deletions, and approximately 33,000 potential post-translational modifications that occur only in one of the four compared species. CAPIH also provides an interactive interface to display the host-HIV-1 protein interaction networks, the presence/absence of orthologous proteins in the model organisms in the networks, the genetic changes that occur in the protein nodes, and the functional domains and potential protein interaction hot sites that may be affected by the genetic changes. The CAPIH interface is freely accessible at http://bioinfo-dbb.nhri.org.tw/capih. CONCLUSION: CAPIH exemplifies that large divergences exist in disease-associated proteins between human and the model animals. Since all of the newly developed medications must be tested in model animals before entering clinical trials, it is advisable that comparative analyses be performed to ensure proper translations of animal-based studies. In the case of AIDS, the host-HIV-1 protein interactions apparently have differed to a great extent among the compared species. An integrated protein network comparison among the four species will probably shed new lights on AIDS studies.

INDELSCAN: a web server for comparative identification of species-specific and non-species-specific insertion/deletion events.

Insertion and deletion (indel) events usually have dramatic effects on genome structure and gene function. Species-specific indels have been demonstrated to be associated with species-unique traits. Currently, indel identifications mainly rely on pair-wise sequence alignments (the 'pair-wise indels'), which suffer lack of discrimination of species specificity and insertion versus deletion. Also, there is no freely accessible web server for genome-wide identification of indels. Therefore, we develop a web server--INDELSCAN--to identify four types of indels using multiple sequence alignments that include sequences from one target, one subject and > or =1 out-group species. The four types of indels identified encompass target species-specific, subject species-specific, non-species-specific and target-subject pair-wise indels. Insertions and deletions are discriminated with reference to out-group sequences. The genomic locations (5'UTR, intron, CDS, 3'UTR and intergenic region) of these indels are also provided for functional analysis. INDELSCAN provides genomic sequences and gene annotations from a wide spectrum of taxa for users to select from, including nine target species (human (Homo sapiens), mouse (Mus musculus), rat (Rattus norvegicus), dog (Canis familiaris), opossum (Monodelphis domestica), chicken (Gallus gallus), zebrafish (Danio rerio), fly (Drosophila melanogaster) and yeast (Saccharomyces cerevisiae) and >35 subject/out-group species, ranging from yeasts to mammals. The server also provides analytic figures and supports indel identification from user-uploaded alignments/annotations. INDELSCAN is freely accessible at http://indelscan.genomics.sinica.edu.tw/IndelScan/.

Inter Partes Review: Patent Killer No More?

Inter Partes review (IPR) can efficiently invalidate drug patents and potentially convey strategic advantages to follow-on drug makers. However, recent changes in the IPR system foretell a tectonic shift in the landscape. Here we summarize these major changes and discuss the implications for the biopharmaceutical community.

China: The Next Pharmacy of the World?

The manufacturing success in China has not yet radiated to the pharmaceutical sector. Recently, China released a policy guideline to foster its follow-on drug industry, revealing its ambition to become a great power in the field. Here, I briefly discuss how this guideline may change the industry landscape, and its associated challenges and opportunities.