Putative Target Sites Research Articles

MiR-370 regulates ISG15 expression and influences IFN-α sensitivity in hepatocellular carcinoma cells.

BACKGROUND: Interferon- (IFN-) is an adjuvant to chemotherapy and radiotherapy for hepatocellular carcinoma (HCC), but some HCC patients do not respond to treatment with IFN-.METHODS: We performed loss-of-function and gain-of-function experiments to examine the role of ISG15 in the IFN- sensitivity of LH86, HLCZ01, SMMC7721, and Huh7 cell lines and tumor samples.RESULTS: The overexpression of ISG15 reduced apoptosis in Huh7 and LH86 cells in the presence of IFN-, whereas the shRNA-mediated knock down of ISG15 expression increased apoptosis in both Huh7 and LH86 cells. We identified a putative miR-370 target site in the 3’-UTR in the ISG15 mRNA, and the level of miR-370 expression in HCC cell lines reflected the level of IFN--induced apoptosis exhibited by each. Both HCC cell lines and tumor samples had significantly lower levels of miR-370 than the control cells and tissues ( 0.05). The overexpression of miR-370 in IFN--treated LH86 and Huh7 cells increased apoptosis and reduced the volume of LH86- and Huh7-derived xenograft tumors in mice treated with IFN- compared with the control tumors.CONCLUSIONS: Our findings suggest that miR-370 functions as an HCC tumor suppressor and regulator of IFN- sensitivity and that miR-370 might be a useful prognostic marker for HCC patients.

Microcystin-LR-Triggered Neuronal Toxicity in Whitefish Does Not Involve MiR124-3p

Microcystin-LR (MC-LR) is a potent hepatotoxin that has also been pointed out of causing neurotoxicity, but the exact mechanisms of action still remain ambiguous and need to be elucidated. Data from studies on mammals show that pathology of astrocyte cells points to perturbations of microRNA signaling. Glial fibrillary acidic protein (GFAP), a neuronal cell/astrocyte-specific protein, and a microRNA-124-3p (MiR124-3p) are among putative triggers and regulators of neuronal cell/astrocyte reactivity. In the present study on whitefish (Coregonus lavaretus), we found that gfap mRNA contains a putative target site for MIR124-3p, to potentially affect its expression changes. qPCR expression study of gfap:MiR124-3p pair in the midbrain of juvenile whitefish, during 28 days of exposure to a repeated subacute dose of MC-LR (100 μg kg−1 body mass), showed marginally significant up-regulation of gfap only on the 7th day of exposure period which suggests neuronal toxicity. During the whole exposure period, neither midbrain nor blood plasma levels of MiR124-3p were changed. Furthermore, double luciferase gene reporter assay confirmed the lack of MiR124-3p involvement in mediating control over gfap mRNA expression. These data show that, although MC-LR may trigger neuronal toxicity in whitefish, this does not involve MiR124-3p in response to the treatment.

Genome-wide analysis of the SPL/miR156 module and its interaction with the AP2/miR172 unit in barley

The SQUAMOSA-promoter binding like (SPL) gene family encodes transcription factors that have been shown in many species to influence plant growth and development, but information about these genes in barley (Hordeum vulgare L.) is limited. This study identified 17 barley SPL genes, within eight distinct groups, that are orthologs of SPL genes described in Arabidopsis, wheat, and rice. Sixteen barley SPLs undergo alternative splicing. Seven SPLs contain a putative miR156 target site and the transcript levels of the miR156-targeted HvSPLs (HvSPL3, 13 and 23) were lower in vegetative than in reproductive phase but this was true also for some SPL genes such as HvSPL6 that were not regulated by miR156. Because SPL gene products regulate miR172, which is also involved in floral development, the expression of miR172 was studied. An antagonistic expression pattern of miR156 and miR172b during the vegetative and the reproductive phases signifies their apparent function in barley growth phase transition. Characterization of a barley mir172 mutant having an abnormal, indeterminate spikelet phenotype suggests the possible feedback role of AP2/miR172 module on HvSPL genes. This is the first comprehensive analysis of the miR156/SPL/miR172 axis in barley that provides a basis to elucidate their roles in various biological processes.

Factors targeting MED12 to drive tumorigenesis?

