miRDeep*: an integrated application tool for miRNA identification from RNA sequencing data

miRDeep and its varieties are widely used to quantify known and novel micro RNA (miRNA) from small RNA sequencing (RNAseq). This article describes miRDeep*, our integrated miRNA identification tool, which is modeled off miRDeep, but the precision of detecting novel miRNAs is improved by introducing new strategies to identify precursor miRNAs. miRDeep* has a user-friendly graphic interface and accepts raw data in FastQ and Sequence Alignment Map (SAM) or the binary equivalent (BAM) format. Known and novel miRNA expression levels, as measured by the number of reads, are displayed in an interface, which shows each RNAseq read relative to the pre-miRNA hairpin. The secondary pre-miRNA structure and read locations for each predicted miRNA are shown and kept in a separate figure file. Moreover, the target genes of known and novel miRNAs are predicted using the TargetScan algorithm, and the targets are ranked according to the confidence score. miRDeep* is an integrated standalone application where sequence alignment, pre-miRNA secondary structure calculation and graphical display are purely Java coded. This application tool can be executed using a normal personal computer with 1.5 GB of memory. Further, we show that miRDeep* outperformed existing miRNA prediction tools using our LNCaP and other small RNAseq datasets. miRDeep* is freely available online at http://www.australianprostatecentre.org/research/software/mirdeep-star

miRPlant : an integrated tool for identification of plant miRNA from RNA sequencing data

Background Small RNA sequencing is commonly used to identify novel miRNAs and to determine their expression levels in plants. There are several miRNA identification tools for animals such as miRDeep, miRDeep2 and miRDeep*. miRDeep-P was developed to identify plant miRNA using miRDeep’s probabilistic model of miRNA biogenesis, but it depends on several third party tools and lacks a user-friendly interface. The objective of our miRPlant program is to predict novel plant miRNA, while providing a user-friendly interface with improved accuracy of prediction. Result We have developed a user-friendly plant miRNA prediction tool called miRPlant. We show using 16 plant miRNA datasets from four different plant species that miRPlant has at least a 10% improvement in accuracy compared to miRDeep-P, which is the most popular plant miRNA prediction tool. Furthermore, miRPlant uses a Graphical User Interface for data input and output, and identified miRNA are shown with all RNAseq reads in a hairpin diagram. Conclusions We have developed miRPlant which extends miRDeep* to various plant species by adopting suitable strategies to identify hairpin excision regions and hairpin structure filtering for plants. miRPlant does not require any third party tools such as mapping or RNA secondary structure prediction tools. miRPlant is also the first plant miRNA prediction tool that dynamically plots miRNA hairpin structure with small reads for identified novel miRNAs. This feature will enable biologists to visualize novel pre-miRNA structure and the location of small RNA reads relative to the hairpin. Moreover, miRPlant can be easily used by biologists with limited bioinformatics skills.

Mapping the regulatory sequences controlling 93 breast cancer-associated miRNA genes leads to the identification of two functional promoters of the Hsa-mir-200b cluster, methylation of which is associated with metastasis or hormone receptor status in advanced breast cancer

MicroRNAs (miRNAs) are small non-coding RNAs of 20 nt in length that are capable of modulating gene expression post-transcriptionally. Although miRNAs have been implicated in cancer, including breast cancer, the regulation of miRNA transcription and the role of defects in this process in cancer is not well understood. In this study we have mapped the promoters of 93 breast cancer-associated miRNAs, and then looked for associations between DNA methylation of 15 of these promoters and miRNA expression in breast cancer cells. The miRNA promoters with clearest association between DNA methylation and expression included a previously described and a novel promoter of the Hsa-mir-200b cluster. The novel promoter of the Hsa-mir-200b cluster, denoted P2, is located 2 kb upstream of the 5′ stemloop and maps within a CpG island. P2 has comparable promoter activity to the previously reported promoter (P1), and is able to drive the expression of miR-200b in its endogenous genomic context. DNA methylation of both P1 and P2 was inversely associated with miR-200b expression in eight out of nine breast cancer cell lines, and in vitro methylation of both promoters repressed their activity in reporter assays. In clinical samples, P1 and P2 were differentially methylated with methylation inversely associated with miR-200b expression. P1 was hypermethylated in metastatic lymph nodes compared with matched primary breast tumours whereas P2 hypermethylation was associated with loss of either oestrogen receptor or progesterone receptor. Hypomethylation of P2 was associated with gain of HER2 and androgen receptor expression. These data suggest an association between miR-200b regulation and breast cancer subtype and a potential use of DNA methylation of miRNA promoters as a component of a suite of breast cancer biomarkers.

