Visualizing RAD51-mediated joint molecules: implications for recombination mechanism and the effect of sequence heterology

The human Shu complex promotes RAD51 activity by modulating RPA dynamics on ssDNA

Templated DNA repair that occurs during homologous recombination and replication stress relies on RAD51. RAD51 activity is positively regulated by BRCA2 and the RAD51 paralogs. The Shu complex is a RAD51 paralog-containing complex consisting of SWSAP1 and SWS1. We demonstrate that SWSAP1-SWS1 binds RAD51, maintains RAD51 filament stability, and enables strand exchange. Using single molecule confocal fluorescence microscopy combined with optical tweezers, we show that SWSAP1-SWS1 decorates RAD51 filaments proficient for homologous recombination. We also find SWSAP1-SWS1 enhances RPA diffusion on ssDNA. Importantly, we show human sgSWSAP1 and sgSWS1 knockout cells are sensitive to pharmacological inhibition of PARP and APE1. Lastly, we identify cancer variants in SWSAP1 that alter SWS1 complex formation. Together, we show that SWSAP1-SWS1 stimulates RAD51-dependent high-fidelity repair and may be an important new cancer therapeutic target.

Modeling the Homologous Recombination Process: Methods, Successes and Challenges

Homologous recombination (HR) is a fundamental process common to all species. HR aims to faithfully repair DNA double strand breaks. HR involves the formation of nucleoprotein filaments on DNA single strands (ssDNA) resected from the break. The nucleoprotein filaments search for homologous regions in the genome and promote strand exchange with the ssDNA homologous region in an unbroken copy of the genome. HR has been the object of intensive studies for decades. Because multi-scale dynamics is a fundamental aspect of this process, studying HR is highly challenging, both experimentally and using computational approaches. Nevertheless, knowledge has built up over the years and has recently progressed at an accelerated pace, borne by increasingly focused investigations using new techniques such as single molecule approaches. Linking this knowledge to the atomic structure of the nucleoprotein filament systems and the succession of unstable, transient intermediate steps that takes place during the HR process remains a challenge; modeling retains a very strong role in bridging the gap between structures that are stable enough to be observed and in exploring transition paths between these structures. However, working on ever-changing long filament systems submitted to kinetic processes is full of pitfalls. This review presents the modeling tools that are used in such studies, their possibilities and limitations, and reviews the advances in the knowledge of the HR process that have been obtained through modeling. Notably, we will emphasize how cooperative behavior in the HR nucleoprotein filament enables modeling to produce reliable information.

MicroRNA biogenesis and activity in plant cell dedifferentiation stimulated by cell wall removal

BACKGROUND

Despite the frequent use of protoplast-to-plant system in in vitro cultures of plants, the molecular mechanisms regulating the first and most limiting stages of this process, i.e., protoplast dedifferentiation and the first divisions leading to the formation of a microcallus, have not been elucidated.

RESULTS

In this study, we investigated the function of miRNAs in the dedifferentiation of A. thaliana mesophyll cells in a process stimulated by the enzymatic removal of the cell wall. Leaf cells, protoplasts and CDPs (cells derived from protoplasts) cultured for 24, 72 and 120 h (first cell division). In protoplasts, a strong decrease in the amount of AGO1 in both the nucleus and the cytoplasm, as well as dicing bodies (DBs), which are considered to be sites of miRNA biogenesis, was shown. However during CDPs division, the amounts of AGO1 and DBs strongly increased. MicroRNA transcriptome studies demonstrated that lower amount of differentially expressed miRNAs are present in protoplasts than in CDPs cultured for 120 h. Then analysis of differentially expressed miRNAs, selected pri-miRNA and mRNA targets were performed.

CONCLUSION

This result indicates that miRNA function is not a major regulation of gene expression in the initial but in later steps of dedifferentiation during CDPs divisions. miRNAs participate in organogenesis, oxidative stress, nutrient deficiencies and cell cycle regulation in protoplasts and CDPs. The important role played by miRNAs in the process of dedifferentiation of mesophyll cells was confirmed by the increased mortality and reduced cell division of CDPs derived from mutants with defective miRNA biogenesis and miR319b expression.