Mediator Complex Subunit 12 (MED12) is part of the transcriptional preinitiation machinery. Mutations of its gene predominantly occur in two types of highly frequent benign tumors, uterine leiomyomas and fibroadenomas of the breast, where they apparently act as driver mutations. Nevertheless, their presence is not restricted to benign tumors having been found at considerable frequencies in uterine leiomyosarcomas, malignant phyllodes tumors, and chronic lymphocytic leukemia also. Most of the mutations are located within exon 2 of the gene but in rare cases the intron 1/exon 2 boundary or exon 1 are affected. As to their type, predominantly single nucleotide exchanges with a hotspot in one codon are found, but small deletions clustering around that hotspot also are not uncommon. These latter deletions are leaving the open reading frame intact. As to the types of mutations, so far no apparent differences between the tumor entities affected have emerged. Interestingly, this pattern with small deletions clustered around the hotspot of single nucleotide exchanges resembles that seen as a result of targeted gene editing. In contrast to other driver mutations the percentage of MED12-mutation positive tumors of independent clonal origin increases with the number of tumors per patient suggesting unknown etiological factors supporting site specific mutagenesis. These factors may act by inducing simultaneous site-specific double strand breaks the erroneous repair of which may lead to corresponding mutations. As inducers of DNA damage and its repair such as foreign nucleic acids of the microbiome displaying sequence homology to the putative target site might play a role. Interestingly, a 16 base pair homology of the hotspot to a putative terminator base-paired hairpin sequence of a Staphylococcus aureus tRNA gene cluster has been noted which might form R-loop like structures with its target sequence thus inducing said changes.

Use of chemical probes to explore the toxicological potential of the K+/Cl− cotransporter (KCC) as a novel insecticide target to control the primary vector of dengue and Zika virus, Aedes aegypti

The majority of commercialized insecticides target the insect nervous system and therefore, neural proteins are well-validated targets for insecticide development. Considering that only a few neural targets are exploited for insecticidal action and the development of insecticide resistance has reduced the efficacy of current insecticidal classes, we sought to test the toxicological potential of the potassium-chloride cotransporter (KCC). In mammals, KCC proteins have seminal roles in shaping GABAergic signaling and inhibitory neurotransmission, thus ion transport through KCC is critical for proper neurotransmission. Therefore, we hypothesized that mosquito KCC represents a putative insecticide target site and that pharmacological inhibition of KCC constructs in Aedes aegypti will be lethal. To test this hypothesis, we developed a robust, cell-based fluorescence assay that enables in vitro characterization of small-molecules against Ae. aegypti KCC and performed a proof-of-concept study employing well characterized mammalian KCC modulators to determine the toxicological potential of Ae. aegypti KCC. The selective inhibitor of mammalian KCC, termed VU0463271, was found to be a potent inhibitor Ae. aegypti KCC and microinjection induced lethality in a concentration-dependent manner to susceptible and pyrethroid resistant strains. Importantly, an analog of VU0463271 was shown to be >40-fold less potent and did not induce toxicity, suggesting that the observed physiological effects and mortality are likely due to KCC inhibition. This proof-of-concept study suggests that Ae. aegypti KCC represents a putative target site for mosquitocide design that can mitigate the current mechanisms of insecticide resistance.

Factors targeting MED12 to drive tumorigenesis?

Mediator Subcomplex 12 (MED12) is part of the transcriptional preinitiation machinery. Mutations of its gene predominantly occur in two types of highly frequent benign tumors, uterine leiomyomas and fibroadenomas of the breast, where they apparently act as driver mutations. Nevertheless, their presence is not restricted to benign tumors having been found at considerable frequencies in uterine leiomyosarcomas, malignant phyllodes tumors, and chronic lymphocytic leukemia also. Most of the mutations are located within exon 2 of the gene but in rare cases the intron 1/exon 2 boundary or exon 1 are affected. As to their type, predominantly single nucleotide exchanges with a hotspot in one codon are found, but small deletions clustering around that hotspot also are not uncommon. According to their presumed classification as gain-of-function mutations, these latter deletions are leaving the open reading frame intact. As to the types of mutations, so far no apparent differences between the tumor entities affected have emerged. Interestingly, this pattern with small deletions clustered around the hotspot of single nucleotide exchanges resembles that seen as a result of targeted gene editing. In contrast to other driver mutations the percentage of MED12-mutation positive tumors of independent clonal origin increases with the number of tumors per patient suggesting unknown etiological factors supporting site specific mutagenesis. These factors may act by inducing simultaneous site-specific double strand breaks the erroneous repair of which may lead to corresponding mutations. As inducers of DNA damage and its repair such as foreign nucleic acids of the microbiome displaying sequence homology to the putative target site might play a role. Interestingly, a 16 base pair homology of the hotspot to a putative terminator base-paired hairpin sequence of a Staphylococcus aureus tRNA gene cluster has been noted which might form R-loop like structures with its target sequence thus inducing said changes.