Study of miRNA polymorphisms and their potential association with breast cancer risk in Australian Caucasian populations

This project established a large and well characterised prospective breast cancer DNA biobank and used this biobank to conduct genetic studies in breast cancer. The thesis presented the results of these high-throughput genotyping studies in two separate Australian Caucasian case-control populations and identified association between three novel genetic variants in microRNA genes and breast cancer risk.

A large-scale analysis of genetic variants within putative miRNA binding sites in prostate cancer

Prostate cancer is the second most common malignancy among men worldwide. Genome-wide association studies have identified 100 risk variants for prostate cancer, which can explain approximately 33% of the familial risk of the disease. We hypothesized that a comprehensive analysis of genetic variations found within the 3' untranslated region of genes predicted to affect miRNA binding (miRSNP) can identify additional prostate cancer risk variants. We investigated the association between 2,169 miRSNPs and prostate cancer risk in a large-scale analysis of 22,301 cases and 22,320 controls of European ancestry from 23 participating studies. Twenty-two miRSNPs were associated (P<2.3×10(-5)) with risk of prostate cancer, 10 of which were within 7 genes previously not mapped by GWAS studies. Further, using miRNA mimics and reporter gene assays, we showed that miR-3162-5p has specific affinity for the KLK3 rs1058205 miRSNP T-allele, whereas miR-370 has greater affinity for the VAMP8 rs1010 miRSNP A-allele, validating their functional role. SIGNIFICANCE Findings from this large association study suggest that a focus on miRSNPs, including functional evaluation, can identify candidate risk loci below currently accepted statistical levels of genome-wide significance. Studies of miRNAs and their interactions with SNPs could provide further insights into the mechanisms of prostate cancer risk.

Genetic association analysis of miRNA SNPs implicates MIR145 in breast cancer susceptibility

Background MicroRNAs (miRNAs) are important small non-coding RNA molecules that regulate gene expression in cellular processes related to the pathogenesis of cancer. Genetic variation in miRNA genes could impact their synthesis and cellular effects and single nucleotide polymorphisms (SNPs) are one example of genetic variants studied in relation to breast cancer. Studies aimed at identifying miRNA SNPs (miR-SNPs) associated with breast malignancies could lead towards further understanding of the disease and to develop clinical applications for early diagnosis and treatment. Methods We genotyped a panel of 24 miR-SNPs using multiplex PCR and chip-based matrix assisted laser desorption ionization time-of-flight (MALDI-TOF) mass spectrometry (MS) analysis in two Caucasian breast cancer case control populations (Primary population: 173 cases and 187 controls and secondary population: 679 cases and 301 controls). Association to breast cancer susceptibility was determined using chi-square (X 2 ) and odds ratio (OR) analysis. Results Statistical analysis showed six miR-SNPs to be non-polymorphic and twelve of our selected miR-SNPs to have no association with breast cancer risk. However, we were able to show association between rs353291 (located in MIR145) and the risk of developing breast cancer in two independent case control cohorts (p = 0.041 and p = 0.023). Conclusions Our study is the first to report an association between a miR-SNP in MIR145 and breast cancer risk in individuals of Caucasian background. This finding requires further validation through genotyping of larger cohorts or in individuals of different ethnicities to determine the potential significance of this finding as well as studies aimed to determine functional significance. Keywords: Association analysis; Breast cancer; microRNA; miR-SNPs; MIR145

Melanoma cell invasiveness is regulated by miR-211 suppression of the BRN2 transcription factor

To identify microRNAs potentially involved in melanomagenesis, we compared microRNA expression profiles between melanoma cell lines and cultured melanocytes. The most differentially expressed microRNA between the normal and tumor cell lines was miR-211. We focused on this pigment-cell-enriched miRNA as it is derived from the microphthalmia-associated transcription factor (MITF)-regulated gene, TRPM1 (melastatin). We find that miR-211 expression is greatly decreased in melanoma cells and melanoblasts compared to melanocytes. Bioinformatic analysis identified a large number of potential targets of miR-211, including POU3F2 (BRN2). Inhibition of miR-211 in normal melanocytes resulted in increased BRN2 protein, indicating that endogenous miR-211 represses BRN2 in differentiated cells. Over-expression of miR-211 in melanoma cell lines changed the invasive potential of the cells in vitro through directly targeting BRN2 translation. We propose a model for the apparent non-overlapping expression levels of BRN2 and MITF in melanoma, mediated by miR-211 expression.