Mismatch sensing by nucleofilament deciphers mechanism of RecA-mediated homologous recombination

Recombinases polymerize along single-stranded DNA (ssDNA) at the end of a broken DNA to form a helical nucleofilament with a periodicity of ∼18 bases. The filament catalyzes the search and checking for homologous sequences and promotes strand exchange with a donor duplex during homologous recombination (HR), the mechanism of which has remained mysterious since its discovery. Here, by inserting mismatched segments into donor duplexes and using single-molecule techniques to catch transient intermediates in HR, we found that, even though 3 base pairs (bp) is still the basic unit, both the homology checking and the strand exchange may proceed in multiple steps at a time, resulting in ∼9-bp large steps on average. More interestingly, the strand exchange is blocked remotely by the mismatched segment, terminating at positions ∼9 bp before the match-mismatch joint. The homology checking and the strand exchange are thus separated in space, with the strand exchange lagging behind. Our data suggest that the strand exchange progresses like a traveling wave in which the donor DNA is incorporated successively into the ssDNA-RecA filament to check homology in ∼9-bp steps in the frontier, followed by a hypothetical transitional segment and then the post-strand-exchanged duplex.

In vitro role of Rad54 in Rad51-ssDNA filament-dependent homology search and synaptic complexes formation

Homologous recombination (HR) uses a homologous template to accurately repair DNA double-strand breaks and stalled replication forks to maintain genome stability. During homology search, Rad51 nucleoprotein filaments probe and interact with dsDNA, forming the synaptic complex that is stabilized on a homologous sequence. Strand intertwining leads to the formation of a displacement-loop (D-loop). In yeast, Rad54 is essential for HR in vivo and required for D-loop formation in vitro, but its exact role remains to be fully elucidated. Using electron microscopy to visualize the DNA-protein complexes, here we find that Rad54 is crucial for Rad51-mediated synaptic complex formation and homology search. The Rad54−K341R ATPase-deficient mutant protein promotes formation of synaptic complexes but not D-loops and leads to the accumulation of stable heterologous associations, suggesting that the Rad54 ATPase is involved in preventing non-productive intermediates. We propose that Rad51/Rad54 form a functional unit operating in homology search, synaptic complex and D-loop formation.

Homologous recombination uses a template to accurately repair DNA double-strand breaks and stalled replication forks to maintain genome stability. Here authors use electron microscopy to investigate the role of Rad54 in homology search and synaptic complex formation.

A stretched conformation of DNA with a biological role?

We have discovered a well-defined extended conformation of double-stranded DNA, which we call Σ-DNA, using laser-tweezers force-spectroscopy experiments. At a transition force corresponding to free energy change ΔG = 1·57 ± 0·12 kcal (mol base pair)-1 60 or 122 base-pair long synthetic GC-rich sequences, when pulled by the 3'-3' strands, undergo a sharp transition to the 1·52 ± 0·04 times longer Σ-DNA. Intriguingly, the same degree of extension is also found in DNA complexes with recombinase proteins, such as bacterial RecA and eukaryotic Rad51. Despite vital importance to all biological organisms for survival, genome maintenance and evolution, the recombination reaction is not yet understood at atomic level. We here propose that the structural distortion represented by Σ-DNA, which is thus physically inherent to the nucleic acid, is related to how recombination proteins mediate recognition of sequence homology and execute strand exchange. Our hypothesis is that a homogeneously stretched DNA undergoes a 'disproportionation' into an inhomogeneous Σ-form consisting of triplets of locally B-like perpendicularly stacked bases. This structure may ensure improved fidelity of base-pair recognition and promote rejection in case of mismatch during homologous recombination reaction. Because a triplet is the length of a gene codon, we speculate that the structural physics of nucleic acids may have biased the evolution of recombinase proteins to exploit triplet base stacks and also the genetic code.

