The conserved 3′ Angio‐domain defines a superfamily of short interspersed nuclear elements (SINEs) in higher plants

Repetitive sequences are ubiquitous components of eukaryotic genomes affecting genome size and evolution as well as gene regulation. Among them, short interspersed nuclear elements (SINEs) are non‐coding retrotransposons usually shorter than 1000 bp. They contain only few short conserved structural motifs, in particular an internal promoter derived from cellular RNAs and a mostly AT‐rich 3′ tail, whereas the remaining regions are highly variable. SINEs emerge and vanish during evolution, and often diversify into numerous families and subfamilies that are usually specific for only a limited number of species. In contrast, at the 3′ end of multiple plant SINEs we detected the highly conserved ‘Angio‐domain’. This 37 bp segment defines the Angio‐SINE superfamily, which encompasses 24 plant SINE families widely distributed across 13 orders within the plant kingdom. We retrieved 28 433 full‐length Angio‐SINE copies from genome assemblies of 46 plant species, frequently located in genes. Compensatory mutations in and adjacent to the Angio‐domain imply selective restraints maintaining its RNA structure. Angio‐SINE families share segmental sequence similarities, indicating a modular evolution with strong Angio‐domain preservation. We suggest that the conserved domain contributes to the evolutionary success of Angio‐SINEs through either structural interactions between SINE RNA and proteins increasing their transpositional efficiency, or by enhancing their accumulation in genes.     

Divergence of 3′ ends as a driver of short interspersed nuclear element (SINE) evolution in the Salicaceae

Short interspersed nuclear elements (SINEs) are small, non‐autonomous and heterogeneous retrotransposons that are widespread in plants. To explore the amplification dynamics and evolutionary history of SINE populations in representative deciduous tree species, we analyzed the genomes of the six following Salicaceae species: Populus deltoides, Populus euphratica, Populus tremula, Populus tremuloides, Populus trichocarpa, and Salix purpurea. We identified 11 Salicaceae SINE families (SaliS‐I to SaliS‐XI), comprising 27 077 full‐length copies. Most of these families harbor segmental similarities, providing evidence for SINE emergence by reshuffling or heterodimerization. We observed two SINE groups, differing in phylogenetic distribution pattern, similarity and 3′ end structure. These groups probably emerged during the ‘salicoid duplication’ (~65 million years ago) in the Salix–Populus progenitor and during the separation of the genus Salix (45–65 million years ago), respectively. In contrast to conserved 5′ start motifs across species and SINE families, the 3′ ends are highly variable in sequence and length. This extraordinary 3′‐end variability results from mutations in the poly(A) tail, which were fixed by subsequent amplificational bursts. We show that the dissemination of newly evolved 3′ ends is accomplished by a displacement of older motifs, leading to various 3′‐end subpopulations within the SaliS families.     

Role of 5′‐ and 3′‐untranslated regions of mRNAs in human diseases

Protein synthesis is often regulated at the level of initiation of translation, making it a critical step. This regulation occurs by both the cis‐regulatory elements, which are located in the 5′‐ and 3′‐UTRs (untranslated regions), and trans‐acting factors. A breakdown in this regulation machinery can perturb cellular metabolism, leading to various physiological abnormalities. The highly structured UTRs, along with features such as GC‐richness, upstream open reading frames and internal ribosome entry sites, significantly influence the rate of translation of mRNAs. In this review, we discuss how changes in the cis‐regulatory sequences of the UTRs, for example, point mutations and truncations, influence expression of specific genes at the level of translation. Such modifications may tilt the physiological balance from healthy to diseased states, resulting in conditions such as hereditary thrombocythaemia, breast cancer, fragile X syndrome, bipolar affective disorder and Alzheimer's disease. This information tends to establish the crucial role of UTRs, perhaps as much as that of coding sequences, in health and disease.     

Extremely slow rate of evolution in the HOX cluster revealed by comparison between Tanzanian and Indonesian coelacanths

