Unveiling the molecular landscape of cognitive aging: insights from polygenic risk scores, DNA methylation, and gene expression.

BACKGROUND: Aging represents a significant risk factor for the occurrence of cerebral small vessel disease, associated with white matter (WM) lesions, and to age-related cognitive alterations, though the precise mechanisms remain largely unknown. This study aimed to investigate the impact of polygenic risk scores (PRS) for WM integrity, together with age-related DNA methylation, and gene expression alterations, on cognitive aging in a cross-sectional healthy aging cohort. The PRSs were calculated using genome-wide association study (GWAS) summary statistics for magnetic resonance imaging (MRI) markers of WM integrity, including WM hyperintensities, fractional anisotropy (FA), and mean diffusivity (MD). These scores were utilized to predict age-related cognitive changes and evaluate their correlation with structural brain changes, which distinguish individuals with higher and lower cognitive scores. To reduce the dimensionality of the data and identify age-related DNA methylation and transcriptomic alterations, Sparse Partial Least Squares-Discriminant Analysis (sPLS-DA) was used. Subsequently, a canonical correlation algorithm was used to integrate the three types of omics data (PRS, DNA methylation, and gene expression data) and identify an individual "omics" signature that distinguishes subjects with varying cognitive profiles. RESULTS: We found a positive association between MD-PRS and long-term memory, as well as a correlation between MD-PRS and structural brain changes, effectively discriminating between individuals with lower and higher memory scores. Furthermore, we observed an enrichment of polygenic signals in genes related to both vascular and non-vascular factors. Age-related alterations in DNA methylation and gene expression indicated dysregulation of critical molecular features and signaling pathways involved in aging and lifespan regulation. The integration of multi-omics data underscored the involvement of synaptic dysfunction, axonal degeneration, microtubule organization, and glycosylation in the process of cognitive aging. CONCLUSIONS: These findings provide valuable insights into the biological mechanisms underlying the association between WM coherence and cognitive aging. Additionally, they highlight how age-associated DNA methylation and gene expression changes contribute to cognitive aging.

Evaluation of the genetic risk for COVID-19 outcomes in COPD and differences among worldwide populations.

BACKGROUND: Populations seem to respond differently to the global pandemic of severe acute respiratory syndrome coronavirus 2. Recent studies show individual variability in both susceptibility and clinical response to COVID-19 infection. People with chronic obstructive pulmonary disease (COPD) constitute one of COVID-19 risk groups, being already associated with a poor prognosis upon infection. This study aims contributing to unveil the underlying reasons for such prognosis in people with COPD and the variability in the response observed across worldwide populations, by looking at the genetic background as a possible answer to COVID-19 infection response heterogeneity. METHODS: SNPs already associated with susceptibility to COVID-19 infection (rs286914 and rs12329760) and severe COVID-19 with respiratory failure (rs657152 and rs11385942) were assessed and their allelic frequencies used to calculate the probability of having multiple risk alleles. This was performed on a Portuguese case-control COPD cohort, previously clinically characterized and genotyped from saliva samples, and also on worldwide populations (European, Spanish, Italian, African, American and Asian), using publicly available frequencies data. A polygenic risk analysis was also conducted on the Portuguese COPD cohort for the two mentioned phenotypes, and also for hospitalization and survival to COVID-19 infection. FINDINGS: No differences in genetic risk for COVID-19 susceptibility, hospitalization, severity or survival were found between people with COPD and the control group (all p-values > 0.01), either considering risk alleles individually, allelic combinations or polygenic risk scores. All populations, even those with European ancestry (Portuguese, Spanish and Italian), showed significant differences from the European population in genetic risk for both COVID-19 susceptibility and severity (all p-values < 0.0001). CONCLUSION: Our results indicate a low genetic contribution for COVID-19 infection predisposition or worse outcomes observed in people with COPD. Also, our study unveiled a high genetic heterogeneity across major world populations for the same alleles, even within European sub-populations, demonstrating the need to build a higher resolution European genetic map, so that differences in the distribution of relevant alleles can be easily accessed and used to better manage diseases, ultimately, safeguarding populations with higher genetic predisposition to such diseases.