Homologous recombination and the repair of DNA double-strand breaks

Homologous recombination enables the cell to access and copy intact DNA sequence information in trans, particularly to repair DNA damage affecting both strands of the double helix. Here, we discuss the DNA transactions and enzymatic activities required for this elegantly orchestrated process in the context of the repair of DNA double-strand breaks in somatic cells. This includes homology search, DNA strand invasion, repair DNA synthesis, and restoration of intact chromosomes. Aspects of DNA topology affecting individual steps are highlighted. Overall, recombination is a dynamic pathway with multiple metastable and reversible intermediates designed to achieve DNA repair with high fidelity.

Akt1 Stimulates Homologous Recombination Repair of DNA Double-Strand Breaks in a Rad51-Dependent Manner

Akt1 is known to promote non-homologous end-joining (NHEJ)-mediated DNA double-strand break (DSB) repair by stimulation of DNA-PKcs. In the present study, we investigated the effect of Akt1 on homologous recombination (HR)-dependent repair of radiation-induced DSBs in non-small cell lung cancer (NSCLC) cells A549 and H460. Akt1-knockdown (Akt1-KD) significantly reduced Rad51 protein level, Rad51 foci formation and its colocalization with γH2AX foci after irradiation. Moreover, Akt1-KD decreased clonogenicity after treatment with Mitomycin C and HR repair, as tested by an HR-reporter assay. Double knockdown of Akt1 and Rad51 did not lead to a further decrease in HR compared to the single knockdown of Rad51. Consequently, Akt1-KD significantly increased the number of residual DSBs after irradiation partially independent of the kinase activity of DNA-PKcs. Likewise, the number of residual BRCA1 foci, indicating unsuccessful HR events, also significantly increased in the irradiated cells after Akt1-KD. Together, the results of the study indicate that Akt1 seems to be a regulatory component in the HR repair of DSBs in a Rad51-dependent manner. Thus, based on this novel role of Akt1 in HR and the previously described role of Akt1 in NHEJ, we propose that targeting Akt1 could be an effective approach to selectively improve the killing of tumor cells by DSB-inducing cytotoxic agents, such as ionizing radiation.

Integrated mRNA and miRNA transcriptome reveal a cross-talk between developing response and hormone signaling for the seed kernels of Siberian apricot

Recently, our transcriptomic analysis has identified some functional genes responsible for oil biosynthesis in developing SASK, yet miRNA-mediated regulation for SASK development and oil accumulation is poorly understood. Here, 3 representative periods of 10, 30 and 60 DAF were selected for sRNA sequencing based on the dynamic patterns of growth tendency and oil content of developing SASK. By miRNA transcriptomic analysis, we characterized 296 known and 44 novel miRNAs in developing SASK, among which 36 known and 6 novel miRNAs respond specifically to developing SASK. Importantly, we performed an integrated analysis of mRNA and miRNA transcriptome as well as qRT-PCR detection to identify some key miRNAs and their targets (miR156-SPL, miR160-ARF18, miR164-NAC1, miR171h-SCL6, miR172-AP2, miR395-AUX22B, miR530-P2C37, miR393h-TIR1/AFB2 and psi-miRn5-SnRK2A) potentially involved in developing response and hormone signaling of SASK. Our results provide new insights into the important regulatory function of cross-talk between development response and hormone signaling for SASK oil accumulation.

PmiRExAt: plant miRNA expression atlas database and web applications