Wolbachia uses host microRNAs to manipulate host gene expression and facilitate colonization of the dengue vector Aedes aegypti

The obligate endosymbiont Wolbachia pipientis is found in a wide range of invertebrates where they are best known for manipulating host reproduction. Recent studies have shown that Wolbachia also can modulate the lifespan of host insects and interfere with the development of human pathogens in mosquito vectors. Despite considerable study, very little is known about the molecular interactions between Wolbachia and its hosts that might mediate these effects. Using microarrays, we show that the microRNA (miRNA) profile of the mosquito, Aedes aegypti, is significantly altered by the wMelPop-CLA strain of W. pipientis. We found that a host miRNA (aae-miR-2940) is induced after Wolbachia infection in both mosquitoes and cell lines. One target of aae-miR-2940 is the Ae. aegypti metalloprotease gene. Interestingly, expression of the target gene was induced after Wolbachia infection, ectopic expression of the miRNA independent of Wolbachia, or transfection of an artificial mimic of the miRNA into mosquito cells. We also confirmed the interaction of aae-miR-2940 with the target sequences using GFP as a reporter gene. Silencing of the metalloprotease gene in both Wolbachia-infected cells and adult mosquitoes led to a significant reduction in Wolbachia density, as did inhibition of the miRNA in cells. These results indicate that manipulation of the mosquito metalloprotease gene via aae-miR-2940 is crucial for efficient maintenance of the endosymbiont. This report shows how Wolbachia alters the host miRNA profile and provides insight into the mechanisms of host manipulation used by this widespread endosymbiont.

Circulating microRNAs involved in multiple sclerosis

Multiple sclerosis (MS) is an immune-mediated, demyelinating and neurodegenerative disease of the central nervous system. After traumatic brain injury, it is the leading cause of neurology disability in young adults. Considerable advances have been made in identifying genes involved in MS but the genetic and phenotypic complexity associated with this disease significantly hinders any progress. A novel class of small RNA molecules, microRNAs (miRNAs) has acquired much attention because they regulate the expression of up to 30% of protein-coding genes and may play a pivotal role in the development of many, if not all, complex diseases. Seven published studies investigated miRNAs from peripheral blood mononuclear cells, CD4+, CD8+ T cell, B lymphocytes, peripheral blood leukocytes, whole blood and brain astrocytes with MS risk. The absence of MS studies investigating plasma miRNA prompted the current investigation of identifying a circulating miRNA signature in MS. We conducted a microarray analysis of over 900 known miRNA transcripts from plasma samples collected from four MS individuals and four sex-aged and ethnicity matched healthy controls. We identified six plasma miRNA (miR-614, miR-572, miR-648, miR-1826, miR-422a and miR-22) that were significantly up-regulated and one plasma miRNA (miR-1979) that was significantly down-regulated in MS individuals. Both miR-422a and miR-22 have previously been implicated in MS. The present study is the first to show a circulating miRNA signature involved in MS that could serve as a potential prognostic and diagnostic biomarker for MS.

A genetic variant located in miR-423 is associated with reduced breast cancer risk

Background/Aim: Since microRNAs (miRNAs) act as translational regulators of multiple genes, single nucleotide polymorphisms (SNP) in them can have potentially wide-ranging effects. Using an association approach, this research examined the effects of the rs6505162 SNP, an A>C polymorphism located in the premiRNA region of miR-423, on breast cancer development. Materials and Methods: Caucasian Australian women with breast cancer and controls matched for age and ethnicity were genotyped for rs6505162 and their genotypic and allelic frequencies analysed for significant differences. Results: Analysis indicated that there were significant differences between the case and control populations (χ 2=6.70, p=0.035), with the CC genotype conferring reduced risk of breast cancer development (odds ratio=0.50 95% confidence interval=0.27-0.92, p=0.03). Conclusion: Further functional research is required to determine the mechanism of action of this SNP on miRNA function.

The association of single nucleotide polymorphisms with endometrial cancer risk and prognosis

