Deleterious alleles in the human genome are on average younger than neutral alleles of the same frequency.

Large-scale population sequencing studies provide a complete picture of human genetic variation within the studied populations. A key challenge is to identify, among the myriad alleles, those variants that have an effect on molecular function, phenotypes, and reproductive fitness. Most non-neutral variation consists of deleterious alleles segregating at low population frequency due to incessant mutation. To date, studies characterizing selection against deleterious alleles have been based on allele frequency (testing for a relative excess of rare alleles) or ratio of polymorphism to divergence (testing for a relative increase in the number of polymorphic alleles). Here, starting from Maruyama's theoretical prediction (Maruyama T (1974), Am J Hum Genet USA 6:669-673) that a (slightly) deleterious allele is, on average, younger than a neutral allele segregating at the same frequency, we devised an approach to characterize selection based on allelic age. Unlike existing methods, it compares sets of neutral and deleterious sequence variants at the same allele frequency. When applied to human sequence data from the Genome of the Netherlands Project, our approach distinguishes low-frequency coding non-synonymous variants from synonymous and non-coding variants at the same allele frequency and discriminates between sets of variants independently predicted to be benign or damaging for protein structure and function. The results confirm the abundance of slightly deleterious coding variation in humans.

Huntingtin with an expanded polyglutamine repeat affects the Jab1-p27(Kip1) pathway.

Expansion of polyglutamine repeats is the cause of at least nine inherited human neurodegenerative disorders, including Huntington's disease (HD). It is widely accepted that deregulation of the transcriptional coactivator CBP by expanded huntingtin (htt) plays an important role in HD molecular pathogenesis. In this study, we report on a novel target of expanded polyglutamine stretches, the transcriptional coactivator Jun activation domain-binding protein 1 (Jab1), which shares DNA-sequence-specific transcription factor targets with CBP. Jab1 also plays a major role in the degradation of the cyclin-dependent-kinase inhibitor and putative transcription cofactor p27(Kip1). We found that Jab1 accumulates in aggregates when co-expressed with either expanded polyglutamine stretches or N-terminal fragments of mutant htt. In addition, the coactivator function of Jab1 was suppressed both by aggregated expanded polyglutamine solely and by mutant htt. Inhibition by mutant htt even preceded the appearance of microscopic aggregation. In an exon 1 HD cell model, we found that endogenous Jab1 could be recruited into aggregates and that this was accompanied by the accumulation of p27(Kip1). Accumulation of p27(Kip1) was also found in brains derived from HD patients. The repression of Jab1 by various mechanisms coupled with an increase of p27(Kip1) at late stages may have important transcriptional effects. In addition, the interference with the Jab1-p27(Kip1) pathway may contribute to the observed lower incidence of cancer in HD patients and may also be relevant for the understanding of the molecular pathogenesis of polyglutamine disorders in general.

Serum matrix metalloproteinase-9 (MMP-9) as a biomarker for monitoring disease progression in Duchenne muscular dystrophy (DMD).

To identify serum biomarkers that allow monitoring of disease progression and treatment effects in Duchenne muscular dystrophy (DMD) patients, levels of matrix metalloproteinase-9 (MMP-9), tissue inhibitors of metalloproteinase-1 (TIMP-1) and osteopontin (OPN) were determined in 63 DMD patients on corticosteroid therapy. These proteins were selected for their role in the pathogenesis of muscular dystrophy. Levels of MMP-9 and TIMP-1 were significantly higher in sera of DMD patients compared to healthy controls, whereas the OPN levels showed no significant difference. MMP-9 levels were also observed to be significantly higher in older, nonambulant patients, compared to ambulant patients. Longitudinal data from a smaller cohort of DMD patients followed up for over 4years showed that MMP-9, but not TIMP-1 increased significantly with age. Hence, MMP-9 is a potential DMD biomarker for disease progression. Future studies have to confirm whether serum MMP-9 levels can be used to monitor therapeutic response.

Normal and mutant HTT interact to affect clinical severity and progression in Huntington disease.

