The Eyes Have It: Regulatory and Structural Changes Both Underlie Cichlid Visual Pigment Diversity

A major goal of evolutionary biology is to unravel the molecular genetic mechanisms that underlie functional diversification and adaptation. We investigated how changes in gene regulation and coding sequence contribute to sensory diversification in two replicate radiations of cichlid fishes. In the clear waters of Lake Malawi, differential opsin expression generates diverse visual systems, with sensitivities extending from the ultraviolet to the red regions of the spectrum. These sensitivities fall into three distinct clusters and are correlated with foraging habits. In the turbid waters of Lake Victoria, visual sensitivity is constrained to longer wavelengths, and opsin expression is correlated with ambient light. In addition to regulatory changes, we found that the opsins coding for the shortest-and longest-wavelength visual pigments have elevated numbers of potentially functional substitutions. Thus, we present a model of sensory evolution in which both molecular genetic mechanisms work in concert. Changes in gene expression generate large shifts in visual pigment sensitivity across the collective opsin spectral range, but changes in coding sequence appear to fine-tune visual pigment sensitivity at the short-and long-wavelength ends of this range, where differential opsin expression can no longer extend visual pigment sensitivity.

Equine cytochrome P450 2B6--genomic identification, expression and functional characterization with ketamine

Ketamine is an anesthetic and analgesic regularly used in veterinary patients. As ketamine is almost always administered in combination with other drugs, interactions between ketamine and other drugs bear the risk of either adverse effects or diminished efficacy. Since cytochrome P450 enzymes (CYPs) play a pivotal role in the phase I metabolism of the majority of all marketed drugs, drug-drug interactions often occur at the active site of these enzymes. CYPs have been thoroughly examined in humans and laboratory animals, but little is known about equine CYPs. The characterization of equine CYPs is essential for a better understanding of drug metabolism in horses. We report annotation, cloning and heterologous expression of the equine CYP2B6 in V79 Chinese hamster fibroblasts. After computational annotation of all CYP2B genes, the coding sequence (CDS) of equine CYP2B6 was amplified by RT-PCR from horse liver total RNA and revealed an amino acid sequence identity of 77% and a similarity of 93.7% to its human ortholog. A non-synonymous variant c.226G>A in exon 2 of the equine CYP2B6 was detected in 97 horses. The mutant A-allele showed an allele frequency of 82%. Two further variants in exon 3 were detected in one and two horses of this group, respectively. Transfected V79 cells were incubated with racemic ketamine and norketamine as probe substrates to determine metabolic activity. The recombinant equine CYP2B6 N-demethylated ketamine to norketamine and produced metabolites of norketamine, such as hydroxylated norketamines and 5,6-dehydronorketamine. V(max) for S-/and R-norketamine formation was 0.49 and 0.45nmol/h/mg cellular protein and K(m) was 3.41 and 2.66μM, respectively. The N-demethylation of S-/R-ketamine was inhibited concentration-dependently with clopidogrel showing an IC(50) of 5.63 and 6.26μM, respectively. The functional importance of the recorded genetic variants remains to be explored. Equine CYP2B6 was determined to be a CYP enzyme involved in ketamine and norketamine metabolism, thus confirming results from inhibition studies with horse liver microsomes. Clopidogrel seems to be a feasible inhibitor for equine CYP2B6. The specificity still needs to be established with other single equine CYPs. Heterologous expression of single equine CYP enzymes opens new possibilities to substantially improve the understanding of drug metabolism and drug interactions in horses.

A Gly98Val mutation in the N-Myc downstream regulated gene 1 (NDRG1) in Alaskan Malamutes with polyneuropathy.

