Compensatory evolution of a precursor messenger RNA secondary structure in the Drosophila melanogaster Adh gene

Epistasis-Driven Evolution of the SARS-CoV-2 Secondary Structure

Epistasis is an evolutionary phenomenon whereby the fitness effect of a mutation depends on the genetic background in which it arises. A key source of epistasis in an RNA molecule is its secondary structure, which contains functionally important topological motifs held together by hydrogen bonds between Watson–Crick (WC) base pairs. Here we study epistasis in the secondary structure of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) by examining properties of derived alleles arising from substitution mutations at ancestral WC base-paired and unpaired (UP) sites in 15 conserved topological motifs across the genome. We uncover fewer derived alleles and lower derived allele frequencies at WC than at UP sites, supporting the hypothesis that modifications to the secondary structure are often deleterious. At WC sites, we also find lower derived allele frequencies for mutations that abolish base pairing than for those that yield G·U “wobbles,” illustrating that weak base pairing can partially preserve the integrity of the secondary structure. Last, we show that WC sites under the strongest epistatic constraint reside in a three-stemmed pseudoknot motif that plays an essential role in programmed ribosomal frameshifting, whereas those under the weakest epistatic constraint are located in 3’ UTR motifs that regulate viral replication and pathogenicity. Our findings demonstrate the importance of epistasis in the evolution of the SARS-CoV-2 secondary structure, as well as highlight putative structural and functional targets of different forms of natural selection.

Adenine·cytosine substitutions are an alternative pathway of compensatory mutation in angiosperm ITS2

Compensatory mutations are crucial for functional RNA because they maintain RNA configuration and thus function. Compensatory mutation has traditionally been considered to be a two-step substitution through the GU-base-pair intermediate. We tested for an alternative AC-mediated compensatory mutation (ACCM). We investigated ACCMs by using a comprehensive sampling of ribosomal internal transcribed spacer 2 (ITS2) from 3934 angiosperm species in 80 genera and 55 families. We predicted ITS2 consensus secondary structures by using LocARNA for structure-based alignment and partitioning paired and unpaired regions. We examined and compared the substitution rates and frequencies among base pairs by using RNA-specific models. Base-pair states of ACCMs were mapped onto the inferred phylogenetic trees to infer their evolution. All types of compensatory mutations involving the AC intermediate were observed, but the most frequent substitutions were with AU or GC pairs, which are part of the AU-AC-GC pathway. Compared with the GU intermediate, AC had a lower frequency and higher mutability. Within the AU-AC-GC pathway, the AU-AC substitution rate was much slower than the AC-GC substitution rate. No consistently higher overall rate was identified for either pathway among all 80 sampled lineages, though compensatory mutations through the AC intermediate averaged about half that through the GU intermediate. These results demonstrate an alternative compensatory mutation between AU and GC that helps address the controversial inference of inferred simultaneous double substitutions.

Association between TAAR6 polymorphisms and airway responsiveness to inhaled corticosteroids in asthmatic patients

BACKGROUND

Genetic polymorphisms may be responsible for the wide variation in response to inhaled corticosteroids in asthmatic patients. We had previously reported that one polymorphism rs7772821, located on the 3'-UTR of trace amine-associated receptor 6 (TAAR6), is significantly associated with percentile changes in the forced expiratory volume in 1 s (%ΔFEV1) after inhaled corticosteroid treatment in asthmatics using a genome-wide association study. The aim of the present study was to validate the association between 15 single-nucleotide polymorphisms (SNPs) on the TAAR6 and airway responsiveness to inhaled corticosteroids in the asthmatics.

METHODS

The %ΔFEV1 induced by 4 weeks' treatment with inhaled fluticasone propionate (1000 μg daily) was measured in 246 asthmatics. The 15 SNPs of TAAR6 were genotyped using a TaqMan assay. An association analysis between %ΔFEV1 and TAAR6 polymorphisms was carried out using a linear regression model controlling for age, sex, smoking status, presence of atopy, and baseline FEV1 as covariates.

RESULTS

Among the 15 SNPs and seven haplotypes of TAAR6, rs7772821 (T>G) on the 3'-UTR showed the strongest correlation with inhaled corticosteroid-induced %ΔFEV1 (Pcorr=0.002 in the codominant model, Pcorr=0.03 in the dominant model, Pcorr=0.01 in the recessive model). The %ΔFEV1 of the rs7772821T>G minor homozygotes (60.77%) was higher than that of patients harboring either the rs7772821 T/G or T/T genotypes (21.32 and 31.60%, respectively).

CONCLUSION

The TAAR6 rs7772821 polymorphism may be one of the important genetic factors for predicting the response to treatment with inhaled corticosteroids in asthmatics.

Alternative splicing of U2AF1 reveals a shared repression mechanism for duplicated exons

The auxiliary factor of U2 small nuclear ribonucleoprotein (U2AF) facilitates branch point (BP) recognition and formation of lariat introns. The gene for the 35-kD subunit of U2AF gives rise to two protein isoforms (termed U2AF35a and U2AF35b) that are encoded by alternatively spliced exons 3 and Ab, respectively. The splicing recognition sequences of exon 3 are less favorable than exon Ab, yet U2AF35a expression is higher than U2AF35b across tissues. We show that U2AF35b repression is facilitated by weak, closely spaced BPs next to a long polypyrimidine tract of exon Ab. Each BP lacked canonical uridines at position -2 relative to the BP adenines, with efficient U2 base-pairing interactions predicted only for shifted registers reminiscent of programmed ribosomal frameshifting. The BP cluster was compensated by interactions involving unpaired cytosines in an upstream, EvoFold-predicted stem loop (termed ESL) that binds FUBP1/2. Exon Ab inclusion correlated with predicted free energies of mutant ESLs, suggesting that the ESL operates as a conserved rheostat between long inverted repeats upstream of each exon. The isoform-specific U2AF35 expression was U2AF65-dependent, required interactions between the U2AF-homology motif (UHM) and the α6 helix of U2AF35, and was fine-tuned by exon Ab/3 variants. Finally, we identify tandem homologous exons regulated by U2AF and show that their preferential responses to U2AF65-related proteins and SRSF3 are associated with unpaired pre-mRNA segments upstream of U2AF-repressed 3′ss. These results provide new insights into tissue-specific subfunctionalization of duplicated exons in vertebrate evolution and expand the repertoire of exon repression mechanisms that control alternative splicing.

Molecular Outcome, Prediction, and Clinical Consequences of Splice Variants in COL1A1, Which Encodes the proα1(I) Chains of Type I Procollagen

Approximately 10%-20% of germline pathogenic variants alter mRNA splicing, with phenotypes often dependent on the stability of the mRNA produced by the mutant allele. To better understand the relationships between genotype, mRNA splicing, and phenotype, we examined clinical and molecular data from 243 probands with osteogenesis imperfecta (OI) representing 145 unique splicing variants within the type I procollagen gene, COL1A1. All individuals with IVSX-1G>A mutations had OI type I because the substitution shifted the splice acceptor site 1 nt downstream and destabilized the mRNA. OI phenotypes were not consistent for any other splice variant identified. We sequenced all cDNA species from cultured dermal fibroblasts from 40 individuals to identify splice outcome and compared those results to splice predictions from Human Splice Finder (HSF), Spliceport (SP), and Automatic Splice Site and Exon Definition Analyses (ASSEDA). Software-based splice predictions were correct in 42%, 55%, and 74% instances for HSF, SP, and ASSEDA, respectively. As molecular diagnostics move increasingly to DNA sequence analysis, the need to understand the effects of splice site variants will increase. These data demonstrate that caution must be exercised when using splice prediction software to predict splice outcome.

