Novel (60%) and recurrent (40%) androgen receptor gene mutations in a series of 59 patients with a 46,XY disorder of sex development.

BACKGROUND Androgen receptor (AR) gene mutations are the most frequent cause of 46,XY disorders of sex development (DSD) and are associated with a variety of phenotypes, ranging from phenotypic women [complete androgen insensitivity syndrome (CAIS)] to milder degrees of undervirilization (partial form or PAIS) or men with only infertility (mild form or MAIS). OBJECTIVE The aim of the study was to characterize the contribution of the AR gene to the molecular cause of 46,XY DSD in a series of Spanish patients. SETTING We studied a series of 133 index patients with 46,XY DSD in whom gonads were differentiated as testes, with phenotypes including varying degrees of undervirilization, and in whom the AR gene was the first candidate for a molecular analysis. METHODS The AR gene was sequenced (exons 1 to 8 with intronic flanking regions) in all patients and in family members of 61% of AR-mutated gene patients. RESULTS AR gene mutations were found in 59 individuals (44.4% of index patients), of whom 46 (78%) were CAIS and 13 (22%) PAIS. Fifty-seven different mutations were found: 21.0% located in exon 1, 15.8% in exons 2 and 3, 57.9% in exons 4-8, and 5.3% intronic. Twenty-three mutations (40.4%) had been previously described and 34 (59.6%) were novel. CONCLUSIONS AR gene mutation is the most frequent cause of 46,XY DSD, with a clearly higher frequency in the complete phenotype. Mutations spread along the whole coding sequence, including exon 1. This series shows that 60% of mutations detected during the period 2002-2009 were novel.

Identification of a novel splice-site mutation in the CYP1A2 gene.

AIMS: To identify the molecular basis for a low CYP1A2 metabolic status, as determined by a caffeine phenotyping test, in a 71-year-old, nonsmoking, Caucasian woman who presented with very high clozapine concentrations despite being administered a standard dose of the drug. METHODS: The nucleotide sequence of the 7 exons, exon-intron boundaries and 5'-flanking region of the CYP1A2 gene was analysed by direct sequencing. RESULTS: Only one heterozygous point mutation was identified in the donor splice site of intron 6 (3534G > A) of CYP1A2. This mutation could cause abnormal RNA splicing and therefore lead to a truncated nonfunctional enzyme. No other carrier of this mutation was identified in a population of 100 unrelated healthy Caucasians. CONCLUSIONS: This is the first report of a splice-site mutation affecting the CYP1A2 gene. This polymorphism is a likely explanation for the low CYP1A2 activity associated with high clozapine concentrations in this patient.

Genomic acute myeloid leukemia-associated inv(16)(p13q22) breakpoints are tightly clustered.

The inv(16) and related t(16;16) are found in 10% of all cases with de novo acute myeloid leukemia. In these rearrangements the core binding factor beta (CBFB) gene on 16q22 is fused to the smooth muscle myosin heavy chain gene (MYH11) on 16p13. To gain insight into the mechanisms causing the inv(16) we have analysed 24 genomic CBFB-MYH11 breakpoints. All breakpoints in CBFB are located in a 15-Kb intron. More than 50% of the sequenced 6.2 Kb of this intron consists of human repetitive elements. Twenty-one of the 24 breakpoints in MYH11 are located in a 370-bp intron. The remaining three breakpoints in MYH11 are located more upstream. The localization of three breakpoints adjacent to a V(D)J recombinase signal sequence in MYH11 suggests a V(D)J recombinase-mediated rearrangement in these cases. V(D)J recombinase-associated characteristics (small nucleotide deletions and insertions of random nucleotides) were detected in six other cases. CBFB and MYH11 duplications were detected in four of six cases tested.

Intronic variants in the NFKB1 gene may influence hearing forecast in patients with unilateral sensorineural hearing loss in Meniere's disease.