The first cases of early-onset progressive polyneuropathy appeared in the Alaskan Malamute population in Norway in the late 1970s. Affected dogs were of both sexes and were ambulatory paraparetic, progressing to non-ambulatory tetraparesis. On neurologic examination, affected dogs displayed predominantly laryngeal paresis, decreased postural reactions, decreased spinal reflexes and muscle atrophy. The disease was considered eradicated through breeding programmes but recently new cases have occurred in the Nordic countries and the USA. The N-myc downstream-regulated gene (NDRG1) is implicated in neuropathies with comparable symptoms or clinical signs both in humans and in Greyhound dogs. This gene was therefore considered a candidate gene for the polyneuropathy in Alaskan Malamutes. The coding sequence of the NDRG1 gene derived from one healthy and one affected Alaskan Malamute revealed a non-synonymous G>T mutation in exon 4 in the affected dog that causes a Gly98Val amino acid substitution. This substitution was categorized to be "probably damaging" to the protein function by PolyPhen2 (score: 1.000). Subsequently, 102 Alaskan Malamutes from the Nordic countries and the USA known to be either affected (n = 22), obligate carriers (n = 7) or healthy (n = 73) were genotyped for the SNP using TaqMan. All affected dogs had the T/T genotype, the obligate carriers had the G/T genotype and the healthy dogs had the G/G genotype except for 13 who had the G/T genotype. A protein alignment showed that residue 98 is conserved in mammals and also that the entire NDRG1 protein is highly conserved (94.7%) in mammals. We conclude that the G>T substitution is most likely the mutation that causes polyneuropathy in Alaskan Malamutes. Our characterization of a novel candidate causative mutation for polyneuropathy offers a new canine model that can provide further insight into pathobiology and therapy of human polyneuropathy. Furthermore, selection against this mutation can now be used to eliminate the disease in Alaskan Malamutes.

Promoter- and RNA polymerase II-dependent hsp-16 gene association with nuclear pores in Caenorhabditis elegans.

Some inducible yeast genes relocate to nuclear pores upon activation, but the general relevance of this phenomenon has remained largely unexplored. Here we show that the bidirectional hsp-16.2/41 promoter interacts with the nuclear pore complex upon activation by heat shock in the nematode Caenorhabditis elegans. Direct pore association was confirmed by both super-resolution microscopy and chromatin immunoprecipitation. The hsp-16.2 promoter was sufficient to mediate perinuclear positioning under basal level conditions of expression, both in integrated transgenes carrying from 1 to 74 copies of the promoter and in a single-copy genomic insertion. Perinuclear localization of the uninduced gene depended on promoter elements essential for induction and required the heat-shock transcription factor HSF-1, RNA polymerase II, and ENY-2, a factor that binds both SAGA and the THO/TREX mRNA export complex. After induction, colocalization with nuclear pores increased significantly at the promoter and along the coding sequence, dependent on the same promoter-associated factors, including active RNA polymerase II, and correlated with nascent transcripts.

THE ROLE OF TYROSINE PHOSPHORYLATION IN THE FUNCTIONS OF THE TUMOR SUPPRESSOR GPRC5A

The retinoic acid inducible G protein coupled receptor family C group 5 type A (GPRC5A) is expressed preferentially in normal lung tissue but its expression is suppressed in the majority of human non-small cell lung cancer cell lines and tissues. This differential expression has led to the idea that GPRC5A is a potential tumor suppressor. This notion was supported by the finding that mice with a deletion of the Gprc5a gene develop spontaneous lung tumors. However, there are various tumor cell lines and tissue samples, including lung, that exhibit higher GPRC5A expression than normal tissues and some reports by other groups that GPRC5A transfection increased cell growth and colony formation. Obviously, GPRC5A has failed to suppress the development of the tumors and the growth of the cell lines where its expression is not suppressed. Since no mutations were detected in the coding sequence of GPRC5A in 20 NSCLC cell lines, it’s possible that GPRC5A acts as a tumor suppressor in the context of some cells but not in others. Alternatively, we raised the hypothesis that the GPRC5A protein may be inactivated by posttranslational modification(s) such as phosphorylation. It is well established that Serine/Threonine phosphorylation of G protein coupled receptors leads to their desensitization and in a few cases Tyrosine phosphorylation of GPCRs has been linked to internalization. Others reported that GPRC5A can undergo tyrosine phosphorylation in the cytoplasmic domain after treatment of normal human mammary epithelial cells (HMECs) with epidermal growth factor (EGF) or Heregulin. This suggested that GPRC5A is a substrate of EGFR. Therefore, we hypothesized that tyrosine phosphorylation of GPRC5A by activation of EGFR signaling may lead to its inactivation. To test this hypothesis, we transfected human embryo kidney (HEK) 293 cells with GPRC5A and EGFR expression vectors and confirmed that GPRC5A can be tyrosine phosphorylated after activation of EGFR by EGF. Further, we found that EGFR and GPRC5A can interact either directly or through other proteins and that inhibition of the EGFR kinase activity decreased the phosphorylation of GPRA5A and the interaction between GPRC5A and EGFR. In c-terminal of GPRC5A, There are four tyrosine residues Y317, Y320, Y347, Y350. We prepared GPRC5A mutants in which all four tyrosine residues had been replaced by phenylalanine (mutant 4F) or each individual Tyr residue was replaced by Phe and found that Y317 is the major site for EGFR mediated phosphorylation in the HEK293T cell line. We also found that EGF can induce GPRC5A internalization both in H1792 transient and stable cell lines. EGF also partially inactivates the suppressive function of GPRC5A on cell invasion activity and anchorage-independent growth ability of H1792 stable cell lines. These finding support our hypothesis that GPRC5A may be inactivated by posttranslational modification- tyrosine phosphorylation.