Intermediate introns in nuclear genes of euglenids - are they a distinct type?

Background

Nuclear genes of euglenids contain two major types of introns: conventional spliceosomal and nonconventional introns. The latter are characterized by variable non-canonical borders, RNA secondary structure that brings intron ends together, and an unknown mechanism of removal. Some researchers also distinguish intermediate introns, which combine features of both types. They form a stable RNA secondary structure and are classified into two subtypes depending on whether they contain one (intermediate/nonconventional subtype) or both (conventional/intermediate subtype) canonical spliceosomal borders. However, it has been also postulated that most introns classified as intermediate could simply be special cases of conventional or nonconventional introns.

Results

Sequences of tubB, hsp90 and gapC genes from six strains of Euglena agilis were obtained. They contain four, six, and two or three introns, respectively (the third intron in the gapC gene is unique for just one strain). Conventional introns were present at three positions: two in the tubB gene (at one position conventional/intermediate introns were also found) and one in the gapC gene. Nonconventional introns are present at ten positions: two in the tubB gene (at one position intermediate/nonconventional introns were also found), six in hsp90 (at four positions intermediate/nonconventional introns were also found), and two in the gapC gene.

Conclusions

Sequence and RNA secondary structure analyses of nonconventional introns confirmed that their most strongly conserved elements are base pairing nucleotides at positions +4, +5 and +6/ -8, −7 and −6 (in most introns CAG/CTG nucleotides were observed). It was also confirmed that the presence of the 5' GT/C end in intermediate/nonconventional introns is not the result of kinship with conventional introns, but is due to evolutionary pressure to preserve the purine at the 5' end. However, an example of a nonconventional intron with GC-AG ends was shown, suggesting the possibility of intron type conversion between nonconventional and conventional. Furthermore, an analysis of conventional introns revealed that the ability to form a stable RNA secondary structure by some introns is probably not a result of their relationship with nonconventional introns. It was also shown that acquisition of new nonconventional introns is an ongoing process and can be observed at the level of a single species. In the recently acquired intron in the gapC gene an extended direct repeats at the intron-exon junctions are present, suggesting that double-strand break repair process could be the source of new nonconventional introns.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0620-5) contains supplementary material, which is available to authorized users.

RNA structure replaces the need for U2AF2 in splicing

RNA secondary structure plays an integral role in catalytic, ribosomal, small nuclear, micro, and transfer RNAs. Discovering a prevalent role for secondary structure in pre-mRNAs has proven more elusive. By utilizing a variety of computational and biochemical approaches, we present evidence for a class of nuclear introns that relies upon secondary structure for correct splicing. These introns are defined by simple repeat expansions of complementary AC and GT dimers that co-occur at opposite boundaries of an intron to form a bridging structure that enforces correct splice site pairing. Remarkably, this class of introns does not require U2AF2, a core component of the spliceosome, for its processing. Phylogenetic analysis suggests that this mechanism was present in the ancestral vertebrate lineage prior to the divergence of tetrapods from teleosts. While largely lost from land dwelling vertebrates, this class of introns is found in 10% of all zebrafish genes.

The role of short RNA loops in recognition of a single-hairpin exon derived from a mammalian-wide interspersed repeat

Splice-site selection is controlled by secondary structure through sequestration or approximation of splicing signals in primary transcripts but the exact role of even the simplest and most prevalent structural motifs in exon recognition remains poorly understood. Here we took advantage of a single-hairpin exon that was activated in a mammalian-wide interspersed repeat (MIR) by a mutation stabilizing a terminal triloop, with splice sites positioned close to each other in a lower stem of the hairpin. We first show that the MIR exon inclusion in mRNA correlated inversely with hairpin stabilities. Employing a systematic manipulation of unpaired regions without altering splice-site configuration, we demonstrate a high correlation between exon inclusion of terminal tri- and tetraloop mutants and matching tri-/tetramers in splicing silencers/enhancers. Loop-specific exon inclusion levels and enhancer/silencer associations were preserved across primate cell lines, in 4 hybrid transcripts and also in the context of a distinct stem, but only if its loop-closing base pairs were shared with the MIR hairpin. Unlike terminal loops, splicing activities of internal loop mutants were predicted by their intramolecular Watson-Crick interactions with the antiparallel strand of the MIR hairpin rather than by frequencies of corresponding trinucleotides in splicing silencers/enhancers. We also show that splicing outcome of oligonucleotides targeting the MIR exon depend on the identity of the triloop adjacent to their antisense target. Finally, we identify proteins regulating MIR exon recognition and reveal a distinct requirement of adjacent exons for C-terminal extensions of Tra2α and Tra2β RNA recognition motifs.

A model of compensatory molecular evolution involving multiple sites in RNA molecules

Consider two sites under compensatory fitness interaction, such as a Watson-Crick base pair in an RNA helix or two interacting residues in a protein. A mutation at any one of these two sites may reduce the fitness of an individual. However, fitness may be restored by the occurrence of a second mutation at the other site. Kimura modeled this process using a two-locus haploid fitness scheme with two alleles at each locus. He predicted that compensatory evolution following this model is very rare unless selection against the deleterious single mutations is weak and linkage between the interacting sites is tight. Here we investigate the question whether the rate of compensatory evolution increases if we take the context of the two directly interacting sites into account. By "context", we mean the effect of neighboring sites in an RNA helix. Interaction between the focal pair of sites under consideration and the context may lead to so-called indirect compensation. Thus, extending Kimura's classical model of compensatory evolution, we study the effects of both direct and indirect compensation on the rate of compensatory evolution. It is shown that the effects of indirect compensation are very strong. We find that recombination does not slow down the rate of compensatory evolution as predicted by the classical model. Instead, compensatory substitutions may be relatively frequent, even if linkage between the focal interacting sites is loose, selection against deleterious mutations is strong, and mutation rate is low. We compare our theoretical results with data on RNA secondary structures from vertebrate introns.

Strong epistatic selection on the RNA secondary structure of HIV