Redox Imbalance and Methylation Disturbances in Early Childhood Obesity.

Obesity is increasing worldwide in prepubertal children, reducing the age of onset of associated comorbidities, including type 2 diabetes. Sulfur-containing amino acids, methionine, cysteine, and their derivatives play important roles in the transmethylation and transsulfuration pathways. Dysregulation of these pathways leads to alterations in the cellular methylation patterns and an imbalanced redox state. Therefore, we tested the hypothesis that one-carbon metabolism is already dysregulated in prepubertal children with obesity. Peripheral blood was collected from 64 children, and the plasma metabolites from transmethylation and transsulfuration pathways were quantified by HPLC. The cohort was stratified by BMI z-scores and HOMA-IR indices into healthy lean (HL), healthy obese (HO), and unhealthy obese (UHO). Fasting insulin levels were higher in the HO group compared to the HL, while the UHO had the highest. All groups presented normal fasting glycemia. Furthermore, high-density lipoprotein (HDL) was lower while triglycerides and lactate levels were higher in the UHO compared to HO subjects. S-adenosylhomocysteine (SAH) and total homocysteine levels were increased in the HO group compared to HL. Additionally, glutathione metabolism was also altered. Free cystine and oxidized glutathione (GSSG) were increased in the HO as compared to HL subjects. Importantly, the adipocyte secretory function was already compromised at this young age. Elevated circulating leptin and decreased adiponectin levels were observed in the UHO as compared to the HO subjects. Some of these alterations were concomitant with alterations in the DNA methylation patterns in the obese group, independent of the impaired insulin levels. In conclusion, our study informs on novel and important metabolic alterations in the transmethylation and the transsulfuration pathways in the early stages of obesity. Moreover, the altered secretory function of the adipocyte very early in life may be relevant in identifying early metabolic markers of disease that may inform on the increased risk for specific future comorbidities in this population.

Plasmodium apicoplast tyrosyl-tRNA synthetase recognizes an unusual, simplified identity set in cognate tRNATyr.

The life cycle of Plasmodium falciparum, the agent responsible for malaria, depends on both cytosolic and apicoplast translation fidelity. Apicoplast aminoacyl-tRNA synthetases (aaRS) are bacterial-like enzymes devoted to organellar tRNA aminoacylation. They are all encoded by the nuclear genome and are translocated into the apicoplast only after cytosolic biosynthesis. Apicoplast aaRSs contain numerous idiosyncratic sequence insertions: An understanding of the roles of these insertions has remained elusive and they hinder efforts to heterologously overexpress these proteins. Moreover, the A/T rich content of the Plasmodium genome leads to A/U rich apicoplast tRNA substrates that display structural plasticity. Here, we focus on the P. falciparum apicoplast tyrosyl-tRNA synthetase (Pf-apiTyrRS) and its cognate tRNATyr substrate (Pf-apitRNATyr). Cloning and expression strategies used to obtain an active and functional recombinant Pf-apiTyrRS are reported. Functional analyses established that only three weak identity elements in the apitRNATyr promote specific recognition by the cognate Pf-apiTyrRS and that positive identity elements usually found in the tRNATyr acceptor stem are excluded from this set. This finding brings to light an unusual behavior for a tRNATyr aminoacylation system and suggests that Pf-apiTyrRS uses primarily negative recognition elements to direct tyrosylation specificity.

Adaptive Mistranslation Accelerates the Evolution of Fluconazole Resistance and Induces Major Genomic and Gene Expression Alterations in Candida albicans.

