Maintenance of pre-mRNA secondary structure by epistatic selection.

Linkage disequilibrium between polymorphisms in a natural population may result from various evolutionary forces, including random genetic drift due to sampling of gametes during reproduction, restricted migration between subpopulations in a subdivided population, or epistatic selection. In this report, we present evidence that the majority of significant linkage disequilibria observed in introns of the alcohol dehydrogenase locus (Adh) of Drosophila pseudoobscura are due to epistatic selection maintaining secondary structure of precursor mRNA (pre-mRNA). Based on phylogenetic-comparative analysis and a likelihood approach, we propose secondary structure models of Adh pre-mRNA for the regions of the adult intron and intron 2 where clustering of linkage disequilibria has been observed. Furthermore, we applied the likelihood ratio test to the phylogenetically predicted secondary structure in intron 1. In contrast to the other two structures, polymorphisms associated with the more conserved stem-loop structure of intron 1 are in low frequency, and linkage disequilibria have not been observed. These findings are qualitatively consistent with a model of compensatory fitness interactions. This model assumes that mutations disrupting pairing in a secondary structural element are individually deleterious if they destabilize a functionally important structure; a second "compensatory" mutation, however, may restabilize the structure and restore fitness.

Identification and characterization of putative transposable DNA elements in solanaceous plants and Caenorhabditis elegans.

Several families of putative transposable elements (TrEs) in both solanaceous plants and Caenorhabditis elegans have been identified by screening the DNA data base for inverted repeated domains present in multiple copies in the genome. The elements are localized within intron and flanking regions of many genes. These elements consist of two inverted repeats flanking sequences ranging from 5 bp to > 500 bp. Identification of multiple elements in which sequence conservation includes both the flanking and internal regions implies that these TrEs are capable of duplicative transposition. Two of the elements were identified in promoter regions of the tomato (Lycoperiscon esculentum) polygalacturonase and potato (Solanum tuberosum) Win1 genes. The element in the polygalacturonase promoter spans a known regulatory region. In both cases, ancestral DNA sequences, which represent potential recombination target sequences prior to insertion of the elements, have been cloned from related species. The sequences of the inverted repeated domains in plants and C. elegans show a high degree of phylogenetic conservation. While frequency of the different elements is variable, some are present in very high copy number. A member of a single C. elegans TrE family is observed approximately once every 20 kb in the genome. The abundance of the described TrEs suggests utility in the genomic analysis of these and related organisms.

Evidence against the exon theory of genes derived from the triose-phosphate isomerase gene.

The exon theory of genes proposes that the introns of protein-encoding nuclear genes are remnants of the DNA spacers between ancient minigenes. The discovery of an intron at a predicted position in the triose-phosphate isomerase (EC 5.3.1.1) gene of Culex mosquitoes has been hailed as an evidential pillar of the theory. We have found that that intron is also present in Aedes mosquitoes, which are closely related to Culex, but not in the phylogenetically more distant Anopheles, nor in the fly Calliphora vicina, nor in the moth Spodoptera littoralis. The presence of this intron in Culex and Aedes is parsimoniously explained as the result of an insertion in a recent common ancestor of these two species rather than as the remnant of an ancient intron. The absence of the intron in 19 species of very diverse organisms requires at least 10 independent evolutionary losses in order to be consistent with the exon theory.

Rat growth hormone gene introns stimulate nucleosome alignment in vitro and in transgenic mice.

Average hepatic expression (mRNA per cell per gene) of a metallothionein-rat growth hormone (rGH) gene with its natural introns was about 15-fold higher than an intronless version when tested in transgenic mice. We examined the idea that intron removal leads to an alteration in chromatin structure that might be responsible for this effect. Using an in vitro chromatin assembly system, we observed that nucleosomes were aligned in a characteristic ordered array over the gene and promoter when all introns were present. Linker histones were necessary for this alignment to occur. In contrast, nucleosome alignment was perturbed in constructs lacking some or all of the introns. A similar disruption of nucleosome alignment was observed when comparing chromatin from livers of transgenic mice carrying rGH transgenes with or without introns. In vitro, sequences at the 3' end of the rGH gene position nucleosomes and facilitate nucleosome alignment upstream; however, nucleosome alignment does not occur on the approximately 3 kb of downstream flanking rat sequence. These observations suggest that signals present in genomic rGH DNA may serve to establish appropriate nucleosome alignment during development and, possibly, to restore nucleosome alignment to the transcribed region after disruption incurred by the passage of an RNA polymerase molecule, thereby facilitating subsequent rounds of transcription.

