Novel Insights Into Rheumatoid Arthritis Through Characterization of Concordant Changes in DNA Methylation and Gene Expression in Synovial Biopsies of Patients With Differing Numbers of Swollen Joints.

In this study, we sought to characterize synovial tissue obtained from individuals with arthralgia and disease-specific auto-antibodies and patients with established rheumatoid arthritis (RA), by applying an integrative multi-omics approach where we investigated differences at the level of DNA methylation and gene expression in relation to disease pathogenesis. We performed concurrent whole-genome bisulphite sequencing and RNA-Sequencing on synovial tissue obtained from the knee and ankle from 4 auto-antibody positive arthralgia patients and thirteen RA patients. Through multi-omics factor analysis we observed that the latent factor explaining the variance in gene expression and DNA methylation was associated with Swollen Joint Count 66 (SJC66), with patients with SJC66 of 9 or more displaying separation from the rest. Interrogating these observed differences revealed activation of the immune response as well as dysregulation of cell adhesion pathways at the level of both DNA methylation and gene expression. We observed differences for 59 genes in particular at the level of both transcript expression and DNA methylation. Our results highlight the utility of genome-wide multi-omics profiling of synovial samples for improved understanding of changes associated with disease spread in arthralgia and RA patients, and point to novel candidate targets for the treatment of the disease.

Single-Step, High-Efficiency CRISPR-Cas9 Genome Editing in Primary Human Disease-Derived Fibroblasts.

Genome editing is a tool that has many applications, including the validation of potential drug targets. However, performing genome editing in low-passage primary human cells with the greatest physiological relevance is notoriously difficult. High editing efficiency is desired because it enables gene knockouts (KO) to be generated in bulk cellular populations and circumvents the problem of having to generate clonal cell isolates. Here, we describe a single-step workflow enabling >90% KO generation in primary human lung fibroblasts via CRISPR ribonucleoprotein delivery in the absence of antibiotic selection or clonal expansion. As proof of concept, we edited two SMAD family members and demonstrated that in response to transforming growth factor beta, SMAD3, but not SMAD2, is critical for deposition of type I collagen in the fibrotic response. The optimization of this workflow can be readily transferred to other primary cell types.

Multiplexed DNA Methylation Analysis of Target Regions Using Microfluidics (Fluidigm).

Whole genome shotgun bisulfite sequencing is a method used to generate genome-wide methylation profiles. There are many available protocols to validate the results of this genome-wide method, but they mostly share the limitation of measuring methylation at a small number of CpG positions in small numbers of samples. We developed a multiplexed DNA methylation analysis protocol, which allows for the simultaneous quantitative measurement of cytosine methylation at single nucleotide resolution in 48 PCR amplicons and 48 samples utilizing the microfluidic system established by Fluidigm. Following bisulfite conversion of 500 ng of the target DNA, a PCR reaction is performed using a 48.48 Access Array, which allows parallel amplification of 48 samples by 48 primer pairs. The products of each reaction are labeled with individual, sample specific tags, pooled in a single library and sequenced using the Illumina MiSeq sequencer. The advantages of this system are: speed, small amount of input material, single nucleotide resolution, high coverage of each locus, low cost of simultaneously assaying multiple CpG loci in multiple DNA samples and high reproducibility.

Strategies for In Vivo Screening and Mitigation of Hepatotoxicity Associated with Antisense Drugs.

Antisense oligonucleotide (ASO) gapmers downregulate gene expression by inducing enzyme-dependent degradation of targeted RNA and represent a promising therapeutic platform for addressing previously undruggable genes. Unfortunately, their therapeutic application, particularly that of the more potent chemistries (e.g., locked-nucleic-acid-containing gapmers), has been hampered by their frequent hepatoxicity, which could be driven by hybridization-mediated interactions. An early de-risking of this liability is a crucial component of developing safe, ASO-based drugs. To rank ASOs based on their effect on the liver, we have developed an acute screen in the mouse that can be applied early in the drug development cycle. A single-dose (3-day) screen with streamlined endpoints (i.e., plasma transaminase levels and liver weights) was observed to be predictive of ASO hepatotoxicity ranking established based on a repeat-dose (15 day) study. Furthermore, to study the underlying mechanisms of liver toxicity, we applied transcriptome profiling and pathway analyses and show that adverse in vivo liver phenotypes correlate with the number of potent, hybridization-mediated off-target effects (OTEs). We propose that a combination of in silico OTE predictions, streamlined in vivo hepatotoxicity screening, and a transcriptome-wide selectivity screen is a valid approach to identifying and progressing safer compounds.