Endometrial cancer is one of the most common female diseases in developed nations and is the most commonly diagnosed gynaecological cancer in Australia. The disease is commonly classified by histology: endometrioid or non-endometrioid endometrial cancer. While non-endometrioid endometrial cancers are accepted to be high-grade, aggressive cancers, endometrioid cancers (comprising 80% of all endometrial cancers diagnosed) generally carry a favourable patient prognosis. However, endometrioid endometrial cancer patients endure significant morbidity due to surgery and radiotherapy used for disease treatment, and patients with recurrent disease have a 5-year survival rate of less than 50%. Genetic analysis of women with endometrial cancer could uncover novel markers associated with disease risk and/or prognosis, which could then be used to identify women at high risk and for the use of specialised treatments. Proteases are widely accepted to play an important role in the development and progression of cancer. This PhD project hypothesised that SNPs from two protease gene families, the matrix metalloproteases (MMPs, including their tissue inhibitors, TIMPs) and the tissue kallikrein-related peptidases (KLKs) would be associated with endometrial cancer susceptibility and/or prognosis. In the first part of this study, optimisation of the genotyping techniques was performed. Results from previously published endometrial cancer genetic association studies were attempted to be validated in a large, multicentre replication set (maximum cases n = 2,888, controls n = 4,483, 3 studies). The rs11224561 progesterone receptor SNP (PGR, A/G) was observed to be associated with increased endometrial cancer risk (per A allele OR 1.31, 95% CI 1.12-1.53; p-trend = 0.001), a result which was initially reported among a Chinese sample set. Previously reported associations for the remaining 8 SNPs investigated for this section of the PhD study were not confirmed, thereby reinforcing the importance of validation of genetic association studies. To examine the effect of SNPs from the MMP and KLK families on endometrial cancer risk, we selected the most significantly associated MMP and KLK SNPs from genome-wide association study analysis (GWAS) to be genotyped in the GWAS replication set (cases n = 4,725, controls n = 9,803, 13 studies). The significance of the MMP24 rs932562 SNP was unchanged after incorporation of the stage 2 samples (Stage 1 per allele OR 1.18, p = 0.002; Combined Stage 1 and 2 OR 1.09, p = 0.002). The rs10426 SNP, located 3' to KLK10 was predicted by bioinformatic analysis to effect miRNA binding. This SNP was observed in the GWAS stage 1 result to exhibit a recessive effect on endometrial cancer risk, a result which was not validated in the stage 2 sample set (Stage 1 OR 1.44, p = 0.007; Combined Stage 1 and 2 OR 1.14, p = 0.08). Investigation of the regions imputed surrounding the MMP, TIMP and KLK genes did not reveal any significant targets for further analysis. Analysis of the case data from the endometrial cancer GWAS to identify genetic variation associated with cancer grade did not reveal SNPs from the MMP, TIMP or KLK genes to be statistically significant. However, the representation of SNPs from the MMP, TIMP and KLK families by the GWAS genotyping platform used in this PhD project was examined and observed to be very low, with the genetic variation of four genes (MMP23A, MMP23B, MMP28 and TIMP1) not captured at all by this technique. This suggests that comprehensive candidate gene association studies will be required to assess the role of SNPs from these genes with endometrial cancer risk and prognosis. Meta-analysis of gene expression microarray datasets curated as part of this PhD study identified a number of MMP, TIMP and KLK genes to display differential expression by endometrial cancer status (MMP2, MMP10, MMP11, MMP13, MMP19, MMP25 and KLK1) and histology (MMP2, MMP11, MMP12, MMP26, MMP28, TIMP2, TIMP3, KLK6, KLK7, KLK11 and KLK12). In light of these findings these genes should be prioritised for future targeted genetic association studies. Two SNPs located 43.5 Mb apart on chromosome 15 were observed from the GWAS analysis to be associated with increased endometrial cancer grade, results that were validated in silico in two independent datasets. One of these SNPs, rs8035725 is located in the 5' untranslated region of a MYC promoter binding protein DENND4A (Stage 1 OR 1.15, p = 9.85 x 10P -5 P, combined Stage 1 and in silico validation OR 1.13, p = 5.24 x 10P -6 P). This SNP has previously been reported to alter the expression of PTPLAD1, a gene involved in the synthesis of very long fatty acid chains and in the Rac1 signaling pathway. Meta-analysis of gene expression microarray data found PTPLAD1 to display increased expression in the aggressive non-endometrioid histology compared with endometrioid endometrial cancer, suggesting that the causal SNP underlying the observed genetic association may influence expression of this gene. Neither rs8035725 nor significant SNPs identified by imputation were predicted bioinformatically to affect transcription factor binding sites, indicating that further studies are required to assess their potential effect on other regulatory elements. The other grade- associated SNP, rs6606792, is located upstream of an inferred pseudogene, ELMO2P1 (Stage 1 OR 1.12, p = 5 x 10P -5 P; combined Stage 1 and in silico validation OR 1.09, p = 3.56 x 10P -5 P). Imputation of the ±1 Mb region surrounding this SNP revealed a cluster of significantly associated variants which are predicted to abolish various transcription factor binding sites, and would be expected to decrease gene expression. ELMO2P1 was not included on the microarray platforms collected for this PhD, and so its expression could not be investigated. However, the high sequence homology of ELMO2P1 with ELMO2, a gene important to cell motility, indicates that ELMO2 could be the parent gene for ELMO2P1 and as such, ELMO2P1 could function to regulate the expression of ELMO2. Increased expression of ELMO2 was seen to be associated with increasing endometrial cancer grade, as well as with aggressive endometrial cancer histological subtypes by microarray meta-analysis. Thus, it is hypothesised that SNPs in linkage disequilibrium with rs6606792 decrease the transcription of ELMO2P1, reducing the regulatory effect of ELMO2P1 on ELMO2 expression. Consequently, ELMO2 expression is increased, cell motility is enhanced leading to an aggressive endometrial cancer phenotype. In summary, these findings have identified several areas of research for further study. The results presented in this thesis provide evidence that a SNP in PGR is associated with risk of developing endometrial cancer. This PhD study also reports two independent loci on chromosome 15 to be associated with increased endometrial cancer grade, and furthermore, genes associated with these SNPs to be differentially expressed according in aggressive subtypes and/or by grade. The studies reported in this thesis support the need for comprehensive SNP association studies on prioritised MMP, TIMP and KLK genes in large sample sets. Until these studies are performed, the role of MMP, TIMP and KLK genetic variation remains unclear. Overall, this PhD study has contributed to the understanding of genetic variation involvement in endometrial cancer susceptibility and prognosis. Importantly, the genetic regions highlighted in this study could lead to the identification of novel gene targets to better understand the biology of endometrial cancer and also aid in the development of therapeutics directed at treating this disease.