SNP rs403212791 in exon 2 of the MTNR1A gene is associated with reproductive seasonality in the Rasa aragonesa sheep breed

The aim of this study was to characterize and identify causative SNPs in the MTNR1A gene responsible for the reproductive seasonality traits in the Rasa aragonesa sheep breed. A total of 290 ewes (155, 84 and 51 mature, young and ewe lambs, respectively) from one flock were controlled from January to August. The following three reproductive seasonality traits were considered: the total days of anoestrus (TDA) and the progesterone cycling months (P4CM); both ovarian function seasonality traits based on blood progesterone levels; and the oestrus cycling months (OCM) based on oestrous detection, which indicate behavioural signs of oestrous. We have sequenced the total coding region plus 733 and 251 bp from the promoter and 3′-UTR regions, respectively, from the gene in 268 ewes. We found 9 and 4 SNPs associated with seasonality traits in the promoter (for TDA and P4CM) and exon 2 (for the three traits), respectively. The SNPs located in the gene promoter modify the putative binding sites for various trans-acting factors. In exon 2, two synonymous SNPs affect RFLP sites, rs406779174/RsaI (for the three traits) and rs430181568/MnlI (for OCM), and they have been related with seasonal reproductive activity in previous association studies with other breeds. SNP rs400830807, which is located in the 3′-UTR, was associated with the three traits, but this did not modify the putative target sites for ovine miRNAs according to in silico predictions. Finally, the SNP rs403212791 (NW_014639035.1: g.15099004G > A), which is also associated with the three seasonality phenotypes, was the most significant SNP detected in this study and was a non-synonymous polymorphism, leading a change from an Arginine to a Cysteine (R336C). Haplotype analyses confirmed the association results and showed that the effects found for the seasonality traits were caused by the SNPs located in exon 2. We have demonstrated that the T allele in the SNP rs403212791 in the MNTR1A gene is associated with a lower TDA and higher P4CM and OCM values in the Rasa Aragonesa breed.

Differentially expressed microRNAs between cattleyak and yak testis

Cattleyak are interspecific hybrids between cattle and yak, exhibiting the same prominent adaptability as yak and much higher performances than yak. However, male infertility of cattleyak resulted from spermatogenic arrest has greatly restricted their effective utilization in yak breeding. In past decades, much work has been done to investigate the mechanisms of spermatogenic arrest, but little is known about the differences of the post-transcriptional regulators between cattleyak and yak, which may contribute to the impaired spermatogenesis. MiRNAs, a class of endogenous non-coding small RNA, were revealed to play crucial roles in regulating gene expression at post-transcriptional level. In the present study, we identified 50 differentially expressed (DE) known miRNAs and 11 novel miRNAs by using Illumina HISeq and bioinformatic analysis. A total of 50 putative target sites for the 13 DE known miRNAs and 30 for the 6 DE novel miRNAs were identified, respectively. GO and KEGG analyses were performed to reveal the functions of target genes for DE miRNAs. In addition, RT-qPCR was performed to validate the expression of the DE miRNAs and its targets. The identification of these miRNAs may provide valuable information for a better understanding of spermatogenic arrest in cattleyak.

The effects of alternative splicing on miRNA binding sites in bladder cancer.