High-throughput small RNA (sRNA) sequencing technology enables an entirely new perspective for plant microRNA (miRNA) research and has immense potential to unravel regulatory networks. Novel insights gained through data mining in publically available rich resource of sRNA data will help in designing biotechnology-based approaches for crop improvement to enhance plant yield and nutritional value. Bioinformatics resources enabling meta-analysis of miRNA expression across multiple plant species are still evolving. Here, we report PmiRExAt, a new online database resource that caters plant miRNA expression atlas. The web-based repository comprises of miRNA expression profile and query tool for 1859 wheat, 2330 rice and 283 maize miRNA. The database interface offers open and easy access to miRNA expression profile and helps in identifying tissue preferential, differential and constitutively expressing miRNAs. A feature enabling expression study of conserved miRNA across multiple species is also implemented. Custom expression analysis feature enables expression analysis of novel miRNA in total 117 datasets. New sRNA dataset can also be uploaded for analysing miRNA expression profiles for 73 plant species. PmiRExAt application program interface, a simple object access protocol web service allows other programmers to remotely invoke the methods written for doing programmatic search operations on PmiRExAt database.Database URL:http://pmirexat.nabi.res.in.

Design, simplified cloning, and in-silico analysis of multisite small interfering RNA-targeting cassettes

Multiple gene silencing is being required to target and tangle metabolic pathways in eukaryotes and researchers have to develop a subtle method for construction of RNA interference (RNAi) cassettes. Although, several vectors have been developed due to different screening and cloning strategies but still some potential limitations remain to be dissolved. Here, we worked out a simple cloning strategy to develop multisite small interfering RNA (siRNA) cassette from different genes by two cloning steps. In this method, effective siRNA sites in the target messenger RNAs (mRNAs) were determined using in silico analysis and consecutively arranged to reduce length of inverted repeats. Here, we used one-step (polymerase chain reaction) PCR by designed long primer sets covering the selected siRNA sites. Rapid screening, cost-effective and shorten procedure are advantages of this method compare to PCR classic cloning. Validity of constructs was confirmed by optimal centroid secondary structures with high stability in plants.

A novel Fanconi anaemia subtype associated with a dominant-negative mutation in RAD51

Fanconi anaemia (FA) is a hereditary disease featuring hypersensitivity to DNA cross-linker-induced chromosomal instability in association with developmental abnormalities, bone marrow failure and a strong predisposition to cancer. A total of 17 FA disease genes have been reported, all of which act in a recessive mode of inheritance. Here we report on a de novo g.41022153G>A; p.Ala293Thr (NM_002875) missense mutation in one allele of the homologous recombination DNA repair gene RAD51 in an FA-like patient. This heterozygous mutation causes a novel FA subtype, 'FA-R', which appears to be the first subtype of FA caused by a dominant-negative mutation. The patient, who features microcephaly and mental retardation, has reached adulthood without the typical bone marrow failure and paediatric cancers. Together with the recent reports on RAD51-associated congenital mirror movement disorders, our results point to an important role for RAD51-mediated homologous recombination in neurodevelopment, in addition to DNA repair and cancer susceptibility.

Mismatch repair and homeologous recombination

DNA mismatch repair influences the outcome of recombination events between diverging DNA sequences. Here we discuss how mismatch repair proteins are active in different homologous recombination subpathways and specific reaction steps, resulting in differential modulation of these recombination events, with a focus on the mechanism of heteroduplex rejection during the inhibition of recombination between slightly diverged (homeologous) DNA sequences.

Integrating multi-scale data on homologous recombination into a new recognition mechanism based on simulations of the RecA-ssDNA/dsDNA structure

RecA protein is the prototypical recombinase. Members of the recombinase family can accurately repair double strand breaks in DNA. They also provide crucial links between pairs of sister chromatids in eukaryotic meiosis. A very broad outline of how these proteins align homologous sequences and promote DNA strand exchange has long been known, as are the crystal structures of the RecA-DNA pre- and postsynaptic complexes; however, little is known about the homology searching conformations and the details of how DNA in bacterial genomes is rapidly searched until homologous alignment is achieved. By integrating a physical model of recognition to new modeling work based on docking exploration and molecular dynamics simulation, we present a detailed structure/function model of homology recognition that reconciles extremely quick searching with the efficient and stringent formation of stable strand exchange products and which is consistent with a vast body of previously unexplained experimental results.