The role of EBV-encoded microRNAs in B cell differentiation

Epstein Barr virus (EBV) is a common γ-herpes virus, infecting approximately 90% of the world‟s population. It is also one of the first known viruses known to be oncogenic, and is associated with a number of tumour types, primarily lymphomas. MicroRNAs (miRNAs) are post-transcriptional regulators of gene expression and many human miRNAs have been associated with the development of malignancies including cancer. EBV was the first human virus identified to express miRNAs and encodes more than 40 miRNAs within its genome. Yet, an understanding of the targets of EBV-miRNAs, and thereby the function of them in pathogenesis remains sadly limited. This study identifies a potential novel target of EBV-miRNAs, MECP2 and characterises the miRNA:mRNA interactions between two previously identified novel targets; Bim and EBF1. In particular, this study focuses upon the interaction between EBF1 and the EBV-miRNA BART11-5p, demonstrating a 151bp region of the EBF1 3‟UTR that is capable of mediating the silencing of luciferase expression by BART11-5p but is not capable of silencing a full length EBF1-3‟UTR luciferase construct. This study provides evidence that EBF1 may be a target of one or more EBV-miRNAs.

DRB2 is required for microRNA biogenesis in Arabidopsis thaliana

Background The Arabidopsis thaliana (Arabidopsis) DOUBLE-STRANDED RNA BINDING (DRB) protein family consists of five members, DRB1 to DRB5. The biogenesis of two developmentally important small RNA (sRNA) species, the microRNAs (miRNAs) and trans-acting small interfering RNAs (tasiRNAs) by DICER-LIKE (DCL) endonucleases requires the assistance of DRB1 and DRB4 respectively. The importance of miRNA-directed target gene expression in plant development is exemplified by the phenotypic consequence of loss of DRB1 activity (drb1 plants). Principal Findings Here we report that the developmental phenotype of the drb235 triple mutant plant is the result of deregulated miRNA biogenesis in the shoot apical meristem (SAM) region. The expression of DRB2, DRB3 and DRB5 in wild-type seedlings is restricted to the SAM region. Small RNA sequencing of the corresponding tissue of drb235 plants revealed altered miRNA accumulation. Approximately half of the miRNAs detected remained at levels equivalent to those of wild-type plants. However, the accumulation of the remaining miRNAs was either elevated or reduced in the triple mutant. Examination of different single and multiple drb mutants revealed a clear association between the loss of DRB2 activity and altered accumulation for both the elevated and reduced miRNA classes. Furthermore, we show that the constitutive over-expression of DRB2 outside of its wild-type expression domain can compensate for the loss of DRB1 activity in drb1 plants. Conclusions/Significance Our results suggest that in the SAM region, DRB2 is both antagonistic and synergistic to the role of DRB1 in miRNA biogenesis, adding an additional layer of gene regulatory complexity in this developmentally important tissue.