In silico and in vitro evaluation of exonic and intronic off-target effects form a critical element of therapeutic ASO gapmer optimization.

With many safety and technical limitations partly mitigated through chemical modifications, antisense oligonucleotides (ASOs) are gaining recognition as therapeutic entities. The increase in potency realized by 'third generation chemistries' may, however, simultaneously increase affinity to unintended targets with partial sequence complementarity. However, putative hybridization-dependent off-target effects (OTEs), a risk historically regarded as low, are not being adequately investigated. Here we show an unexpectedly high OTEs confirmation rate during screening of fully phosphorothioated (PS)-LNA gapmer ASOs designed against the BACH1 transcript. We demonstrate in vitro mRNA and protein knockdown of off-targets with a wide range of mismatch (MM) and gap patterns. Furthermore, with RNase H1 activity residing within the nucleus, hybridization predicted against intronic regions of pre-mRNAs was tested and confirmed. This dramatically increased ASO-binding landscape together with relatively high potency of such interactions translates into a considerable safety concern. We show here that with base pairing-driven target recognition it is possible to predict the putative off-targets and address the liability during lead design and optimization phases. Moreover, in silico analysis performed against both primary as well as spliced transcripts will be invaluable in elucidating the mechanism behind the hepatoxicity observed with some LNA-modified gapmers.

The zinc transporter ZIP12 regulates the pulmonary vascular response to chronic hypoxia.

The typical response of the adult mammalian pulmonary circulation to a low oxygen environment is vasoconstriction and structural remodelling of pulmonary arterioles, leading to chronic elevation of pulmonary artery pressure (pulmonary hypertension) and right ventricular hypertrophy. Some mammals, however, exhibit genetic resistance to hypoxia-induced pulmonary hypertension. We used a congenic breeding program and comparative genomics to exploit this variation in the rat and identified the gene Slc39a12 as a major regulator of hypoxia-induced pulmonary vascular remodelling. Slc39a12 encodes the zinc transporter ZIP12. Here we report that ZIP12 expression is increased in many cell types, including endothelial, smooth muscle and interstitial cells, in the remodelled pulmonary arterioles of rats, cows and humans susceptible to hypoxia-induced pulmonary hypertension. We show that ZIP12 expression in pulmonary vascular smooth muscle cells is hypoxia dependent and that targeted inhibition of ZIP12 inhibits the rise in intracellular labile zinc in hypoxia-exposed pulmonary vascular smooth muscle cells and their proliferation in culture. We demonstrate that genetic disruption of ZIP12 expression attenuates the development of pulmonary hypertension in rats housed in a hypoxic atmosphere. This new and unexpected insight into the fundamental role of a zinc transporter in mammalian pulmonary vascular homeostasis suggests a new drug target for the pharmacological management of pulmonary hypertension.

Genetic analysis of the cardiac methylome at single nucleotide resolution in a model of human cardiovascular disease.

Epigenetic marks such as cytosine methylation are important determinants of cellular and whole-body phenotypes. However, the extent of, and reasons for inter-individual differences in cytosine methylation, and their association with phenotypic variation are poorly characterised. Here we present the first genome-wide study of cytosine methylation at single-nucleotide resolution in an animal model of human disease. We used whole-genome bisulfite sequencing in the spontaneously hypertensive rat (SHR), a model of cardiovascular disease, and the Brown Norway (BN) control strain, to define the genetic architecture of cytosine methylation in the mammalian heart and to test for association between methylation and pathophysiological phenotypes. Analysis of 10.6 million CpG dinucleotides identified 77,088 CpGs that were differentially methylated between the strains. In F1 hybrids we found 38,152 CpGs showing allele-specific methylation and 145 regions with parent-of-origin effects on methylation. Cis-linkage explained almost 60% of inter-strain variation in methylation at a subset of loci tested for linkage in a panel of recombinant inbred (RI) strains. Methylation analysis in isolated cardiomyocytes showed that in the majority of cases methylation differences in cardiomyocytes and non-cardiomyocytes were strain-dependent, confirming a strong genetic component for cytosine methylation. We observed preferential nucleotide usage associated with increased and decreased methylation that is remarkably conserved across species, suggesting a common mechanism for germline control of inter-individual variation in CpG methylation. In the RI strain panel, we found significant correlation of CpG methylation and levels of serum chromogranin B (CgB), a proposed biomarker of heart failure, which is evidence for a link between germline DNA sequence variation, CpG methylation differences and pathophysiological phenotypes in the SHR strain. Together, these results will stimulate further investigation of the molecular basis of locally regulated variation in CpG methylation and provide a starting point for understanding the relationship between the genetic control of CpG methylation and disease phenotypes.