Meniere's disease is an episodic vestibular syndrome associated with sensorineural hearing loss (SNHL) and tinnitus. Patients with MD have an elevated prevalence of several autoimmune diseases (rheumatoid arthritis, systemic lupus erythematosus, ankylosing spondylitis and psoriasis), which suggests a shared autoimmune background. Functional variants of several genes involved in the NF-κB pathway, such as REL, TNFAIP3, NFKB1 and TNIP1, have been associated with two or more immune-mediated diseases and allelic variations in the TLR10 gene may influence bilateral affectation and clinical course in MD. We have genotyped 716 cases of MD and 1628 controls by using the ImmunoChip, a high-density genotyping array containing 186 autoimmune loci, to explore the association of immune system related-loci with sporadic MD. Although no single nucleotide polymorphism (SNP) reached a genome-wide significant association (p<10(-8)), we selected allelic variants in the NF-kB pathway for further analyses to evaluate the impact of these SNPs in the clinical outcome of MD in our cohort. None of the selected SNPs increased susceptibility for MD in patients with uni or bilateral SNHL. However, two potential regulatory variants in the NFKB1 gene (rs3774937 and rs4648011) were associated with a faster hearing loss progression in patients with unilateral SNHL. So, individuals with unilateral MD carrying the C allele in rs3774937 or G allele in rs4648011 had a shorter mean time to reach hearing stage 3 (>40 dB HL) (log-rank test, corrected p values were p = 0.009 for rs3774937 and p = 0.003 for rs4648011, respectively). No variants influenced hearing in bilateral MD. Our data support that the allelic variants rs3774937 and rs4648011 can modify hearing outcome in patients with MD and unilateral SNHL.

CYP3A5 and ABCB1 genes and hypertension.

Hypertension is the first single modifiable cause of disease burden worldwide. Genes encoding proteins that are involved in the metabolism (CYP3A5) and transport (ABCB1) of drugs and hormones might contribute to blood pressure control in humans. Indeed, recent data have suggested that CYP3A5 and ABCB1 gene polymorphisms are associated with blood pressure in the rat as well as in humans. Interestingly, the effects of these genes on blood pressure appear to be modified by dietary salt intake. This review summarizes what is known regarding the relationships of the ABCB1 and CYP3A5 genes with blood pressure, and discusses the potential underlying mechanisms of the association. If the role of these genes in blood pressure control is confirmed in other populations and other ethnic groups, these findings would point toward a new pathway for blood pressure control in humans.

Systematics of snow voles (Chionomys, Arvicolinae) revisited.

To elucidate the evolutionary history of snow voles, genus Chionomys, we studied the phylogeography of Chionomysnivalis across its range and investigated its relationships with two congeneric species, Chionomysgud and Chionomysroberti, using independent molecular markers. Analyses were based on mitochondrial (~940 bp cyt b) and Y-chromosomal (~2020 bp from three introns) genetic variation. Our data provide conclusive evidence for a Caucasian and Middle Eastern origin for the three species and a subsequent westward expansion of C.nivalis. In addition, we discuss the taxonomic status of the genus Chionomys in relation to the genus Microtus.

Splicing and the evolution of proteins in mammals.

It is often supposed that a protein's rate of evolution and its amino acid content are determined by the function and anatomy of the protein. Here we examine an alternative possibility, namely that the requirement to specify in the unprocessed RNA, in the vicinity of intron-exon boundaries, information necessary for removal of introns (e.g., exonic splice enhancers) affects both amino acid usage and rates of protein evolution. We find that the majority of amino acids show skewed usage near intron-exon boundaries, and that differences in the trends for the 2-fold and 4-fold blocks of both arginine and leucine show this to be owing to effects mediated at the nucleotide level. More specifically, there is a robust relationship between the extent to which an amino acid is preferred/avoided near boundaries and its enrichment/paucity in splice enhancers. As might then be expected, the rate of evolution is lowest near intron-exon boundaries, at least in part owing to splice enhancers, such that domains flanking intron-exon junctions evolve on average at under half the rate of exon centres from the same gene. In contrast, the rate of evolution of intronless retrogenes is highest near the domains where intron-exon junctions previously resided. The proportion of sequence near intron-exon boundaries is one of the stronger predictors of a protein's rate of evolution in mammals yet described. We conclude that after intron insertion selection favours modification of amino acid content near intron-exon junctions, so as to enable efficient intron removal, these changes then being subject to strong purifying selection even if nonoptimal for protein function. Thus there exists a strong force operating on protein evolution in mammals that is not explained directly in terms of the biology of the protein.

Comparison of splice sites in mammals and chicken