The advance of tomato disease-related microRNAs

Tomato is a model plant for studying plant-pathogen interactions. As regulatory factors, microRNAs (miRNAs) have been widely identified and play crucial roles in tomato-pathogen interactions, including host defense and pathogen counter-defense. Here, the review summarizes the discoveries and highlights of miRNAs in tomato diseases. Roles of artificial miRNAs in disease resistance are further discussed. Hence, a better understanding of the contribution of miRNAs in tomato disease will shed light on strategies in enhancing tomato-pathogen resistance.

Identification of Cold-Responsive miRNAs and Their Target Genes in Nitrogen-Fixing Nodules of Soybean

As a warm climate species, soybean is highly sensitive to chilling temperatures. Exposure to chilling temperatures causes a significant reduction in the nitrogen fixation rate in soybean plants and subsequent yield loss. However, the molecular basis for the sensitivity of soybean to chilling is poorly understood. In this study, we identified cold-responsive miRNAs in nitrogen-fixing nodules of soybean. Upon chilling, the expression of gma-miR397a, gma-miR166u and gma-miR171p was greatly upregulated, whereas the expression of gma-miR169c, gma-miR159b, gma-miR319a/b and gma-miR5559 was significantly decreased. The target genes of these miRNAs were predicted and validated using 5' complementary DNA ends (5'-RACE) experiments, and qPCR analysis identified putative genes targeted by the cold-responsive miRNAs in response to chilling temperatures. Taken together, our results reveal that miRNAs may be involved in the protective mechanism against chilling injury in mature nodules of soybean.

Taking it one step at a time in homologous recombination repair

The individual steps in the process of homologous recombination are particularly amenable to analysis by single-molecule imaging and manipulation experiments. Over the past 20 years these have provided a wealth of new information on the DNA transactions that make up this vital process. Exciting progress in developing new tools and techniques to analyze more complex components, dynamic reaction steps and molecular coordination continues at a rapid pace. Here we highlight recent results and indicate some emerging techniques likely to produce the next stage of advanced insight into homologous recombination. In this and related fields the future is bright.

Swi5-Sfr1 protein stimulates Rad51-mediated DNA strand exchange reaction through organization of DNA bases in the presynaptic filament

The Swi5-Sfr1 heterodimer protein stimulates the Rad51-promoted DNA strand exchange reaction, a crucial step in homologous recombination. To clarify how this accessory protein acts on the strand exchange reaction, we have analyzed how the structure of the primary reaction intermediate, the Rad51/single-stranded DNA (ssDNA) complex filament formed in the presence of ATP, is affected by Swi5-Sfr1. Using flow linear dichroism spectroscopy, we observe that the nucleobases of the ssDNA are more perpendicularly aligned to the filament axis in the presence of Swi5-Sfr1, whereas the bases are more randomly oriented in the absence of Swi5-Sfr1. When using a modified version of the natural protein where the N-terminal part of Sfr1 is deleted, which has no affinity for DNA but maintained ability to stimulate the strand exchange reaction, we still observe the improved perpendicular DNA base orientation. This indicates that Swi5-Sfr1 exerts its activating effect through interaction with the Rad51 filament mainly and not with the DNA. We propose that the role of a coplanar alignment of nucleobases induced by Swi5-Sfr1 in the presynaptic Rad51/ssDNA complex is to facilitate the critical matching with an invading double-stranded DNA, hence stimulating the strand exchange reaction.

Combined optical and topographic imaging reveals different arrangements of human RAD54 with presynaptic and postsynaptic RAD51-DNA filaments

