Expressional dynamics of minisatellite 33.15 tagged spermatozoal transcriptome in Bubalus bubalis
© Srivastava et al; licensee BioMed Central Ltd. 2009
Received: 03 February 2009
Accepted: 07 July 2009
Published: 07 July 2009
Transcriptionally quiescent spermatozoa have been established to be a repository of mRNA coding for several functionally essential cellular proteins. This entourage of mRNA is envisaged to be involved in post-fertilization and early embryogenesis. Minisatellites tagged with mRNA transcripts have been implicated with gene organization, regulation and function. However, the organization and expression of the minisatellite tagged transcript diversity, particularly in spermatozoa, remains unclear.
In the present study, we identified and characterized 12 mRNA transcripts from the spermatozoa of water buffalo Bubalus bubalis employing minisatellite associated sequence amplification (MASA) and a consensus sequence of 33.15 repeat loci. Of these 33.15 tagged transcripts, only one was found to be homologous to Bovine steroid 21-hydroxylase (P-450-c21) gene. Other ten transcripts showed significant similarity with various mRNAs or chromosomal contigs across the species. The remaining one construed to be novel since this was unreported in the database (NCBI GenBank). All these uncharacterized and known transcripts showed highest expression in testis and spermatozoa compared to that in somatic tissues and ovary. Of these 12 mRNA transcripts, 4 showed differential expression in the forebrain and hindbrain of buffalo. Moreover, genes corresponding to all the 33.15 tagged spermatozoal transcripts were found to be conserved across 13 other species analyzed.
Our results show MASA as an important tool to capture mRNA transcript diversity tagged with minisatellites in the spermatozoa. Comprehensive characterization of these transcripts is envisaged to augment our understanding on the genes involved in testicular functions and sustenance of a viable paternal genome during pre- and post- fertilization events and early stages of development. Prospects of this approach in genome analysis in general and comparative genomics in particular are highlighted.
Ejaculated spermatozoa represent terminally differentiated cells in which transcription and/or translation of nuclear encoded mRNAs are considered to be unlikely. Therefore, the paternal genome is the only cargo carried by the spermatozoa to the ooplasm. Discovery of several transcription factors, soluble signaling molecules and structures delivered by spermatozoan into the zygotic cytoplasm upon fertilization have changed this perception [1–3]. Despite the transcriptionally inactive state, spermatozoa retain an entourage of mRNA transcripts encoding transcription factors and proteins involved in signal transduction, cell proliferation, chromatin condensation, regulation of sperm motility, capacitation and acrosome reaction [2–7]. Delivery of such spermatozoal transcripts to ooplasm is envisaged to have their potential significance during fertilization, embryogenesis and morphogenesis. Approximately 3000–5000 mRNA transcripts have been suggested to be present in the spermatozoa [2–7]. However, these remain to be characterized for their organization, expression and association with different regulatory elements and repetitive sequences.
Repetitive DNA sequences are dynamic components of the genome encompassing transposable elements, major satellites, minisatellites and microsatellites [8, 9]. Most of these repeats are found in the non-coding regions of the genomes while a small fraction is retained within the transcriptome [10–12] and participate in gene regulation through transcription, translation or gene silencing [13–15]. Surprisingly, organization of the repeats within the transcriptomes of any mammalian species, particularly in the spermatozoa has remained undeciphered. To explore the organization of such repeats and the repeat tagged genes, we undertook the transcriptome analysis of water buffalo Bubalus bubalis, an important player in the agriculture, dairy and meat industries in the Indian sub-continent.
Minisatellite 33.15 (5' CACCTCTCCACCTGCC 3') originating from the human myoglobin gene (7q35-q36) has been studied in a number of species [16–18]. This repeat has also been found to be associated with the heterochromatic sequences of the human Y chromosome . In our previous in-silico study, we demonstrated presence of minisatellite 33.15 within the transcriptomes of several eukaryotes. Following this, we isolated and characterized several known and novel mRNA transcripts tagged with the consensus of 33.15 repeat loci from different somatic tissues and gonads of water buffalo Bubalus bubalis .
Owing to the envisaged participation of the spermatozoal mRNA transcripts during early zygotic and embryonic development, we studied the spermatozoal transcriptome of water buffalo Bubalus bubalis tagged with minisatellite 33.15 employing Minisatellite associated sequence amplification (MASA). These mRNA transcripts were further characterized for their sequence organization, homology status, variation of gene expression, copy number and evolutionary status of the corresponding genes.
Differential distribution of the consensus sequence of 33.15 repeat loci within the spermatozoal and somatic transcriptomes
Detailed analysis of the spermatozoal mRNA transcripts tagged with consensus of 33.15 repeat loci in buffalo Bubalus bubalis and their homology or similarity with genes across the species#
Size (In bp)
Description of Blast search
Accession no. of similar genes
Gene length (bp)
Nucleotide position of 33.6 tagged transcripts
1. Bos taurus chromosome 1 reference assembly (based on Btau 4.0); featuring in this part similar to uncharacterized protein C3orf44
2. Homo sapiens cDNA FLJ30354 fis, clone BRACE2007682
3. Homo sapiens chromosome 19 clone CTD-2292J18
4. Homo sapiens tetratricopeptide repeat protein 12 (TTC12) gene, alternatively spliced
1. Bos taurus chromosome 1 reference assembly (based on Btau 4.0)
2. Bovine steroid 21-hydroxylase gene (P-450-c21) gene
3. Ovis aries cytochrome P-450 steroid 21 hydroxylase (P-450-c21) gene
4. Sheep cytochrome protein (P450C21C) mRNA
5. Bovine microsomal cytochrome P-450(C21) mRNA
6. Bos taurus cytochrome P450, subfamily XXI (CYP21A2), mRNA
1. Homo sapiens BAC clone RP11-498O22 from 2
1. Bos taurus chromosome 5 reference assembly (based on Btau 4.0)
1. Bos taurus chromosome 9 reference assembly (based on Btau 4.0)
2. Human DNA sequence from clone RP1-230L10 on chromosome 6 Contains part of a novel gene
1. Bos taurus chromosome 29, reference assembly (based on Btau_4.0),
2. Bos taurus chromosome 16, reference assembly (based on Btau_4.0),
3. Bos taurus Y Chr NOVECTOR CH240-223B23 (Children's Hospital Oakland Research Institute Bovine BAC Library (male
4. Bos taurus molybdenum cofactor synthesis 3 (MOCS3), mRNA
5. Bos taurus T cell receptor gamma cluster 2 (TCRG2) gene
6. Bos taurus clone p68.3 MHC class I antigen gene, 5' UTR and partial cds
1. Bos taurus chromosome 2, reference assembly (based on Btau_4.0)
2. Bos taurus chromosome 11, reference assembly
2. Bos taurus chromosome X, reference assembly
3. Bos taurus ER degradation enhancer, mannosidase alpha-like 2 (EDEM2), mRNA
4. Bos taurus similar to mutant UDP-N-acetylglucosamine 2-epimerase/N-acetylman-nosamine kinase (GNE), mRNA
5. Rhinolophus ferrumequinum clone VMRC7-44G7
6.Felis catus clone RP86-177C24
7. Equus caballus subclone 355H03 3.3 MHC class I antigen gene
1. Bos taurus chromosome 10, reference assembly (based on Btau_4.0)
2. Homo sapiens chromosome 15 clone RP11-358M11 map 15q21.3
1. Bos taurus chromosome 2, reference assembly (based on Btau_4.0)
2. Bos taurus Y chr BAC CH240-96I24 (Children's Hospital Oakland Research Institute Bovine BAC Library (male)
3. Bos taurus Fanconi anemia, complementation group M, mRNA
4. Bos taurus similar to zinc finger protein 568 (ZNF793), mRNA
5. Bos taurus similar to peroxisomal long-chain acyl-coA thioesterase (LOC516537), mRNA
1. Bos taurus chromosome 12, reference assembly (based on Btau_4.0
2. Bos taurus BAC CH240-10G15 (Children's Hospital Oakland Research Institute Bovine BAC Library (male)
3. Bos taurus cDNA clone IMAGE:8284348 & CR4-NOT transcription complex, subunit 8
4. O. aries KII-9 gene for hair type II keratin intermediate filament
5. Bos taurus similar to Lysophosphatidic acid receptor Edg-7 (LPA receptor 3) (LPA-3) (EDG7), mRNA
6. Bos taurus similar to Myotubularin-related protein 2, mRNA
1. Bos taurus chromosome 16, reference assembly (based on Btau_4.0)
Conserved 33.15 tagged spermatozoal genes are present as single to double copy in buffalo genome
Relative quantitative expression and copy number status of the spermatozoal genes/gene fragments tagged with the consensus of 33.15 repeat loci
Relative Expression in Different Tissues and Spermatozoa (In folds)
Spermatozoa of Four Different Animals
Status Per Haploid
Differential expression profiles of 33.15 tagged mRNA transcripts
The primer specificity for each individual gene was established using five fold dilution series of the template cDNA (Additional file 1C–D). Straight standard curve with a slope = (-3.4) to (-3.6) and a single peak in the dissociation protocol established the maximum specificity and efficiency of the primers. For more accuracy, the expression study was conducted using three different dilutions of cDNA. The relative quantitation of expression so obtained was substantiated further using somatic/gonadal/spermatozoal cDNA from five additional animals.
Minisatellites have been implicated with gene regulation, chromosomal fragile sites and genome imprinting [15, 21, 22] but the biological significance of their association with the coding regions remains largely unresolved. In the present study, we established the association of 33.15 repeats with the spermatozoal transcriptome of water buffalo, Bubalus bubalis and observed highest expression of these transcripts in testis and/or spermatozoa indicating their crucial roles in male gametogenesis.
Potential implications of association of repeats with the spermatozoal transcriptome
Our previous study showed distribution of 33.15 sequences in the flanking regions and within the mRNA transcripts of many structural, functional and regulatory genes in buffalo and six other species . Previously, various 33.15 tagged mRNA transcripts were uncovered from somatic tissues and gonads of buffalo. Here, a total of 12 different mRNA transcripts, representing one known and others mostly uncharacterized ones from spermatozoa were identified (Figure 1 and Table 1) complementing our earlier study. Notably, none of the transcripts was shared between somatic/gonadal tissues and spermatozoa. The differential transcript profile may be a reflection of their diverse functions in somatic tissues, gonads (testis/ovary), and spermatozoa. Alternatively, they may have specific roles to play during different stages of development. These transcripts could not be picked up in somatic tissues/gonad and vice-versa due to either polymorphic nature of SSRs or their much lower numbers. A comparison of the mRNA transcripts uncovered in the present study with those reported earlier  confirmed that these were the different ones.
Thus, present study can be used as a basis for evaluation of other repeats to unveil their combined significance within and adjacent to the coding regions. Spermatozoa have come a long way from being regarded as an artifact to a repository of a variety of RNAs. The initial controversy which shrouded the existence of spermatozoal RNA pool has now been cleared by the identification of several transcripts in different species including human [2–7].
Many signaling molecules and transcription factors have been identified which are transported into the zygotic cytoplasm upon fertilization by the spermatozoa [2–7]. However, several other spermatozoal transcripts remain to be characterized. Present study seems to be the first one, unveiling existence of the minisatellite tagged transcripts in the buffalo spermatozoa suggesting their predominant involvements in the pre- and post-fertilization events. It becomes more significant, when these transcripts showed faithful conservation across thirteen different species suggesting their testicular functions in much broader number of higher eukaryotes.
Prospects of repeat tagged transcripts carrying highest expression in the testis and spermatozoa
All the 33.15 tagged genes/gene fragments showing highest expression in the testis and spermatozoa suggested their strong participations in the testicular development and spermatogenesis. Further, uniform expression of steroid 21-hydroxylase gene suggested its consistent requirement in all the tissues except spermatozoa. The steroid 21-hydroxylase; P450 (C21) belongs to the cytochrome P-450 super-gene family and plays a crucial role in the synthesis of steroid hormones such as cortisol and aldosterone. Mutations in this gene have been implicated in causing Adrenal Hyperplasia .
This gene has been reported to show maximum expression in adrenal glands and minimum in liver and testis. However, in our study, this gene showed maximum expression in spermatozoa thus, adding a new dimension towards its unidentified functions in the sperm development and maintenance. Expression of P450 (C21) gene could not be confirmed in the adrenal gland of buffalo due to unavailability of the tissue sample. Since the spermatozoal mRNA is delivered to the ooplasm, and has its potential roles in the chromatin repackaging, early zygotic and embryonic development and post syngamy events [2, 3, 5, 23], the observed highest expression of all the 33.15 tagged mRNA transcripts in spermatozoa suggest their important roles towards specific functions of paternal genome sustenance.
Present work demonstrated the association of consensus sequence of minisatellite 33.15 with buffalo spermatozoal transcriptome. Maximum expression of all the 33.15-tagged mRNA transcripts in testis and spermatozoa suggested their involvement in various testicular functions and fertilization. Following this approach, other repeat tagged mRNA transcripts may be studied within the buffalo genome and across the species. A comparative mRNA fingerprint based on several minisatellites resulting in detailed expression profile would allow the identification of common and uncommon transcripts in somatic tissues and gonads. Indeed, this would help understanding the more focused roles of repetitive DNA sequences in the paternal genome during and post fertilization events and different stages of embryonic development.
Sperm purification and RNA isolation
Fresh ejaculated semen samples of buffaloes were obtained from the local dairy farm following strictly the guidelines of the Institute's Bio-Ethical and Biosafety Committees. All the sperm samples were procured from the buffalo in-vitro fertilization (IVF) center where utmost care was taken to avoid the diploid cell contaminations in the sperm preparations using highly sensitive methodologies and modern techniques. Samples were subjected to percoll gradient method to select only motile sperms as described earlier . Sperm cells were washed in sperm wash buffer (0.15 mM NaCl and 10 mM EDTA) and pellets were processed for RNA extraction using TRIzol reagent (Invitrogen, CA). After re-suspension of the cells in TRIzol reagent, cells were incubated at 60°C for 30 min as described earlier , vortexing every 10 min for lysis following manufacturer's instructions. Each RNA sample was treated with RNase free DNase I (Ambion, USA), and precipitated with ammonium acetate and ethanol following standard protocols [23–25].
Quantification of RNA was done with spectrophotometer (Amersham Life Sciences). A fraction of each RNA sample (5 ng/μl) was tested for residual DNA contamination by PCR using primers against β-actin (ACTB; GenBank accession no. DQ661647) and Protamine-1 (PRM1; GenBank accession no. NM_174156) following standard procedures [23, 25] (Additional file 2). These exon-specific primers encompassed one intron, thus, distinguishing cDNA from the genomic sequence. The PCR was conducted using initial denaturation at 95°C for 5 min; 35 cycles of denaturation at 94°C for 30 sec, primer annealing at 55°C for PRM1 and 60°C for ACTB for 30 sec, extension at 72°C for 1 min; and final extension at 72°C for 10 min.
Isolation of genomic DNA, total RNA and cDNA synthesis
Blood and tissue samples of both the sexes of water buffalo were collected from local slaughterhouse, following guidelines of the Institute's Ethical and Bio-safety Committees. Details of the genomic DNA isolation from buffalo and other species used in this study for cross hybridization have been reported earlier [20, 26]. Total RNA was also isolated from Testis, Ovary, Spleen, Liver, Lung, Kidney, Heart and Brain using TRIzol reagent (Invitrogen, CA) as per the standard procedure [23–25]. Following quantification of RNA, the cDNA synthesis was conducted using commercially available cDNA archive kit (ABI, USA). The quality of cDNA so prepared was confirmed by PCR amplification using β-actin and PRM1 primers. The absence of somatic and other non-sperm cells in each sperm RNA preparation was confirmed by PCR reactions using each cDNA sample, and common leukocyte antigen (CD45) and epithelial E-cadherin (CDH1) gene markers (Additional files 2 &3) at 55°C annealing temperature.
Minisatellite associated sequence amplification (MASA)
For conducting minisatellite associated sequence amplification (MASA), a 16 base long oligo (5' CACCTCTCCACCTGCC 3') based on consensus of 33.15 repeat loci, was purchased from Microsynth GmbH (Balgach, Switzerland). MASA reactions were performed using cDNA templates from the spermatozoa of different animals following standard procedures [20, 26]. The resultant amplicons were resolved on 20 cm long 2.5% (w/v) agarose gel using 0.5× TBE buffer.
Cloning, sequencing and characterization of MASA uncovered amplicons
From the MASA reactions, a total of 72 amplicons were uncovered (Table 1). These amplicons resolved on the agarose gel were sliced; DNA eluted (Qiagen Gel Extraction kit, Germany) and processed independently for cloning into pGEMT-easy vector (Promega, USA). A total of 20–25 recombinant clones corresponding to each amplicon were characterized by restriction digestion and slot blot hybridization using labeled buffalo genomic DNA . Sequences of all the clones generated from every single amplicon were independently subjected to ClustalW alignment to ascertain interclonal variation. Finally, a total of 12 sequences were deposited in the NCBI GenBank and their accession numbers are given in the table 1.
Database search was conducted to determine homology of these sequences independently with other entries in the GenBank using default NCBI server http://blast.ncbi.nlm.nih.gov/Blast.cgi having megablast "highly similar" and blastn "somewhat similar" algorithms. Each sequence was first subjected to blast search across the databases of Nucleotide collection (nr/nt), Reference mRNA sequences (refseq_mRNA), NCBI genomes (chromosome) and Expressed sequence tags (ESTs) of all the organisms, and thereafter individually to those of cattle, sheep, goat and human.
Cross-hybridization of the uncovered genes/gene fragments
Cross hybridization studies were conducted using DNA from across the species. DNA was extracted from the peripheral blood of buffalo Bubalus bubalis, cattle Bos indicus, sheep Ovis aries, goat Capra hircus, human Homo sapiens, Pigeon Columba livia, pig Sus scrofa, Baboon Papio hamadryas, Bonnet monkey Macaca radiata, Langur Presbytis entellus, Rhesus monkey Macaca mulatta, Lion Panthera leo, Tiger Tigris tigris following standard protocols [25, 26]. Lion and Tiger blood samples were procured with due approval of the competent authorities of the States and Union Government of India. Hybridization of genomic DNA from different sources using recombinant cloned probes was conducted following standard procedures [25, 27].
RNA slot blot analysis, RT-PCR and Southern Blotting
For RNA slot blot analysis, approximately 2 μg of total RNA from different tissues of buffalo in 100 μl of 2 × SSC was slot blotted onto the nylon membrane (Minifold Apparatus, Schleicher & Schuell, Germany) and UV fixed. For positive control, 5 ng of recombinant plasmid, each, was included in the blot(s). Hybridizations were performed under high stringent conditions using standard procedure [25, 26]. Individual probes for each fragment was labeled with [32P] α-dCTP using rediprime™ II kit (Amersham Pharmacia Biotech, USA). In order to confirm the RNA slot blot, internal primers were designed from each fragment (Additional file 4) and RT-PCR was conducted using cDNA templates from different tissues following standard thermal profiles (Additional file 4). The products were resolved on the agarose gel and transferred onto the nylon membrane followed by hybridization with [32P] α-dCTP labeled recombinant clones corresponding to each uncovered fragment using standard procedures [25–27]. Bubaline derived β-actin gene probe and bacterial genomic DNA were used as positive and negative controls, respectively.
Relative Expression Studies Using Real Time PCR
For expression quantification of different genes/fragments (JS1–JS12), SYBR green assays were conducted using Real Time PCR (Sequence Detection System, 7500, ABI) for individual fragments using equal amount of cDNA from all the tissues and spermatozoa. The GAPDH and ACTB genes were used as endogenous controls for relative expression. The primers for assessing relative expression of each of the transcripts (JS1–JS12, GAPDH &ACTB) and the copy numbers of their corresponding genes in the haploid genome were designed using Primer Express 3.0 software (Applied Biosystems, USA). The primer details have been given in additional file 4. Standard curves for all the genes were composed of five points covering 10 log orders of magnitude using a serial dilution of a purified cDNA template. Slopes of these standard curves were obtained between -3.3 and -3.6, reflecting amplification kinetics ranging from 80% to 100%.
Melting-curve analysis was performed for SYBR Green amplifications by plotting fluorescence intensity in a graphic model. To validate the specificity of the qRT-PCR assays, a single melting temperature peak representing a single amplicon (SYBR Green amplifications) was confirmed. Thereafter, the quantitative RT-PCR assay was performed individually for JS1–JS12, GAPDH and ACTB genes (Additional file 4) using SYBR Green chemistry containing 5 μl of 2× Power SYBR Green Master Mix (Applied Biosystems), 0.3 μM of forward and reverse primers each, and 0.25 U of AmpErase UNG (Applied Biosystems). All the runs included a triplicate of the melting curve for each sample and from three to six negative controls (absence of target DNA) that were partitioned to detect potential contamination during preparation of the plate. Thermal cycling conditions used were: 50°C for 2 min (incubation for the AmpErase UNG) and then the first denaturation step at 95°C for 10 min, followed by 40 cycles of 95°C for 15 sec and 60°C for 1 min. Further, each experiment was repeated three times at different concentrations of the template to ensure consistency of the results. Expression level of the genes was calculated using the formula: expression status = (1+E)-ΔΔCt, where E is the efficiency of the PCR and ΔΔCt is the cycle threshold normalized first with the endogenous control GAPDH (Ct sample - Ct GAPDH = ΔCt) and then with the calibrator sample (ΔCt Sample - ΔCt Calibrator = ΔΔCt) [23, 25, 26].
Copy number calculation of the 33.15 tagged genes using Real Time PCR
The copy number status of the genes/gene fragments corresponding to each 33.15 tagged transcripts was calculated based on absolute quantification assay using SYBR green chemistry. The primers standardized and checked for the quantitative expression for JS1–JS12 were employed also for copy number calculation (Additional file 4). The gene copy number of individual genes/gene fragments corresponding to the transcripts tagged with minisatellite 33.15 was calculated using standard curve method. A ten fold dilution series of known copies of each recombinant clone (for individual JS1–JS12) was used to plot standard curve with slope -3.3 to -3.4 substantiating maximum efficiency of the real time PCR reaction. Extrapolation of standard curves for individual recombinant clones uncovered copy number status of each 33.15 tagged gene in the haploid buffalo genome.
This work was supported by a DBT Grant No. BT/PR8476/AAQ/01/315/2006 to SA and a core grant from the Department of Biotechnology, Govt. of India to the National Institute of Immunology, New Delhi. SP acknowledges the Senior Research Fellowship from the Council of Scientific and Industrial Research, New Delhi. Equipment donation from the Alexander Von Humboldt Foundation, Bonn, Germany is gratefully acknowledged. We thank Shri Khem Singh Negi for the technical assistance.
- Saunders CM, Larman GM, Parrington J, Cox LJ, Royse J, Blayney LM, Swann K, Lai FA: PLC zeta: a sperm-specific trigger of Ca(2+) oscillations in eggs and embryo development. Development. 2002, 129 (15): 3533-3544.PubMedGoogle Scholar
- Krawetz SA: Paternal contribution: new insights and future challenges. Nat Rev Genet. 2005, 6: 633-642. 10.1038/nrg1654.View ArticlePubMedGoogle Scholar
- Miller D, Ostermeier GC, Krawetz SA: The controversy, potential and roles of spermatozoal RNA. Trends Mol Med. 2005, 11 (4): 156-163. 10.1016/j.molmed.2005.02.006.View ArticlePubMedGoogle Scholar
- Wykes SM, Visscher DW, Krawetz SA: Haploid transcripts persist in mature human spermatozoa. Mol Hum Reprod. 1997, 3 (1): 15-19. 10.1093/molehr/3.1.15.View ArticlePubMedGoogle Scholar
- Miller D: Analysis and significance of messenger RNA in human ejaculated spermatozoa. Mol Reprod Dev. 2000, 56: 259-264. 10.1002/(SICI)1098-2795(200006)56:2+<259::AID-MRD10>3.0.CO;2-R.View ArticlePubMedGoogle Scholar
- Lambard S, Galeraud-Denis I, Martin R, Levy G, Chocat A, Carreau S: Analysis and significance of mRNA in human ejaculated sperm from normozoospermic donors: relationship to sperm motility and capacitation. Mol Hum Reprod. 2004, 10 (7): 535-541. 10.1093/molehr/gah064.View ArticlePubMedGoogle Scholar
- Ostermeier GC, Goodrich RJ, Moldenhauer JS, Diamond MP, Krawetz SA: A suite of novel human spermatozoal RNAs. J Androl. 2005, 26 (1): 70-74.PubMedGoogle Scholar
- Charlesworth B, Sniegowski P, Stephan W: The evolutionary dynamics of repetitive DNA in eukaryotes. Nature. 1994, 371: 215-220. 10.1038/371215a0.View ArticlePubMedGoogle Scholar
- Jeffereys AJ, Royle NJ, Wilson V, Wong Z: Spontaneous mutation rates to new length alleles at tandem-repetitive hyper-variable loci in human DNA. Nature. 1998, 332 (6161): 278-281. 10.1038/332278a0.View ArticleGoogle Scholar
- Bennett P: Microsatellites. J Clin Pathol. 2000, 53: 177-183.View ArticleGoogle Scholar
- Jasinska A, Krzyzosiak WJ: Repetitive sequences that shape the human transcriptome. FEBS letters. 2004, 567: 136-141. 10.1016/j.febslet.2004.03.109.View ArticlePubMedGoogle Scholar
- Borstnik B, Pumpernik D: Tandem repeats in protein coding regions of primate genes. Genome Res. 2002, 12: 909-915. 10.1101/gr.138802.PubMed CentralView ArticlePubMedGoogle Scholar
- Rocha EP, Matic I, Taddei F: Over-representation of repeats in stress response genes: a strategy to increase versatility under stressful conditions?. Nucl Acids Res. 2002, 30: 1886-1894. 10.1093/nar/30.9.1886.PubMed CentralView ArticlePubMedGoogle Scholar
- Li YC, Korol AB, Fahima T, Nevo E: Microsatellites within genes: structure, function and evolution. Mol Biol Evol. 2004, 21: 991-1007. 10.1093/molbev/msh073.View ArticlePubMedGoogle Scholar
- Vergnaud G, Denoeud F: Minisatellites: Mutability and Genome architecture. Genome Res. 2000, 10: 899-907. 10.1101/gr.10.7.899.View ArticlePubMedGoogle Scholar
- Ali S, Wallace RB: Intrinsic polymorphism of variable number tandem repeat loci in the human genome. Nucl Acids Res. 1988, 16: 8487-8496. 10.1093/nar/16.17.8487.PubMed CentralView ArticlePubMedGoogle Scholar
- Jeffereys AJ, Wilson V Thein SL: Hypervariable 'minisatellite' regions in human DNA. Nature. 1985, 314 (6006): 67-73. 10.1038/314067a0.View ArticleGoogle Scholar
- Weitzel JN, Hows JM, Jeffreys AJ, Min GL, Goldman JM: Use of a hypervariable minisatellite DNA probe (33.15) for evaluating engraftment two or more years after bone marrow transplantation for aplastic anaemia. Br J Haematol. 1998, 70 (1): 91-97. 10.1111/j.1365-2141.1988.tb02439.x.View ArticleGoogle Scholar
- Bashamboo A, Ali S: Minisatellite associated sequence amplification (MASA) of the hypervariable repeat marker 33.15 reveals a male specific band in humans. Mol Cell Probes. 2001, 15: 89-92. 10.1006/mcpr.2001.0350.View ArticlePubMedGoogle Scholar
- Srivastava J, Premi S, Pathak D, Ahsan Z, Tiwari M, Garg LC, Ali S: Transcriptional Status of Known and Novel Genes Tagged with Consensus of 33.15 Repeat Loci Employing Minisatellite Associated Sequence Amplification (MASA) and Real Time PCR in Water Buffalo Bubalus bubalis. DNA Cell Biol. 2006, 25 (1): 31-48. 10.1089/dna.2006.25.31.View ArticlePubMedGoogle Scholar
- Toth G, Gaspari Z, Jurka J: Microsatellites in different eukaryotic genomes: survey and analysis. Genome Res. 2000, 10: 967-981. 10.1101/gr.10.7.967.PubMed CentralView ArticlePubMedGoogle Scholar
- Katti MV, Ranjekar PK, Gupta VS: Differential distribution of simple sequence repeats in eukaryotic genome sequences. Mol Biol Evol. 2001, 18 (7): 1161-1167.View ArticlePubMedGoogle Scholar
- Lalancette C, Thibault C, Bachand L, Caron N, Bissonnette N: Transcriptome Analysis of Bull Semen with Extreme Nonreturn Rate: Use of Suppression-Subtractive Hybridization to Identify Functional Markers for Fertility. Biol Reprod. 2008, 78: 618-635. 10.1095/biolreprod.106.059030.View ArticlePubMedGoogle Scholar
- Miller D, Tang PZ, Skinner C, Lilford R: Differential RNA fingerprinting as a tool in the analysis of spermatozoal gene expression. Hum Reprod. 1994, 9 (5): 864-869.PubMedGoogle Scholar
- Srivastava J, Premi S, Kumar S, Parwez I, Ali S: Characterization of Smoc-1 uncovers two transcript variants showing differential tissue and age specific expression in Bubalus bubalis. BMC Genomics. 2007, 8: 436-10.1186/1471-2164-8-436.PubMed CentralView ArticlePubMedGoogle Scholar
- Srivastava J, Premi S, Kumar S, Ali S: Organization and differential expression of the GACA/GATA tagged somatic and spermatozoal transcriptomes in buffalo Bubalus bubalis. BMC Genomics. 2008, 9: 132-10.1186/1471-2164-9-132.PubMed CentralView ArticlePubMedGoogle Scholar
- Srivastava J, Premi S, Garg LC, Ali S: Organizational and Expressional Uniqueness of a Testis Specific mRNA Transcript of Proto-oncogene c-kit Receptor in Water Buffalo Bubalus bubalis. DNA Cell Biol. 2006, 25 (9): 501-513. 10.1089/dna.2006.25.501.View ArticlePubMedGoogle 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.