We have carried out an initial analysis of the dynamics of the recent evolution of the splice-sites sequences on a large collection of human, rodent (mouse and rat), and chicken introns. Our results indicate that the sequences of splice sites are largely homogeneous within tetrapoda. We have also found that orthologous splice signals between human and rodents and within rodents are more conserved than unrelated splice sites, but the additional conservation can be explained mostly by background intron conservation. In contrast, additional conservation over background is detectable in orthologous mammalian and chicken splice sites. Our results also indicate that the U2 and U12 intron classes seem to have evolved independently since the split of mammals and birds; we have not been able to find a convincing case of interconversion between these two classes in our collections of orthologous introns. Similarly, we have not found a single case of switching between AT-AC and GT-AG subtypes within U12 introns, suggesting that this event has been a rare occurrence in recent evolutionary times. Switching between GT-AG and the noncanonical GC-AG U2 subtypes, on the contrary, does not appear to be unusual; in particular, T to C mutations appear to be relatively well tolerated in GT-AG introns with very strong donor sites.

Improving gene annotation using peptide mass spectrometry

Annotation of protein-coding genes is a key goal of genome sequencing projects. In spite of tremendous recent advances in computational gene finding, comprehensive annotation remains a challenge. Peptide mass spectrometry is a powerful tool for researching the dynamic proteome and suggests an attractive approach to discover and validate protein-coding genes. We present algorithms to construct and efficiently search spectra against a genomic database, with no prior knowledge of encoded proteins. By searching a corpus of 18.5 million tandem mass spectra (MS/MS) from human proteomic samples, we validate 39,000 exons and 11,000 introns at the level of translation. We present translation-level evidence for novel or extended exons in 16 genes, confirm translation of 224 hypothetical proteins, and discover or confirm over 40 alternative splicing events. Polymorphisms are efficiently encoded in our database, allowing us to observe variant alleles for 308 coding SNPs. Finally, we demonstrate the use of mass spectrometry to improve automated gene prediction, adding 800 correct exons to our predictions using a simple rescoring strategy. Our results demonstrate that proteomic profiling should play a role in any genome sequencing project.

Gene finding in the chicken genome

Background: Despite the continuous production of genome sequence for a number of organisms,reliable, comprehensive, and cost effective gene prediction remains problematic. This is particularlytrue for genomes for which there is not a large collection of known gene sequences, such as therecently published chicken genome. We used the chicken sequence to test comparative andhomology-based gene-finding methods followed by experimental validation as an effective genomeannotation method.Results: We performed experimental evaluation by RT-PCR of three different computational genefinders, Ensembl, SGP2 and TWINSCAN, applied to the chicken genome. A Venn diagram wascomputed and each component of it was evaluated. The results showed that de novo comparativemethods can identify up to about 700 chicken genes with no previous evidence of expression, andcan correctly extend about 40% of homology-based predictions at the 5' end.Conclusions: De novo comparative gene prediction followed by experimental verification iseffective at enhancing the annotation of the newly sequenced genomes provided by standardhomology-based methods.

The quest for a message: budding yeast, a model organism to study the control of pre-mRNA splicing

Removal of introns during pre-mRNA splicing is a critical process in gene expression, and understanding its control at both single-gene and genomic levels is one of the great challenges in Biology. Splicing takes place in a dynamic, large ribonucleoprotein complex known as the spliceosome. Combining Genetics and Biochemistry, Saccharomyces cerevisiae provides insights into its mechanisms, including its regulation by RNA-protein interactions. Recent genome-wide analyses indicate that regulated splicing is broad and biologically relevant even in organisms with a relatively simple intronic structure, such as yeast. Furthermore, the possibility of coordination in splicing regulation at genomic level is becoming clear in this model organism. This should provide a valuable system to approach the complex problem of the role of regulated splicing in genomic expression.

Deciphering 3'ss selection in the yeast genome reveals an RNA thermosensor that mediates alternative splicing

Poor understanding of the spliceosomal mechanisms to select intronic 3' ends (3'ss) is a major obstacle to deciphering eukaryotic genomes. Here, we discern the rules for global 3'ss selection in yeast. We show that, in contrast to the uniformity of yeast splicing, the spliceosome uses all available 3'ss within a distance window from the intronic branch site (BS), and that in 70% of all possible 3'ss this is likely to be mediated by pre-mRNA structures. Our results reveal that one of these RNA folds acts as an RNA thermosensor, modulating alternative splicing in response to heat shock by controlling alternate 3'ss availability. Thus, our data point to a deeper role for the pre-mRNA in the control of its own fate, and to a simple mechanism for some alternative splicing.