Analysis of mutations in $\beta$-tubulin genes affecting tubulin stability and assembly

Cmd4 is a colcemid-sensitive CHO cell line that is temperature sensitive for growth and expresses an altered $\beta$-tubulin, $\beta\sb1$. One revertant of this cell line, D2, exhibits a further alteration in $\beta\sb1$ resulting in an acidic shift in its isoelectric point and a decrease in its molecular weight to 40 kD, as measured by two dimensional gel electrophoresis. This $\beta$-tubulin variant has been shown to be assembly-defective and unstable. Characterization of the mutant $\beta\sb1$ in D2 by high pressure liquid chromatography (HPLC) revealed the loss of methionine containing tryptic peptides 7,8,9, and 10. Southern analysis of the genomic DNA digested with several different restriction enzymes resulted in the appearance of new restriction fragments 250 base pairs shorter than the corresponding fragments from the wild-type $\beta\sb1$-tubulin gene. Northern analysis on mRNA from D2 revealed two new message products that also differed by 250 bases from the corresponding wild type $\beta$-tubulin transcripts. To precisely define the region of the alteration, cloning and sequencing of the mutant and wild type genomic $\beta$-tubulin genes were conducted. A size-selected EcoRI genomic library was prepared using the Stratagene lambda Zap II phage cloning system. Using subclones of CHO $\beta$-tubulin cDNA as probes, a 2.5 kb wild type clone and a 2.3 kb mutant clone were identified from this library. Each of these was shown to contain a portion of the gene extending from intron 3 through the end of the coding sequence in exon 4 and into the 3$\sp\prime$ untranslated region on the basis of alignment with the published human $\beta$-tubulin sequence. Sequencing of the mutant 2.3 kb clone revealed that the mutation is due to a 246 base pair internal deletion in exon 4 (base pair 756-1001) that encodes amino acids 253-334. This deletion results in the loss of a putative binding site for GTP which could potentially explain the phenotype of this mutant $\beta$-tubulin. Also sequence comparison of the 3$\sp\prime$ untranslated region between different species revealed the conservation of 200 base pairs with 78% homology. It is proposed that this region could play an important role in the regulation of $\beta$-tubulin gene expression. ^

Mechanistic aspects of mRNA targeting for nonsense-mediated mRNA decay in human cells