OBJECTIVE: Huntington disease (HD) is an autosomal dominant neurodegenerative disorder caused by a CAG repeat expansion in the HD gene (HTT). We aimed to assess whether interaction between CAG repeat sizes in the mutant and normal allele could affect disease severity and progression. METHODS: Using linear regression and mixed-effects models, the influence of mutant and normal CAG repeat sizes interaction was assessed on 1) age at onset in 921 patients with HD, 2) clinical severity and progression in 512 of these patients with follow-up data available, and 3) basal ganglia volume on magnetic resonance images in 16 premanifest HD mutation carriers. RESULTS: Normal and mutant CAG repeat sizes interacted to influence 1) age at onset (p = 0.001), 2) severity or progression of motor, cognitive, and functional, but not behavioral, symptoms in patients with HD (all p < 0.05), and 3) in premanifest subjects, basal ganglia volumes (p < 0.05). In subjects with mutant CAG expansions in the low range, increasing size of the normal repeat correlated with more severe symptoms and pathology, whereas for those subjects with expansions in the high range, increasing size of the normal repeat correlated with less severe symptoms and pathology. CONCLUSIONS: Increasing CAG repeat size in normal HTT diminishes the association between mutant CAG repeat size and disease severity and progression in Huntington disease. The underlying mechanism may involve interaction of the polyglutamine domains of normal and mutant huntingtin (fragments) and needs further elucidation. These findings may have predictive value and are essential for the design and interpretation of future therapeutic trials.

Relative power and sample size analysis on gene expression profiling data.

BACKGROUND: With the increasing number of expression profiling technologies, researchers today are confronted with choosing the technology that has sufficient power with minimal sample size, in order to reduce cost and time. These depend on data variability, partly determined by sample type, preparation and processing. Objective measures that help experimental design, given own pilot data, are thus fundamental. RESULTS: Relative power and sample size analysis were performed on two distinct data sets. The first set consisted of Affymetrix array data derived from a nutrigenomics experiment in which weak, intermediate and strong PPARalpha agonists were administered to wild-type and PPARalpha-null mice. Our analysis confirms the hierarchy of PPARalpha-activating compounds previously reported and the general idea that larger effect sizes positively contribute to the average power of the experiment. A simulation experiment was performed that mimicked the effect sizes seen in the first data set. The relative power was predicted but the estimates were slightly conservative. The second, more challenging, data set describes a microarray platform comparison study using hippocampal deltaC-doublecortin-like kinase transgenic mice that were compared to wild-type mice, which was combined with results from Solexa/Illumina deep sequencing runs. As expected, the choice of technology greatly influences the performance of the experiment. Solexa/Illumina deep sequencing has the highest overall power followed by the microarray platforms Agilent and Affymetrix. Interestingly, Solexa/Illumina deep sequencing displays comparable power across all intensity ranges, in contrast with microarray platforms that have decreased power in the low intensity range due to background noise. This means that deep sequencing technology is especially more powerful in detecting differences in the low intensity range, compared to microarray platforms. CONCLUSION: Power and sample size analysis based on pilot data give valuable information on the performance of the experiment and can thereby guide further decisions on experimental design. Solexa/Illumina deep sequencing is the technology of choice if interest lies in genes expressed in the low-intensity range. Researchers can get guidance on experimental design using our approach on their own pilot data implemented as a BioConductor package, SSPA http://bioconductor.org/packages/release/bioc/html/SSPA.html.

Limb-girdle muscular dystrophy in the Netherlands: gene defect identified in half the families.

Pheno- and genotype correlation is attempted in a Dutch cross-sectional study on limb- girdle muscular dystrophy. Sarcoglycans, caveolin-3, calpain-3, and dysferlin were analyzed on muscle tissue. Mutation analysis of the calpain-3, caveolin-3, and fukutin-related protein gene was executed in successive order for all samples. In 51% of all families a classifying diagnosis was made. Several new mutations in LGMD2A, B, and C patients have been found in this population.

Therapeutic modulation of DMD splicing by blocking exonic splicing enhancer sites with antisense oligonucleotides.

Antisense oligonucleotides (AONs) can be used to correct the disrupted reading frame of Duchenne muscular dystophy patients (DMD). We have a collection of 121 AONs, of which 79 are effective in inducing the specific skipping of 38 out of the 79 different DMD exons. All AONs are located within exons and were hypothesized to act by steric hindrance of serine-arginine rich (SR) protein binding to exonic splicing enhancer (ESE) sites. Indeed, retrospective in silico analysis of effective versus ineffective AONs revealed that the efficacy of AONs is correlated to the presence of putative ESE sites (as predicted by the ESEfinder and RESCUE-ESE software). ESE predicting software programs are thus valuable tools for the optimization of exon-internal antisense target sequences.

