Analysis of the Genome and Transcriptome of Cryptococcus neoformans var. grubii Reveals Complex RNA Expression and Microevolution Leading to Virulence Attenuation

Cryptococcus neoformans is a pathogenic basidiomycetous yeast responsible for more than 600,000 deaths each year. It occurs as two serotypes (A and D) representing two varieties (i.e. grubii and neoformans, respectively). Here, we sequenced the genome and performed an RNA-Seq-based analysis of the C. neoformans var. grubii transcriptome structure. We determined the chromosomal locations, analyzed the sequence/structural features of the centromeres, and identified origins of replication. The genome was annotated based on automated and manual curation. More than 40,000 introns populating more than 99% of the expressed genes were identified. Although most of these introns are located in the coding DNA sequences (CDS), over 2,000 introns in the untranslated regions (UTRs) were also identified. Poly(A)-containing reads were employed to locate the polyadenylation sites of more than 80% of the genes. Examination of the sequences around these sites revealed a new poly(A)-site-associated motif (AUGHAH). In addition, 1,197 miscRNAs were identified. These miscRNAs can be spliced and/or polyadenylated, but do not appear to have obvious coding capacities. Finally, this genome sequence enabled a comparative analysis of strain H99 variants obtained after laboratory passage. The spectrum of mutations identified provides insights into the genetics underlying the micro-evolution of a laboratory strain, and identifies mutations involved in stress responses, mating efficiency, and virulence.

Rrp1B gene polymorphism (1307T > C) in metastatic progression of breast cancer

Rrp1B (ribosomal RNA processing1 homolog B) is a novel candidate metastasis modifier gene in breast cancer. Functional gene assays demonstrated that a physical and functional interaction existing between Rrp1b and metastasis modifier gene SIPA1 causes reduction in the tumor growth and metastatic potential. Ectopic expression of Rrp1B modulates various metastasis predictive extra cellular matrix (ECM) genes associated with tumor suppression. The aim of this study is to determine the functional significance of single nucleotide polymorphism (SNP) in human Rrp1B gene (1307 T > C; rs9306160) with breast cancer development and progression. The study consists of 493 breast cancer cases recruited from Nizam's Institute of Medical Sciences, Hyderabad, and 558 age-matched healthy female controls from rural and urban areas. Genomic DNA was isolated by non-enzymatic method. Genotyping was done by amplification refractory mutation system (ARMS-PCR) method. Genotypes were reconfirmed by sequencing and results were analyzed statistically. We have performed Insilco analysis to know the RNA secondary structure by using online tool m fold. The TT genotype and T allele frequencies of Rrp1B1307 T > C polymorphism were significantly elevated in breast cancer (chi (2); p = < 0.008) cases compared to controls under different genetic models. The presence of T allele had conferred 1.75-fold risk for breast cancer development (OR = 1.75; 95 % CI = 1.15-2.67). The frequency of TT genotype of Rrp1b 1307T > C polymorphism was significantly elevated in obese patients (chi (2); p = 0.008) and patients with advanced disease (chi (2); p = 0.01) and with increased tumor size (chi (2); p = 0.01). Moreover, elevated frequency of T allele was also associated with positive lymph node status (chi (2); p = 0.04) and Her2 negative receptor status (chi (2); p = 0.006). Presence of Rrp1b1307TT genotype and T allele confer strong risk for breast cancer development and progression.

The beta hairpin structure within ribosomal protein S5 mediates interplay between domains II and IV and regulates HCV IRES function