The nonsense-mediated mRNA decay (NMD) pathway is best known as a translation-coupled quality control system that recognizes and degrades aberrant mRNAs with ORF-truncating premature termination codons (PTCs), but a more general role of NMD in posttranscriptional regulation of gene expression is indicated by transcriptome-wide mRNA profilings that identified a plethora of physiological mRNAs as NMD substrates. We try to decipher the mechanism of mRNA targeting to the NMD pathway in human cells. Recruitment of the conserved RNA-binding helicase UPF1 to target mRNAs has been reported to occur through interaction with release factors at terminating ribosomes, but evidence for translation-independent interaction of UPF1 with the 3’ untranslated region (UTR) of mRNAs has also been reported. We have transcriptome-wide determined the UPF1 binding sites by individual-nucleotide resolution UV crosslinking and immunoprecipitation (iCLIP) in human cells, untreated or after inhibiting translation. We detected a strongly enriched association of UPF1 with 3’ UTRs in undisturbed, translationally active cells. After translation inhibition, a significant increase in UPF1 binding to coding sequence (CDS) was observed, indicating that UPF1 binds RNA before translation and gets displaced from the CDS by translating ribosomes. This suggests that the decision to trigger NMD occurs after association of UPF1 with mRNA, presumably through activation of RNA-bound UPF1 by aberrant translation termination. In a second recent study, we re-visited the reported restriction of NMD in mammals to the ‘pioneer round of translation’, i.e. to cap-binding complex (CBC)-bound mRNAs. The limitation of mammalian NMD to early rounds of translation would indicate a – from an evolutionary perspective – unexpected mechanistic difference to NMD in yeast and plants, where PTC-containing mRNAs seem to be available to NMD at each round of translation. In contrast to previous reports, our comparison of decay kinetics of two NMD reporter genes in mRNA fractions bound to either CBC or the eukaryotic initiation factor 4E (eIF4E) in human cells revealed that NMD destabilizes eIF4E-bound transcripts as efficiently as those associated with CBC. These results corroborate an emerging unified model for NMD substrate recognition, according to which NMD can ensue at every aberrant translation termination event.

Deletion in the EVC2 gene causes chondrodysplastic dwarfism in Tyrolean Grey cattle.

During the summer of 2013 seven Italian Tyrolean Grey calves were born with abnormally short limbs. Detailed clinical and pathological examination revealed similarities to chondrodysplastic dwarfism. Pedigree analysis showed a common founder, assuming autosomal monogenic recessive transmission of the defective allele. A positional cloning approach combining genome wide association and homozygosity mapping identified a single 1.6 Mb genomic region on BTA 6 that was associated with the disease. Whole genome re-sequencing of an affected calf revealed a single candidate causal mutation in the Ellis van Creveld syndrome 2 (EVC2) gene. This gene is known to be associated with chondrodysplastic dwarfism in Japanese Brown cattle, and dwarfism, abnormal nails and teeth, and dysostosis in humans with Ellis-van Creveld syndrome. Sanger sequencing confirmed the presence of a 2 bp deletion in exon 19 (c.2993_2994ACdel) that led to a premature stop codon in the coding sequence of bovine EVC2, and was concordant with the recessive pattern of inheritance in affected and carrier animals. This loss of function mutation confirms the important role of EVC2 in bone development. Genetic testing can now be used to eliminate this form of chondrodysplastic dwarfism from Tyrolean Grey cattle.

The genomic substrate for adaptive radiation in African cichlid fish

Cichlid fishes are famous for large, diverse and replicated adaptive radiations in the Great Lakes of East Africa. To understand the molecular mechanisms underlying cichlid phenotypic diversity, we sequenced the genomes and transcriptomes of five lineages of African cichlids: the Nile tilapia (Oreochromis niloticus), an ancestral lineage with low diversity; and four members of the East African lineage: Neolamprologus brichardi/pulcher (older radiation, Lake Tanganyika), Metriaclima zebra (recent radiation, Lake Malawi), Pundamilia nyererei (very recent radiation, Lake Victoria), and Astatotilapia burtoni (riverine species around Lake Tanganyika). We found an excess of gene duplications in the East African lineage compared to tilapia and other teleosts, an abundance of non-coding element divergence, accelerated coding sequence evolution, expression divergence associated with transposable element insertions, and regulation by novel microRNAs. In addition, we analysed sequence data from sixty individuals representing six closely related species from Lake Victoria, and show genome-wide diversifying selection on coding and regulatory variants, some of which were recruited from ancient polymorphisms. We conclude that a number of molecular mechanisms shaped East African cichlid genomes, and that amassing of standing variation during periods of relaxed purifying selection may have been important in facilitating subsequent evolutionary diversification.