Eukaryotic organisms have developed a variety of mechanisms to regulate translation post-transcriptionally, including but not limited to the use of miRNA silencing in many species. One method of post-transcriptional regulation is through miRNAs that bind to the 3′ UTRs to regulate mRNA abundance and influence protein expression. Therefore, the diversity of mRNA 3′ UTRs mediating miRNA binding sites influence miRNA-mediated regulation. Alternative polyadenylation, by shortening mRNA isoforms, increases the diversity of 3′ UTRs; moreover, short mRNA isoforms elude miRNA-medicated repression. Because no current prediction methods for putative miRNA target sites consider whether or not 1) splicing-informed miRNA binding sites and/or 2) the use of 3′ UTRs provide higher resolution or functionality, we sought to identify not only the genome-wide impact of using exons in mRNA 3′ UTRs but also their functional connection to miRNA regulation and clinical outcomes in cancer. With a genome-wide expression of mRNA and miRNA quantified by 395 bladder cancer cases from The Cancer Genome Atlas (TCGA), we 1) demonstrate the diversity of 3′ UTRs affecting miRNA efficiency and 2) identify a set of genes clinically associated with mRNA expression in bladder cancer. Knowledge of 3′ UTR diversity will not only be a useful addition to current miRNA target prediction algorithms but also enhance the clinical utility of mRNA isoforms in the expression of mRNA in cancer. Thus, variability among cancer patient’s variability in molecular signatures based on these exon usage events in 3′ UTR along with miRNAs in bladder cancer may lead to better prognostic/treatment strategies for improved precision medicine.

Post-transcriptional Regulation of PCSK9 by miR-191, miR-222, and miR-224.

Since proprotein convertase subtilisin kexin 9 (PCSK9) discovery, a gene involved in LDL metabolism regulation and cardiovascular diseases (CVD), many therapeutic strategies have been introduced for direct targeting of PCSK9. The main goal of these strategies has been to reduce PCSK9 protein level either by application of antibodies or inhibition of its production. In this study, we have tried to discover microRNAs (miRNAs) which can target, and hence regulate, PCSK9 expression. Using bioinformatics tools, we selected three microRNAs with binding sites on 3′-UTR of PCSK9. The expression level of these miRNAs was examined in three different cell lines using real-time RT-PCR. We observed a reciprocal expression pattern between expression level of miR-191, miR-222, and miR-224 with that of PCSK9. Accordingly, the expression levels were highest in Huh7 cells which expressed the lowest level of PCSK9, compared to HepG2 and A549 cell lines. PCSK9 mRNA level also showed a significant decline in HepG2 cells transfected with the vectors overexpressing the aforementioned miRNAs. Furthermore, the miRNAs target sites were cloned in psiCHECK-2 vector, and a direct interaction of the miRNAs and the PCSK9 3′-UTR putative target sites was investigated by means of luciferase assay. Our findings revealed that miR-191, miR-222, and miR-224 can directly interact with PCSK9 3′-UTR and regulate its expression. In conclusion, our data introduces a role for miRNAs to regulate PCSK9 expression.

MicroRNA-30a suppresses non-small-cell lung cancer by targeting Myb-related protein B.

Non-small-cell lung cancer (NSCLC) is one of the leading causes of cancer mortality worldwide. A growing body of evidence indicates that microRNA (miR) have important and diverse roles in the proliferation, apoptosis and metastasis of human cancer cells. In the present study, the molecular regulation mechanism of miR-30a and its potential target, Myb-related protein B (MYBL2) was investigated in NSCLC. Reverse transcription-quantitative polymerase chain reaction results showed that miR-30a was significantly downregulated in NSCLC tissues compared with adjacent normal tissues (P<0.05). MYBL2 has a putative miR-30a target site in its 3'untranslated region according to previous data, prediction databases and TargetScan software. In the present study, a negative correlation was demonstrated between miR-30a and MYBL2 expression in NSCLC. Direct interaction between miR-30a and MYBL2 was also confirmed via a dual-luciferase reporter assay. miR-30a overexpression inhibited the growth of A549 and H460 cells via MTT and bromodeoxyuridine incorporation assays, whereas miR-30a downregulation promoted cell proliferation. In addition, miR-30a overexpression not only increased cell apoptosis and induced cell cycle arrest in A549 and H460 cell lines, but also attenuated tumor growth, and mRNA and protein expression levels of MYBL2. The present findings suggest that miR-30a may suppress NSCLC by targeting MYBL2.