Natural polymorphisms in Tap2 influence negative selection and CD4∶CD8 lineage commitment in the rat.

Genetic variation in the major histocompatibility complex (MHC) affects CD4∶CD8 lineage commitment and MHC expression. However, the contribution of specific genes in this gene-dense region has not yet been resolved. Nor has it been established whether the same genes regulate MHC expression and T cell selection. Here, we assessed the impact of natural genetic variation on MHC expression and CD4∶CD8 lineage commitment using two genetic models in the rat. First, we mapped Quantitative Trait Loci (QTLs) associated with variation in MHC class I and II protein expression and the CD4∶CD8 T cell ratio in outbred Heterogeneous Stock rats. We identified 10 QTLs across the genome and found that QTLs for the individual traits colocalized within a region spanning the MHC. To identify the genes underlying these overlapping QTLs, we generated a large panel of MHC-recombinant congenic strains, and refined the QTLs to two adjacent intervals of ∼0.25 Mb in the MHC-I and II regions, respectively. An interaction between these intervals affected MHC class I expression as well as negative selection and lineage commitment of CD8 single-positive (SP) thymocytes. We mapped this effect to the transporter associated with antigen processing 2 (Tap2) in the MHC-II region and the classical MHC class I gene(s) (RT1-A) in the MHC-I region. This interaction was revealed by a recombination between RT1-A and Tap2, which occurred in 0.2% of the rats. Variants of Tap2 have previously been shown to influence the antigenicity of MHC class I molecules by altering the MHC class I ligandome. Our results show that a restricted peptide repertoire on MHC class I molecules leads to reduced negative selection of CD8SP cells. To our knowledge, this is the first study showing how a recombination between natural alleles of genes in the MHC influences lineage commitment of T cells.

Genome sequencing reveals loci under artificial selection that underlie disease phenotypes in the laboratory rat.

Large numbers of inbred laboratory rat strains have been developed for a range of complex disease phenotypes. To gain insights into the evolutionary pressures underlying selection for these phenotypes, we sequenced the genomes of 27 rat strains, including 11 models of hypertension, diabetes, and insulin resistance, along with their respective control strains. Altogether, we identified more than 13 million single-nucleotide variants, indels, and structural variants across these rat strains. Analysis of strain-specific selective sweeps and gene clusters implicated genes and pathways involved in cation transport, angiotensin production, and regulators of oxidative stress in the development of cardiovascular disease phenotypes in rats. Many of the rat loci that we identified overlap with previously mapped loci for related traits in humans, indicating the presence of shared pathways underlying these phenotypes in rats and humans. These data represent a step change in resources available for evolutionary analysis of complex traits in disease models.

Isolation and high-throughput sequencing of two closely linked epistatic hypertension susceptibility loci with a panel of bicongenic strains.

Interactions or epistasis between genetic factors may contribute to "missing heritability." While linkage analyses detect epistasis, defining the limits of the interacting segments poses a significant challenge especially when the interactions are between loci in close proximity. The goal of the present study was to isolate two such epistatic blood pressure (BP) loci on rat chromosome 5. A panel of S.LEW bicongenic strains along with the corresponding monocongenic strains was constructed. BP of each set comprising of one bicongenic and two corresponding monocongenic strains were determined along with the parental Salt-sensitive (S) strain. Epistasis was observed in one out of four sets of congenic strains, wherein systolic blood pressures (SBP) of the two monocongenic strains S.LEW(5)x6Bx9x5a and S.LEW(5)x6Bx9x5b were comparable to that of S, but the SBP of the bicongenic strain S.LEW(5)x6Bx9x5 (157 ± 4.3 mmHg) was significantly lower than that of S (196 ± 6.8 mmHg, P < 0.001). A two-way ANOVA indicated significant interactions between the LEW alleles at the two loci. The interacting loci were 2.02 Mb apart and located within genomic segments spanning 7.77 and 4.18 Mb containing 7,360 and 2,753 candidate variants, respectively. The current study demonstrates definitive evidence for epistasis and provides genetic tools for further dissection of the isolated epistatic BP loci.