A key question in evolutionary genomics is how populations navigate the adaptive landscape in the presence of epistasis, or interactions among loci. This problem can be directly addressed by studying the evolution of RNA secondary structures, for which there is constraint to maintain pairing between Watson-Crick (WC) sites. Replacement of a nucleotide at one site of a WC pair reduces fitness by disrupting binding, which can be restored via a compensatory replacement at the interacting site. Here, I present the first genome-scale analysis of epistasis on the RNA secondary structure of human immunodeficiency virus type 1 (HIV-1). Comparison of polymorphism frequencies at ancestrally conserved sites reveals that selection against replacements is ∼2.7 times stronger at WC than at non-WC sites, such that nearly 50% of constraint can be attributed to epistasis. However, almost all epistatic constraint is due to selection against conversions of WC pairs to unpaired (UP) nucleotides, whereas conversions to GU wobbles are only slightly deleterious. This disparity is also evident in pairs with second-site compensatory replacements; conversions from UP nucleotides to WC pairs increase median fitness by ∼4.2%, whereas conversions from GU wobbles to WC pairs only increase median fitness by ∼0.3%. Moreover, second-site replacements that convert UP nucleotides to GU wobbles also increase median fitness by ∼4%, indicating that such replacements are nearly as compensatory as those that restore WC pairing. Thus, WC peaks of the HIV-1 epistatic adaptive landscape are connected by high GU ridges, enabling the viral population to rapidly explore distant peaks without traversing deep UP valleys.

Epistasis is an evolutionary process in which the effect of a nucleotide at one site in the genome is dependent on the presence or absence of particular nucleotides at other sites in the genome. One of the simplest types of epistasis occurs between Watson-Crick (WC) nucleotides in RNA secondary structures, which are under constraint to maintain base-pairing. In this study, I examine the effects of mutations at WC sites in the RNA secondary structure of HIV-1. I show that while epistasis plays a major role in the evolution of the HIV-1 secondary structure, different types of mutations have variable effects on fitness. Therefore, by favoring certain mutational trajectories, HIV-1 can evolve rapidly despite strong epistatic constraint on its RNA secondary structure.

Distribution of conventional and nonconventional introns in tubulin (α and β) genes of euglenids

The nuclear genomes of euglenids contain three types of introns: conventional spliceosomal introns, nonconventional introns for which a splicing mechanism is unknown (variable noncanonical borders, RNA secondary structure bringing together intron ends), and so-called intermediate introns, which combine features of conventional and nonconventional introns. Analysis of two genes, tubA and tubB, from 20 species of euglenids reveals contrasting distribution patterns of conventional and nonconventional introns--positions of conventional introns are conserved, whereas those of the nonconventional ones are unique to individual species or small groups of closely related taxa. Moreover, in the group of phototrophic euglenids, 11 events of conventional intron loss versus 15 events of nonconventional intron gain were identified. A comparison of all nonconventional intron sequences highlighted the most conserved elements in their sequence and secondary structure. Our results led us to put forward two hypotheses. 1) The first one posits that mutational changes in intron sequence could lead to a change in their excision mechanism--intermediate introns would then be a transitional form between the conventional and nonconventional introns. 2) The second hypothesis concerns the origin of nonconventional introns--because of the presence of inverted repeats near their ends, insertion of MITE-like transposon elements is proposed as a possible source of new introns.

Coevolution in RNA molecules driven by selective constraints: evidence from 5S rRNA

Understanding intra-molecular coevolution helps to elucidate various structural and functional constraints acting on molecules and might have practical applications in predicting molecular structure and interactions. In this study, we used 5S rRNA as a template to investigate how selective constraints have shaped the RNA evolution. We have observed the nonrandom occurrence of paired differences along the phylogenetic trees, the high rate of compensatory evolution, and the high TIR scores (the ratio of the numbers of terminal to intermediate states), all of which indicate that significant positive selection has driven the evolution of 5S rRNA. We found three mechanisms of compensatory evolution: Watson-Crick interaction (the primary one), complex interactions between multiple sites within a stem, and interplay of stems and loops. Coevolutionary interactions between sites were observed to be highly dependent on the structural and functional environment in which they occurred. Coevolution occurred mostly in those sites closest to loops or bulges within structurally or functionally important helices, which may be under weaker selective constraints than other stem positions. Breaking these pairs would directly increase the size of the adjoining loop or bulge, causing a partial or total structural rearrangement. In conclusion, our results indicate that sequence coevolution is a direct result of maintaining optimal structural and functional integrity.

Association between WDR21A polymorphisms and airway responsiveness to inhaled corticosteroids in asthmatic patients

OBJECTIVE

Genetic polymorphism is partially responsible for the wide variation in the response of moderate-to-severe asthmatic patients to inhaled corticosteroids. The goal of the study was to examine polymorphisms in WDR21A, which encodes a putative glucocorticoid receptor (GR)-interacting protein, for their possible association with corticosteroid responsiveness.

METHODS

The change in forced expiratory volume in 1 s [FEV(1) (ΔFEV(1))] induced by 4 weeks of inhaled treatment with fluticasone propionate (1000 µg daily) was measured in 230 asthmatic patients. Fifteen single nucleotide polymorphisms (SNPs) of WDR21A were genotyped using a TaqMan assay, and 11 SNPs were used for further analysis. WDR21A transcripts were analyzed for variant splicing using reverse transcriptase-PCR. The WDR21A protein structure was predicted using a template-based modeling method and docked to a GR using Zdock.

RESULTS

Of the 11 SNPs and three haplotypes of WDR21A analyzed, only the intronic SNP -10073G>C appeared to affect ΔFEV(1). The ΔFEV(1) of the -10073C/C homozygous genotype was twice that of the -10073G/G and -10073C/G genotypes (P(corr)=0.04 in recessive model). No splicing variant of WDR21A was observed, regardless of genotype. The predicted WDR21A protein structure was similar to the Gβ(1) protein structure (template modeling-score=0.93).

CONCLUSION

The minor allele -10073C of WDR21A may induce a good response to inhaled corticosteroids possibly through competition with the Gβ(1) proteins for binding to GRs.

Evidence for widespread association of mammalian splicing and conserved long-range RNA structures

Pre-mRNA structure impacts many cellular processes, including splicing in genes associated with disease. The contemporary paradigm of RNA structure prediction is biased toward secondary structures that occur within short ranges of pre-mRNA, although long-range base-pairings are known to be at least as important. Recently, we developed an efficient method for detecting conserved RNA structures on the genome-wide scale, one that does not require multiple sequence alignments and works equally well for the detection of local and long-range base-pairings. Using an enhanced method that detects base-pairings at all possible combinations of splice sites within each gene, we now report RNA structures that could be involved in the regulation of splicing in mammals. Statistically, we demonstrate strong association between the occurrence of conserved RNA structures and alternative splicing, where local RNA structures are generally more frequent at alternative donor splice sites, while long-range structures are more associated with weak alternative acceptor splice sites. As an example, we validated the RNA structure in the human SF1 gene using minigenes in the HEK293 cell line. Point mutations that disrupted the base-pairing of two complementary boxes between exons 9 and 10 of this gene altered the splicing pattern, while the compensatory mutations that reestablished the base-pairing reverted splicing to that of the wild-type. There is statistical evidence for a Dscam-like class of mammalian genes, in which mutually exclusive RNA structures control mutually exclusive alternative splicing. In sum, we propose that long-range base-pairings carry an important, yet unconsidered part of the splicing code, and that, even by modest estimates, there must be thousands of such potentially regulatory structures conserved throughout the evolutionary history of mammals.

DNApol-ϵ gene is indispensable for the survival and growth of Drosophila melanogaster

