Prediction of single-nucleotide substitutions that result in exon skipping: identification of a splicing silencer in BRCA1 exon 6.

Missense, nonsense, and translationally silent mutations can inactivate genes by altering the inclusion of mutant exons in mRNA, but their overall frequency among disease-causing exonic substitutions is unknown. Here, we have tested missense and silent mutations deposited in the BRCA1 mutation databases of unclassified variants for their effects on exon inclusion. Analysis of 21 BRCA1 variants using minigene assays revealed a single exon-skipping mutation c.231G>T. Comprehensive mutagenesis of an adjacent 12-nt segment showed that this silent mutation resulted in a higher level of exon skipping than the 35 other single-nucleotide substitutions. Exon inclusion levels of mutant constructs correlated significantly with predicted splicing enhancers/silencers, prompting the development of two online utilities freely available at http://www.dbass.org.uk. EX-SKIP quickly estimates which allele is more susceptible to exon skipping, whereas HOT-SKIP examines all possible mutations at each exon position and identifies candidate exon-skipping positions/substitutions. We demonstrate that the distribution of exon-skipping and disease-associated substitutions previously identified in coding regions was biased toward top-ranking HOT-SKIP mutations. Finally, we show that proteins 9G8, SC35, SF2/ASF, Tra2, and hnRNP A1 were associated with significant alterations of BRCA1 exon 6 inclusion in the mRNA. Together, these results facilitate prediction of exonic substitutions that reduce exon inclusion in mature transcripts.

Proviruses selected for high and stable expression of transduced genes accumulate in broadly transcribed genome areas.

Retroviruses and retrovirus-derived vectors integrate nonrandomly into the genomes of host cells with specific preferences for transcribed genes, gene-rich regions, and CpG islands. However, the genomic features that influence the transcriptional activities of integrated retroviruses or retroviral vectors are poorly understood. We report here the cloning and characterization of avian sarcoma virus integration sites from chicken tumors. Growing progressively, dependent on high and stable expression of the transduced v-src oncogene, these tumors represent clonal expansions of cells bearing transcriptionally active replication-defective proviruses. Therefore, integration sites in our study distinguished genomic loci favorable for the expression of integrated retroviruses and gene transfer vectors. Analysis of integration sites from avian sarcoma virus-induced tumors showed strikingly nonrandom distribution, with proviruses found prevalently within or close to transcription units, particularly in genes broadly expressed in multiple tissues but not in tissue-specifically expressed genes. We infer that proviruses integrated in these genomic areas efficiently avoid transcriptional silencing and remain active for a long time during the growth of tumors. Defining the differences between unselected retroviral integration sites and sites selected for long-terminal-repeat-driven gene expression is relevant for retrovirus-mediated gene transfer and has ramifications for gene therapy.

The chicken Z chromosome is enriched for genes with preferential expression in ovarian somatic cells.

Theory predicts that sexually antagonistic mutations will be over- or under-represented on the X and Z chromosomes, depending on their average dominance coefficients. However, as little is known about the dominance coefficients for new mutations, the effect of sexually antagonistic selection is difficult to predict. To elucidate the role of sexually antagonistic selection in the evolution of Z chromosome gene content in chicken, we analyzed publicly available microarray data from several somatic tissues as well as somatic and germ cells of the ovary. We found that the Z chromosome is enriched for genes showing preferential expression in ovarian somatic cells, but not for genes with preferential expression in primary oocytes or non-sex-specific somatic tissues. Our results suggest that sexual antagonism leads to a higher abundance of female-benefit alleles on the Z chromosome. No bias toward Z-linkage for oocyte-enriched genes can be explained by lower intensity of sexually antagonistic selection in ovarian germ cells compared to ovarian somatic cells. An alternative explanation would be that meiotic Z chromosome inactivation hinders accumulation of oocyte-expressed genes on the Z chromosome. Our results are consistent with findings in mammals and indicate that recessive rather than dominant sexually antagonistic mutations shape the gene content of the X and Z chromosomes.

Ab initio prediction of mutation-induced cryptic splice-site activation and exon skipping.