Alternative splicing and transcriptome profiling of experimental autoimmune encephalomyelitis using genome-wide exon arrays.

BACKGROUND: Multiple Sclerosis (MS) is a chronic inflammatory disease causing demyelination and nerve loss in the central nervous system. Experimental autoimmune encephalomyelitis (EAE) is an animal model of MS that is widely used to investigate complex pathogenic mechanisms. Transcriptional control through isoform selection and mRNA levels determines pathway activation and ultimately susceptibility to disease. METHODOLOGY/PRINCIPAL FINDINGS: We have studied the role of alternative splicing and differential expression in lymph node cells from EAE-susceptible Dark Agouti (DA) and EAE-resistant Piebald Virol Glaxo.AV1 (PVG) inbred rat strains using Affymetrix Gene Chip Rat Exon 1.0 ST Arrays. Comparing the two strains, we identified 11 differentially spliced and 206 differentially expressed genes at day 7 post-immunization, as well as 9 differentially spliced and 144 differentially expressed genes upon autoantigen re-stimulation. Functional clustering and pathway analysis implicate genes for glycosylation, lymphocyte activation, potassium channel activity and cellular differentiation in EAE susceptibility. CONCLUSIONS/SIGNIFICANCE: Our results demonstrate that alternative splicing occurs during complex disease and may govern EAE susceptibility. Additionally, transcriptome analysis not only identified previously defined EAE pathways regulating the immune system, but also novel mechanisms. Furthermore, several identified genes overlap known quantitative trait loci, providing novel causative candidate targets governing EAE.

Angiotensin-converting enzyme is a modifier of hypertensive end organ damage.

Severe forms of hypertension are characterized by high blood pressure combined with end organ damage. Through the development and refinement of a transgenic rat model of malignant hypertension incorporating the mouse renin gene, we previously identified a quantitative trait locus on chromosome 10, which affects malignant hypertension severity and morbidity. We next generated an inducible malignant hypertensive model where the timing, severity, and duration of hypertension was placed under the control of the researcher, allowing development of and recovery from end organ damage to be investigated. We have now generated novel consomic Lewis and Fischer rat strains with inducible hypertension and additional strains that are reciprocally congenic for the refined chromosome 10 quantitative trait locus. We have captured a modifier of end organ damage within the congenic region and, using a range of bioinformatic, biochemical and molecular biological techniques, have identified angiotensin-converting enzyme as the modifier of hypertension-induced tissue microvascular injury. Reciprocal differences between angiotensin-converting enzyme and the anti-inflammatory tetrapeptide, N-acetyl-Ser-Asp-Lys-Pro in the kidney, a tissue susceptible to end organ damage, suggest a mechanism for the amelioration of hypertension-dependent damage.

Comparison of dorsal root ganglion gene expression in rat models of traumatic and HIV-associated neuropathic pain.

To elucidate the mechanisms underlying peripheral neuropathic pain in the context of HIV infection and antiretroviral therapy, we measured gene expression in dorsal root ganglia (DRG) of rats subjected to systemic treatment with the anti-retroviral agent, ddC (Zalcitabine) and concomitant delivery of HIV-gp120 to the rat sciatic nerve. L4 and L5 DRGs were collected at day 14 (time of peak behavioural change) and changes in gene expression were measured using Affymetrix whole genome rat arrays. Conventional analysis of this data set and Gene Set Enrichment Analysis (GSEA) was performed to discover biological processes altered in this model. Transcripts associated with G protein coupled receptor signalling and cell adhesion were enriched in the treated animals, while ribosomal proteins and proteasome pathways were associated with gene down-regulation. To identify genes that are directly relevant to neuropathic mechanical hypersensitivity, as opposed to epiphenomena associated with other aspects of the response to a sciatic nerve lesion, we compared the gp120+ddC-evoked gene expression with that observed in a model of traumatic neuropathic pain (L5 spinal nerve transection), where hypersensitivity to a static mechanical stimulus is also observed. We identified 39 genes/expressed sequence tags that are differentially expressed in the same direction in both models. Most of these have not previously been implicated in mechanical hypersensitivity and may represent novel targets for therapeutic intervention. As an external control, the RNA expression of three genes was examined by RT-PCR, while the protein levels of two were studied using western blot analysis.