Coelacanths are known as "living fossils" because their morphology has changed very little from that in the fossil record. To elucidate why coelacanths have evolved so slowly is thus of primary importance in evolutionary biology. In the present study, we determined the entire sequence of the HOX cluster of the Tanzanian coelacanth (Latimeria chalumnae) and compared it with that of the Indonesian coelacanth (L. menadoensis), which was available in the literature. The most intriguing result was the extremely small genetic divergence between the two coelacanths. The synonymous divergence of the HOX coding region between the two coelacanths was estimated to be 0.07%, which is ~11-fold smaller than that of human-chimp. When we applied the estimated divergence time of the two coelacanths of 6 million years ago (MYA) and 30 MYA, which were proposed in independent mitochondrial DNA analyses, the synonymous substitution rate of the coelacanth HOX cluster was estimated to be ~11-fold and 56-fold smaller than that of human-chimp, respectively. Thus, the present study implies that the reduction of the nucleotide substitution rate in coelacanth HOX genes may account for the conservation of coelacanth morphology during evolution.     

Characterization and evolution of 5' and 3' untranslated regions in eukaryotes

Untranslated regions (UTRs) in eukaryotes play a significant role in the regulation of translation and mRNA half-life, as well as interacting with specific RNA-binding proteins. However, UTRs receive less attention than more crucial elements such as genes, and the basic structural and evolutionary characteristics of UTRs of different species, and the relationship between these UTRs and the genome size and species gene number is not well understood. To address these questions, we performed a comparative analysis of 5' and 3' untranslated regions of different species by analyzing the basic characteristics of 244,976 UTRs from three eukaryote kingdoms (Plantae, Fungi, and Protista). The results showed that the UTR lengths and SSR frequencies in UTRs increased significantly with increasing species gene number while the length and G+C content in 5' UTRs and different types of repetitive sequences in 3' UTRs increased with the increase of genome size. We also found that the sequence length of 5' UTRs was significantly positively correlated with the presence of transposons and SSRs while the sequence length of 3' UTRs was significantly positively correlated with the presence of tandem repeat sequences. These results suggested that evolution of species complexity from lower organisms to higher organisms is accompanied by an increase in the regulatory complexity of UTRs, mediated by increasing UTR length, increasing G+C content of 5' UTRs, and insertion and expansion of repetitive sequences.     

Identification of single nucleotide polymorphisms in the bovine follicle‐stimulating hormone receptor and effects of genotypes on superovulatory response traits

In dairy cows, there is evidence that failure to respond to superovulation protocols is a heritable trait. In women, genotyping for the p.N680S single nucleotide polymorphism (SNP) in the follicle‐stimulating hormone receptor (FSHR) gene may help identify poor responders before ovarian stimulation is initiated. Our objectives were to identify SNPs in the coding region of the bovine FSHR gene and to investigate the effect of FSHR genotypes on superovulatory response in Holstein cattle. Sequencing of FSHR exons 1–10 revealed seven SNPs. Three were non‐synonymous mutations (c.337C>G, c.871A>G and c.1973C>G). SNP c.337C>G encodes for a proline‐to‐alanine (p.Pro113Ala) amino acid replacement in the extracellular ligand‐binding domain of the receptor. PCR‐RFLP analyses showed that homozygous GG Holstein cows present a higher percentage of viable embryos, whereas GG and CG animals have less unfertilised oocytes. SNP c.871A>G results in an isoleucine‐to‐valine (p.Ile291Val) modification, and homozygous AA animals present lower embryo yield after superovulatory treatments. SNP c.1973C>G corresponds to a threonine‐to‐serine (p.The658Ser) modification in the intracellular carboxyl‐terminal domain of the FSHR protein, and homozygous GG Holstein cows were associated with a lower embryo yield and a higher percentage of unfertilised oocytes. Our results suggest that specific alleles of the bovine FSHR gene are associated with variations in embryo yield and in the number of unfertilised oocytes.     