Essential genome transactions, such as homologous recombination, are achieved by concerted and dynamic interactions of multiple protein components with DNA. Which proteins do what and how, will be reflected in their relative arrangements. However, obtaining high-resolution structural information on the variable arrangements of these complex assemblies is a challenge. Here we demonstrate the versatility of a combined total internal reflection fluorescence and scanning force microscope (TIRF-SFM) to pinpoint fluorescently labeled human homologous recombination protein RAD54 interacting with presynaptic (ssDNA) and postsynaptic (dsDNA) human recombinase RAD51 nucleoprotein filaments. Labeled proteins were localized by superresolution imaging on complex structures in the SFM image with high spatial accuracy. We observed some RAD54 at RAD51 filament ends, as expected. More commonly, RAD54 interspersed along RAD51-DNA filaments. RAD54 promotes RAD51-mediated DNA strand exchange and has been described to both stabilize and destabilize RAD51-DNA filaments. The different architectural arrangements we observe for RAD54 with RAD51-DNA filaments may reflect the diverse roles of this protein in homologous recombination.

Caffeine suppresses homologous recombination through interference with RAD51-mediated joint molecule formation

Caffeine is a widely used inhibitor of the protein kinases that play a central role in the DNA damage response. We used chemical inhibitors and genetically deficient mouse embryonic stem cell lines to study the role of DNA damage response in stable integration of the transfected DNA and found that caffeine rapidly, efficiently and reversibly inhibited homologous integration of the transfected DNA as measured by several homologous recombination-mediated gene-targeting assays. Biochemical and structural biology experiments revealed that caffeine interfered with a pivotal step in homologous recombination, homologous joint molecule formation, through increasing interactions of the RAD51 nucleoprotein filament with non-homologous DNA. Our results suggest that recombination pathways dependent on extensive homology search are caffeine-sensitive and stress the importance of considering direct checkpoint-independent mechanisms in the interpretation of the effects of caffeine on DNA repair.

Genome-wide identification and characterization of microRNA genes and their targets in flax (Linum usitatissimum): Characterization of flax miRNA genes

MicroRNAs (miRNAs) are small (20-24 nucleotide long) endogenous regulatory RNAs that play important roles in plant growth and development. They regulate gene expression at the post-transcriptional level by translational repression or target degradation and gene silencing. In this study, we identified 116 conserved miRNAs belonging to 23 families from the flax (Linum usitatissimum L.) genome using a computational approach. The precursor miRNAs varied in length; while most of the mature miRNAs were 21 nucleotide long, intergenic and showed conserved signatures of RNA polymerase II transcripts in their upstream regions. Promoter region analysis of the flax miRNA genes indicated prevalence of MYB transcription factor binding sites. Four miRNA gene clusters containing members of three phylogenetic groups were identified. Further, 142 target genes were predicted for these miRNAs and most of these represent transcriptional regulators. The miRNA encoding genes were expressed in diverse tissues as determined by digital expression analysis as well as real-time PCR. The expression of fourteen miRNAs and nine target genes was independently validated using the quantitative reverse transcription PCR (qRT-PCR). This study suggests that a large number of conserved plant miRNAs are also found in flax and these may play important roles in growth and development of flax.

Reconstitution of recombination-associated DNA synthesis with human proteins

The repair of DNA breaks by homologous recombination is a high-fidelity process, necessary for the maintenance of genome integrity. Thus, DNA synthesis associated with recombinational repair must be largely error-free. In this report, we show that human DNA polymerase delta (δ) is capable of robust DNA synthesis at RAD51-mediated recombination intermediates dependent on the processivity clamp PCNA. Translesion synthesis polymerase eta (η) also extends these substrates, albeit far less processively. The single-stranded DNA binding protein RPA facilitates recombination-mediated DNA synthesis by increasing the efficiency of primer utilization, preventing polymerase stalling at specific sequence contexts, and overcoming polymerase stalling caused by topological constraint allowing the transition to a migrating D-loop. Our results support a model whereby the high-fidelity replicative DNA polymerase δ performs recombination-associated DNA synthesis, with translesion synthesis polymerases providing a supportive role as in normal replication.

Probing Rad51-DNA interactions by changing DNA twist