Long Noncoding RNA MEG3 Suppresses Glioma Cell Proliferation, Migration, and Invasion by Acting as a Competing Endogenous RNA of miR-19a.

Glioma, with varying malignancy grades and histological subtypes, is the most common primary brain tumor in adults. Long noncoding RNAs (lncRNAs) are non-protein-coding transcripts and have been proven to play an important role in tumorigenesis. Our study aims to elucidate the combined effect of lncRNA maternally expressed gene 3 (MEG3) and microRNA-19a (miR-19a) in human glioma U87 and U251 cell lines. Real-time PCR revealed that MEG3 was downregulated and miR-19a was upregulated in malignant glioma tissues and cell lines. Bioinformatics analyses (TargetScan, miRanda, and starBase V2.0) showed that phosphatase and tensin homolog (PTEN) is a target of miR-19a with complementary binding sites in the 3′-UTR. As expected, luciferase results verified the putative target site and also revealed the complementary binding between miR-19a and MEG3. miR-19a represses the expression of PTEN and promotes glioma cell proliferation, migration, and invasion. However, MEG3 could directly bind to miR-19a and effectively act as a competing endogenous RNA (ceRNA) for miR-19a to suppress tumorigenesis. Our study is the first to demonstrate that lncRNA MEG3 suppresses glioma cell proliferation, migration, and invasion by acting as a ceRNA of miR-19a, which provides a novel insight about the pathogenesis of glioma.

Aurora Kinase B, a novel regulator of TERF1 binding and telomeric integrity.

AURKB (Aurora Kinase B) is a serine/threonine kinase better known for its role at the mitotic kinetochore during chromosome segregation. Here, we demonstrate that AURKB localizes to the telomeres in mouse embryonic stem cells, where it interacts with the essential telomere protein TERF1. Loss of AURKB function affects TERF1 telomere binding and results in aberrant telomere structure. In vitro kinase experiments successfully identified Serine 404 on TERF1 as a putative AURKB target site. Importantly, in vivo overexpression of S404-TERF1 mutants results in fragile telomere formation. These findings demonstrate that AURKB is an important regulator of telomere structural integrity.

Conjecture Regarding Posttranslational Modifications to the Arabidopsis Type I Proton-Pumping Pyrophosphatase (AVP1).

Agbiotechnology uses genetic engineering to improve the output and value of crops. Altering the expression of the plant Type I Proton-pumping Pyrophosphatase (H+-PPase) has already proven to be a useful tool to enhance crop productivity. Despite the effective use of this gene in translational research, information regarding the intracellular localization and functional plasticity of the pump remain largely enigmatic. Using computer modeling several putative phosphorylation, ubiquitination and sumoylation target sites were identified that may regulate Arabidopsis H+-PPase (AVP1- Arabidopsis Vacuolar Proton-pump 1) subcellular trafficking and activity. These putative regulatory sites will direct future research that specifically addresses the partitioning and transport characteristics of this pump. We posit that fine-tuning H+-PPases activity and cellular distribution will facilitate rationale strategies for further genetic improvements in crop productivity.

The highly pathogenic H5N1 influenza A virus down-regulated several cellular MicroRNAs which target viral genome.

Higher and prolonged viral replication is critical for the increased pathogenesis of the highly pathogenic avian influenza (HPAI) subtype of H5N1 influenza A virus (IAV) over the lowly pathogenic H1N1 IAV strain. Recent studies highlighted the considerable roles of cellular miRNAs in host defence against viral infection. In this report, using a 3′UTR reporter system, we identified several putative miRNA target sites buried in the H5N1 virus genome. We found two miRNAs, miR‐584‐5p and miR‐1249, that matched with the PB2 binding sequence. Moreover, we showed that these miRNAs dramatically down‐regulated PB2 expression, and inhibited replication of H5N1 and H1N1 IAVs in A549 cells. Intriguingly, these miRNAs expression was differently regulated in A549 cells infected with the H5N1 and H1N1 viruses. Furthermore, transfection of miR‐1249 inhibitor enhanced the PB2 expression and promoted the replication of H5N1 and H1N1 IAVs. These results suggest that H5N1 virus may have evolved a mechanism to escape host‐mediated inhibition of viral replication through down‐regulation of cellular miRNAs, which target its viral genome.