Structure, inheritance, and transcriptional effects of Pce1, an insertional element within Phanerochaete chrysosporium lignin peroxidase gene lipI.

A 1747-bp insertion within a lignin peroxidase allele of Phanerochaete chrysosporium BKM-F-1767 is described. Pce1, the element, lies immediately adjacent to the fourth intron of lip12. Southern blots reveal the presence of Pce1-homologous sequences in other P. chrysosporium strains. Transposon-like features include inverted terminal repeats and a dinucleotide (TA) target duplication. Atypical of transposons, Pce1 is present at very low copy numbers (one to five copies), and conserved transposase motifs are lacking. The mutation transcriptionally inactivates lip12 and is inherited in a 1:1 Mendelian fashion among haploid progeny. Thus, Pce1 is a transposon-like element that may play a significant role in generating ligninolytic variation in certain P. chrysosporium strains.

Considerations on the folding topology and evolutionary origin of cadherin domains.

Cell-cell adhesion in zonula adherens and desmosomal junctions is mediated by cadherins, and recent crystal structures of the first domain from murine N-cadherin provide a plausible molecular basis for this adhesive action. A structure-based sequence analysis of this adhesive domain indicates that its fold is common to all extracellular cadherin domains. The cadherin folding topology is also shown to be similar to immunoglobulin-like domains and to other Greek-key beta-sandwich structures, as diverse as domains from plant cytochromes, bacterial cellulases, and eukaryotic transcription factors. Sequence similarities between cadherins and these other molecules are very low, however, and intron patterns are also different. On balance, independent origins for a favorable folding topology seem more likely than evolutionary divergence from an ancestor common to cadherins and immunoglobulins.

kappa-Opioid receptor in humans: cDNA and genomic cloning, chromosomal assignment, functional expression, pharmacology, and expression pattern in the central nervous system.

Using the mouse delta-opioid receptor cDNA as a probe, we have isolated genomic clones encoding the human mu- and kappa-opioid receptor genes. Their organization appears similar to that of the human delta receptor gene, with exon-intron boundaries located after putative transmembrane domains 1 and 4. The kappa gene was mapped at position q11-12 in human chromosome 8. A full-length cDNA encoding the human kappa-opioid receptor has been isolated. The cloned receptor expressed in COS cells presents a typical kappa 1 pharmacological profile and is negatively coupled to adenylate cyclase. The expression of kappa-opioid receptor mRNA in human brain, as estimated by reverse transcription-polymerase chain reaction, is consistent with the involvement of kappa-opioid receptors in pain perception, neuroendocrine physiology, affective behavior, and cognition. In situ hybridization studies performed on human fetal spinal cord demonstrate the presence of the transcript specifically in lamina II of the dorsal horn. Some divergences in structural, pharmacological, and anatomical properties are noted between the cloned human and rodent receptors.

Persistence of immunoglobulin heavy chain/c-myc recombination-positive lymphocyte clones in the blood of human immunodeficiency virus-infected homosexual men.

We studied blood lymphocytes of human immunodeficiency virus (HIV)-seropositive and -negative homosexual men for the presence of T(8;14) translocations that recombine c-myc and immunoglobulin heavy-chain (IgH) mu/IgH alpha switch regions. Clones with T(8;14) translocations were detected in 10.5% (12/114) of the HIV-positive and in 2.0% of the 99 uninfected patients. The majority of recombinations were found at a single time point only. Four patients, however, harbored multiple (up to four) and persistent (up to 9 years) translocation-positive cell clones. No correlation between the presence of these aberrant lymphocytes and a later lymphoma could be established. The exon 1/intron 1 region of the recombined c-myc was investigated for the presence of point mutations and these were found in the nonpersistent clones. Additional alterations detected in these clones included duplications and a deletion in the c-myc gene. The pattern of base substitution indicates that they were introduced after the translocation event.

Pax-6 and lens-specific transcription of the chicken delta 1-crystallin gene.