The Arabidopsis thaliana double-stranded RNA binding protein DRB1 directs guide strand selection from microRNA duplexes

In Arabidopsis thaliana (Arabidopsis), DICER-LIKE1 (DCL1) functions together with the double-stranded RNA binding protein (dsRBP), DRB1, to process microRNAs (miRNAs) from their precursor transcripts prior to their transfer to the RNA-induced silencing complex (RISC). miRNA-loaded RISC directs RNA silencing of cognate mRNAs via ARGONAUTE1 (AGO1)-catalyzed cleavage. Short interefering RNAs (siRNAs) are processed from viral-derived or transgene-encoded molecules of doublestranded RNA (dsRNA) by the DCL/dsRBP partnership, DCL4/DRB4, and are also loaded to AGO1-catalyzed RISC for cleavage of complementary mRNAs. Here, we use an artificial miRNA (amiRNA) technology, transiently expressed in Nicotiana benthamiana, to produce a series of amiRNA duplexes with differing intermolecular thermostabilities at the 5′ end of duplex strands. Analyses of amiRNA duplex strand accumulation and target transcript expression revealed that strand selection (amiRNA and amiRNA*) is directed by asymmetric thermostability of the duplex termini. The duplex strand possessing a lower 59 thermostability was preferentially retained by RISC to guide mRNA cleavage of the corresponding target transgene. In addition, analysis of endogenous miRNA duplex strand accumulation in Arabidopsis drb1 and drb2345 mutant plants revealed that DRB1 dictates strand selection, presumably by directional loading of the miRNA duplex onto RISC for passenger strand degradation. Bioinformatic and Northern blot analyses of DCL4/DRB4-dependent small RNAs (miRNAs and siRNAs) revealed that small RNAs produced by this DCL/dsRBP combination do not conform to the same terminal thermostability rules as those governing DCL1/DRB1-processed miRNAs. This suggests that small RNA processing in the DCL1/DRB1-directed miRNA and DCL4/DRB4-directed sRNA biogenesis pathways operates via different mechanisms.

DRB2, DRB3 and DRB5 function in a non-canonical microRNA pathway in Arabidopsis thaliana

DOUBLE-STRANDED RNA BIN DIN G (DRB) proteins have been functionally characterized in viruses, prokaryotes and eukaryotes and are involved in all aspects of RNA biology. Arabidopsis thaliana (Arabidopsis) encodes five closely related DRB proteins, DRB1 to DRB5. DRB1 and DRB4 are required by DICER-LIKE (DCL) proteins DCL1 and DCL4 to accurately and efficiently process structurally distinct double-stranded RNA (dsRNA) precursor substrates in the microRNA (miRNA) and trans-acting small-interfering RNA (tasiRNA) biogenesis pathways respectively. We recently reported that DRB2 is also involved in the biogenesis of specific miRNA subsets. Furthermore, the severity of the developmental phenotype displayed by the drb235 triple mutant plant, compared with those expressed by either drb2, drb3 and drb5 single mutants, or double mutant combinations thereof, indicates that DRB3 and DRB5 function in the same non-canonical miRNA pathway as DRB2. Through the use of our artificial miRNA (amiRNA) plant expression vector, pBlueGreen 2,3 we demonstrate here that unlike DRB2, DRB3 and DRB5 are not involved in the dsRNA processing stages of the miRNA biogenesis pathway, but are required to mediate RNA silencing of target genes of DRB2-associated miRNA s. © 2012 Landes Bioscience.