Transcriptome-wide identification of host genes targeted by tomato spotted wilt virus-derived small interfering RNAs

RNA silencing mechanism functions as a major defense against invading viruses. The caveat in the RNA silencing mechanism is that the effector small interfering RNAs (siRNAs) act on any RNA transcripts with sequence complementarity irrespective of target’s origin. A subset of highly expressed viral small interfering RNAs (vsiRNAs) derived from the tomato spotted wilt virus (TSWV; Tospovirus: Bunyaviridae) genome was analyzed for their propensity to downregulate the tomato transcriptome. A total of 11898 putative target sites on tomato transcripts were found to exhibit a propensity for down regulation by TSWV-derived vsiRNAs. In total, 2450 unique vsiRNAs were found to have potential cross-reacting capability with the tomato transcriptome. VsiRNAs were found to potentially target a gamut of host genes involved in basal cellular activities including enzymes, transcription factors, membrane transporters, and cytoskeletal proteins. KEGG pathway annotation of targets revealed that the vsiRNAs were mapped to secondary metabolite biosynthesis, amino acids, starch and sucrose metabolism, and carbon and purine metabolism. Transcripts for protein processing, hormone signalling, and plant-pathogen interactions were the most likely targets from the genetic, environmental information processing, and organismal systems, respectively. qRT-PCR validation of target gene expression showed that none of the selected transcripts from tomato cv. Marglobe showed up regulation, and all were down regulated even upto 20 folds (high affinity glucose transporter). However, the expression levels of transcripts from cv. Red Defender revealed differential regulation as three among the target transcripts showed up regulation (Cc-nbs-lrr, resistance protein, AP2-like ethylene-responsive transcription factor, and heat stress transcription factor A3). Accumulation of tomato target mRNAs of corresponding length was proved in both tomato cultivars using 5′ RACE analysis. The TSWV-tomato interaction at the sRNA interface points to the ability of tomato cultivars to overcome vsiRNA-mediated targeting of NBS-LRR class R genes. These results suggest the prevalence of vsiRNA-induced RNA silencing of host transcriptome, and the interactome scenario is the first report on the interaction between tospovirus genome-derived siRNAs and tomato transcripts, and provide a deeper understanding of the role of vsiRNAs in pathogenicity and in perturbing host machinery.

Overexpression of miRNA-221 promotes cell proliferation by targeting the apoptotic protease activating factor-1 and indicates a poor prognosis in ovarian cancer.

MicroRNAs are a class of small non-coding, endogenous RNAs involved in cancer development and progression. MicroRNA-221 (mir-221) has been reported to have both an oncogenic and tumor-suppressive role in human tumors, but the role of miR-221 in ovarian cancer is poorly understood. In the present study, the expression levels of miR-221 and the apoptosis protease activating factor 1 (APAF1) protein in 63 samples of ovarian cancer tissues and the cell lines, IOSE25, A2780, OVCAR3, SKOV3 and 3AO were detected by quantitative reverse-transcription polymerase chain reaction (qRT-PCR) and western blot analysis, respectively. Cell proliferation was measured using Cell Counting kit-8 (CCK-8); cell migration and invasion were detected using a Transwell assay; cell apoptosis was evaluated by flow cytometry and Hoechst staining, and a luciferase assay was performed to verify a putative target site of miR-221 in the 3′-UTR of APAF1 mRNA. Expression of miR-221 was upregulated in ovarian cancer tissues. Patients with increased miR-221 expression levels had a reduced disease-free survival (P=0.0014) and overall survival (P=0.0058) compared with those with low miR-221 expression. Transfection of SKOV3 and A2780 cell lines with miR-221 inhibitor induced APAF1 protein expression, suppressed cell proliferation and migration and promoted tumor cell apoptosis. In conclusion, the APAF1 gene was confirmed as a direct target of miR-221 and overexpression of APAF1 suppressed ovarian cancer cell proliferation and induced cell apoptosis in vitro. These findings indicate that miR-221-APAF1 should be studied further as a potential new diagnostic or prognostic biomarker for ovarian cancer.