Common pathological mechanisms in mouse models for muscular dystrophies.

Duchenne/Becker and limb-girdle muscular dystrophies share clinical symptoms like muscle weakness and wasting but differ in clinical presentation and severity. To get a closer view on the differentiating molecular events responsible for the muscular dystrophies, we have carried out a comparative gene expression profiling of hindlimb muscles of the following mouse models: dystrophin-deficient (mdx, mdx(3cv)), sarcoglycan-deficient (Sgca null, Sgcb null, Sgcg null, Sgcd null), dysferlin-deficient (Dysf null, SJL(Dysf)), sarcospan-deficient (Sspn null), and wild-type (C57Bl/6, C57Bl/10) mice. The expression profiles clearly discriminated between severely affected (dystrophinopathies and sarcoglycanopathies) and mildly or nonaffected models (dysferlinopathies, sarcospan-deficiency, wild-type). Dystrophin-deficient and sarcoglycan-deficient profiles were remarkably similar, sharing inflammatory and structural remodeling processes. These processes were also ongoing in dysferlin-deficient animals, albeit at lower levels, in agreement with the later age of onset of this muscular dystrophy. The inflammatory proteins Spp1 and S100a9 were up-regulated in all models, including sarcospan-deficient mice, which points, for the first time, at a subtle phenotype for Sspn null mice. In conclusion, we identified biomarker genes for which expression correlates with the severity of the disease, which can be used for monitoring disease progression. This comparative study is an integrating step toward the development of an expression profiling-based diagnostic approach for muscular dystrophies in humans.

Muscle regeneration in dystrophin-deficient mdx mice studied by gene expression profiling.

BACKGROUND: Duchenne muscular dystrophy (DMD), caused by mutations in the dystrophin gene, is lethal. In contrast, dystrophin-deficient mdx mice recover due to effective regeneration of affected muscle tissue. To characterize the molecular processes associated with regeneration, we compared gene expression levels in hindlimb muscle tissue of mdx and control mice at 9 timepoints, ranging from 1-20 weeks of age. RESULTS: Out of 7776 genes, 1735 were differentially expressed between mdx and control muscle at at least one timepoint (p < 0.05 after Bonferroni correction). We found that genes coding for components of the dystrophin-associated glycoprotein complex are generally downregulated in the mdx mouse. Based on functional characteristics such as membrane localization, signal transduction, and transcriptional activation, 166 differentially expressed genes with possible functions in regeneration were analyzed in more detail. The majority of these genes peak at the age of 8 weeks, where the regeneration activity is maximal. The following pathways are activated, as shown by upregulation of multiple members per signalling pathway: the Notch-Delta pathway that plays a role in the activation of satellite cells, and the Bmp15 and Neuregulin 3 signalling pathways that may regulate proliferation and differentiation of satellite cells. In DMD patients, only few of the identified regeneration-associated genes were found activated, indicating less efficient regeneration processes in humans. CONCLUSION: Based on the observed expression profiles, we describe a model for muscle regeneration in mdx mice, which may provide new leads for development of DMD therapies based on the improvement of muscle regeneration efficacy.

Deletion of the Caenorhabditis elegans homologues of the CLN3 gene, involved in human juvenile neuronal ceroid lipofuscinosis, causes a mild progeric phenotype.