Regulated erroneous protein translation (adaptive mistranslation) increases proteome diversity and produces advantageous phenotypic variability in the human pathogen Candida albicans. It also increases fitness in the presence of fluconazole, but the underlying molecular mechanism is not understood. To address this question, we evolved hypermistranslating and wild-type strains in the absence and presence of fluconazole and compared their fluconazole tolerance and resistance trajectories during evolution. The data show that mistranslation increases tolerance and accelerates the acquisition of resistance to fluconazole. Genome sequencing, array-based comparative genome analysis, and gene expression profiling revealed that during the course of evolution in fluconazole, the range of mutational and gene deregulation differences was distinctively different and broader in the hypermistranslating strain, including multiple chromosome duplications, partial chromosome deletions, and polyploidy. Especially, the increased accumulation of loss-of-heterozygosity events, aneuploidy, translational and cell surface modifications, and differences in drug efflux seem to mediate more rapid drug resistance acquisition under mistranslation. Our observations support a pivotal role for adaptive mistranslation in the evolution of drug resistance in C. albicans. IMPORTANCE Infectious diseases caused by drug-resistant fungi are an increasing threat to public health because of the high mortality rates and high costs associated with treatment. Thus, understanding of the molecular mechanisms of drug resistance is of crucial interest for the medical community. Here we investigated the role of regulated protein mistranslation, a characteristic mechanism used by C. albicans to diversify its proteome, in the evolution of fluconazole resistance. Such codon ambiguity is usually considered highly deleterious, yet recent studies found that mistranslation can boost adaptation in stressful environments. Our data reveal that CUG ambiguity diversifies the genome in multiple ways and that the full spectrum of drug resistance mechanisms in C. albicans goes beyond the traditional pathways that either regulate drug efflux or alter the interactions of drugs with their targets. The present work opens new avenues to understand the molecular and genetic basis of microbial drug resistance.

The Fungus Candida albicans Tolerates Ambiguity at Multiple Codons.

The ascomycete Candida albicans is a normal resident of the gastrointestinal tract of humans and other warm-blooded animals. It occurs in a broad range of body sites and has high capacity to survive and proliferate in adverse environments with drastic changes in oxygen, carbon dioxide, pH, osmolarity, nutrients, and temperature. Its biology is unique due to flexible reassignment of the leucine CUG codon to serine and synthesis of statistical proteins. Under standard growth conditions, CUG sites incorporate leucine (3% of the times) and serine (97% of the times) on a proteome wide scale, but leucine incorporation fluctuates in response to environmental stressors and can be artificially increased up to 98%. In order to determine whether such flexibility also exists at other codons, we have constructed several serine tRNAs that decode various non-cognate codons. Expression of these tRNAs had minor effects on fitness, but growth of the mistranslating strains at different temperatures, in medium with different pH and nutrients composition was often enhanced relatively to the wild type (WT) strain, supporting our previous data on adaptive roles of CUG ambiguity in variable growth conditions. Parallel evolution of the recombinant strains (100 generations) followed by full genome resequencing identified various strain specific single nucleotide polymorphisms (SNP) and one SNP in the deneddylase (JAB1) gene in all strains. Since JAB1 is a subunit of the COP9 signalosome complex, which interacts with cullin (Cdc53p) to mediate degradation of a variety of cellular proteins, our data suggest that neddylation plays a key role in tolerance and adaptation to codon ambiguity in C. albicans.

TRNA mutations that affect decoding fidelity deregulate development and the proteostasis network in zebrafish.

Mutations in genes that encode tRNAs, aminoacyl-tRNA syntheases, tRNA modifying enzymes and other tRNA interacting partners are associated with neuropathies, cancer, type-II diabetes and hearing loss, but how these mutations cause disease is unclear. We have hypothesized that levels of tRNA decoding error (mistranslation) that do not fully impair embryonic development can accelerate cell degeneration through proteome instability and saturation of the proteostasis network. To test this hypothesis we have induced mistranslation in zebrafish embryos using mutant tRNAs that misincorporate Serine (Ser) at various non-cognate codon sites. Embryo viability was affected and malformations were observed, but a significant proportion of embryos survived by activating the unfolded protein response (UPR), the ubiquitin proteasome pathway (UPP) and downregulating protein biosynthesis. Accumulation of reactive oxygen species (ROS), mitochondrial and nuclear DNA damage and disruption of the mitochondrial network, were also observed, suggesting that mistranslation had a strong negative impact on protein synthesis rate, ER and mitochondrial homeostasis. We postulate that mistranslation promotes gradual cellular degeneration and disease through protein aggregation, mitochondrial dysfunction and genome instability.