Differentiated evolutionary rates in alternative exons and the implications for splicing regulation

Background: Alternatively spliced exons play an important role in the diversification of gene function in most metazoans and are highly regulated by conserved motifs in exons and introns. Two contradicting properties have been associated to evolutionary conserved alternative exons: higher sequence conservation and higher rate of non-synonymous substitutions, relative to constitutive exons. In order to clarify this issue, we have performed an analysis of the evolution of alternative and constitutive exons, using a large set of protein coding exons conserved between human and mouse and taking into account the conservation of the transcript exonic structure. Further, we have also defined a measure of the variation of the arrangement of exonic splicing enhancers (ESE-conservation score) to study the evolution of splicing regulatory sequences. We have used this measure to correlate the changes in the arrangement of ESEs with the divergence of exon and intron sequences. Results: We find evidence for a relation between the lack of conservation of the exonic structure and the weakening of the sequence evolutionary constraints in alternative and constitutive exons. Exons in transcripts with non-conserved exonic structures have higher synonymous (dS) and non-synonymous (dN) substitution rates than exons in conserved structures. Moreover, alternative exons in transcripts with non-conserved exonic structure are the least constrained in sequence evolution, and at high EST-inclusion levels they are found to be very similar to constitutive exons, whereas alternative exons in transcripts with conserved exonic structure have a dS significantly lower than average at all EST-inclusion levels. We also find higher conservation in the arrangement of ESEs in constitutive exons compared to alternative ones. Additionally, the sequence conservation at flanking introns remains constant for constitutive exons at all ESE-conservation values, but increases for alternative exons at high ESE-conservation values. Conclusion: We conclude that most of the differences in dN observed between alternative and constitutive exons can be explained by the conservation of the transcript exonic structure. Low dS values are more characteristic of alternative exons with conserved exonic structure, but not of those with non-conserved exonic structure. Additionally, constitutive exons are characterized by a higher conservation in the arrangement of ESEs, and alternative exons with an ESE-conservation similar to that of constitutive exons are characterized by a conservation of the flanking intron sequences higher than average, indicating the presence of more intronic regulatory signals.

Genome organization and gene expression shape the transposable element distribution in the Drosophila melanogaster euchromatin.

The distribution of transposable elements (TEs) in a genome reflects a balance between insertion rate and selection against new insertions. Understanding the distribution of TEs therefore provides insights into the forces shaping the organization of genomes. Past research has shown that TEs tend to accumulate in genomic regions with low gene density and low recombination rate. However, little is known about the factors modulating insertion rates across the genome and their evolutionary significance. One candidate factor is gene expression, which has been suggested to increase local insertion rate by rendering DNA more accessible. We test this hypothesis by comparing the TE density around germline- and soma-expressed genes in the euchromatin of Drosophila melanogaster. Because only insertions that occur in the germline are transmitted to the next generation, we predicted a higher density of TEs around germline-expressed genes than soma-expressed genes. We show that the rate of TE insertions is greater near germline- than soma-expressed genes. However, this effect is partly offset by stronger selection for genome compactness (against excess noncoding DNA) on germline-expressed genes. We also demonstrate that the local genome organization in clusters of coexpressed genes plays a fundamental role in the genomic distribution of TEs. Our analysis shows that-in addition to recombination rate-the distribution of TEs is shaped by the interaction of gene expression and genome organization. The important role of selection for compactness sheds a new light on the role of TEs in genome evolution. Instead of making genomes grow passively, TEs are controlled by the forces shaping genome compactness, most likely linked to the efficiency of gene expression or its complexity and possibly their interaction with mechanisms of TE silencing.

Assessment of transcript reconstruction methods for RNA-seq.

We evaluated 25 protocol variants of 14 independent computational methods for exon identification, transcript reconstruction and expression-level quantification from RNA-seq data. Our results show that most algorithms are able to identify discrete transcript components with high success rates but that assembly of complete isoform structures poses a major challenge even when all constituent elements are identified. Expression-level estimates also varied widely across methods, even when based on similar transcript models. Consequently, the complexity of higher eukaryotic genomes imposes severe limitations on transcript recall and splice product discrimination that are likely to remain limiting factors for the analysis of current-generation RNA-seq data.