The CLN3 gene is involved in juvenile neuronal ceroid lipofuscinosis (JNCL), or Batten-Spielmeyer-Vogt disease, a severe hereditary neurodegenerative lysosomal storage disorder characterized by progressive disease pathology, with loss of vision as the first symptom. Another characteristic of JNCL is the lysosomal accumulation of autofluorescent lipopigments, forming fingerprint storage patterns visible by electron microscopy. The function of the CLN3 protein is still unknown, although the evolutionarily conserved CLN3 protein is being functionally analysed using different experimental models. We have explored the potential of the nematode Caenorhabditis elegans as a model for Batten disease in order to bridge the gap between the unicellular yeast and very complex mouse JNCL models. C. elegans has three genes homologous to CLN3, for each of which deletion mutants were isolated. Cln-3.1 deletion mutants have a decreased lifespan, and cln-3.2 deletion mutants a decreased brood size. However, the neuronal or movement defects and aberrant lipopigment distribution or accumulation observed in JNCL were not found in the worms. To detect possible redundancy, single deletion mutants were crossed to obtain double and triple mutants, which were viable but showed no JNCL-specific defects. The cln-3 triple mutants show a more prominent decrease in lifespan and brood size, the latter most conspicuously at the end of the egg-laying period, suggesting premature ageing. To focus our functional analysis we examined the C. elegans cln-3 expression patterns, using promoter-GFP (green fluorescent protein) gene fusions. Fluorescence patterns suggest cln-3.1 expression in the intestine, cln-3.2 expression in the hypoderm, and cln-3.3 expression in intestinal muscle, male-specific posterior muscle and hypoderm. Further life stage- and tissue-specific analysis of the processes causing the phenotype of the cln-3 triple mutants may provide more information about the function of the cln-3 protein and contribute to a better understanding of the basic processes affected in Batten disease patients.

Comparative analysis of antisense oligonucleotide analogs for targeted DMD exon 46 skipping in muscle cells.

As small molecule drugs for Duchenne muscular dystrophy (DMD), antisense oligonucleotides (AONs) have been shown to restore the disrupted reading frame of DMD transcripts by inducing specific exon skipping. This allows the synthesis of largely functional Becker muscular dystrophy (BMD)-like dystrophins and potential conversion of severe DMD into milder BMD phenotypes. Thus far we have used 2'-O-methyl phosphorothioate (2OMePS) AONs. Here, we assessed the skipping efficiencies of different AON analogs containing morpholino-phosphorodiamidate, locked nucleic acid (LNA) or peptide nucleic acid (PNA) backbones. In contrast to PNAs and morpholinos, LNAs have not yet been tested as splice modulators. Compared to the most effective 2OMePS AON directed at exon 46, the LNA induced higher skipping levels in myotubes from a human control (85 versus 20%) and an exon 45 deletion DMD patient (98 versus 75%). The morpholino-induced skipping levels were only 5-6%, whereas the PNA appeared to be ineffective. Further comparative analysis of LNA and 2OMePS AONs containing up to three mismatches revealed that LNAs, while inducing higher skipping efficiencies, show much less sequence specificity. This limitation increases the risk of adverse effects elsewhere in the human genome. Awaiting further improvements in oligochemistry, we thus consider 2OMePS AONs currently the most favorable compounds, at least for targeted DMD exon 46 skipping.

Genomic imbalances in mental retardation.

INTRODUCTION: It has been estimated that cytogenetically visible rearrangements are present in approximately 1% of newborns. These chromosomal changes can cause a wide range of deleterious developmental effects, including mental retardation (MR). It is assumed that many other cases exist where the cause is a submicroscopic deletion or duplication. To facilitate the detection of such cases, different techniques have been developed, which have differing efficiency as to the number of loci and patients that can be tested. METHODS: We implemented multiplex amplifiable probe hybridisation (MAPH) to test areas known to be rearranged in MR patients (for example, subtelomeric/pericentromeric regions and those affected in microdeletion syndromes) and to look for new regions that might be related to MR. RESULTS: In this study, over 30 000 screens for duplications and deletions were carried out; 162 different loci tested in each of 188 developmentally delayed patients. The analysis resulted in the detection of 19 rearrangements, of which approximately 65% would not have been detected by conventional cytogenetic analysis. A significant fraction (46%) of the rearrangements found were interstitial, despite the fact that only a limited number of these loci have so far been tested. DISCUSSION: Our results strengthen the arguments for whole genome screening within this population, as it can be assumed that many more interstitial rearrangements would be detected. The strengths of MAPH for this analysis are the simplicity, the high throughput potential, and the high resolution of analysis. This combination should help in the future identification of the specific genes that are responsible for MR.

Fluorescent labelling of cRNA for microarray applications.