MiMiR--an integrated platform for microarray data sharing, mining and analysis.

BACKGROUND: Despite considerable efforts within the microarray community for standardising data format, content and description, microarray technologies present major challenges in managing, sharing, analysing and re-using the large amount of data generated locally or internationally. Additionally, it is recognised that inconsistent and low quality experimental annotation in public data repositories significantly compromises the re-use of microarray data for meta-analysis. MiMiR, the Microarray data Mining Resource was designed to tackle some of these limitations and challenges. Here we present new software components and enhancements to the original infrastructure that increase accessibility, utility and opportunities for large scale mining of experimental and clinical data. RESULTS: A user friendly Online Annotation Tool allows researchers to submit detailed experimental information via the web at the time of data generation rather than at the time of publication. This ensures the easy access and high accuracy of meta-data collected. Experiments are programmatically built in the MiMiR database from the submitted information and details are systematically curated and further annotated by a team of trained annotators using a new Curation and Annotation Tool. Clinical information can be annotated and coded with a clinical Data Mapping Tool within an appropriate ethical framework. Users can visualise experimental annotation, assess data quality, download and share data via a web-based experiment browser called MiMiR Online. All requests to access data in MiMiR are routed through a sophisticated middleware security layer thereby allowing secure data access and sharing amongst MiMiR registered users prior to publication. Data in MiMiR can be mined and analysed using the integrated EMAAS open source analysis web portal or via export of data and meta-data into Rosetta Resolver data analysis package. CONCLUSION: The new MiMiR suite of software enables systematic and effective capture of extensive experimental and clinical information with the highest MIAME score, and secure data sharing prior to publication. MiMiR currently contains more than 150 experiments corresponding to over 3000 hybridisations and supports the Microarray Centre's large microarray user community and two international consortia. The MiMiR flexible and scalable hardware and software architecture enables secure warehousing of thousands of datasets, including clinical studies, from microarray and potentially other -omics technologies.

Mouse MAELSTROM: the link between meiotic silencing of unsynapsed chromatin and microRNA pathway?

Meiotic silencing of unsynapsed chromatin (MSUC) is a key mechanism in spermatogenesis and a model system to study the dynamics of gene silencing. Here we show that MAEL, the ortholog of Drosophila's high mobility group box protein Maelstrom, is associated not only with the silenced XY body, but also with unsynapsed autosomes. Characterization of MAEL revealed that it interacts directly with the chromatin remodeler SNF5/INI1 and chromatin-associated protein SIN3B, which we also find localized to the XY body. This is the first time that a chromatin remodeler has been shown to associate with whole chromosomes. In addition, we show that MAEL is a component of the mouse meiotic nuage and its haploid cell counterpart, the chromatoid body. This is a site of accumulation of RNA and RNA processing enzymes, including proteins involved in the microRNA (miRNA) pathway. Furthermore, in the nuage, MAEL is present in a complex with germ cell specific MVH, an RNA helicase and Argonaute family members, MILI and MIWI. The presence of MAEL in these critical compartments of male germ cells and its interactions provide a link suggesting the involvement of the miRNA pathway in MSUC.

Dazl binds in vivo to specific transcripts and can regulate the pre-meiotic translation of Mvh in germ cells.

Gametogenesis is a complex process subject to strict controls at both levels of transcription and translation. Members of a family of conserved RNA-binding proteins encoded by the DAZ genes are required for the translational regulation of gene expression essential for this process. Although loss of DAZ family genes is associated with infertility in several organisms including humans, the identity of the transcripts regulated in vivo is unknown. Using a combination of immunoprecipitation and microarray analysis, we have identified a number of mRNAs that are bound by the murine Dazl protein both in vivo and in vitro. Sequence analysis shows that these transcripts contain binding sites for Dazl, which have been conserved during evolution between human, rat and mouse. We have focussed on mouse vasa homologue (Mvh), a gene that is essential for male gametogenesis, and show that Dazl stimulates translation via the Mvh 3'-UTR. Finally, we show that germ cells of Dazl null mice contain reduced levels of Mvh protein, indicating that Dazl-mediated regulation of Mvh translation is crucial for mammalian spermatogenesis.