Based on deletion and complementation mapping and DNA sequencing, a new recessive fully penetrant mutation (DNApol-ϵpl10R), causing prolonged larval life and larval/early pupal lethality, is identified as the first mutant allele of the DNApol-ϵ (CG6768) gene of Drosophila melanogaster. A same-sense base pair substitution in exon 1 of the DNApol-ϵ gene is associated with retention of the first intron and significant reduction in DNApol-ϵ transcripts in DNApol-ϵpl10R homozygotes. Homozygous mutant larvae show small imaginal discs with fewer cells and reduced polyteny in salivary glands, presumably because of the compromised DNA polymerase function following exhaustion of the maternal contribution. Extremely small and rare DNApol-ϵpl10R homozygous somatic clones in DNApol-ϵpl10R/+imaginal discs confirm their poor mitotic activity. The DNApol-ϵpl10R homozygotes, like those expressing DNApol-ϵ-RNAi transgene, show high sensitivity to DNA damaging agents. The first mutant allele of the DNApol-ϵ gene will facilitate functional characterization of this enzyme in the genetically tractable Drosophila model.

Mitochondrial-nuclear interactions and accelerated compensatory evolution: evidence from the primate cytochrome C oxidase complex

Accelerated rates of mitochondrial protein evolution have been proposed to reflect Darwinian coadaptation for efficient energy production for mammalian flight and brain activity. However, several features of mammalian mtDNA (absence of recombination, small effective population size, and high mutation rate) promote genome degradation through the accumulation of weakly deleterious mutations. Here, we present evidence for "compensatory" adaptive substitutions in nuclear DNA- (nDNA) encoded mitochondrial proteins to prevent fitness decline in primate mitochondrial protein complexes. We show that high mutation rate and small effective population size, key features of primate mitochondrial genomes, can accelerate compensatory adaptive evolution in nDNA-encoded genes. We combine phylogenetic information and the 3D structure of the cytochrome c oxidase (COX) complex to test for accelerated compensatory changes among interacting sites. Physical interactions among mtDNA- and nDNA-encoded components are critical in COX evolution; amino acids in close physical proximity in the 3D structure show a strong tendency for correlated evolution among lineages. Only nuclear-encoded components of COX show evidence for positive selection and adaptive nDNA-encoded changes tend to follow mtDNA-encoded amino acid changes at nearby sites in the 3D structure. This bias in the temporal order of substitutions supports compensatory weak selection as a major factor in accelerated primate COX evolution.

The role of the effective population size in compensatory evolution