Microarrays of oligonucleotide expression libraries can be hybridised with either cDNA, generated from mRNA during reverse transcription, or cRNA, generated in an Eberwine mRNA amplification procedure. While methods for fluorescent labelling of cDNA have been thoroughly investigated, methods for cRNA labelling have not. To this purpose, we developed an aminoallyl-UTP (aa-UTP) driven cRNA labelling protocol and compared it in expression profiling studies using spotted 7.5 K 65mer murine oligonucleotide arrays with labelling via direct incorporation of Cy-UTPs. The presence of dimethylsulfoxide during coupling of aa-modified cRNA with N-hydroxysuccinimide-modified, fluorescent Cy dyes greatly enhanced the labelling efficiency, as analysed by spectrophotometry and fluorescent hybridisation signals. Indirect labelling using aa-UTP resulted in 2- to 3-fold higher degrees of labelling and fluorescent signals than labelling by direct incorporation of Cy-UTP. By variation of the aa-UTP:UTP ratio, a clear optimal degree of labelling was found (1 dye per 20-25 nt). Incorporation of more label increased Cy3 signal but lowered Cy5 fluorescence. This effect is probably due to quenching, which is more prominent for Cy5 than for Cy3. In conclusion, the currently developed method is an efficient, robust and inexpensive technique for fluorescent labelling of cRNA and allows sensitive detection of gene expression profiles on oligonucleotide microarrays.

The Human Genome Project and the future of diagnostics, treatment and prevention.

The Human Genome Project, the mapping of our 30,000-50,000 genes and the sequencing of all of our DNA, will have major impact on biomedical research and the whole of therapeutic and preventive health care. The tracing of genetic diseases to their molecular causes is rapidly expanding diagnostic and preventive options. The increased insights into molecular pathways, gained from high-throughput 'functional genomics', using DNA-chip and protein-chip approaches and specially designed animal model systems, will open great prospects for pharmacological and genetic therapies. Powerful bioinformatics and biostatistics will further improve our pattern recognition and accelerate progress. A rapidly expanding area of high expectations is that of 'pharmacogenomics': the design of more effective drugs with lower toxicity through tailoring of drug treatment to individual, genetically determined differences in drug metabolism. Not only will this decrease the cost of health care through reduction of adverse drug reactions, but a better stratification of populations will also provide more statistical power farther upstream in drug trials. However, the optimal benefits from the current explosion of 'data mining' will only be realized when the basic data are made and kept publicly accessible, while at the same time safeguarding the protection of intellectual property arising from downstream inventions. This is one of the goals of HUGO, the international Human Genome Organization, established 13 years ago to assist coordination of data acquisition and exchange and societal implementation of the genome project. Additional points of attention in this historic endeavour are the prevention of stigmatization and discrimination and the safeguarding of a worldwide balance in the contribution by--and benefits to--different populations, while respecting the diversity in cultures and traditions.

FcgammaRI (CD64) contributes substantially to severity of arthritis, hypersensitivity responses, and protection from bacterial infection.

The high-affinity receptor for IgG, FcgammaRI, shares its capacity to bind IgG2a immune complexes (IgG2a-IC) with the low-affinity receptor FcgammaRIII and complement factors, hampering the definition of its biological role. Moreover, in vivo, FcgammaRI is occupied by monomeric IgG2a, reducing its accessibility to newly formed IgG2a-IC. By using a variety of FcgammaR(-/-) mice, we demonstrate that in the absence of FcgammaRI, the IgG2a-IC-induced cellular processes of phagocytosis, cytokine release, cellular cytotoxicity, and antigen presentation are impaired. FcgammaRI(-/-) mice showed impaired hypersensitivity responses, strongly reduced cartilage destruction in an arthritis model, and impaired protection from a bacterial infection. We conclude that FcgammaRI contributes substantially to a variety of IgG2a-IC-dependent immune functions and immunopathological responses.

A partial BRCA1 sequence homology mapping to 4q28.

Using a BRCA1 cDNA probe in Southern analysis, we detected a sequence of 348 bp on 4q28 that is homologous to the 3' end of BRCA1. A 28-kb sequence contig has been assembled spanning the homologous region, which we designated BRCA1-h. An open reading frame was identified encoding a sequence of 82 amino acids; 22 of the last 23 amino acids are identical to the last 23 residues of BRCA1. BLAST-searches, RT-PCR and RACE-experiments have been unable to provide evidence that BRCA1-h is part of an expressed gene.