The RAD51 and DMC1 homoeologous genes of bread wheat: cloning, molecular characterization and expression analysis
© Mayes et al; licensee BioMed Central Ltd. 2010
Received: 25 May 2010
Accepted: 29 September 2010
Published: 29 September 2010
Meiotic recombination in eukaryotes requires two homologues of the E. coli RecA proteins: Rad51 and Dmc1. Both proteins play important roles in the binding of single stranded DNA, homology search, strand invasion and strand exchange. Meiotic recombination has been well studied in Arabidopsis, rice, maize and the orthologues of RAD51 and DMC1 have been characterized. However genetic analysis of the RAD51 and DMC1 genes in bread wheat has been hampered due to the absence of complete sequence information and because of the existence of multiple copies of each gene in the hexaploid wheat genome.
In this study we have identified that TaRAD51 and TaDMC1 homoeologues are located on group 7 and group 5 chromosomes of hexaploid wheat, respectively. Comparative sequence analysis of cDNA derived from the TaRAD51 and TaDMC1 homoeologues revealed limited sequence divergence at both the nucleotide and the amino acid level. Indeed, comparisons between the predicted amino acid sequences of TaRAD51 and TaDMC1 and those of other eukaryotes reveal a high degree of evolutionary conservation. Despite the high degree of sequence conservation at the nucleotide level, genome-specific primers for cDNAs of TaRAD51 and TaDMC1 were developed to evaluate expression patterns of individual homoeologues during meiosis. QRT-PCR analysis showed that expression of the TaRAD51 and TaDMC1 cDNA homoeologues was largely restricted to meiotic tissue, with elevated levels observed during the stages of prophase I when meiotic recombination occurs. All three homoeologues of both strand-exchange proteins (TaRAD51 and TaDMC1) are expressed in wheat.
Bread wheat contains three expressed copies of each of the TaRAD51 and TaDMC1 homoeologues. While differences were detected between the three cDNA homoeologues of TaRAD51 as well as the three homoeologues of TaDMC1, it is unlikely that the predicted amino acid substitutions would have an effect on the protein structure, based on our three-dimensional structure prediction analyses. There are differences in the levels of expression of the three homoeologues of TaRAD51 and TaDMC1 as determined by QRT-PCR and if these differences are reflected at the protein level, bread wheat may be more dependent upon a particular homoeologue to achieve full fertility than all three equally.
Meiosis is central to all sexually reproducing living organisms and is a highly conserved process. During meiosis a single round of DNA replication is followed by two rounds of cell division resulting in the formation of four daughter cells each having half the chromosome number of the parent. Recombination at meiosis is pivotal not only for creating genetic diversity through crossing-over but also has a mechanistic role ensuring physical connections in the form of chiasmata that allow the chromosomes to segregate accurately during the first division of meiosis. Bacterial RecA was the first gene to be cloned from E. coli that was shown to play an important role in recombination [1, 2]. Meiotic recombination in eukaryotic cells generally requires two homologues of the E. coli RecA protein - Rad51 and Dmc1 - that were first identified in Saccharomyces cerevisiae[3, 4]. The Rad51 protein forms a nucleofilament with ATP-dependant strand-exchange activity, which has been shown to have roles in meiosis, homologous recombination and DSB repair . Dmc1 also forms a nucleofilament in cooperation with Rad51 and is involved in generating hybrid joint molecules specifically in meiosis . The RAD51 gene exists as a single copy in tomato and Arabidopsis, while maize and Physcomitrella patens have two copies each [6–11]. The DMC1 gene sequences have been reported only in a few plant species. While the Arabidopsis genome contains one copy of the DMC1 gene , rice has two copies, OsDMC1A and OsDMC1B, which are located on chromosomes 12 and 11, respectively, most likely produced by chromosomal duplication .
Common wheat or bread wheat is an allohexaploid composed of three genomes, each consisting of seven chromosomes (AABBDD; 2n = 6x = 42). Information on the genes and gene homoeologues of RAD51 and DMC1 in bread wheat is limited, primarily due to the added complexity of polyploidy. Given that RAD51 and DMC1 gene orthologues are reported to be involved in meiotic recombination in other cereal species and also the model plant Arabidopsis thaliana, we have investigated their orthologues in bread wheat. A previous study on meiotic recombination genes in wheat reported the existence of only one copy of TaDMC1 gene [GenBank accession number EU915561] and two copies of TaRAD51 gene - TaRAD51A1 [GenBank accession number EU915557] and TaRAD51A2 [GenBank accession number EU915558] . We report here for the first time the existence of three homoeologous copies of the TaRAD51 and TaDMC1 genes, their genomic location and expression patterns in inflorescence and vegetative tissues in bread wheat. These observations in the Chinese Spring wheat do not completely agree with the results reported earlier  and these apparent differences are discussed in this paper.
Plant materials and DNA isolation
Wild type hexaploid wheat (Triticum aestivum c.v Highbury) and Nullisomic-Tetrasomic (NT) wheat stocks (Triticum aestivum c.v Chinese Spring) were kindly provided by Prof John Snape and Dr Steve Reader at the John Innes Centre (Norwich, UK). The NT stocks, originally developed by Sears , were used for chromosomal assignment of the putative homoeologous genes. Lines Nullisomic for genome A chromosomes were tetrasomics either for B or D chromosomes, lines Nullisomic for genome B chromosomes were tetrasomics either for A or D chromosomes and Nullisomic for genome D chromosomes were tetrasomics either for A or B chromosomes.
Seeds of wheat were first germinated on moistened filter papers in Petri dishes incubated at 20°C in a growth chamber. When coleoptile length had reached ~1 cm, seedlings were transplanted to 13 cm pots containing John Innes No. 3 compost and grown under standard glasshouse conditions. Total genomic DNA was extracted from young leaf tissue (2 - 4 weeks post emergence) using the GenElute™ Plant Genomic DNA Miniprep kit (Sigma Aldrich) according to the manufacturer's instructions.
RNA extraction and first strand cDNA synthesis
Total RNA was isolated from the following tissues: root tips (RT, 14 days after germination), mature leaf (L) and immature inflorescence material (I). In addition, meiotically staged anthers were collected for RNA extraction from anther (microspore) tissues. Stages included are: pre-meiosis (PM), leptotene to pachytene (LP), diplotene to anaphase I (DA), telophase I to telophase II (TT), tetrads (T) and immature pollen (IP). Extractions used Trizol reagent (Gibco-BRL) according to the manufacturer's instructions. The isolated RNA was treated with RQ1 RNase-Free DNase (Promega) to remove any residual genomic DNA according to the manufacturers' instructions and first strand cDNA was synthesized using Superscript III Reverse Transcriptase (Invitrogen) using standard methods.
Isolation of full-length cDNA sequences of wheat RAD51 and DMC1 genes
List of primers used in the present study
Primer sequence (5'-3')
OsRAD51 Exon 1-2-F
Cloning of full-length cDNA of wheat RAD51
OsRAD51 Exon 1-2-R
OsRAD51 Exon 8-9-F
OsRAD51 Exon 8-9-R
OsRAD51 Exon 1-7-F
OsRAD51 Exon 1-7-R
OsRAD51 Exon 2-9-F
OsRAD51 Exon 2-9-R
Cloning of full-length cDNA of wheat DMC1
Not I(dT) anchor
TaRAD51 Exon 8-F
Amplification of intronic regions in TaRAD51 &TaDMC1 genes
450 & 500
TaRAD51 &TaDMC1 genome-specific primers based on intronic regions for localization on wheat genome
TaRAD51 cDNA full-length primers including UTR's
TaRAD51 genome-specific primers based on exons for QRT-PCR analyses
TaDMC1 cDNA full-length primers including UTR's
TaDMC1 genome-specific primers based on exons QRT-PCR analyses
QRT-PCR control primers
Generation of Genome-specific primers for TaRAD51 and TaDMC1 using intron sequences
OsRAD51A1 [GenBank accession number ABO8O261] and OsDMC1A [GenBank accession number AB079873] cDNA sequences were used to predict the exon/intron boundaries of TaRAD51 and TaDMC1, respectively. Once the boundaries were identified, exon-anchored primers that were expected to amplify products spanning two or three introns were designed (Table 1). Multiple alignment of sequences obtained from the TaRAD51 and TaDMC1 cDNA exon-anchored primers on Chinese Spring genomic DNA identified three distinct groups of sequences for each gene with respect to the intronic sequences in bread wheat (data not shown). Putative primers that were expected to amplify a particular homoeologue were designed based on the sequence polymorphisms detected and were coupled with a conserved exon-anchored primer in a PCR reaction containing genomic DNA from Chinese Spring and the appropriate nulli-tetrasomic stocks as templates.
Generation of Genome-specific primers for TaRAD51 and TaDMC1 sequences at the exon level
Highbury spikelet cDNA was amplified using primers that generate full-length copies of the TaRAD51 and TaDMC1 genes. A major amplification band, approximately 1.5 kb was obtained with both sets of primers. In the absence of differentially sized fragments, genome-specific primers for TaRAD51 and TaDMC1 cDNA homoeologues were designed based on sequence polymorphisms identified at the 3' UTRs for TaRAD51 cDNA homoeologues and at the 5' UTRs for the TaDMC1 cDNA homoeologues and were paired with a conserved gene primer (Table 1). The genome specificity of the primers was tested using appropriate nulli-tetrasomic stocks as templates.
Comparative amino acid sequence analysis and 3D protein modelling
Comparative amino acid sequence analysis was undertaken using the ClustalW tool of European Bioinformatics Institute (http://www.ebi.ac.uk/Tools/clustalw2). The amino acid sequences were processed using the Swiss model workspace (version 8.05, 2009) . The computed outputs were then used to predict the 3D protein structures using Swiss-Pdb Viewer (DeepView; Version 4.0, 2009) .
SIFT (Sorting Intolerant From Tolerant) analysis
The potential impact of amino-acid changes on the final protein structure were assessed with the SIFT program . SIFT is a program which predicts whether an amino acid substitution in the protein affects its protein function or not. SIFT assumes that important amino acids will be conserved in the protein family and so changes at well-conserved positions tend to be predicted as deleterious. SIFT score ranges from 0 to 1. The amino acid substitution is predicted deleterious if the score is ≤ 0.05 and tolerated if the score is > 0.05. This approach can only give an indication of potential changes to function through mutation.
The alignment of nucleotide sequences and subsequent selection of best-fit molecular evolution model was performed using Molecular Evolutionary Genetics Analysis (MEGA) software (version 5.0) . The construction of a phylogenetic tree was performed under maximum likelihood (ML) using RAxML version 7.2.6  on the CIPRES portal  with the following parameters: nucleotide search using the GTR + G model, random seed for initial parsimony inference, rapid bootstrapping with automatic bootstopping function (frequency criterion) followed by the final ML search.
Quantitative Real-Time PCR
Harvesting and staging of meiotic anthers for QRT-PCR was performed from immature wheat spikes . Harvested material was immediately frozen and stored at -80°C. Eight independent QRT-PCR amplifications were made for the following tissues: PM, LP, DA, TT, T, IP, L and RT. QRT-PCR was performed in triplicates using wheat Tubulin (Table 1) as an internal standard in a total volume of 12 μl reaction that contained 6 μl of Brilliant® SYBR® Green (Qiagen), 50 ng cDNA and 1 μl of an equal mix of forward and reverse gene specific primers (Table 1). The reaction conditions for QRT-PCR were as follows: 15 min at 95°C followed by 45 cycles of 30 s at 95°C, 30 s at gene specific annealing temperature and 30 s at 72°C. Primer sequences and annealing temperatures are given in Table 1.
Q-PCR amplification reactions were performed in the LightCycler® 480 instrument (Roche Applied Science, UK) using Brilliant® SYBR® Green (Qiagen) according to the manufacturer's instructions.
The cDNA sequences of the TaRAD51 and TaDMC1 homoeologues were deposited in GenBank at NCBI [GenBank accession numbers FJ594479, FJ594480 and FJ594481 for Ta cRAD51-7A, Ta cRAD51-7B and Ta cRAD51-7D homoeologues, respectively and FJ594476, FJ594477 and FJ594478 for Ta cDMC1-5A, Ta cDMC1-5B and Ta cDMC1-5D homoeologues, respectively)].
Results and Discussion
Isolation and genome localization of RAD51 and DMC1 gene homoeologous sequences of hexaploid wheat
TaRAD51 and TaDMC1 cDNA homoeologues are highly conserved at the amino acid level
Phylogenetic analysis of TaRAD51 and TaDMC1 cDNA homoeologues reveals evolutionary conservation
This study also found that there was only one copy of TaRAD51 per haploid genome of wheat in contrast to two copies of the RAD51 gene reported in rice (OsRAD51A1/A2; rice chromosomes 11 and 12) and maize (ZmRAD51A1/A2; maize chromosomes 7 and 8) [8, 13]. The presence of two RAD51 copies in rice and maize and only one copy of RAD51 per genome in wheat may argue for gene loss in the ancestors of wheat. This conclusion is supported by the phylogenetic tree (Figure 6). The genomes of rice and maize contain both products of a predicted earlier gene duplication before divergence of the two species i.e RAD51A1 and RAD51A2, while all three of the RAD51 genes in wheat appear to be homologues of the RAD51A2 gene, suggesting that the RAD51A1 lineage has been lost in bread wheat. There are also two versions of the DMC1 gene in rice and maize and only one in the haploid wheat genome. This may have a similar origin, but the greater conservation of DMC1 in these species makes an assessment more difficult and the presence of two copies may simply reflect a local duplication of R11 and R12 in rice and the fact that maize is an ancient tetraploid. In both case, it is worth noting that DMC1 homologues are more highly conserved than the RAD51 homologues, arguing for a higher level of selective pressure preventing mutation or a more recent duplication event.
All homoeologues of TaRAD51 and TaDMC1 genes were highly expressed during prophase I of meiosis
QRT-PCR analysis of TaDMC1 cDNA homoeologues indicated that all three homoeologues were expressed in both vegetative and meiotic tissues (Figure 7b). However higher expression levels were found in meiotic stages (Figure 7b, stages PM to TT) compared to non-meiotic tissues and we observed an increase in expression of all three TaDMC1 cDNA homoeologues during the Tetrad stage. There were differences in the expression levels among all three cDNA homoeologues of TaDMC1 in the observed tissues, which suggests that they could contribute to differing extents in meiotic recombination. Both TaDMC1-5B and TaDMC1-5D homoeologues were highly expressed in all meiotic stages, although not the Tetrad stage, compared to the transcript level of TaDMC1-5A suggesting they could be the versions with the greatest role in meiosis.
Overall, we have reported the cloning and characterisation of the hexaploid wheat versions of the major strand-exchange proteins DMC1 and RAD51. Primer tests have been developed which are specific to individual homoeologues and these have been used to characterise the levels of genome-specific expression in a range of tissues through Q-PCR. The same primers could be used to screen mutation and deletion stocks of wheat to detect loss-of-function mutations, allowing in planta characterisation of the effect of these mutations on meiosis.
Double Strand Breaks
Rapid Amplification of cDNA Ends
Reverse Transcription Polymerase Chain Reaction.
The authors gratefully acknowledge Prof John Snape for Highbury wheat lines, Dr Steve Reader for Nullisomic-Tetrasomic wheat stocks, Dr. Hongying Li for helping with Q-PCR experiments and Dr Jason Able for providing valuable suggestions for meiotic staging of wheat anthers and for useful discussion during the project. UKD was supported by a University of Nottingham Fees PhD Scholarship. The work was partly supported by a New Lecturer grant to SM.
- Clark AJ, Margulies AD: Isolation and characterization of recombination-deficient mutants of Escherichia coli K12. Proc Natl Acad Sci USA. 1965, 53: 451-459. 10.1073/pnas.53.2.451.PubMed CentralPubMedView ArticleGoogle Scholar
- Lusetti , Cox MM: The bacterial RecA protein and the recombinational DNA repair of stalled replication forks. Annu Rev Biochem. 2002, 71: 71-100. 10.1146/annurev.biochem.71.083101.133940.PubMedView ArticleGoogle Scholar
- Bishop DK, Park D, Xu L, Kleckner N: DMC1: a meiosis-specific yeast homolog of E. coli recA required for recombination, synaptonemal complex formation, and cell cycle progression. Cell. 1992, 69: 439-456. 10.1016/0092-8674(92)90446-J.PubMedView ArticleGoogle Scholar
- Shinohara A, Ogawa H, Ogawa T: Rad51 protein involved in repair and recombination in S. cerevisiae is a RecA-like protein. Cell. 1992, 69: 457-470. 10.1016/0092-8674(92)90447-K.PubMedView ArticleGoogle Scholar
- Symington LS: Role of RAD52 epistasis group genes in homologous recombination and double-strand break repair. Microbiol Mol Biol Rev. 2002, 66: 630-670. 10.1128/MMBR.66.4.630-670.2002.PubMed CentralPubMedView ArticleGoogle Scholar
- Stassen NY, Logsdon JM, Vora GV, Offenberg HH, Palmer JD, Zolan ME: Isolation and characterization of rad51 orthologs from Coprinus cinereus and Lycopersicon esculentum, and phylogenetic analysis of eukaryotic RecA homologs. Curr Genet. 1997, 31: 144-157. 10.1007/s002940050189.PubMedView ArticleGoogle Scholar
- Doutriaux MP, Couteau F, Bergounioux C, White C: Isolation and characterisation of the RAD51 and DMC1 homologs from Arabidopsis thaliana. Mol Gen Genet. 1998, 257: 283-291. 10.1007/s004380050649.PubMedView ArticleGoogle Scholar
- Franklin AE, McElver J, Sunjevaric I, Rothstein R, Bowen B, Cande WZ: Three-dimension microscopy of the RAD51 recombination protein during meiotic prophase. The Plant Cell. 1999, 11: 809-824. 10.2307/3870816.PubMed CentralPubMedView ArticleGoogle Scholar
- Ayora S, Piruat JI, Luna R, Reiss B, Russo VEA, Aguilera A, Alonso JC: Characterization of two highly similar Rad51 homologs of Physcomitrella patens. J Mol Biol. 2002, 316: 35-49. 10.1006/jmbi.2001.5336.PubMedView ArticleGoogle Scholar
- Markmann-Mulisch U, Hadi MZ, Koepchen K, Alonso JC, Russo VEA, Schell J, Reiss B: The organization of Physcomitrella patens RAD51 genes is unique among eukaryotic organisms. Proc Natl Acad Sci USA. 2002, 99: 2959-2964. 10.1073/pnas.032668199.PubMed CentralPubMedView ArticleGoogle Scholar
- Li W, Changbin C, Ullrich MM, Ljudmilla T, Elmon S, Ma H, Reiss R: The Arabidopsis AtRAD51 gene is dispensable for vegetative development but required for meiosis. Proc Natl Acad Sci USA. 2004, 101: 10596-10601. 10.1073/pnas.0404110101.PubMed CentralPubMedView ArticleGoogle Scholar
- Klimyuk VI, Jones JD: AtDMC1, the Arabidopsis homologue of the yeast DMC1 gene: characterization, transposon-induced allelic variation and meiosis-assosciated expression. Plant J. 1997, 11: 1-14. 10.1046/j.1365-313X.1997.11010001.x.PubMedView ArticleGoogle Scholar
- Kathiresan A, Khush GS, Bennett J: Two rice DMC1 genes are differentially expressed during meiosis and during haploid and diploid mitosis. Sex Plant Reprod. 2002, 14: 257-267. 10.1007/s00497-001-0113-5.View ArticleGoogle Scholar
- Khoo KHP, Jolly HR, Able JA: The RAD51 gene family in bread wheat is highly conserved across eukaryotes, with RAD51A1 upregulated during early meiosis. Funct Plant Biol. 2008, 35: 1267-1277. 10.1071/FP08203.View ArticleGoogle Scholar
- Sears ER: Nullisomic-tetrasomic combinations in hexaploid wheat. Chromosome manipulation and plant Genetics. Edited by: Riley R, Lewis KR. 1966, Oliver & Boyd, Edinburgh, 29-45.View ArticleGoogle Scholar
- Bordoli L, Kiefer F, Arnold K, Benkert P, Battey J, Schwede T: Protein structure homology modelling using SWISS-MODEL Workspace. Nat Protoc. 2009, 4: 1-10.1038/nprot.2008.197.PubMedView ArticleGoogle Scholar
- Kaplan Warren, Littlejohn TG: Swiss-PDB Viewer (Deep View). Brief Bioinform. 2001, 2: 195-197. 10.1093/bib/2.2.195.PubMedView ArticleGoogle Scholar
- Pauline CN, Henikoff S: SIFT: predicting amino acid changes that affect protein function. Nucl Acid Res. 2003, 31: 3812-3814. 10.1093/nar/gkg509.View ArticleGoogle Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol and Evol. 2007, 24: 1596-1599. 10.1093/molbev/msm092.View ArticleGoogle Scholar
- Stamatakis A, Ludwig T, Meier H: Raxml-iii: A fast program for maximum likelihood-based inference of large phylogenetic trees. Bioinformatics. 2005, 21: 456-463. 10.1093/bioinformatics/bti191.PubMedView ArticleGoogle Scholar
- Miller MA, Holder MT, Vos R, Midford PE, Liebowitz T, Chan L, Hoover P, Warnow T: The CIPRES Portals. --- Either ISSN or Journal title must be supplied.. [http://www.phylo.org/sub_sections/portal]
- Crismani W, Baumann U, Sutton T, Shirley N, Webster T, Spangenberg G, Langridge P, Able J: Microarray expression analysis of meiosis and microsporogenesis in hexaploid bread wheat. BMC Genomics. 2006, 7-Google Scholar
- Lin Z, Kong H, Nei M, Ma H: Origins and evolution of the recA/RAD51 gene family: Evidence for ancient gene duplication and endosymbiotic gene transfer. Proc Natl Acad Sci USA. 2006, 103: 10328-10333. 10.1073/pnas.0604232103.PubMed CentralPubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.