Two mutations in the 5′‐flanking region of the KITLG gene are associated with litter size of dairy goats

In this study, Xinong Saanen (SN) and Guanzhong (GZ) dairy goat breeds were used to detect single nucleotide polymorphisms (SNPs) in the 5′‐flanking region of the KITLG gene by DNA sequencing and primer‐introduced restriction analysis–polymerase chain reaction. Two novel SNPs (g.13090G>T and g.13664C>A) were identified (GenBank Accession no. KM658964). Furthermore, g.13090G>T and g.13664C>A loci were closely linked in SN and GZ breeds (r² > 0.33). Association analysis results showed that g.13090G>T and g.13664C>A SNPs significantly affected litter size (P < 0.05). The litter size of individuals with the combined genotype GG/CC from both dairy goat breeds was greater than that of individuals with TT/AA in average parity (P < 0.05). Known biochemical and physiological functions, along with our results, indicated that GG/CC could be used in marker‐assisted selection to choose individuals with greater litter size from both breeds. These results extend the spectrum of genetic variation in the caprine KITLG gene and may contribute to genetic resources and breeding of goats.     

Individual variation of human S1P1 coding sequence leads to heterogeneity in receptor function and drug interactions

Sphingosine 1-phosphate receptor 1 (S1P₁), an abundantly-expressed G protein-coupled receptor which regulates key vascular and immune responses, is a therapeutic target in autoimmune diseases. Fingolimod/Gilenya (FTY720), an oral medication for relapsing-remitting multiple sclerosis, targets S1P₁ receptors on immune and neural cells to suppress neuroinflammation. However, suppression of endothelial S1P₁ receptors is associated with cardiac and vascular adverse effects. Here we report the genetic variations of the S1P₁ coding region from exon sequencing of >12,000 individuals and their functional consequences. We conducted functional analyses of 14 nonsynonymous single nucleotide polymorphisms (SNPs) of the S1PR1 gene. One SNP mutant (Arg¹²⁰ to Pro) failed to transmit sphingosine 1-phosphate (S1P)-induced intracellular signals such as calcium increase and activation of p44/42 MAPK and Akt. Two other mutants (Ile⁴⁵ to Thr and Gly³⁰⁵ to Cys) showed normal intracellular signals but impaired S1P-induced endocytosis, which made the receptor resistant to FTY720-induced degradation. Another SNP mutant (Arg¹³ to Gly) demonstrated protection from coronary artery disease in a high cardiovascular risk population. Individuals with this mutation showed a significantly lower percentage of multi-vessel coronary obstruction in a risk factor-matched case-control study. This study suggests that individual genetic variations of S1P₁ can influence receptor function and, therefore, infer differential disease risks and interaction with S1P₁-targeted therapeutics.     

Common genetic variants of the β2‐adrenergic receptor affect its translational efficiency and are associated with human longevity