Mutations that affect splicing of precursor messenger RNAs play a major role in the development of hereditary diseases. Most splicing mutations have been found to eliminate GT or AG dinucleotides that define the 5' and 3' ends of introns, leading to exon skipping or cryptic splice-site activation. Although accurate description of the mis-spliced transcripts is critical for predicting phenotypic consequences of these alterations, their exact nature in affected individuals cannot often be determined experimentally. Using a comprehensive collection of exons that sustained cryptic splice-site activation or were skipped as a result of splice-site mutations, we have developed a multivariate logistic discrimination procedure that distinguishes the two aberrant splicing outcomes from DNA sequences. The new algorithm was validated using an independent sample of exons and implemented as a free online utility termed CRYP-SKIP (http://www.dbass.org.uk/cryp-skip/). The web application takes up one or more mutated alleles, each consisting of one exon and flanking intronic sequences, and provides a list of important predictor variables and their values, the overall probability of activating cryptic splice vs exon skipping, and the location and intrinsic strength of predicted cryptic splice sites in the input sequence. These results will facilitate phenotypic prediction of splicing mutations and provide further insights into splicing enhancer and silencer elements and their relative importance for splice-site selection in vivo.

Mouse consomic strains: exploiting genetic divergence between Mus m. musculus and Mus m. domesticus subspecies.

Consomic (chromosome substitution) strains (CSs) represent the most recent addition to the mouse genetic resources aimed to genetically analyze complex trait loci (QTLs). In this study, we report the development of a set of 28 mouse intersubspecific CSs. In each CS, we replaced a single chromosome of the C57BL/6J (B6) inbred strain (mostly Mus m. domesticus) with its homolog from the PWD/Ph inbred strain of the Mus m. musculus subspecies. These two progenitor subspecies diverged less than 1 million years ago and accumulated a large number of genetic differences that constitute a rich resource of genetic variation for QTL analyses. Altogether, the 18 consomic, nine subconsomic, and one conplastic strain covered all 19 autosomes, X and Y sex chromosomes, and mitochondrial DNA. Most CSs had significantly lower reproductive fitness compared with the progenitor strains. CSs homosomic for chromosomes 10 and 11, and the C57BL/6J-Chr X males, failed to reproduce and were substituted by less affected subconsomics carrying either a proximal, central, or distal part of the respective chromosome. A genome-wide scan of 965 DNA markers revealed 99.87% purity of the B6 genetic background. Thirty-three nonsynonymous substitutions were uncovered in the protein-coding regions of the mitochondrial DNA of the B6.PWD-mt conplastic strain. A pilot-phenotyping experiment project revealed a high number of variations among B6.PWD consomics.

Development of a human mitochondrial oligonucleotide microarray (h-MitoArray) and gene expression analysis of fibroblast cell lines from 13 patients with isolated F1Fo ATP synthase deficiency.

BACKGROUND: To strengthen research and differential diagnostics of mitochondrial disorders, we constructed and validated an oligonucleotide microarray (h-MitoArray) allowing expression analysis of 1632 human genes involved in mitochondrial biology, cell cycle regulation, signal transduction and apoptosis. Using h-MitoArray we analyzed gene expression profiles in 9 control and 13 fibroblast cell lines from patients with F1Fo ATP synthase deficiency consisting of 2 patients with mt9205deltaTA microdeletion and a genetically heterogeneous group of 11 patients with not yet characterized nuclear defects. Analysing gene expression profiles, we attempted to classify patients into expected defect specific subgroups, and subsequently reveal group specific compensatory changes, identify potential phenotype causing pathways and define candidate disease causing genes. RESULTS: Molecular studies, in combination with unsupervised clustering methods, defined three subgroups of patient cell lines--M group with mtDNA mutation and N1 and N2 groups with nuclear defect. Comparison of expression profiles and functional annotation, gene enrichment and pathway analyses of differentially expressed genes revealed in the M group a transcription profile suggestive of synchronized suppression of mitochondrial biogenesis and G1/S arrest. The N1 group showed elevated expression of complex I and reduced expression of complexes III, V, and V-type ATP synthase subunit genes, reduced expression of genes involved in phosphorylation dependent signaling along MAPK, Jak-STAT, JNK, and p38 MAP kinase pathways, signs of activated apoptosis and oxidative stress resembling phenotype of premature senescent fibroblasts. No specific functionally meaningful changes, except of signs of activated apoptosis, were detected in the N2 group. Evaluation of individual gene expression profiles confirmed already known ATP6/ATP8 defect in patients from the M group and indicated several candidate disease causing genes for nuclear defects. CONCLUSION: Our analysis showed that deficiency in the ATP synthase protein complex amount is generally accompanied by only minor changes in expression of ATP synthase related genes. It also suggested that the site (mtDNA vs nuclear DNA) and the severity (ATP synthase content) of the underlying defect have diverse effects on cellular gene expression phenotypes, which warrants further investigation of cell cycle regulatory and signal transduction pathways in other OXPHOS disorders and related pharmacological models.

Construction and characterization of a genomic BAC library for the Mus m. musculus mouse subspecies (PWD/Ph inbred strain).

BACKGROUND: The genome of classical laboratory strains of mice is an artificial mosaic of genomes originated from several mouse subspecies with predominant representation (>90%) of the Mus m. domesticus component. Mice of another subspecies, East European/Asian Mus m. musculus, can interbreed with the classical laboratory strains to generate hybrids with unprecedented phenotypic and genotypic variations. To study these variations in depth we prepared the first genomic large insert BAC library from an inbred strain derived purely from the Mus m. musculus-subspecies. The library will be used to seek and characterize genomic sequences controlling specific monogenic and polygenic complex traits, including modifiers of dominant and recessive mutations. RESULTS: A representative mouse genomic BAC library was derived from a female mouse of the PWD/Ph inbred strain of Mus m. musculus subspecies. The library consists of 144,768 primary clones from which 97% contain an insert of 120 kb average size. The library represents an equivalent of 6.7 x mouse haploid genome, as estimated from the total number of clones carrying genomic DNA inserts and from the average insert size. The clones were arrayed in duplicates onto eight high-density membranes that were screened with seven single-copy gene probes. The individual probes identified four to eleven positive clones, corresponding to 6.9-fold coverage of the mouse genome. Eighty-seven BAC-ends of PWD/Ph clones were sequenced, edited, and aligned with mouse C57BL/6J (B6) genome. Seventy-three BAC-ends displayed unique hits on B6 genome and their alignment revealed 0.92 single nucleotide polymorphisms (SNPs) per 100 bp. Insertions and deletions represented 0.3% of the BAC end sequences. CONCLUSION: Analysis of the novel genomic library for the PWD/Ph inbred strain demonstrated coverage of almost seven mouse genome equivalents and a capability to recover clones for specific regions of PWD/Ph genome. The single nucleotide polymorphism between the strains PWD/Ph and C57BL/6J was 0.92/100 bp, a value significantly higher than between classical laboratory strains. The library will serve as a resource for dissecting the phenotypic and genotypic variations between mice of the Mus m. musculus subspecies and classical laboratory mouse strains.

Global transcriptome analysis of the C57BL/6J mouse testis by SAGE: evidence for nonrandom gene order.

BACKGROUND: We generated the gene expression profile of the total testis from the adult C57BL/6J male mice using serial analysis of gene expression (SAGE). Two high-quality SAGE libraries containing a total of 76 854 tags were constructed. An extensive bioinformatic analysis and comparison of SAGE transcriptomes of the total testis, testicular somatic cells and other mouse tissues was performed and the theory of male-biased gene accumulation on the X chromosome was tested. RESULTS: We sorted out 829 genes predominantly expressed from the germinal part and 944 genes from the somatic part of the testis. The genes preferentially and specifically expressed in total testis and testicular somatic cells were identified by comparing the testis SAGE transcriptomes to the available transcriptomes of seven non-testis tissues. We uncovered chromosomal clusters of adjacent genes with preferential expression in total testis and testicular somatic cells by a genome-wide search and found that the clusters encompassed a significantly higher number of genes than expected by chance. We observed a significant 3.2-fold enrichment of the proportion of X-linked genes specific for testicular somatic cells, while the proportions of X-linked genes specific for total testis and for other tissues were comparable. In contrast to the tissue-specific genes, an under-representation of X-linked genes in the total testis transcriptome but not in the transcriptomes of testicular somatic cells and other tissues was detected. CONCLUSION: Our results provide new evidence in favor of the theory of male-biased genes accumulation on the X chromosome in testicular somatic cells and indicate the opposite action of the meiotic X-inactivation in testicular germ cells.

Genetic analysis of X-linked hybrid sterility in the house mouse.

Hybrid sterility is a common postzygotic reproductive isolation mechanism that appears in the early stages of speciation of various organisms. Mus musculus musculus and Mus musculus domesticus represent two recently separated mouse subspecies particularly suitable for genetic studies of hybrid sterility. Here we show that the introgression of Chr X of M. m. musculus origin (PWD/Ph inbred strain, henceforth PWD) into the genetic background of the C57BL/6J (henceforth B6) inbred strain (predominantly of M. m. domesticus origin) causes male sterility. The X-linked hybrid sterility is associated with reduced testes weight, lower sperm count, and morphological abnormalities of sperm heads. The analysis of recombinant Chr Xs in sterile and fertile males as well as quantitative trait locus (QTL) analysis of several fertility parameters revealed an oligogenic nature of the X-linked hybrid sterility. The Hstx1 locus responsible for male sterility was mapped near DXMit119 in the central part of Chr X. To ensure full sterility, the PWD allele of Hstx1 has to be supported with the PWD allelic form of loci in at least one proximal and/or one distal region of Chr X. Mapping and cloning of Hstx1 and other genes responsible for sterility of B6-X PWD Y B6 males could help to elucidate the special role of Chr X in hybrid sterility and consequently in speciation.

The Mouse SAGE Site: database of public mouse SAGE libraries.

The Mouse SAGE Site is a web-based database of all available public libraries generated by the Serial Analysis of Gene Expression (SAGE) from various mouse tissues and cell lines. The database contains mouse SAGE libraries organized in a uniform way and provides web-based tools for browsing, comparing and searching SAGE data with reliable tag-to-gene identification. A modified approach based on the SAGEmap database is used for reliable tag identification. The Mouse SAGE Site is maintained on an ongoing basis at the Institute of Molecular Genetics, Academy of Sciences of the Czech Republic and is accessible at the internet address http://mouse.biomed.cas.cz/sage/.