A mutation in hairless dogs implicates FOXI3 in ectodermal development

Mexican and Peruvian hairless dogs and Chinese crested dogs are characterized by missing hair and teeth, a phenotype termed canine ectodermal dysplasia (CED). CED is inherited as a monogenic autosomal semidominant trait. With genomewide association analysis we mapped the CED mutation to a 102-kilo-base pair interval on chromosome 17. The associated interval contains a previously uncharacterized member of the forkhead box transcription factor family (FOXI3), which is specifically expressed in developing hair and teeth. Mutation analysis revealed a frameshift mutation within the FOXI3 coding sequence in hairless dogs. Thus, we have identified FOXI3 as a regulator of ectodermal development.

The gene for histone RNA hairpin binding protein is located on human chromosome 4 and encodes a novel type of RNA binding protein.

The hairpin structure at the 3' end of animal histone mRNAs controls histone RNA 3' processing, nucleocytoplasmic transport, translation and stability of histone mRNA. Functionally overlapping, if not identical, proteins binding to the histone RNA hairpin have been identified in nuclear and polysomal extracts. Our own results indicated that these hairpin binding proteins (HBPs) bind their target RNA as monomers and that the resulting ribonucleoprotein complexes are extremely stable. These features prompted us to select for HBP-encoding human cDNAs by RNA-mediated three-hybrid selection in Saccharomyces cerevesiae. Whole cell extract from one selected clone contained a Gal4 fusion protein that interacted with histone hairpin RNA in a sequence- and structure-specific manner similar to a fraction enriched for bovine HBP, indicating that the cDNA encoded HBP. DNA sequence analysis revealed that the coding sequence did not contain any known RNA binding motifs. The HBP gene is composed of eight exons covering 19.5 kb on the short arm of chromosome 4. Translation of the HBP open reading frame in vitro produced a 43 kDa protein with RNA binding specificity identical to murine or bovine HBP. In addition, recombinant HBP expressed in S. cerevisiae was functional in histone pre-mRNA processing, confirming that we have indeed identified the human HBP gene.

Whole-Genome Sequencing of a Canine Family Trio Reveals a FAM83G Variant Associated with Hereditary Footpad Hyperkeratosis

Over 250 Mendelian traits and disorders, caused by rare alleles have been mapped in the canine genome. Although each disease is rare in the dog as a species, they are collectively common and have major impact on canine health. With SNP-based genotyping arrays, genome-wide association studies (GWAS) have proven to be a powerful method to map the genomic region of interest when 10-20 cases and 10-20 controls are available. However, to identify the genetic variant in associated regions, fine-mapping and targeted re-sequencing is required. Here we present a new approach using whole-genome sequencing (WGS) of a family trio without prior GWAS. As a proof-of-concept, we chose an autosomal recessive disease known as hereditary footpad hyperkeratosis (HFH) in Kromfohrl änder dogs. To our knowledge, this is the first time this family trio WGS-approach, has successfully been used to identify a genetic variant that perfectly segregates with a canine disorder. The sequencing of three Kromfohrl änder dogs from a family trio (an affected offspring and both its healthy parents) resulted in an average genome coverage of 9.2X per individual. After applying stringent filtering criteria for candidate causative coding variants, 527 single nucleotide variants (SNVs) and 15 indels were found to be homozygous in the affected offspring and heterozygous in the parents. Using the computer software packages ANNOVAR and SIFT to functionally annotate coding sequence differences and to predict their functional effect, resulted in seven candidate variants located in six different genes. Of these, only FAM83G:c155G>C (p.R52P) was found to be concordant in eight additional cases and 16 healthy Kromfohrl änder dogs.

A transposable element insertion in APOB causes cholesterol deficiency in Holstein cattle