In eukaryotes, Rad51 protein is responsible for the recombinational repair of double-strand DNA breaks. Rad51 monomers cooperatively assemble on exonuclease-processed broken ends forming helical nucleo-protein filaments that can pair with homologous regions of sister chromatids. Homologous pairing allows the broken ends to be reunited in a complex but error-free repair process. Rad51 protein has ATPase activity but its role is poorly understood, as homologous pairing is independent of adenosine triphosphate (ATP) hydrolysis. Here we use magnetic tweezers and electron microscopy to investigate how changes of DNA twist affect the structure of Rad51-DNA complexes and how ATP hydrolysis participates in this process. We show that Rad51 protein can bind to double-stranded DNA in two different modes depending on the enforced DNA twist. The stretching mode is observed when DNA is unwound towards a helical repeat of 18.6 bp/turn, whereas a non-stretching mode is observed when DNA molecules are not permitted to change their native helical repeat. We also show that the two forms of complexes are interconvertible and that by enforcing changes of DNA twist one can induce transitions between the two forms. Our observations permit a better understanding of the role of ATP hydrolysis in Rad51-mediated homologous pairing and strand exchange.

Mechanistic insight from chaos: how RecA mediates DNA strand exchange

Cleverly designed single-molecule FRET experiments reported in this issue of Structure by Ragunathan et al. coax RecA to reveal some of its secrets. Observing individual events identifies intermediate steps and provides clues for how to drive strand exchange forward.

Expression profile of small RNAs in Acacia mangium secondary xylem tissue with contrasting lignin content - potential regulatory sequences in monolignol biosynthetic pathway

BACKGROUND

Lignin, after cellulose, is the second most abundant biopolymer accounting for approximately 15-35% of the dry weight of wood. As an important component during wood formation, lignin is indispensable for plant structure and defense. However, it is an undesirable component in the pulp and paper industry. Removal of lignin from cellulose is costly and environmentally hazardous process. Tremendous efforts have been devoted to understand the role of enzymes and genes in controlling the amount and composition of lignin to be deposited in the cell wall. However, studies on the impact of downregulation and overexpression of monolignol biosynthesis genes in model species on lignin content, plant fitness and viability have been inconsistent. Recently, non-coding RNAs have been discovered to play an important role in regulating the entire monolignol biosynthesis pathway. As small RNAs have critical functions in various biological process during wood formation, small RNA profiling is an important tool for the identification of complete set of differentially expressed small RNAs between low lignin and high lignin secondary xylem.

RESULTS

In line with this, we have generated two small RNAs libraries from samples with contrasting lignin content using Illumina GAII sequencer. About 10 million sequence reads were obtained in secondary xylem of Am48 with high lignin content (41%) and a corresponding 14 million sequence reads were obtained in secondary xylem of Am54 with low lignin content (21%). Our results suggested that A. mangium small RNAs are composed of a set of 12 highly conserved miRNAs families found in plant miRNAs database, 82 novel miRNAs and a large proportion of non-conserved small RNAs with low expression levels. The predicted target genes of those differentially expressed conserved and non-conserved miRNAs include transcription factors associated with regulation of the lignin biosynthetic pathway genes. Some of these small RNAs play an important role in epigenetic silencing. Differential expression of the small RNAs between secondary xylem tissues with contrasting lignin content suggests that a cascade of miRNAs play an interconnected role in regulating the lignin biosynthetic pathway in Acacia species.

CONCLUSIONS

Our study critically demonstrated the roles of small RNAs during secondary wall formation. Comparison of the expression pattern of small RNAs between secondary xylem tissues with contrasting lignin content strongly indicated that small RNAs play a key regulatory role during lignin biosynthesis. Our analyses suggest an evolutionary mechanism for miRNA targets on the basis of the length of their 5' and 3' UTRs and their cellular roles. The results obtained can be used to better understand the roles of small RNAs during lignin biosynthesis and for the development of gene constructs for silencing of specific genes involved in monolignol biosynthesis with minimal effect on plant fitness and viability. For the first time, small RNAs were proven to play an important regulatory role during lignin biosynthesis in A. mangium.