Translation initiation in Hepatitis C Virus (HCV) is mediated by Internal Ribosome Entry Site (IRES), which is independent of cap-structure and uses a limited number of canonical initiation factors. During translation initiation IRES-40S complex formation depends on high affinity interaction of IRES with ribosomal proteins. Earlier, it has been shown that ribosomal protein S5 (RPS5) interacts with HCV IRES. Here, we have extensively characterized the HCV IRES-RPS5 interaction and demonstrated its role in IRES function. Computational modelling and RNA-protein interaction studies demonstrated that the beta hairpin structure within RPS5 is critically required for the binding with domains II and IV. Mutations disrupting IRES-RPS5 interaction drastically reduced the 80S complex formation and the corresponding IRES activity. Computational analysis and UV cross-linking experiments using various IRES-mutants revealed interplay between domains II and IV mediated by RPS5. In addition, present study demonstrated that RPS5 interaction is unique to HCV IRES and is not involved in 40S-3 ` UTR interaction. Further, partial silencing of RPS5 resulted in preferential inhibition of HCV RNA translation. However, global translation was marginally affected by partial silencing of RPS5. Taken together, results provide novel molecular insights into IRES-RPS5 interaction and unravel its functional significance in mediating internal initiation of translation.

Is the cellular initiation of translation an exclusive property of the initiator tRNAs?

Translation of mRNAs is the primary function of the ribosomal machinery. Although cells allow for a certain level of translational errors/mistranslation (which may well be a strategic need), maintenance of the fidelity of translation is vital for the cellular function and fitness. The P-site bound initiator tRNA selects the start codon in an mRNA and specifies the reading frame. A direct P-site binding of the initiator tRNA is a function of its special structural features, ribosomal elements, and the initiation factors. A highly conserved feature of the 3 consecutive G:C base pairs (3GC pairs) in the anticodon stem of the initiator tRNAs is vital in directing it to the P-site. Mutations in the 3GC pairs diminish/abolish initiation under normal physiological conditions. Using molecular genetics approaches, we have identified conditions that allow initiation with the mutant tRNAs in Escherichia coli. During our studies, we have uncovered a novel phenomenon of in vivo initiation by elongator tRNAs. Here, we recapitulate how the cellular abundance of the initiator tRNA, and nucleoside modifications in rRNA are connected with the tRNA selection in the P-site. We then discuss our recent finding of how a conserved feature in the mRNA, the Shine-Dalgarno sequence, influences tRNA selection in the P-site.

Identification of Sex and Female’s Reproductive Stage in Commercial Fish Species through the Quantification of Ribosomal Transcripts in Gonads

The estimation of maturity and sex of fish stocks in European waters is a requirement of the EU Data Collection Framework as part of the policy to improve fisheries management. On the other hand, research on fish biology is increasingly focused in molecular approaches, researchers needing correct identification of fish sex and reproductive stage without necessarily having in house the histological know-how necessary for the task. Taking advantage of the differential gene transcription occurring during fish sex differentiation and gametogenesis, the utility of 5S ribosomal RNA (5S rRNA) and General transcription factor IIIA (gtf3a) in the molecular identification of sex and gametogenic stage was tested in different economically-relevant fish species from the Bay of Biscay. Gonads of 9 fish species (, Atlantic, Atlantic-chub and horse mackerel, blue whiting, bogue, European anchovy, hake and pilchard and megrim), collected from local commercial fishing vessels were histologically sexed and 5S and 18S rRNA concentrations were quantified by capillary electrophoresis to calculate a 5S/18S rRNA index. Degenerate primers permitted cloning and sequencing of gtf3a fragments in 7 of the studied species. 5S rRNA and gtf3a transcript levels, together with 5S/18S rRNA index, distinguished clearly ovaries from testis in all of the studied species. The values were always higher in females than in males. 5S/18S rRNA index values in females were always highest when fish were captured in early phases of ovary development whilst, in later vitellogenic stages, the values decreased significantly. In megrim and European anchovy, where gonads in different oogenesis stages were obtained, the 5S/18S rRNA index identified clearly gametogenic stage. This approach, to the sexing and the quantitative non-subjective identification of the maturity stage of female fish, could have multiple applications in the study of fish stock dynamics, fish reproduction and fecundity and fish biology in general.

Deep RNA sequencing at single base-pair resolution reveals high complexity of the rice transcriptome

Understanding the dynamics of eukaryotic transcriptome is essential for studying the complexity of transcriptional regulation and its impact on phenotype. However, comprehensive studies of transcriptomes at single base resolution are rare, even for modern organisms, and lacking for rice. Here, we present the first transcriptome atlas for eight organs of cultivated rice. Using high-throughput paired-end RNA-seq, we unambiguously detected transcripts expressing at an extremely low level, as well as a substantial number of novel transcripts, exons, and untranslated regions. An analysis of alternative splicing in the rice transcriptome revealed that alternative cis-splicing occurred in similar to 33% of all rice genes. This is far more than previously reported. In addition, we also identified 234 putative chimeric transcripts that seem to be produced by trans-splicing, indicating that transcript fusion events are more common than expected. In-depth analysis revealed a multitude of fusion transcripts that might be by-products of alternative splicing. Validation and chimeric transcript structural analysis provided evidence that some of these transcripts are likely to be functional in the cell. Taken together, our data provide extensive evidence that transcriptional regulation in rice is vastly more complex than previously believed.

Identification of a Giardia krr1 homolog gene and the secondarily anucleolate condition of Giaridia lamblia

Giaridia lamblia was long considered to be one of the most primitive eukaryotes and to lie close to the transition between prokaryotes and eukaryotes, but several supporting features, such as lack of mitochondrion and Golgi, have been challenged recently. It was also reported previously that G. lamblia lacked nucleolus, which is the site of pre-rRNA processing and ribosomal assembling in the other eukaryotic cells. Here, we report the identification of the yeast homolog gene, krr1, in the anucleolate eukaryote, G. lamblia. The krr1 gene, encoding one of the pre-rRNA processing proteins in yeast, is actively transcribed in G. lamblia. The deduced protein sequence of G. lamblia krr1 is highly similar to yeast KRR1p that contains a single-KH domain. Our database searches indicated that krr1 genes actually present in diverse eukaryotes and also seem to present in Archaea. However, only the eukaryotic homologs, including that of G. lamblia, have the single-KH domain, which contains the conserved motif KR(K)R. Fibrillarin, another important pre-rRNA processing protein has also been identified previously in G. lamblia. Moreover, our database search shows that nearly half of the other nucleolus-localized protein genes of eukaryotic cells also have their homologs in Giardia. Therefore, we suggest that a common mechanism of pre-RNA processing may operate in the anucleolate eukaryote G. lamblia and in the other eukaryotes and that like the case of "lack of mitochondrion," "lack of nucleolus" may not be a primitive feature, but a secondarily evolutionary condition of the parasite.

Phylogeny of east Asian bufonids inferred from mitochondrial DNA sequences (Anura : Amphibia)

We investigated the relationships of Asian bufonids using partial sequences of mitochondrial DNA genes. Twenty-six samples representing 14 species of Bufo from China and Vietnam and 2 species of Torrentophryne from China were examined. Three samples of Bufo viridis from Armenia and Georgia were also sequenced to make a comparison to its sibling tetraploid species B. danatensis. Bufo americanus, from Canada, was used as the outgroup. Sequences from the 12S ribosomal RNA, 16S ribosomal RNA, cytochrome b, and the control region were analyzed using parsimony. East Asian bufonids were grouped into two major clades. One clade included B. andrewsi, B. bankorensis, B. gargarizans, B. tibetanus, B. tuberculatus, its sister clade B. cryptotympanicus, and the 2 species of Torrentophryne. The second clade consisted of B. galeatus, B. himalayanus, B. melanostictus, and a new species from Vietnam. The placement of three taxa (B. raddei B. viridis, and its sister species, B. danatensis) was problematic. The genus Torrentophryne should be synonymized with Bufo to remove paraphyly. Because B. raddei does not belong to the clade that includes B. viridis and B. danatensis, it was removed from the viridis species group. The species status of B bankorensis from Taiwan is evaluated. (C) 2000 Academic Press.

Origin and evolution of the RIG-I like RNA helicase gene family

Background: The DExD/H domain containing RNA helicases such as retinoic acid-inducible gene I (RIG-I) and melanoma differentiation-associated gene 5 (MDA5) are key cytosolic pattern recognition receptors (PRRs) for detecting nucleotide pathogen associated molecular patterns (PAMPs) of invading viruses. The RIG-I and MDA5 proteins differentially recognise conserved PAMPs in double stranded or single stranded viral RNA molecules, leading to activation of the interferon system in vertebrates. They share three core protein domains including a RNA helicase domain near the C terminus (HELICc), one or more caspase activation and recruitment domains (CARDs) and an ATP dependent DExD/H domain. The RIG-I/MDA5 directed interferon response is negatively regulated by laboratory of genetics and physiology 2 (LGP2) and is believed to be controlled by the mitochondria antiviral signalling protein (MAVS), a CARD containing protein associated with mitochondria. Results: The DExD/H containing RNA helicases including RIG-I, MDA5 and LGP2 were analysed in silico in a wide spectrum of invertebrate and vertebrate genomes. The gene synteny of MDA5 and LGP2 is well conserved among vertebrates whilst conservation of the gene synteny of RIG-I is less apparent. Invertebrate homologues had a closer phylogenetic relationship with the vertebrate RIG-Is than the MDA5/LGP2 molecules, suggesting the RIG-I homologues may have emerged earlier in evolution, possibly prior to the appearance of vertebrates. Our data suggest that the RIG-I like helicases possibly originated from three distinct genes coding for the core domains including the HELICc, CARD and ATP dependent DExD/H domains through gene fusion and gene/domain duplication. Furthermore, presence of domains similar to a prokaryotic DNA restriction enzyme III domain (Res III), and a zinc finger domain of transcription factor (TF) IIS have been detected by bioinformatic analysis. Conclusion: The RIG-I/MDA5 viral surveillance system is conserved in vertebrates. The RIG-I like helicase family appears to have evolved from a common ancestor that originated from genes encoding different core functional domains. Diversification of core functional domains might be fundamental to their functional divergence in terms of recognition of different viral PAMPs.

Variation patterns of the mitochondrial 16S rRNA gene with secondary structure constraints and their application to phylogeny of cyprinine fishes (Teleostei : Cypriniformes)

The mitochondrial 16S ribosomal RNA (rRNA) gene sequences from 93 cyprinid fishes were examined to reconstruct the phylogenetic relationships within the diverse and economically important subfamily Cyprininae. Within the subfamily a biased nucleotide composition (A > T, C > G) was observed in the loop regions of the gene, and in stem regions apparent selective pressures of base pairing showed a bias in favor of G over C and T over A. The bias may be associated with transition-transversion bias. Rates of nucleotide substitution were lower in stems than in loops. Analysis of compensatory substitutions across these taxa demonstrates 68% covariation in the gene and a logical weighting factor to account for dependence in mutations for phylogenetic inference should be 0.66. Comparisons of varied stem-loop weighting schemes indicate that the down-weightings for stem regions could improve the phylogenetic analysis and the degree of non-independence of stem substitutions was not as important as expected. Bayesian inference under four models of nucleotide substitution indicated that likelihood-based phylogenetic analyses were more effective in improving the phylogenetic performance than was weighted parsimony analysis. In Bayesian analyses, the resolution of phylogenies under the 16-state models for paired regions, incorporating GTR + G + I models for unpaired regions was better than those under other models. The subfamily Cyprininae was resolved as a monophyletic group, as well as tribe Labein and several genera. However, the monophyly of the currently recognized tribes, such as Schizothoracin, Barbin, Cyprinion + Onychostoma lineages, and some genera was rejected. Furthermore, comparisons of the parsimony and Bayesian analyses and results of variable length bootstrap analysis indicates that the mitochondrial 16S rRNA gene should contain important character variation to recover well-supported phylogeny of cyprinid taxa whose divergences occurred within the recent 8 MY, but could not provide resolution power for deep phylogenies spanning 10-19 MYA. (c) 2008 Published by Elsevier Inc.

Molecular phylogenetics of the family Cyprinidae (Actinopterygii : Cypriniformes) as evidenced by sequence variation in the first intron of S7 ribosomal protein-coding gene: Further evidence from a nuclear gene of the systematic chaos in the family

The family Cyprinidae is the largest freshwater fish group in the world, including over 200 genera and 2100 species. The phylogenetic relationships of major clades within this family are simply poorly understood, largely because of the overwhelming diversity of the group; however, several investigators have advanced different hypotheses of relationships that pre- and post-date the use of shared-derived characters as advocated through phylogenetic systematics. As expected, most previous investigations used morphological characters. Recently, mitochondrial DNA (mtDNA) sequences and combined morphological and mtDNA investigations have been used to explore and advance our understanding of species relationships and test monophyletic groupings. Limitations of these studies include limited taxon sampling and a strict reliance upon maternally inherited mtDNA variation. The present study is the first endeavor to recover the phylogenetic relationships of the 12 previously recognized monophyletic subfamilies within the Cyprinidae using newly sequenced nuclear DNA (nDNA) for over 50 species representing members of the different previously hypothesized subfamily and family groupings within the Cyprinidae and from other cypriniform families as outgroup taxa. Hypothesized phylogenetic relationships are constructed using maximum parsimony and Basyesian analyses of 1042 sites, of which 971 sites were variable and 790 were phylogenetically informative. Using other appropriate cypriniform taxa of the families Catostomidae (Myxocyprinus asiaticus), Gyrinocheilidae (Gyrinocheilus aymonieri), and Balitoridae (Nemacheilus sp. and Beaufortia kweichotvensis) as outgroups, the Cyprinidae is resolved as a monophyletic group. Within the family the genera Raiamas, Barilius, Danio, and Rasbora, representing many of the tropical cyprinids, represent basal members of the family. All other species can be classified into variably supported and resolved monophyletic lineages, depending upon analysis, that are consistent with or correspond to Barbini and Leuciscini. The Barbini includes taxa traditionally aligned with the subfamily Cyprininae sensu previous morphological revisionary studies by Howes (Barbinae, Labeoninae, Cyprininae and Schizothoracinae). The Leuciscini includes six other subfamilies that are mainly divided into three separate lineages. The relationships among genera and subfamilies are discussed as well as the possible origins of major lineages. (c) 2008 Published by Elsevier Inc.

The complete mitochondrial genome of the helmet catfish Cranoglanis bouderius (Silurifonnes : Cranoglanididae) and the phylogeny of otophysan fishes

The complete sequence of the 16,539 nucleotide mitochondrial genome from the single species of the catfish family Cranoglanididae, the helmet catfish Cranoglanis bouderius, was determined using the long and accurate polymerase chain reaction (LA PCR) method. The nucleotide sequences of C. bouderius mitochondrial DNA have been compared with those of three other catfish species in the same order. The contents of the C. bouderius mitochondrial genome are 13 protein-coding genes, two ribosomal RNA and 22 transfer RNA genes, and a non-coding control region, the gene order of which is identical to that observed in most other vertebrates. Phylogenetic analyses for 13 otophysan fishes were performed using Bayesian method based on the concatenated mtDNA protein-coding gene sequence and the individual protein-coding gene sequence data set. The competing otophysan topologies were then tested by using the approximately unbiased test, the Kishino-Hasegawa test, and the Shimodaira-Hasegawa test. The results show that the grouping ((((Characifonnes, Gymnotiformes), Siluriformes), Cyprinifionnes), outgroup) is the most likely but there is no significant difference between this one and the other alternative hypotheses. In addition, the phylogenetic placement of the family Cranoglanididae among siluriform families was also discussed. (c) 2006 Elsevier B.V. All rights reserved.

Efficient RNA interference in zebrafish embryos using siRNA synthesized with SP6 RNA polymerase

Double-stranded RNA (dsRNA) has been shown to be a useful tool for silencing genes in zebrafish (Danio rerio), while the blocking specificity of dsRNA is still of major concern for application. It was reported that siRNA (small interfering RNA) prepared by endoribonuclease digestion (esiRNA) could efficiently silence endogenous gene expression in mammalian embryos. To test whether esiRNA could work in zebrafish, we utilized Escherichia coli RNaseIII to digest dsRNA of zebrafish no tail (ntl), a mesoderm determinant in zebrafish and found that esi-ntl could lead to developmental defects, however, the effective dose was so close to the toxic dose that esi-ntl often led to non-specific developmental defects. Consequently, we utilized SP6 RNA polymerase to produce si-ntl, siRNA designed against ntl, by in vitro transcription. By injecting in vitro synthesized si-ntl into zebrafish zygotes, we obtained specific phenocopies of reported mutants of ntl. We achieved up to a 59%no tail phenotype when the injection concentration was as high as 4 mu g/mu L. Quantitative reverse transcription-polymerase chain reaction (RT-PCR) and whole-mount in situ hybridization analysis showed that si-ntl could largely and specifically reduce mRNA levels of the ntl gene. As a result, our data indicate that esiRNA is unable to cause specific developmental defects in zebrafish, while siRNA should be an alternative for downregulation of specific gene expression in zebrafish in cases where RNAi techniques are applied to zebrafish reverse genetics.

Involvement of p90 ribosomal S6 kinase in termination of cell cycle arrest during development of Artemia-encysted embryos

Artemia has evolved a unique developmental pattern of encysted embryos to cope with various environmental threats. Cell divisions totally cease during the preemergence developmental stage from gastrula to prenauplius. The molecular mechanism of this, however, remains unknown. Our study focuses on the involvement of p90 ribosomal S6 kinase (RSK), a family of serine/threonine kinase-mediating signal transduction downstream of mitogen-activated protein kinase cascades, in the termination of cell cycle arrest during the post-embryonic development of Artemia-encysted gastrula. With immunochemistry, morphology, and cell cycle analysis, the identified Artemia RSK was established to be specifically activated during the post-embryonic and early larval developmental stages when arrested cells of encysted embryos resumed mitoses. In vivo knockdown of RSK activity by RNA interference, kinase inhibition, and antibody neutralization consistently induced defective larvae with distinct gaps between the exoskeleton and internal tissues. In these abnormal individuals, mitoses were detected to be largely inhibited in the affected regions. These results display the requirement of RSK activity during Artemia development and suggest its role in termination of cell cycle (G(2)/M phase) arrest and promotion of mitogenesis. Our findings may, thus, provide insights into the regulation of cell division during Artemia post-embryonic development and reveal further aspects of RSK functions.