Reversion of a fungal genetic code alteration links proteome instability with genomic and phenotypic diversification.

Many fungi restructured their proteomes through incorporation of serine (Ser) at thousands of protein sites coded by the leucine (Leu) CUG codon. How these fungi survived this potentially lethal genetic code alteration and its relevance for their biology are not understood. Interestingly, the human pathogen Candida albicans maintains variable Ser and Leu incorporation levels at CUG sites, suggesting that this atypical codon assignment flexibility provided an effective mechanism to alter the genetic code. To test this hypothesis, we have engineered C. albicans strains to misincorporate increasing levels of Leu at protein CUG sites. Tolerance to the misincorporations was very high, and one strain accommodated the complete reversion of CUG identity from Ser back to Leu. Increasing levels of Leu misincorporation decreased growth rate, but production of phenotypic diversity on a phenotypic array probing various metabolic networks, drug resistance, and host immune cell responses was impressive. Genome resequencing revealed an increasing number of genotype changes at polymorphic sites compared with the control strain, and 80% of Leu misincorporation resulted in complete loss of heterozygosity in a large region of chromosome V. The data unveil unanticipated links between gene translational fidelity, proteome instability and variability, genome diversification, and adaptive phenotypic diversity. They also explain the high heterozygosity of the C. albicans genome and open the door to produce microorganisms with genetic code alterations for basic and applied research.

Molecular reconstruction of a fungal genetic code alteration.

Fungi of the CTG clade translate the Leu CUG codon as Ser. This genetic code alteration is the only eukaryotic sense-to-sense codon reassignment known to date, is mediated by an ambiguous serine tRNA (tRNACAG(Ser)), exposes unanticipated flexibility of the genetic code and raises major questions about its selection and fixation in this fungal lineage. In particular, the origin of the tRNACAG(Ser) and the evolutionary mechanism of CUG reassignment from Leu to Ser remain poorly understood. In this study, we have traced the origin of the tDNACAG(Ser) gene and studied critical mutations in the tRNACAG(Ser) anticodon-loop that modulated CUG reassignment. Our data show that the tRNACAG(Ser) emerged from insertion of an adenosine in the middle position of the 5'-CGA-3'anticodon of a tRNACGA(Ser) ancestor, producing the 5'-CAG-3' anticodon of the tRNACAG(Ser), without altering its aminoacylation properties. This mutation initiated CUG reassignment while two additional mutations in the anticodon-loop resolved a structural conflict produced by incorporation of the Leu 5'-CAG-3'anticodon in the anticodon-arm of a tRNA(Ser). Expression of the mutant tRNACAG(Ser) in yeast showed that it cannot be expressed at physiological levels and we postulate that such downregulation was essential to maintain Ser misincorporation at sub-lethal levels during the initial stages of CUG reassignment. We demonstrate here that such low level CUG ambiguity is advantageous in specific ecological niches and we propose that misreading tRNAs are targeted for degradation by an unidentified tRNA quality control pathway.

Dre-miR-2188 targets Nrp2a and mediates proper intersegmental vessel development in zebrafish embryos.