The abundance of delta-crystallin in the chicken eye lens provides an advantageous marker for tissue-specific gene expression during cellular differentiation. The lens-specific expression of the delta 1-crystallin gene is governed by an enhancer in the third intron, which binds a positive (delta EF2) and negative (delta EF1) factor in its core region. Here we show by DNase I footprinting, electrophoretic mobility-shift assays, and cotransfection experiments with the delta 1-promoter/enhancer fused to the chloramphenicol acetyltransferase reporter gene that the delta 1-crystallin enhancer has two adjacent functional Pax-6 binding sites. We also demonstrate by DNase I footprinting that the delta EF1 site can bind the transcription factor USF, raising the possibility that USF may cooperate with Pax-6 in activation of the chicken delta 1- and alpha A-crystallin genes. These data, coupled with our recent demonstration that Pax-6 activates the alpha A-crystallin gene, suggest that Pax-6 may have been used extensively throughout evolution to recruit and express crystallin genes in the lens.

Isolation of a differentially regulated splicing isoform of human NF-E2.

The transcription factor NF-E2 (nuclear factor erythroid 2), interacting via DNA motifs within regulatory regions of several hematopoietic genes, is thought to mediate the enhancer activity of the globin locus control regions. By screening a human fetal liver cDNA library with probes derived from mouse NF-E2, we have isolated a splicing variant of the NF-E2 gene (fNF-E2) that differs in the 5' untranslated region from the previously reported cDNA (aNF-E2). The fNF-E2 isoform is transcribed from an alternative promoter located in the 3' end of the first intron and joined by alternative splicing to the second and third exons, which are shared by both RNA isoforms. Although the two forms produce the same protein, they are expressed in different ratios during development. fNF-E2 is more abundant in the fetal liver and less abundant in the adult bone marrow compared to the previously described form. Their distribution apparently follows the differential expression of fetal and adult hemoglobins.

Modulation of 5' splice site choice in pre-messenger RNA by two distinct steps.

Ser/Arg-rich proteins (SR proteins) are essential splicing factors that commit pre-messenger RNAs to splicing and also modulate 5' splice site choice in the presence or absence of functional U1 small nuclear ribonucleoproteins (snRNPs). Here, we perturbed the U1 snRNP in HeLa cell nuclear extract by detaching the U1-specific A protein using a 2'-O-methyl oligonucleotide (L2) complementary to its binding site in U1 RNA. In this extract, the standard adenovirus substrate is spliced normally, but excess amounts of SR proteins do not exclusively switch splicing from the normal 5' splice site to a proximal site (site 125 within the adenovirus intron), suggesting that modulation of 5' splice site choice exerted by SR proteins requires integrity of the U1 snRNP. The observation that splicing does not necessarily follow U1 binding indicates that interactions between the U1 snRNP and components assembled on the 3' splice site via SR proteins may also be critical for 5' splice site selection. Accordingly, we found that SR proteins promote the binding of the U2 snRNP to the branch site and stabilize the complex formed on a 3'-half substrate in the presence or absence of functional U1 snRNPs. A novel U2/U6/3'-half substrate crosslink was also detected and promoted by SR proteins. Our results suggest that SR proteins in collaboration with the U1 snRNP function in two distinct steps to modulate 5' splice site selection.

Fate of a redundant gamma-globin gene in the atelid clade of New World monkeys: implications concerning fetal globin gene expression.

Conclusive evidence was provided that gamma 1, the upstream of the two linked simian gamma-globin loci (5'-gamma 1-gamma 2-3'), is a pseudogene in a major group of New World monkeys. Sequence analysis of PCR-amplified genomic fragments of predicted sizes revealed that all extant genera of the platyrrhine family Atelidae [Lagothrix (woolly monkeys), Brachyteles (woolly spider monkeys), Ateles (spider monkeys), and Alouatta (howler monkeys)] share a large deletion that removed most of exon 2, all of intron 2 and exon 3, and much of the 3' flanking sequence of gamma 1. The fact that two functional gamma-globin genes were not present in early ancestors of the Atelidae (and that gamma 1 was the dispensible gene) suggests that for much or even all of their evolution, platyrrhines have had gamma 2 as the primary fetally expressed gamma-globin gene, in contrast to catarrhines (e.g., humans and chimpanzees) that have gamma 1 as the primary fetally expressed gamma-globin gene. Results from promoter sequences further suggest that all three platyrrhine families (Atelidae, Cebidae, and Pitheciidae) have gamma 2 rather than gamma 1 as their primary fetally expressed gamma-globin gene. The implications of this suggestion were explored in terms of how gene redundancy, regulatory mutations, and distance of each gamma-globin gene from the locus control region were possibly involved in the acquisition and maintenance of fetal, rather than embryonic, expression.