Transcriptome Analysis of Leaves, Flowers and Fruits Perisperm of Coffea arabica L. Reveals the Differential Expression of Genes Involved in Raffinose Biosynthesis.

Coffea arabica L. is an important crop in several developing countries. Despite its economic importance, minimal transcriptome data are available for fruit tissues, especially during fruit development where several compounds related to coffee quality are produced. To understand the molecular aspects related to coffee fruit and grain development, we report a large-scale transcriptome analysis of leaf, flower and perisperm fruit tissue development. Illumina sequencing yielded 41,881,572 high-quality filtered reads. De novo assembly generated 65,364 unigenes with an average length of 1,264 bp. A total of 24,548 unigenes were annotated as protein coding genes, including 12,560 full-length sequences. In the annotation process, we identified nine candidate genes related to the biosynthesis of raffinose family oligossacarides (RFOs). These sugars confer osmoprotection and are accumulated during initial fruit development. Four genes from this pathway had their transcriptional pattern validated by quantitative reverse transcription polymerase chain reaction (qRT-PCR). Furthermore, we identified ~24,000 putative target sites for microRNAs (miRNAs) and 134 putative transcriptionally active transposable elements (TE) sequences in our dataset. This C. arabica transcriptomic atlas provides an important step for identifying candidate genes related to several coffee metabolic pathways, especially those related to fruit chemical composition and therefore beverage quality. Our results are the starting point for enhancing our knowledge about the coffee genes that are transcribed during the flowering and initial fruit development stages.

Prediction of human miRNA target genes using computationally reconstructed ancestral mammalian sequences.

MicroRNAs (miRNA) are short single-stranded RNA molecules derived from hairpin-forming precursors that play a crucial role as post-transcriptional regulators in eukaryotes and viruses. In the past years, many microRNA target genes (MTGs) have been identified experimentally. However, because of the high costs of experimental approaches, target genes databases remain incomplete. Although several target prediction programs have been developed in the recent years to identify MTGs in silico, their specificity and sensitivity remain low. Here, we propose a new approach called MirAncesTar, which uses ancestral genome reconstruction to boost the accuracy of existing MTGs prediction tools for human miRNAs. For each miRNA and each putative human target UTR, our algorithm makes uses of existing prediction tools to identify putative target sites in the human UTR, as well as in its mammalian orthologs and inferred ancestral sequences. It then evaluates evidence in support of selective pressure to maintain target site counts (rather than sequences), accounting for the possibility of target site turnover. It finally integrates this measure with several simpler ones using a logistic regression predictor. MirAncesTar improves the accuracy of existing MTG predictors by 26% to 157%. Source code and prediction results for human miRNAs, as well as supporting evolutionary data are available at http://cs.mcgill.ca/∼blanchem/mirancestar.

Genome-Wide Identification, Characterization and Expression Analysis of the TCP Gene Family in Prunus mume.

TCP proteins, belonging to a plant-specific transcription factors family, are known to have great functions in plant development, especially flower and leaf development. However, there is little information about this gene family in Prunus mume, which is widely cultivated in China as an ornamental and fruit tree. Here a genome-wide analysis of TCP genes was performed to explore their evolution in P. mume. Nineteen PmTCPs were identified and three of them contained putative miR319 target sites. Phylogenetic and comprehensive bioinformatics analyses of these genes revealed that different types of TCP genes had undergone different evolutionary processes and the genes in the same clade had similar chromosomal location, gene structure, and conserved domains. Expression analysis of these PmTCPs indicated that there were diverse expression patterns among different clades. Most TCP genes were predominantly expressed in flower, leaf, and stem, and showed high expression levels in the different stages of flower bud differentiation, especially in petal formation stage and gametophyte development. Genes in TCP-P subfamily had main roles in both flower development and gametophyte development. The CIN genes in double petal cultivars might have key roles in the formation of petal, while they were correlated with gametophyte development in the single petal cultivar. The CYC/TB1 type genes were highly detected in the formation of petal and pistil. The less-complex flower types of P. mume might result from the fact that there were only two CYC type genes present in P. mume and a lack of CYC2 genes to control the identity of flower types. These results lay the foundation for further study on the functions of TCP genes during flower development.