The impact of the effective population size (N_e) on the efficacy of selection has been the focus of many theoretical and empirical studies over the recent years. Yet, the effect of N_e on evolution under epistatic fitness interactions is not well understood. In this study, we compare selective constraints at independently evolving (unpaired) and coevolving (paired) sites in orthologous transfer RNAs (tRNA molecules for vertebrate and drosophilid species pairs of different N_e. We show that patterns of nucleotide variation for the two classes of sites are explained well by Kimura's one- and two-locus models of sequence evolution under mutational pressure. We find that constraints in orthologous tRNAs increase with increasing N_e of the investigated species pair. Thereby, the effect of N_e on the efficacy of selection is stronger at unpaired sites than at paired sites. Furthermore, we identify a “core” set of tRNAs with high structural similarity to tRNAs from all major kingdoms of life and a “peripheral” set with lower similarity. We observe that tRNAs in the former set are subject to higher constraints and less prone to the effect of N_e, whereas constraints in tRNAs of the latter set show a large influence of N_e. Finally, we are able to demonstrate that constraints are relaxed in X-linked drosophilid tRNAs compared with autosomal tRNAs and suggest that N_e is responsible for this difference. The observed effects of N_e are consistent with the hypothesis that evolution of most tRNAs is governed by slightly to moderately deleterious mutations (i.e., |N_es| ≤ 5).

Association analysis of RGS7BP gene polymorphisms with aspirin intolerance in asthmatic patients

BACKGROUND

Signal-regulated palmitoylation of RGS7BP(regulator of G-protein-signaling 7-binding protein) initiates the activation of G-protein-coupled receptors (GPCRs), including muscarinic receptors, which contribute to the development of asthma and its subphenotypes.

OBJECTIVE

To determine the association of RGS7BP gene polymorphisms with the development of aspirin-exacerbated respiratory disease (AERD).

METHODS

We evaluated the association of RGS7BP gene polymorphisms with response to oral aspirin challenge and with responsiveness to methacholine challenge. RGS7BP messenger RNA splice variants in peripheral blood platelets from patients with different single-nucleotide polymorphisms were analyzed by reverse-transcription polymerase chain reaction.

RESULTS

Logistic regression analysis of RGS7BP gene polymorphisms in patients with AERD (n = 102) and aspirin-tolerant asthma (n = 429) revealed that a haplotype of block 3 consisting of rare alleles +98092 C>G, +98853 C>T, and +104450 T>G of the RGS7BP gene was associated with AERD. Multiple linear regression analysis showed that asthmatic patients carrying ht2/ht2 in block 3 were more responsive to aspirin challenge than those not carrying ht2 (P = .008 in a codominant model). The log-transformed provocation concentration that caused a decrease in forced expiratory volume in 1 second of 20% for methacholine was significantly dependent on the BL3-ht2 haplotype. No significant differences in platelet expression of different RGS7BP messenger RNA splice variants were detected between those with and without the BL3-ht2 haplotype.

CONCLUSION

BL3-ht2 of RGS7BP may be an important genetic variant associated with AERD. The haplotype of block 3 may play a protective role against aspirin hypersensitivity in asthma, perhaps by altering the responsiveness of muscarinic receptors.

Selective constraints in conserved folded RNAs of drosophilid and hominid genomes

Small noncoding RNAs as well as folded RNA structures in genic regions are crucial for many cellular processes. They are involved in posttranscriptional gene regulation (microRNAs), RNA modification (small nucleolar RNAs), regulation of splicing, correct localization of proteins, and many other processes. In most cases, a distinct secondary structure of the molecule is necessary for its correct function. Hence, selection should act to retain the structure of the molecule, although the underlying sequence is allowed to vary. Here, we present the first genome-wide estimates of selective constraints in folded RNA molecules in the nuclear genomes of drosophilids and hominids. In comparison to putatively neutrally evolving sites, we observe substantially reduced rates of substitutions at paired and unpaired sites of folded molecules. We estimated evolutionary constraints to be in the ranges of (0.974,0.991) and (0.895,1.000) for paired nucleotides in drosophilids and hominids, respectively. These values are significantly higher than for constraints at nonsynonymous sites of protein-coding genes in both genera. Nonetheless, valleys of only moderately reduced fitness (s ≈ 10(-4)) are sufficient to generate the observed fraction of nucleotide changes that are removed by purifying selection. In addition, a comparison of selective coefficients between drosophilids and hominids revealed significantly higher constraints in drosophilids, which can be attributed to the difference in long-term effective population size between these two groups of species. This difference is particularly apparent at the independently evolving (unpaired) sites.

Association between colony-stimulating factor 1 receptor gene polymorphisms and asthma risk

Colony-stimulating factor 1 receptor (CSF1R) is expressed in monocytes/macrophages and dendritic cells. These cells play important roles in the innate immune response, which is regarded as an important aspect of asthma development. Genetic alterations in the CSF1R gene may contribute to the development of asthma. We investigated whether CSF1R gene polymorphisms were associated with the risk of asthma. Through direct DNA sequencing of the CSF1R gene, we identified 28 single nucleotide polymorphisms (SNPs) and genotyped them in 303 normal controls and 498 asthmatic patients. Expression of CSF1R protein and mRNA were measured on CD14-positive monocytes and neutrophils in peripheral blood of asthmatic patients using flow cytometry and real-time PCR. Among the 28 polymorphisms, two intronic polymorphism (+20511C>T and +22693T>C) were associated with the risk of asthma by logistic regression analysis. The frequencies of the minor allele at CSF1R +20511C>T and +22693T>C were higher in asthmatic subjects than in normal controls (4.6 vs. 7.7%, p = 0.001 in co-dominant and dominant models; 16.4 vs. 25.8%, p = 0.0006 in a recessive model). CSF1R mRNA levels in neutrophils of the asthmatic patients having the +22693CC allele were higher than in those having the +22693TT allele (p = 0.026). Asthmatic patients with the +22693CC allele also showed significantly higher CSF1R expression on CD14-positive monocytes and neutrophils than did those with the +22693TT allele (p = 0.045 and p = 0.044). The +20511C>T SNP had no association with CSF1R mRNA or protein expression. In conclusion, the minor allele at CSF1R +22693T>C may have a susceptibility effect in the development of asthma, via increased CSF1R protein and mRNA expression in inflammatory cells.

Compensatory evolution in mitochondrial tRNAs navigates valleys of low fitness

A long-standing controversy in evolutionary biology is whether or not evolving lineages can cross valleys on the fitness landscape that correspond to low-fitness genotypes, which can eventually enable them to reach isolated fitness peaks. Here we study the fitness landscapes traversed by switches between different AU and GC Watson-Crick nucleotide pairs at complementary sites of mitochondrial transfer RNA stem regions in 83 mammalian species. We find that such Watson-Crick switches occur 30-40 times more slowly than pairs of neutral substitutions, and that alleles corresponding to GU and AC non-Watson-Crick intermediate states segregate within human populations at low frequencies, similar to those of non-synonymous alleles. Substitutions leading to a Watson-Crick switch are strongly correlated, especially in mitochondrial tRNAs encoded on the GT-nucleotide-rich strand of the mitochondrial genome. Using these data we estimate that a typical Watson-Crick switch involves crossing a fitness valley of a depth of about 10(-3) or even about 10(-2), with AC intermediates being slightly more deleterious than GU intermediates. This compensatory evolution must proceed through rare intermediate variants that never reach fixation. The ubiquitous nature of compensatory evolution in mammalian mitochondrial tRNAs and other molecules implies that simultaneous fixation of two alleles that are individually deleterious may be a common phenomenon at the molecular level.

Experimentally increased codon bias in the Drosophila Adh gene leads to an increase in larval, but not adult, alcohol dehydrogenase activity

Although most amino acids can be encoded by more than one codon, the synonymous codons are not used with equal frequency. This phenomenon is known as codon bias and appears to be a universal feature of genomes. The translational selection hypothesis posits that the use of optimal codons, which match the most abundant species of isoaccepting tRNAs, results in increased translational efficiency and accuracy. Previous work demonstrated that the experimental reduction of codon bias in the Drosophila alcohol dehydrogenase (Adh) gene led to a significant decrease in ADH protein expression. In this study we performed the converse experiment: we replaced seven suboptimal leucine codons that occur naturally in the Drosophila melanogaster Adh gene with the optimal codon. We then compared the in vivo ADH activities imparted by the wild-type and mutant alleles. The introduction of optimal leucine codons led to an increase in ADH activity in third-instar larvae. In adult flies, however, the introduction of optimal codons led to a decrease in ADH activity. There is no evidence that other selectively constrained features of the Adh gene, or its rate of transcription, were altered by the synonymous replacements. These results are consistent with translational selection for codon bias being stronger in the larval stage and suggest that there may be a selective conflict over optimal codon usage between different developmental stages.

Role of RNA structure in regulating pre-mRNA splicing

Pre-mRNA splicing involves removing non-coding introns from RNA transcripts. It is carried out by the spliceosome, along with other auxiliary factors. In general, research in splicing has focused on the sequences within the pre-mRNA, without considering the structures that these sequences might form. We propose that the role of RNA structure deserves more consideration when thinking about splicing mechanisms. RNA structures can inhibit or aid binding of spliceosomal components to the pre-mRNA, or can increase splicing efficiency by bringing important sequences into close proximity. Recent reports have identified proteins and small molecules that can regulate splicing by modulating RNA structures, thereby expanding our knowledge of the mechanisms used to regulate splicing.

Patterns of DNA-sequence divergence between Drosophila miranda and D. pseudoobscura

Contrary to the classical view, a large amount of non-coding DNA seems to be selectively constrained in Drosophila and other species. Here, using Drosophila miranda BAC sequences and the Drosophila pseudoobscura genome sequence, we aligned coding and non-coding sequences between D. pseudoobscura and D. miranda, and investigated their patterns of evolution. We found two patterns that have previously been observed in comparisons between Drosophila melanogaster and its relatives. First, there is a negative correlation between intron divergence and intron length, suggesting that longer non-coding sequences may contain more regulatory elements than shorter sequences. Our other main finding is a negative correlation between the rate of non-synonymous substitutions (d(N)) and codon usage bias (F(op)), showing that fast-evolving genes have a lower codon usage bias, consistent with strong positive selection interfering with weak selection for codon usage.

Epistatic interactions modulate the evolution of mammalian mitochondrial respiratory complex components

BACKGROUND

The deleterious effect of a mutation can be reverted by a second-site interacting residue. This is an epistatic compensatory process explaining why mutations that are deleterious in some species are tolerated in phylogenetically related lineages, rendering evident that those mutations are, by all means, only deleterious in the species-specific context. Although an extensive and refined theoretical framework on compensatory evolution does exist, the supporting evidence remains limited, especially for protein models. In this current study, we focused on the molecular mechanism underlying the epistatic compensatory process in mammalian mitochondrial OXPHOS proteins using a combination of in-depth structural and sequence analyses.

RESULTS

Modeled human structures were used in this study to predict the structural impairment and recovery of deleterious mutations alone and combined with an interacting compensatory partner, respectively. In two cases, COI and COIII, intramolecular interactions between spatially linked residues restore the folding pattern impaired by the deleterious mutation. In a third case, intermolecular contact between mitochondrial CYB and nuclear CYT1 encoded components of the cytochrome bc1 complex are likely to restore protein binding. Moreover, we observed different modes of compensatory evolution that have resulted in either a quasi-simultaneous occurrence of a mutation and corresponding compensatory partner, or in independent occurrences of mutations in distinct lineages that were always preceded by the compensatory site.

CONCLUSION

Epistatic interactions between individual replacements involving deleterious mutations seems to follow a parsimonious model of evolution in which genomes hold pre-compensating states that subsequently tolerate deleterious mutations. This phenomenon is likely to have been constraining the variability at coevolving sites and shaping the interaction between the mitochondrial and the nuclear genome.

Modulation of alternative splicing by long-range RNA structures in Drosophila

Accurate and efficient recognition of splice sites during pre-mRNA splicing is essential for proper transcriptome expression. Splice site usage can be modulated by secondary structures, but it is unclear if this type of modulation is commonly used or occurs to a significant degree with secondary structures forming over long distances. Using phlyogenetic comparisons of intronic sequences among 12 Drosophila genomes, we elucidated a group of 202 highly conserved pairs of sequences, each at least nine nucleotides long, capable of forming stable stem structures. This set was highly enriched in alternatively spliced introns and introns with weak acceptor sites and long introns, and most occurred over long distances (>150 nucleotides). Experimentally, we analyzed the splicing of several of these introns using mini-genes in Drosophila S2 cells. Wild-type splicing patterns were changed by mutations that opened the stem structure, and restored by compensatory mutations that re-established the base-pairing potential, demonstrating that these secondary structures were indeed implicated in the splice site choice. Mechanistically, the RNA structures masked splice sites, brought together distant splice sites and/or looped out introns. Thus, base-pairing interactions within introns, even those occurring over long distances, are more frequent modulators of alternative splicing than is currently assumed.

Discovery of a long inverted repeat in human POTE genes

POTE gene family is tightly related to prostate, ovary, testis and placenta cancers. We recently identified an intronic long inverted repeat (LIR) in some members of the POTE gene family. Due to the capacity of inducing gene amplification, the POTE intronic LIRs may be involved in over-expression of the POTE genes. Our study aimed to understand the origin of the LIR in primates. We collected the LIR and its flanking sequences within rhesus monkey, chimpanzee and human genomes. The rhesus monkey genome only has half-sized LIRs (lack one repeat copy), whereas the human and chimpanzee genomes contain both full-sized and half-sized LIRs. Phylogenetic tree indicates that the LIR is formed after divergence of rhesus monkey and the common ancestor of human and chimpanzee. The POTE genes containing a full-sized LIR were amplified in the human genome.

Acquisition of inverted GSTM exons by an intron of primate GSTM5 gene

The human GSTM gene family is composed of five gene members, GSTM1-5, and plays an important role in detoxification. In this study, the human GSTM5 gene was found to have a long inverted repeat (LIR) in intron 5. The LIR is able to form a stem-loop structure with a 31-bp stem and a 9-nt loop. The intronic LIR was also identified in other primates but not in non-primates. The human and chimpanzee LIRs had undergone compensating mutations that make the stem loop more stable, suggesting a functional role for the LIR. Sequence homology showed that the LIR was actually a part of inverted exons acquired by the intron. Results of phylogenetic analysis indicate that the inverted exons were derived from exon 5 of GSTM4 and exon 5 of GSTM1. The intronic LIR and inverted GSTM exons can probably introduce complexity in the expression of GSTM gene family.

A study on genomic distribution and sequence features of human long inverted repeats reveals species-specific intronic inverted repeats

The inverted repeats present in a genome play dual roles. They can induce genomic instability and, on the other hand, regulate gene expression. In the present study, we report the distribution and sequence features of recombinogenic long inverted repeats (LIRs) that are capable of forming stable stem-loops or palindromes within the human genome. A total of 2551 LIRs were identified, and 37% of them were located in long introns (largely > 10 kb) of genes. Their distribution appears to be random in introns and is not restrictive, even for regions near intron-exon boundaries. Almost half of them comprise TG/CA-rich repeats, inversely arranged Alu repeats and MADE1 mariners. The remaining LIRs are mostly unique in their sequence features. Comparative studies of human, chimpanzee, rhesus monkey and mouse orthologous genes reveal that human genes have more recombinogenic LIRs than other orthologs, and over 80% are human-specific. The human genes associated with the human-specific LIRs are involved in the pathways of cell communication, development and the nervous system, as based on significantly over-represented Gene Ontology terms. The functional pathways related to the development and functions of the nervous system are not enriched in chimpanzee and mouse orthologs. The findings of the present study provide insight into the role of intronic LIRs in gene regulation and primate speciation.

Analyzing the evolution of RNA secondary structures in vertebrate introns using Kimura's model of compensatory fitness interactions

Previous studies have shown that splicing efficiency, and thus maturation of pre-mRNA, depends on the correct folding of the RNA molecule into a secondary or higher order structure. When disrupted by a mutation, aberrant folding may result in a lower splicing efficiency. However, the structure can be restored by a second, compensatory mutation. Here, we present a logistic regression approach to analyze the evolutionary dynamics of RNA secondary structures. We apply our approach to a set of computationally predicted RNA secondary structures in vertebrate introns. Our results are consistent with the hypothesis of a negative influence of the physical distance between pairing nucleotides on the occurrence of covariations, as predicted by Kimura's model of compensatory evolution. We also confirm the hypothesis that longer local secondary structure elements (helices) can accommodate a larger number of covariations, wobbles, and mismatches. Furthermore, we find that wobbles and mismatches are more frequent in the middle of a helix, whereas covariations occur preferentially at the helix ends. The GC content is a major determinant of this pattern.

Longer first introns are a general property of eukaryotic gene structure

While many properties of eukaryotic gene structure are well characterized, differences in the form and function of introns that occur at different positions within a transcript are less well understood. In particular, the dynamics of intron length variation with respect to intron position has received relatively little attention. This study analyzes all available data on intron lengths in GenBank and finds a significant trend of increased length in first introns throughout a wide range of species. This trend was found to be even stronger when using high-confidence gene annotation data for three model organisms (Arabidopsis thaliana, Caenorhabditis elegans, and Drosophila melanogaster) which show that the first intron in the 5' UTR is--on average--significantly longer than all downstream introns within a gene. A partial explanation for increased first intron length in A. thaliana is suggested by the increased frequency of certain motifs that are present in first introns. The phenomenon of longer first introns can potentially be used to improve gene prediction software and also to detect errors in existing gene annotations.

Determining the relationship of gene expression and global mRNA stability in Drosophila melanogaster and Escherichia coli using linear models

There are several sequence-dependent factors regulating gene expression. Some of them have been extensively studied, among the most prominent are GC content and codon usage bias. Other factors hypothesized to have an impact on gene expression are gene length and the thermodynamic stability of mRNA secondary structure. In this work, we analyzed two different microarray datasets of Drosophila melanogaster gene expression and one dataset of Escherichia coli. To investigate the relationship between gene expression, codon usage bias and GC content of first, second and third codon position, gene length and mRNA stability we employed a multiple regression analysis using a comprehensive linear model. It is shown that codon usage bias and GC content of the first, second and third codon position show a significant influence on gene expression, whereas no significant effect of mRNA secondary structure stability is observed.

Conflicting selection pressures on synonymous codon use in yeast suggest selection on mRNA secondary structures

BACKGROUND

Eukaryotic mRNAs often contain secondary structures in their untranslated regions that are involved in expression regulation. Whether secondary structures in the protein coding regions are of functional importance remains unclear: laboratory studies suggest stable secondary structures within the protein coding sequence interfere with translation, while several bioinformatic studies indicate stable mRNA structures are more frequent than expected.

RESULTS

In contrast to several studies testing for unexpected structural stabilities, I directly compare the selective constraint of sites that differ in their structural importance. I.e. for each nucleotide, I identify whether it is paired with another nucleotide, or unpaired, in the predicted secondary structure. I assume paired sites are more important for the predicted secondary structure than unpaired sites. I look at protein coding yeast sequences and use optimal codons and synonymous substitutions to test for structural constraints. As expected under selection for secondary structures, paired sites experience higher constraint than unpaired sites, i.e. significantly lower numbers of conserved optimal codons and consistently lower numbers of synonymous substitutions. This is true for structures predicted by different algorithms.

CONCLUSION

The results of this study are consistent with purifying selection on mRNA secondary structures in yeast protein coding sequences and suggest their biological importance. One should be aware, however, that accuracy of structure prediction is unknown for mRNAs and interrelated selective forces may contribute as well. Note that if selection pressures alternative to translational selection affect synonymous (and optimal) codon use, this may lead to under- or over-estimates of selective strength on optimal codon use depending on strength and direction of translational selection.

Conserved secondary structures in Aspergillus

Background

Recent evidence suggests that the number and variety of functional RNAs (ncRNAs as well as cis-acting RNA elements within mRNAs ) is much higher than previously thought; thus, the ability to computationally predict and analyze RNAs has taken on new importance. We have computationally studied the secondary structures in an alignment of six Aspergillus genomes. Little is known about the RNAs present in this set of fungi, and this diverse set of genomes has an optimal level of sequence conservation for observing the correlated evolution of base-pairs seen in RNAs.

Methodology/Principal Findings

We report the results of a whole-genome search for evolutionarily conserved secondary structures, as well as the results of clustering these predicted secondary structures by structural similarity. We find a total of 7450 predicted secondary structures, including a new predicted ∼60 bp long hairpin motif found primarily inside introns. We find no evidence for microRNAs. Different types of genomic regions are over-represented in different classes of predicted secondary structures. Exons contain the longest motifs (primarily long, branched hairpins), 5′ UTRs primarily contain groupings of short hairpins located near the start codon, and 3′ UTRs contain very little secondary structure compared to other regions. There is a large concentration of short hairpins just inside the boundaries of exons. The density of predicted intronic RNAs increases with the length of introns, and the density of predicted secondary structures within mRNA coding regions increases with the number of introns in a gene.

Conclusions/Sigificance

There are many conserved, high-confidence RNAs of unknown function in these Aspergillus genomes, as well as interesting spatial distributions of predicted secondary structures. This study increases our knowledge of secondary structure in these aspergillus organisms.

Correlation between the secondary structure of pre-mRNA introns and the efficiency of splicing in Saccharomyces cerevisiae

BACKGROUND

Secondary structure interactions within introns have been shown to be essential for efficient splicing of several yeast genes. The nature of these base-pairing interactions and their effect on splicing efficiency were most extensively studied in ribosomal protein gene RPS17B (previously known as RP51B). It was determined that complementary pairing between two sequence segments located downstream of the 5' splice site and upstream of the branchpoint sequence promotes efficient splicing of the RPS17B pre-mRNA, presumably by shortening the branchpoint distance. However, no attempts were made to compute a shortened, 'structural' branchpoint distance and thus the functional relationship between this distance and the splicing efficiency remains unknown.

RESULTS

In this paper we use computational RNA secondary structure prediction to analyze the secondary structure of the RPS17B intron. We show that it is necessary to consider suboptimal structure predictions and to compute the structural branchpoint distances in order to explain previously published splicing efficiency results. Our study reveals that there is a tight correlation between this distance and splicing efficiency levels of intron mutants described in the literature. We experimentally test this correlation on additional RPS17B mutants and intron mutants within two other yeast genes.

CONCLUSION

The proposed model of secondary structure requirements for efficient splicing is the first attempt to specify the functional relationship between pre-mRNA secondary structure and splicing. Our findings provide further insights into the role of pre-mRNA secondary structure in gene splicing in yeast and also offer basis for improvement of computational methods for splice site identification and gene-finding.

Positive and negative selection on noncoding DNA in Drosophila simulans

There is now a wealth of evidence that some of the most important regions of the genome are found outside those that encode proteins, and noncoding regions of the genome have been shown to be subject to substantial levels of selective constraint, particularly in Drosophila. Recent work has suggested that these regions may also have been subject to the action of positive selection, with large fractions of noncoding divergence having been driven to fixation by adaptive evolution. However, this work has focused on Drosophila melanogaster, which is thought to have experienced a reduction in effective population size (N(e)), and thus a reduction in the efficacy of selection, compared with its closest relative Drosophila simulans. Here, we examine patterns of evolution at several classes of noncoding DNA in D. simulans and find that all noncoding DNA is subject to the action of negative selection, indicated by reduced levels of polymorphism and divergence and a skew in the frequency spectrum toward rare variants. We find that the signature of negative selection on noncoding DNA and nonsynonymous sites is obscured to some extent by purifying selection acting on preferred to unpreferred synonymous codon mutations. We investigate the extent to which divergence in noncoding DNA is inferred to be the product of positive selection and to what extent these inferences depend on selection on synonymous sites and demography. Based on patterns of polymorphism and divergence for different classes of synonymous substitution, we find the divergence excess inferred in noncoding DNA and nonsynonymous sites in the D. simulans lineage difficult to reconcile with demographic explanations.

An Ancient Repeat Sequence in the ATP Synthase β-Subunit Gene of Forcipulate Sea Stars

A novel repeat sequence with a conserved secondary structure is described from two nonadjacent introns of the ATP synthase beta-subunit gene in sea stars of the order Forcipulatida (Echinodermata: Asteroidea). The repeat is present in both introns of all forcipulate sea stars examined, which suggests that it is an ancient feature of this gene (with an approximate age of 200 Mya). Both stem and loop regions show high levels of sequence constraint when compared to flanking nonrepetitive intronic regions. The repeat was also detected in (1) the family Pterasteridae, order Velatida and (2) the family Korethrasteridae, order Velatida. The repeat was not detected in (1) the family Echinasteridae, order Spinulosida, (2) the family Astropectinidae, order Paxillosida, (3) the family Solasteridae, order Velatida, or (4) the family Goniasteridae, order Valvatida. The repeat lacks similarity to published sequences in unrestricted GenBank searches, and there are no significant open reading frames in the repeat or in the flanking intron sequences. Comparison via parametric bootstrapping to a published phylogeny based on 4.2 kb of nuclear and mitochondrial sequence for a subset of these species allowed the null hypothesis of a congruent phylogeny to be rejected for each repeat, when compared separately to the published phylogeny. In contrast, the flanking nonrepetitive sequences in each intron yielded separate phylogenies that were each congruent with the published phylogeny. In four species, the repeat in one or both introns has apparently experienced gene conversion. The two introns also show a correlated pattern of nucleotide substitutions, even after excluding the putative cases of gene conversion.

RNA structure is a key regulatory element in pathological ATM and CFTR pseudoexon inclusion events

Genomic variations deep in the intronic regions of pre-mRNA molecules are increasingly reported to affect splicing events. However, there is no general explanation why apparently similar variations may have either no effect on splicing or cause significant splicing alterations. In this work we have examined the structural architecture of pseudoexons previously described in ATM and CFTR patients. The ATM case derives from the deletion of a repressor element and is characterized by an aberrant 5′ss selection despite the presence of better alternatives. The CFTR pseudoexon instead derives from the creation of a new 5′ss that is used while a nearby pre-existing donor-like sequence is never selected. Our results indicate that RNA structure is a major splicing regulatory factor in both cases. Furthermore, manipulation of the original RNA structures can lead to pseudoexon inclusion following the exposure of unused 5′ss already present in their wild-type intronic sequences and prevented to be recognized because of their location in RNA stem structures. Our data show that intrinsic structural features of introns must be taken into account to understand the mechanism of pseudoexon activation in genetic diseases. Our observations may help to improve diagnostics prediction programmes and eventual therapeutic targeting.

The iStem, a long-range RNA secondary structure element required for efficient exon inclusion in the Drosophila Dscam pre-mRNA

The Drosophila Dscam gene encodes 38,016 different proteins, due to alternative splicing of 95 of its 115 exons, that function in axon guidance and innate immunity. The alternative exons are organized into four clusters, and the exons within each cluster are spliced in a mutually exclusive manner. Here we describe an evolutionarily conserved RNA secondary structure we call the Inclusion Stem (iStem) that is required for efficient inclusion of all 12 variable exons in the exon 4 cluster. Although the iStem governs inclusion or exclusion of the entire exon 4 cluster, it does not play a significant role in determining which variable exon is selected. Thus, the iStem is a novel type of regulatory element that simultaneously controls the splicing of multiple alternative exons.

Patterns of intron sequence evolution in Drosophila are dependent upon length and GC content

BACKGROUND

Introns comprise a large fraction of eukaryotic genomes, yet little is known about their functional significance. Regulatory elements have been mapped to some introns, though these are believed to account for only a small fraction of genome wide intronic DNA. No consistent patterns have emerged from studies that have investigated general levels of evolutionary constraint in introns.

RESULTS

We examine the relationship between intron length and levels of evolutionary constraint by analyzing inter-specific divergence at 225 intron fragments in Drosophila melanogaster and Drosophila simulans, sampled from a broad distribution of intron lengths. We document a strongly negative correlation between intron length and divergence. Interestingly, we also find that divergence in introns is negatively correlated with GC content. This relationship does not account for the correlation between intron length and divergence, however, and may simply reflect local variation in mutational rates or biases.

CONCLUSION

Short introns make up only a small fraction of total intronic DNA in the genome. Our finding that long introns evolve more slowly than average implies that, while the majority of introns in the Drosophila genome may experience little or no selective constraint, most intronic DNA in the genome is likely to be evolving under considerable constraint. Our results suggest that functional elements may be ubiquitous within longer introns and that these introns may have a more general role in regulating gene expression than previously appreciated. Our finding that GC content and divergence are negatively correlated in introns has important implications for the interpretation of the correlation between divergence and levels of codon bias observed in Drosophila.

Exonization of AluYa5 in the human ACE gene requires mutations in both 3' and 5' splice sites and is facilitated by a conserved splicing enhancer

Ancient Alu elements have been shown to be included in mature transcripts by point mutations that improve their 5′ or 3′ splice sites. We have examined requirements for exonization of a younger, disease-associated AluYa5 in intron 16 of the human ACE gene. A single G>C transversion in position −3 of the new Alu exon was insufficient for Alu exonization and a significant inclusion in mRNA was only observed when improving several potential splice donor sites in the presence of 3′ CAG. Since complete Alu exonization was not achieved by optimizing traditional splicing signals, including the branch site, we tested whether auxiliary elements in AluYa5 were required for constitutive inclusion. Exonization was promoted by a SELEX-predicted heptamer in Alu consensus sequence 222–228 and point mutations in highly conserved nucleotides of this heptamer decreased Alu inclusion. In addition, we show that Alu exonization was facilitated by a subset of serine/arginine-rich (SR) proteins through activation of the optimized 3′ splice site. Finally, the haplotype- and allele-specific ACE minigenes generated similar splicing patterns in both ACE-expressing and non-expressing cells, suggesting that previously reported allelic association with plasma ACE activity and cardiovascular disease is not attributable to differential splicing of introns 16 and 17.

Functional analysis of Drosophila melanogaster gene regulatory sequences by transgene coplacement

The function of putative regulatory sequences identified by comparative genomics can be elucidated only through experimentation. Here the effectiveness of using heterologous gene constructs and transgene coplacement to characterize regulatory sequence function is demonstrated. This method shows that a sequence in the Adh 3'-untranslated region negatively regulates expression, independent of gene or chromosomal context.

Insertion/deletion and nucleotide polymorphism data reveal constraints in Drosophila melanogaster introns and intergenic regions

Our study of nucleotide sequence and insertion/deletion polymorphism in Drosophila melanogaster noncoding DNA provides evidence for selective pressures in both intergenic regions and introns (of the large size class). Intronic and intergenic sequences show a similar polymorphic deletion bias. Insertions have smaller sizes and higher frequencies than deletions, supporting the hypothesis that insertions are selected to compensate for the loss of DNA caused by deletion bias. Analysis of a simple model of selective constraints suggests that the blocks of functional elements located in intergenic sequences are on average larger than those in introns, while the length distribution of relatively unconstrained sequences interspaced between these blocks is similar in intronic and intergenic regions.

Effects of population structure and sex on association between serotonin receptors and Drosophila heart rate

As a first step toward population and quantitative genetic analysis of neurotransmitter receptors in Drosophila melanogaster, we describe the parameters of nucleotide variation in three serotonin receptors and their association with pupal heart rate. Thirteen kilobases of DNA including the complete coding regions of 5-HT1A, 5-HT1B, and 5-HT2 were sequenced in 216 highly inbred lines extracted from two North American populations in California and North Carolina. Nucleotide and amino acid polymorphism is in the normal range for Drosophila genes and proteins, and linkage disequilibrium decays rapidly such that haplotype blocks are typically only a few SNPs long. However, intron 1 of 5-HT1A consists of two haplotypes that are at significantly different frequencies in the two populations. Neither this region of the gene nor any of the common amino acid polymorphisms in the three loci associate with either heart rate or heart rate variability. A cluster of SNPs in intron 2 of 5-HT1A, including a triallelic site, do show a highly significant interaction between genotype, sex, and population. While it is likely that a combination of weak, complex selection pressures and population structure has helped shape variation in the serotonin receptors of Drosophila, much larger sampling strategies than are currently adopted in evolutionary genetics will be required to disentangle these effects.

The origin of recent introns: transposons?

A recent study makes considerable progress towards answering the question of how genes acquire introns by identifying numerous recently gained introns in nematodes.

The long-standing question of how genes acquire introns has provoked much debate. A recent study makes considerable progress by identifying numerous recently gained introns in nematodes - although it remains difficult to distinguish definitively between models of intron gain.