Cholesterol deficiency, a new autosomal recessive inherited genetic defect in Holstein cattle, has been recently reported to have an influence on the rearing success of calves. The affected animals show unresponsive diarrhea accompanied by hypocholesterolemia and usually die within the first weeks or months of life. Here, we show that whole genome sequencing combined with the knowledge about the pedigree and inbreeding status of a livestock population facilitates the identification of the causative mutation. We resequenced the entire genomes of an affected calf and a healthy partially inbred male carrying one copy of the critical 2.24-Mb chromosome 11 segment in its ancestral state and one copy of the same segment with the cholesterol deficiency mutation. We detected a single structural variant, homozygous in the affected case and heterozygous in the non-affected carrier male. The genetic makeup of this key animal provides extremely strong support for the causality of this mutation. The mutation represents a 1.3kb insertion of a transposable LTR element (ERV2-1) in the coding sequence of the APOB gene, which leads to truncated transcripts and aberrant splicing. This finding was further supported by RNA sequencing of the liver transcriptome of an affected calf. The encoded apolipoprotein B is an essential apolipoprotein on chylomicrons and low-density lipoproteins, and therefore, the mutation represents a loss of function mutation similar to autosomal recessive inherited familial hypobetalipoproteinemia-1 (FHBL1) in humans. Our findings provide a direct gene test to improve selection against this deleterious mutation in Holstein cattle.

Expression of a bacterial serine acetyltransferase in transgenic potato plants leads to increased levels of cysteine and glutathione

The coding sequence of the wild-type, cys-sensitive, cysE gene from Escherichia coli, which encodes an enzyme of the cysteine biosynthetic pathway, namely serine acetyltransferase (SAT, EC 2.3.1.30), was introduced into the genome of potato plants under the control of the cauliflower mosaic virus 35S promoter. In order to target the protein into the chloroplast, cysE was translationally fused to the 5′-signal sequence of rbcS from Arabidopsis thaliana. Transgenic plants showed a high accumulation of the cysE mRNA. The chloroplastic localisation of the E. coli SAT protein was demonstrated by determination of enzymatic activities in enriched organelle fractions. Crude leaf extracts of these plants exhibited up to 20-fold higher SAT activity than those prepared from wild-type plants. The transgenic potato plants expressing the E. coli gene showed not only increased levels of enzyme activity but also exhibited elevated levels of cysteine and glutathione in leaves. Both were up to twofold higher than in control plants. However, the thiol content in tubers of transgenic lines was unaffected. The alterations observed in leaf tissue had no effect on the expression of O-acetylserine(thiol)-lyase, the enzyme which converts O-acetylserine, the product of SAT, to cysteine. Only a minor effect on its enzymatic activity was observed. In conclusion, the results presented here demonstrate the importance of SAT in plant cysteine biosynthesis and show that production of cysteine and related sulfur-containing compounds can be enhanced by metabolic engineering.

Cloning of a full-length cDNA encoding bovine interleukin 4 by the polymerase chain reaction

A human interleukin 4 (hIL-4)-encoding cDNA (hIL4) probe was used to screen a bovine genomic library, and three clones containing sequences with homology to the human and mouse IL4 cDNAs were isolated. Sequence information obtained from one of these genomic clones was used to design an oligodeoxyribonucleotide primer corresponding to the transcription start point region for use in the polymerase chain reaction (PCR). The PCR-RACE protocol, designed for the rapid amplification of cDNA ends, was successfully used to generate a full-length bovine IL4 (bIL4) cDNA clone from polyadenylated RNA isolated from concanavalin A-stimulated bovine lymph node cells. The bIL4 cDNA is 570 bp in length and contains an open reading frame of 405 nucleotides (nt), coding for a 15.1-kDa precursor of 135 amino acids (aa), which should be reduced to 12.6 kDa for unglycosylated bIL4 after cleavage of a putative hydrophobic leader sequence of 24 aa. The aa sequence contains one possible Asn-linked glycosylation site. Bovine IL4 is shorter than mouse (mIL4) and hIL4, because of a 51-nt deletion in the coding region. Comparison of the overall nt and deduced aa sequences shows a greater homology of bIL4 with hIL4 than with mIL4. This homology is not evenly distributed, however, with the nt sequences 5' and 3' of the coding region showing a much greater homology between all three species than the coding sequence.