BACKGROUND: MicroRNAs (miRNAs) are a class of small RNAs that are implicated in the control of eukaryotic gene expression by binding to the 3'UTR of target mRNAs. Several algorithms have been developed for miRNA target prediction however, experimental validation is still essential for the correct identification of miRNA targets. We have recently predicted that Neuropilin2a (Nrp2a), a vascular endothelial growth factor receptor which is essential for normal developmental angiogenesis in zebrafish, is a dre-miR-2188 target. METHODOLOGY: Here we show that dre-miR-2188 targets the 3'-untranslated region (3'UTR) of Nrp2a mRNA and is implicated in proper intersegmental vessel development in vivo. Over expression of miR-2188 in zebrafish embryos down regulates Nrp2a expression and results in intersegmental vessel disruption, while its silencing increases Nrp2a expression and intersegmental vessel sprouting. An in vivo GFP sensor assay based on a fusion between the GFP coding region and the Nrp2a 3'UTR confirms that miR-2188 binds to the 3'UTR of Nrp2a and inhibits protein translation. CONCLUSIONS: We demonstrate that miR-2188 targets Nrp2a and affects intersegmental vessel development in zebrafish embryos.

Low level genome mistranslations deregulate the transcriptome and translatome and generate proteotoxic stress in yeast.

BACKGROUND: Organisms use highly accurate molecular processes to transcribe their genes and a variety of mRNA quality control and ribosome proofreading mechanisms to maintain intact the fidelity of genetic information flow. Despite this, low level gene translational errors induced by mutations and environmental factors cause neurodegeneration and premature death in mice and mitochondrial disorders in humans. Paradoxically, such errors can generate advantageous phenotypic diversity in fungi and bacteria through poorly understood molecular processes. RESULTS: In order to clarify the biological relevance of gene translational errors we have engineered codon misreading in yeast and used profiling of total and polysome-associated mRNAs, molecular and biochemical tools to characterize the recombinant cells. We demonstrate here that gene translational errors, which have negligible impact on yeast growth rate down-regulate protein synthesis, activate the unfolded protein response and environmental stress response pathways, and down-regulate chaperones linked to ribosomes. CONCLUSIONS: We provide the first global view of transcriptional and post-transcriptional responses to global gene translational errors and we postulate that they cause gradual cell degeneration through synergistic effects of overloading protein quality control systems and deregulation of protein synthesis, but generate adaptive phenotypes in unicellular organisms through activation of stress cross-protection. We conclude that these genome wide gene translational infidelities can be degenerative or adaptive depending on cellular context and physiological condition.

Ethanol exposure induces upregulation of specific microRNAs in zebrafish embryos.

Prenatal exposure to ethanol leads to a myriad of developmental disorders known as fetal alcohol spectrum disorder, often characterized by growth and mental retardation, central nervous system damage, and specific craniofacial dysmorphic features. The mechanisms of ethanol toxicity are not fully understood, but exposure during development affects the expression of several genes involved in cell cycle control, apoptosis, and transcriptional regulation. MicroRNAs (miRNAs) are implicated in some of these processes, however, it is not yet clear if they are involved in ethanol-induced toxicity. In order to clarify this question, we have exposed zebrafish embryos to ethanol and evaluated whether a miRNA deregulation signature could be obtained. Zebrafish embryos were exposed to 1 and 1.5% of ethanol from 4 h postfertilization (hpf) to 24 hpf. The miRNA expression profiles obtained reveal significant miRNA deregulation and show that both ethanol concentrations upregulate miR-153a, miR-725, miR-30d, let-7k, miR-100, miR-738, and miR-732. Putative gene targets of deregulated miRNAs are involved in cell cycle control, apoptosis, and transcription, which are the main processes affected by ethanol toxicity. The conservation of affected mechanisms among vertebrates leads us to postulate that similar miRNA deregulation occurs in humans, highlighting a relevant role of miRNAs in ethanol toxicology.

Species-specific codon context rules unveil non-neutrality effects of synonymous mutations.

BACKGROUND: Codon pair usage (codon context) is a species specific gene primary structure feature whose evolutionary and functional roles are poorly understood. The data available show that codon-context has direct impact on both translation accuracy and efficiency, but one does not yet understand how it affects these two translation variables or whether context biases shape gene evolution. METHODOLOGIES/PRINCIPAL FINDINGS: Here we study codon-context biases using a set of 72 orthologous highly conserved genes from bacteria, archaea, fungi and high eukaryotes to identify 7 distinct groups of codon context rules. We show that synonymous mutations, i.e., neutral mutations that occur in synonymous codons of codon-pairs, are selected to maintain context biases and that non-synonymous mutations, i.e., non-neutral mutations that alter protein amino acid sequences, are also under selective pressure to preserve codon-context biases. CONCLUSIONS: Since in vivo studies provide evidence for a role of codon context on decoding fidelity in E. coli and for decoding efficiency in mammalian cells, our data support the hypothesis that, like codon usage, codon context modulates the evolution of gene primary structure and fine tunes the structure of open reading frames for high genome translational fidelity and efficiency in the 3 domains of life.

The genetic code of the fungal CTG clade.

Genetic code alterations discovered over the last 40 years in bacteria and eukaryotes invalidate the hypothesis that the code is universal and frozen. Mitochondria of various yeast species translate the UGA stop codon as tryptophan (Trp) and leucine (Leu) CUN codons (N = any nucleotide) as threonine (Thr) and fungal CTG clade species reassigned Leu CUG codons to serine and translate them ambiguously in their cytoplasms. This unique sense-to-sense genetic code alteration is mediated by a Ser-tRNA containing a Leu 5'-CAG-3'anticodon (ser-tRNA(CAG)), which is recognized and charged with Ser (~97%) by the seryl-tRNA synthetase (SerRS) and with Leu (~3%) by the leucyl-tRNA synthetase (LeuRS). This unusual tRNA appeared 272 ± 25 million years ago and had a profound impact on the evolution of the CTG clade species. Here, we review the most recent results and concepts arising from the study of this codon reassignment and we highlight how its study is changing our views of the evolution of the genetic code.

Expression variability of co-regulated genes differentiates Saccharomyces cerevisiae strains.

BACKGROUND: Saccharomyces cerevisiae (Baker's yeast) is found in diverse ecological niches and is characterized by high adaptive potential under challenging environments. In spite of recent advances on the study of yeast genome diversity, little is known about the underlying gene expression plasticity. In order to shed new light onto this biological question, we have compared transcriptome profiles of five environmental isolates, clinical and laboratorial strains at different time points of fermentation in synthetic must medium, during exponential and stationary growth phases. RESULTS: Our data unveiled diversity in both intensity and timing of gene expression. Genes involved in glucose metabolism and in the stress response elicited during fermentation were among the most variable. This gene expression diversity increased at the onset of stationary phase (diauxic shift). Environmental isolates showed lower average transcript abundance of genes involved in the stress response, assimilation of nitrogen and vitamins, and sulphur metabolism, than other strains. Nitrogen metabolism genes showed significant variation in expression among the environmental isolates. CONCLUSIONS: Wild type yeast strains respond differentially to the stress imposed by nutrient depletion, ethanol accumulation and cell density increase, during fermentation of glucose in synthetic must medium. Our results support previous data showing that gene expression variability is a source of phenotypic diversity among closely related organisms.

Development of the genetic code: insights from a fungal codon reassignment.

The high conservation of the genetic code and its fundamental role in genome decoding suggest that its evolution is highly restricted or even frozen. However, various prokaryotic and eukaryotic genetic code alterations, several alternative tRNA-dependent amino acid biosynthesis pathways, regulation of tRNA decoding by diverse nucleoside modifications and recent in vivo incorporation of non-natural amino acids into prokaryotic and eukaryotic proteins, show that the code evolves and is surprisingly flexible. The cellular mechanisms and the proteome buffering capacity that support such evolutionary processes remain unclear. Here we explore the hypothesis that codon misreading and reassignment played fundamental roles in the development of the genetic code and we show how a fungal codon reassignment is enlightening its evolution.

Genetic code ambiguity: an unexpected source of proteome innovation and phenotypic diversity.

Translation of the genome into the proteome is a highly accurate biological process. However, the molecular mechanisms involved in protein synthesis are not error free and downstream protein quality control systems are needed to counteract the negative effects of translational errors (mistranslation) on proteome and cell homeostasis. This plus human and mice diseases caused by translational error generalized the idea that codon ambiguity is detrimental to life. Here we depart from this classical view of deleterious translational error and highlight how codon ambiguity can play important roles in the evolution of novel proteins. We also explain how tRNA mischarging can be relevant for the synthesis of functional proteomes, how codon ambiguity generates phenotypic and genetic diversity and how advantageous phenotypes can be selected, fixed, and inherited. A brief introduction to the molecular nature of translational error is provided; however, detailed information on the mechanistic aspects of mistranslation or comprehensive literature reviews of this topic should be obtained elsewhere.

The yeast PNC1 longevity gene is up-regulated by mRNA mistranslation.

Translation fidelity is critical for protein synthesis and to ensure correct cell functioning. Mutations in the protein synthesis machinery or environmental factors that increase synthesis of mistranslated proteins result in cell death and degeneration and are associated with neurodegenerative diseases, cancer and with an increasing number of mitochondrial disorders. Remarkably, mRNA mistranslation plays critical roles in the evolution of the genetic code, can be beneficial under stress conditions in yeast and in Escherichia coli and is an important source of peptides for MHC class I complex in dendritic cells. Despite this, its biology has been overlooked over the years due to technical difficulties in its detection and quantification. In order to shed new light on the biological relevance of mistranslation we have generated codon misreading in Saccharomyces cerevisiae using drugs and tRNA engineering methodologies. Surprisingly, such mistranslation up-regulated the longevity gene PNC1. Similar results were also obtained in cells grown in the presence of amino acid analogues that promote protein misfolding. The overall data showed that PNC1 is a biomarker of mRNA mistranslation and protein misfolding and that PNC1-GFP fusions can be used to monitor these two important biological phenomena in vivo in an easy manner, thus opening new avenues to understand their biological relevance.

Parallel DNA pyrosequencing unveils new zebrafish microRNAs.

BACKGROUND: MicroRNAs (miRNAs) are a new class of small RNAs of approximately 22 nucleotides in length that control eukaryotic gene expression by fine tuning mRNA translation. They regulate a wide variety of biological processes, namely developmental timing, cell differentiation, cell proliferation, immune response and infection. For this reason, their identification is essential to understand eukaryotic biology. Their small size, low abundance and high instability complicated early identification, however cloning/Sanger sequencing and new generation genome sequencing approaches overcame most technical hurdles and are being used for rapid miRNA identification in many eukaryotes. RESULTS: We have applied 454 DNA pyrosequencing technology to miRNA discovery in zebrafish (Danio rerio). For this, a series of cDNA libraries were prepared from miRNAs isolated at different embryonic time points and from fully developed organs. Each cDNA library was tagged with specific sequences and was sequenced using the Roche FLX genome sequencer. This approach retrieved 90% of the 192 miRNAs previously identified by cloning/Sanger sequencing and bioinformatics. Twenty five novel miRNAs were predicted, 107 miRNA star sequences and also 41 candidate miRNA targets were identified. A miRNA expression profile built on the basis of pyrosequencing read numbers showed high expression of most miRNAs throughout zebrafish development and identified tissue specific miRNAs. CONCLUSION: This study increases the number of zebrafish miRNAs from 192 to 217 and demonstrates that a single DNA mini-chip pyrosequencing run is effective in miRNA identification in zebrafish. This methodology also produced sufficient information to elucidate miRNA expression patterns during development and in differentiated organs. Moreover, some zebrafish miRNA star sequences were more abundant than their corresponding miRNAs, suggesting a functional role for the former in gene expression control in this vertebrate model organism.