RNA-Seq expression profile of genes related to neurodegenerative disorders affecting the human retina

Human neurodegenerative diseases, such as Parkinson’s disease (PD) and the neuromuscular disorders called dystroglycanopathies (DGPs), cause retinal impairments. We have used RNA-Seq technology to catalog all known genes linked to PD and DGPs expressed in the human retina and quantitate their mRNA levels in terms of FPKM. We have also characterized their expression profiles in the retina by determining their exonic, intronic and exon-intron junction expression levels, as well as the alternative splicing pattern of particular genes. We believe these data could pave the way toward understanding the molecular bases of sight deficiencies associated with neurodegenerative disorders.

Genetic association of CD247 (CD3ζ) with SLE in a large-scale multiethnic study

A classic T-cell phenotype in systemic lupus erythematosus (SLE) is the downregulation and replacement of the CD3ζ chain that alters T-cell receptor signaling. However, genetic associations with SLE in the human CD247 locus that encodes CD3ζ are not well established and require replication in independent cohorts. Our aim was therefore to examine, localize and validate CD247-SLE association in a large multiethnic population. We typed 44 contiguous CD247 single-nucleotide polymorphisms (SNPs) in 8922 SLE patients and 8077 controls from four ethnically distinct populations. The strongest associations were found in the Asian population (11 SNPs in intron 1, 4.99 × 10(-4) < P < 4.15 × 10(-2)), where we further identified a five-marker haplotype (rs12141731-rs2949655-rs16859085-rs12144621-rs858554; G-G-A-G-A; P(hap) = 2.12 × 10(-5)) that exceeded the most associated single SNP rs858554 (minor allele frequency in controls = 13%; P = 4.99 × 10(-4), odds ratio = 1.32) in significance. Imputation and subsequent association analysis showed evidence of association (P < 0.05) at 27 additional SNPs within intron 1. Cross-ethnic meta-analysis, assuming an additive genetic model adjusted for population proportions, showed five SNPs with significant P-values (1.40 × 10(-3) < P< 3.97 × 10(-2)), with one (rs704848) remaining significant after Bonferroni correction (P(meta) = 2.66 × 10(-2)). Our study independently confirms and extends the association of SLE with CD247, which is shared by various autoimmune disorders and supports a common T-cell-mediated mechanism.

Next-generation sequencing of iron-metabolism related genes in Portuguese patients with iron overload: novel pathogenic genetic variants

Objective: In Southern European countries up to one-third of the patients with hereditary hemochromatosis (HH) do not present the common HFE risk genotype. In order to investigate the molecular basis of these cases we have designed a gene panel for rapid and simultaneous analysis of 6 HH-related genes (HFE, TFR2, HJV, HAMP, SLC40A1 and FTL) by next-generation sequencing (NGS). Materials and Methods: Eighty-eight iron overload Portuguese patients, negative for the common HFE mutations, were analysed. A TruSeq Custom Amplicon kit (TSCA, by Illumina) was designed in order to generate 97 amplicons covering exons, intron/exon junctions and UTRs of the mentioned genes with a cumulative target sequence of 12115bp. Amplicons were sequenced in the MiSeq instrument (IIlumina) using 250bp paired-end reads. Sequences were aligned against human genome reference hg19 using alignment and variant caller algorithms in the MiSeq reporter software. Novel variants were validated by Sanger sequencing and their pathogenic significance were assessed by in silico studies. Results: We found a total of 55 different genetic variants. These include novel pathogenic missense and splicing variants (in HFE and TFR2), a very rare variant in IRE of FTL, a variant that originates a novel translation initiation codon in the HAMP gene, among others. Conclusion: The merging of TSCA methodology and NGS technology appears to be an appropriate tool for simultaneous and fast analysis of HH-related genes in a large number of samples. However, establishing the clinical relevance of NGS-detected variants for HH development remains a hard-working task, requiring further functional studies.