Two novel proteins recruited by synaptonemal complex protein 1 (SYCP1) are at the centre of meiosis.

Completion of meiosis in mammals depends on the formation of the synaptonemal complex, a tripartite structure that physically links homologous chromosomes during prophase I. Several components of the synaptonemal complex are known, including constituents of the cohesin core, the axial/lateral element and the transverse filaments. No protein has previously been identified as an exclusive component of the central element. Mutations in some synaptonemal-complex proteins results in impaired meiosis. In humans, cases of male infertility have been associated with failure to build the synaptonemal complex. To search for new components of the meiotic machinery, we have used data from microarray expression profiling and found two proteins localising solely to the central element of the mammalian synaptonemal complex. These new proteins, SYCE1 and CESC1, interact with the transverse filament protein SYCP1, and their localisation to the central element appears to depend on recruitment by SYCP1. This suggests a role for SYCE1 and CESC1 in synaptonemal-complex assembly, and perhaps also stability and recombination.

Triple-target microarray experiments: a novel experimental strategy.

BACKGROUND: High-throughput, parallel gene expression analysis by means of microarray technology has become a widely used technique in recent years. There are currently two main dye-labelling strategies for microarray studies based on custom-spotted cDNA or oligonucleotides arrays: (I) Dye-labelling of a single target sample with a particular dye, followed by subsequent hybridisation to a single microarray slide, (II) Dye-labelling of two different target samples with two different dyes, followed by subsequent co-hybridisation to a single microarray slide. The two dyes most frequently used for either method are Cy3 and Cy5. We propose and evaluate a novel experiment set-up utilising three differently labelled targets co-hybridised to one microarray slide. In addition to Cy3 and Cy5, this incorporates Alexa 594 as a third dye-label. We evaluate this approach in line with current data processing and analysis techniques for microarrays, and run separate analyses on Alexa 594 used in single-target, dual-target and the intended triple-target experiment set-ups (a total of 18 microarray slides). We follow this by pointing out practical applications and suitable analysis methods, and conclude that triple-target microarray experiments can add value to microarray research by reducing material costs for arrays and related processes, and by increasing the number of options for pragmatic experiment design. RESULTS: The addition of Alexa 594 as a dye-label for an additional--third--target sample works within the framework of more commonplace Cy5/Cy3 labelled target sample combinations. Standard normalisation methods are still applicable, and the resulting data can be expected to allow identification of expression differences in a biological experiment, given sufficient levels of biological replication (as is necessary for most microarray experiments). CONCLUSION: The use of three dye-labelled target samples can be a valuable addition to the standard repertoire of microarray experiment designs. The method enables direct comparison between two experimental populations as well as measuring these two populations in relation to a third reference sample, allowing comparisons within the slide and across slides. These benefits are only offset by the added level of consideration required in the experimental design and data processing of a triple-target study design. Common methods for data processing and analysis are still applicable, but there is scope for the development of custom models for triple-target data. In summary, we do not consider the triple-target approach to be a new standard, but a valuable addition to the existing microarray study toolkit.

Expression profiling of the developing testis in wild-type and Dazl knockout mice.

Genetic understanding of male-factor infertility requires knowledge of gene expression patterns associated with normal germ cell differentiation. The mouse is one of the best models of mammalian fertility due to its well-characterized genetics and the existence of many infertile mutants both naturally occurring and experimentally induced. We used cDNA microarrays firstly to investigate normal gene expression in the wild-type (wt) testis and secondly to gain a better insight into the effect of the disruption of the Dazl gene on spermatogenesis. We constructed a cDNA microarray from a subtracted and normalized adult testis library and focused on six developmental time-points during the initial synchronous wave of spermatogenesis. The results suggest that in the wild-type testis, 89.5% of genes on our chip change expression dramatically during the time-course. To identify patterns in the gene-expression data, a k-means clustering algorithm and principal component analysis were used. In the Dazl knockout testes, the majority of genes remain at baseline levels of expression, because absence of Dazl has a severe effect on cell-types present in the testis. Although in the prepubescent Dazl-null mice the final point reached in germ cell development is the leptotene-zygotene stage, the microarray results suggest that lack of Dazl expression has a detectable effect on the mRNA complement of germ cells as early as day 5 when only type A spermatogonia are present. Mol. Reprod. Dev. 67: 26-54, 2004.