β‐adrenoceptors are the common pharmacological targets for the treatment of cardiovascular diseases and asthma. Genetic modifications of β‐adrenergic system in engineered mice affect their lifespan. Here, we tested whether genes encoding for key components of the β‐adrenergic signaling pathway are associated with human longevity. We performed a 10‐year follow‐up study of the Chinese longitudinal healthy longevity survey. The Han Chinese population in this study consisted of 963 long‐lived and 1028 geography‐matched young individuals. Sixteen SNPs from ADRB1, ADRB2, ADCY5, ADCY6, and MAPK1 were selected and genotyped. Two SNPs, rs1042718 (C/A) and rs1042719 (G/C), of ADRB2 in linkage disequilibrium (D' = 1.0; r² = 0.67) were found to be associated with enhanced longevity in men in two geographically isolated populations. Bonferroni‐corrected P‐values in a combined analysis were 0.00053–0.010. Men with haplotype A‐C showed an increased probability to become centenarians (the frequency of A‐C in long‐lived and young individuals are 0.332 and 0.250, respectively, OR = 1.49, CI 95% = 1.17–1.88, P = 0.0007), in contrast to those with haplotype C‐G (the frequency of C‐G in long‐lived and young individuals are 0.523 and 0.635, respectively, OR = 0.63, CI 95% = 0.51–0.78, P = 0.000018). The permuted P‐values were 0.00005 and 0.0009, respectively. ADRB2 encodes the β2‐adrenergic receptor; the haplotype A‐C markedly reduced its translational efficiency compared with C‐G (P = 0.002) in transfected HEK293 cells. Thus, our data indicate that enhanced production of β2‐adrenergic receptors caused by genetic variants is inversely associated with human lifespan.     

Alanine‐(87)‐threonine polymorphism impairs signaling and internalization of the human P2Y11 receptor,when co‐expressed with the P2Y1 receptor

The P2Y₁₁ nucleotide receptor detects high extracellular ATP concentrations. Mutations of the human P2RY₁₁ gene can play a role in brain autoimmune responses, and the P2Y₁₁ receptor alanine‐87‐threonine (A87T) polymorphism has been suggested to affect immune‐system functions. We investigated receptor functionality of the P2Y₁₁A87T mutant using HEK293 and 1321N1 astrocytoma cells. In HEK293 cells, the P2Y₁₁ receptor agonist 3′‐O‐(4‐benzoylbenzoyl)adenosine 5′‐triphosphate (BzATP) was completely inactive in evoking intracellular calcium release while the potency of ATP was reduced. ATP was also less potent in triggering cAMP generation. However, 1321N1 astrocytoma cells, which lack any endogenous P2Y₁ receptors, did not display a reduction. Only when 1321N1 cells were co‐transfected with P2Y₁₁A87T and P2Y₁ receptors, the calcium responses to the P2Y₁₁ receptor‐specific agonist BzATP were reduced. It is already known that P2Y₁ and P2Y₁₁ receptors interact. We thus conclude that the physiological impact of A87T mutation of the P2Y₁₁ receptor derives from detrimental effects on P2Y₁–P2Y₁₁ receptor interaction. We additionally investigated alanine‐87‐serine and alanine‐87‐tyrosine P2Y₁₁ receptor mutants. Both mutations rescue the response to BzATP in HEK293 cells, thus ruling out polarity of amino acid‐87 to be the molecular basis for altered receptor characteristics. We further found that the P2Y₁₁A87T receptor shows complete loss of nucleotide‐induced internalization in HEK293 cells. Thus, we demonstrate impaired signaling of the P2Y₁₁ A87T‐mutated receptors when co‐operating with P2Y₁ receptors.

     

A targeted genotyping approach enhances identification of variants in taste receptor and appetite/reward genes of potential functional importance for obesity‐related porcine traits

Taste receptors (TASRs) and appetite and reward (AR) mechanisms influence eating behaviour, which in turn affects food intake and risk of obesity. In a previous study, we used next generation sequencing to identify potentially functional mutations in TASR and AR genes and found indications for genetic associations between identified variants and growth and fat deposition in a subgroup of animals (n = 38) from the UNIK resource pig population. This population was created for studying obesity and obesity‐related diseases. In the present study we validated results from our previous study by investigating genetic associations between 24 selected single nucleotide variants in TASR and AR gene variants and 35 phenotypes describing obesity and metabolism in the entire UNIK population (n = 564). Fifteen variants showed significant association with specific obesity‐related phenotypes after Bonferroni correction. Six of the 15 genes, namely SIM1, FOS, TAS2R4, TAS2R9, MCHR2 and LEPR, showed good correlation between known biological function and associated phenotype. We verified a genetic association between potentially functional variants in TASR/AR genes and growth/obesity and conclude that the combination of identification of potentially functional variants by next generation sequencing followed by targeted genotyping and association studies is a powerful and cost‐effective approach for increasing the power of genetic association studies.     

Prevalence of the fibrinogen β‐chain,angiotensin‐converting enzyme and plasminogen activator inhibitor‐1 polymorphisms in Costa Rican young adults with thrombotic disease

Thrombotic disease is a multifactorial condition that involves both classical and genetic risk factors. We studied the association between the classical risk factors of hypertension and smoking, and polymorphisms on the genes of the angiotensin‐converting enzyme (ACE), the β‐chain of fibrinogen (FG), and the plasminogen activator inhibitor‐1 (PAI‐1) in patients with venous and arterial thrombosis. The present investigation is a retrospective case–control study. A total of 340 participants were analyzed, including 162 patients and 178 healthy controls. Hypertension and smoking showed a significant association with thrombotic disease (p < 0.05) but FG level was found significant risk factor only for the venous thrombosis (VT) group (p < 0.04). Significant differences between thrombotic groups were found for the studied polymorphisms of PAI‐1 (p < 0.0014), but for both FG β‐chain gene polymorphisms, none of the molecular analyses showed a positive sample for any mutating allele (p > 0.05). For the ACE polymorphism, the I allele present a protective effect in the general thrombotic group. This is one of the first reports in a Latin‐American population dealing with these molecular markers and thrombotic diseases. Copyright © 2010 John Wiley & Sons, Ltd.     

Erythroid 5‐aminolevulinate synthase mediates the upregulation of membrane band 3 protein expression by iron

Iron deficiency leads to abnormal expression and function of band 3 protein in erythrocytes, but the underlying mechanisms remain elusive. The mRNA of erythroid‐specific 5‐aminolevulinate synthase (eALAS) contains an iron response element and the eALAS protein is an important mediator of iron utilization by erythrocytes. In this study, we investigated the effect of short hairpin RNA (shRNA) mediated silencing of eALAS on the expression of band 3 protein induced by iron. By real‐time RT‐PCR and Western blot we showed that at mRNA and protein level iron‐induced expression of band 3 protein was lower in eALAS‐shRNA transfected K562 cells than in control cells. Of note, the lowest expression was detected in K562 cells cultured in iron deficiency condition (p < 0.01). Thus either iron deficiency or depletion of eALAS could suppress the expression of erythroid band 3 protein. These results demonstrated for the first time that iron and the iron‐regulatory system regulate the expression of the erythrocyte membrane proteins. Copyright © 2010 John Wiley & Sons, Ltd.     

Enantioselective σ1 receptor binding and biotransformation of the spirocyclic PET tracer 1′‐benzyl‐3‐(3‐fluoropropyl)‐3H‐spiro[[2]benzofuran‐1,4′‐piperidine]

It was shown that racemic (±)‐ 2 [1′‐benzyl‐3‐(3‐fluoropropyl)‐3H‐spiro[[2]benzofuran‐1,4′‐piperidine], WMS‐1813 ] represents a promising positron emission tomography (PET) tracer for the investigation of centrally located σ₁ receptors. To study the pharmacological activity of the enantiomers of 2 , a preparative HPLC separation of (R)‐2 and (S)‐2 was performed. The absolute configuration of the enantiomers was determined by CD‐spectroscopy together with theoretical calculations of the CD‐spectrum of a model compound. In receptor binding studies with the radioligand [³H]‐(+)‐pentazocine, (S)‐2 was thrice more potent than its (R)‐configured enantiomer (R)‐2 . The metabolic degradation of the more potent (S)‐enantiomer was considerably slower than the metabolism of (R)‐2 . The structures of the main metabolites of both enantiomers were elucidated by determination of the exact mass using an Orbitrap‐LC‐MS system. These experiments showed a stereoselective biotransformation of the enantiomers of 2 . Chirality, 2011. © 2010 Wiley‐Liss, Inc.