Biased gene conversion and GC-content evolution in the coding sequences of reptiles and vertebrates.

Mammalian and avian genomes are characterized by a substantial spatial heterogeneity of GC-content, which is often interpreted as reflecting the effect of local GC-biased gene conversion (gBGC), a meiotic repair bias that favors G and C over A and T alleles in high-recombining genomic regions. Surprisingly, the first fully sequenced nonavian sauropsid (i.e., reptile), the green anole Anolis carolinensis, revealed a highly homogeneous genomic GC-content landscape, suggesting the possibility that gBGC might not be at work in this lineage. Here, we analyze GC-content evolution at third-codon positions (GC3) in 44 vertebrates species, including eight newly sequenced transcriptomes, with a specific focus on nonavian sauropsids. We report that reptiles, including the green anole, have a genome-wide distribution of GC3 similar to that of mammals and birds, and we infer a strong GC3-heterogeneity to be already present in the tetrapod ancestor. We further show that the dynamic of coding sequence GC-content is largely governed by karyotypic features in vertebrates, notably in the green anole, in agreement with the gBGC hypothesis. The discrepancy between third-codon positions and noncoding DNA regarding GC-content dynamics in the green anole could not be explained by the activity of transposable elements or selection on codon usage. This analysis highlights the unique value of third-codon positions as an insertion/deletion-free marker of nucleotide substitution biases that ultimately affect the evolution of proteins.

Evolution of vitellogenin genes: comparative analysis of the nucleotide sequences downstream of the transcription initiation site of four Xenopus laevis and one chicken gene.

Electron microscopic analysis of heteroduplexes between the most distantly related Xenopus vitellogenin genes (A genes X B genes) has revealed the distribution of homologous regions that have been preferentially conserved after the duplication events that gave rise to the multigene family in Xenopus laevis. DNA sequence analysis was limited to the region downstream of the transcription initiation site of the Xenopus genes A1, B1 and B2 and a comparison with the Xenopus A2 and the major chicken vitellogenin gene is presented. Within the coding regions of the first three exons, nucleotide substitutions resulting in amino acid changes accumulate at a rate similar to that observed in globin genes. This suggests that the duplication event which led to the formation of the A and B ancestral genes in Xenopus laevis occurred about 150 million years ago. Homologous exons of the A1-A2 and B1-B2 gene pairs, which formed about 30 million years ago, show a quite similar sequence divergence. In contrast, A1-A2 homologous introns seem to have evolved much faster than their B1-B2 counterparts.

Isolation of a fragment homologous to the rp49 constitutive gene of Drosophila in the Neotropical malaria vector Anopheles aquasalis (Diptera: Culicidae)

The constitutive ribosomal gene rp49 is frequently used as an endogenous control in Drosophila gene expression experiments. Using the degenerate primer PCR technique we have cloned a fragment homologous to this gene in Anopheles aquasalis Curry, a Neotropical vector of malaria. In addition, based on this first sequence, a new primer was designed, which allowed the isolation of fragments of rp49 in two other species, Aedes aegypti (Linnaeus) and Culex quinquefasciatus Say, suggesting that it could be used to clone fragments of this gene in a number of other mosquito species. Primers were also designed to specifically amplify rp49 cDNA fragments in An. aquasalis and Ae. aegypti, showing that rp49 could be used as a good constitutive control in gene expression studies of these and other vectorially important mosquito species.

Mapping of a gene coding for a major late structural polypeptide on the vaccinia virus genome.

Cell-free translation of total RNA isolated from vaccinia virus-infected cells late in infection results in a complex mixture of polypeptides. A monospecific antibody directed against one of the major structural proteins of the virus particle immunoprecipitated a single polypeptide with a molecular weight of 11,000 (11K) from this mixture. Immunoprecipitation was therefore used to identify the structural polypeptide among the in vitro translation products of RNA purified by hybridization selection to restriction fragments of the vaccinia virus genome. This allowed us to map the mRNA coding for the 11K polypeptide to the extreme left-hand end of the HindIII E fragment. Detailed transcriptional mapping of this region of the genome by nuclease S1 analysis revealed the presence of a late RNA transcribed from the rightward-reading strand. Its 5' end mapped at ca. 130 base pairs to the left of the HindIII site at the junction between the HindIII F and E fragments. The map position of this RNA coincided precisely with the map position of the late message coding for the 11K polypeptide.

SGP-1: prediction and validation of homologous genes based on sequence alignments

Conventional methods of gene prediction rely on the recognition of DNA-sequence signals, the coding potential or the comparison of a genomic sequence with a cDNA, EST, or protein database. Reasons for limited accuracy in many circumstances are species-specific training and the incompleteness of reference databases. Lately, comparative genome analysis has attracted increasing attention. Several analysis tools that are based on human/mouse comparisons are already available. Here, we present a program for the prediction of protein-coding genes, termed SGP-1 (Syntenic Gene Prediction), which is based on the similarity of homologous genomic sequences. In contrast to most existing tools, the accuracy of SGP-1 depends little on species-specific properties such as codon usage or the nucleotide distribution. SGP-1 may therefore be applied to nonstandard model organisms in vertebrates as well as in plants, without the need for extensive parameter training. In addition to predicting genes in large-scale genomic sequences, the program may be useful to validate gene structure annotations from databases. To this end, SGP-1 output also contains comparisons between predicted and annotated gene structures in HTML format. The program can be accessed via a Web server at http://soft.ice.mpg.de/sgp-1. The source code, written in ANSI C, is available on request from the authors.

One hundred base pairs of 5' flanking sequence of a vaccinia virus late gene are sufficient to temporally regulate late transcription.

A vaccinia virus late gene coding for a major structural polypeptide of 11 kDa was sequenced. Although the 5' flanking gene region is very A+T rich, it shows little homology either to the corresponding region of vaccinia early genes or to consensus sequences characteristic of most eukaryotic genes. Three DNA fragments (100, 200, and 500 base pairs, respectively), derived from the flanking region and including the late gene mRNA start site, were inserted into the coding sequence of the vaccinia virus thymidine kinase (TK) early gene by homologous in vivo recombination. Recombinants were selected on the basis of their TK- phenotype. Cells were infected with the recombinant viruses and RNA was isolated at 1-hr intervals. Transcripts initiating either from the TK early promoter, or from the late gene promoter at its authentic position, or from the translocated late gene promoters within the early gene were detected by nuclease S1 mapping. Early after infection, only transcripts from the TK early promoter were detected. Later in infection, however, transcripts were also initiated from the translocated late promoters. This RNA appeared at the same time and in similar quantities as the RNA from the late promoter at its authentic position. No quantitative differences in promoter efficiency between the 100-, 200-, and 500-base-pair insertions were observed. We conclude that all necessary signals for correct regulation of late-gene expression reside within only 100 base pairs of 5' flanking sequence.

Vitellogenesis and the vitellogenin gene family.

Vitellogenin is synthesized under estrogen control in the liver, extensively modified, transported to the ovary, and there processed to the yolk proteins lipovitellin and phosvitin. In the frog Xenopus laevis there are at least four distinct but related vitellogenin genes. The two genes A1 and A2 have a 95 percent sequence homology in their messenger RNA coding regions, and contain 33 introns that interrupt the coding region (exons) at homologous positions. Sequences and lengths of analogous introns differ, and many introns contain repetitive DNA elements. The introns in these two genes that have apparently arisen by duplication have diverged extensively by events that include deletions, insertions, and probably duplications. Rapid evolutionary change involving rearrangements and the presence of repeated DNA suggests that the bulk of the sequences within introns may not have any specific function.

High-level transgene expression by homologous recombination-mediated gene transfer.

Gene transfer and expression in eukaryotes is often limited by a number of stably maintained gene copies and by epigenetic silencing effects. Silencing may be limited by the use of epigenetic regulatory sequences such as matrix attachment regions (MAR). Here, we show that successive transfections of MAR-containing vectors allow a synergistic increase of transgene expression. This finding is partly explained by an increased entry into the cell nuclei and genomic integration of the DNA, an effect that requires both the MAR element and iterative transfections. Fluorescence in situ hybridization analysis often showed single integration events, indicating that DNAs introduced in successive transfections could recombine. High expression was also linked to the cell division cycle, so that nuclear transport of the DNA occurs when homologous recombination is most active. Use of cells deficient in either non-homologous end-joining or homologous recombination suggested that efficient integration and expression may require homologous recombination-based genomic integration of MAR-containing plasmids and the lack of epigenetic silencing events associated with tandem gene copies. We conclude that MAR elements may promote homologous recombination, and that cells and vectors can be engineered to take advantage of this property to mediate highly efficient gene transfer and expression.

Sequence homologies within the 5' end region of the estrogen-controlled vitellogenin gene in Xenopus and chicken.

In oviparous vertebrates vitellogenin, the precursor of the major yolk proteins, is synthesized in the liver of mature females under the control of estrogen. We have established the organization and primary structure of the 5' end region of the Xenopus laevis vitellogenin A2 gene and of the major chicken vitellogenin gene. The first three homologous exons have exactly the same length in both species, namely 53, 21 and 152 nucleotides, and present an overall sequence homology of 60%. In both species, the 5'-non-coding region of the vitellogenin mRNA measures only 13 nucleotides, nine of which are conserved. In contrast, the corresponding introns of the Xenopus and the chicken vitellogenin gene show no significant sequence homology. Within the 500 nucleotides preceding the 5' end of the genes, at least six blocks with sequence homologies of greater than 70% were detected. It remains to be demonstrated which of these conserved sequences, if any, are involved in the hormone-regulated expression of the vitellogenin genes.

Evolution of the P-type II ATPase gene family in the fungi and presence of structural genomic changes among isolates of Glomus intraradices.

BACKGROUND: The P-type II ATPase gene family encodes proteins with an important role in adaptation of the cell to variation in external K+, Ca2+ and Na2+ concentrations. The presence of P-type II gene subfamilies that are specific for certain kingdoms has been reported but was sometimes contradicted by discovery of previously unknown homologous sequences in newly sequenced genomes. Members of this gene family have been sampled in all of the fungal phyla except the arbuscular mycorrhizal fungi (AMF; phylum Glomeromycota), which are known to play a key-role in terrestrial ecosystems and to be genetically highly variable within populations. Here we used highly degenerate primers on AMF genomic DNA to increase the sampling of fungal P-Type II ATPases and to test previous predictions about their evolution. In parallel, homologous sequences of the P-type II ATPases have been used to determine the nature and amount of polymorphism that is present at these loci among isolates of Glomus intraradices harvested from the same field. RESULTS: In this study, four P-type II ATPase sub-families have been isolated from three AMF species. We show that, contrary to previous predictions, P-type IIC ATPases are present in all basal fungal taxa. Additionally, P-Type IIE ATPases should no longer be considered as exclusive to the Ascomycota and the Basidiomycota, since we also demonstrate their presence in the Zygomycota. Finally, a comparison of homologous sequences encoding P-type IID ATPases showed unexpectedly that indel mutations among coding regions, as well as specific gene duplications occur among AMF individuals within the same field. CONCLUSION: On the basis of these results we suggest that the diversification of P-Type IIC and E ATPases followed the diversification of the extant fungal phyla with independent events of gene gains and losses. Consistent with recent findings on the human genome, but at a much smaller geographic scale, we provided evidence that structural genomic changes, such as exonic indel mutations and gene duplications are less rare than previously thought and that these also occur within fungal populations.

Color differences among feral pigeons (Columba livia) are not attributable to sequence variation in the coding region of the melanocortin-1 receptor gene (MC1R)

Genetic variation at the melanocortin-1 receptor (MC1R) gene is correlated with melanin color variation in many birds. Feral pigeons (Columba livia) show two major melanin-based colorations: a red coloration due to pheomelanic pigment and a black coloration due to eumelanic pigment. Furthermore, within each color type, feral pigeons display continuous variation in the amount of melanin pigment present in the feathers, with individuals varying from pure white to a full dark melanic color. Coloration is highly heritable and it has been suggested that it is under natural or sexual selection, or both. Our objective was to investigate whether MC1R allelic variants are associated with plumage color in feral pigeons.We sequenced 888 bp of the coding sequence of MC1R among pigeons varying both in the type, eumelanin or pheomelanin, and the amount of melanin in their feathers. We detected 10 non-synonymous substitutions and 2 synonymous substitution but none of them were associated with a plumage type. It remains possible that non-synonymous substitutions that influence coloration are present in the short MC1R fragment that we did not sequence but this seems unlikely because we analyzed the entire functionally important region of the gene.Our results show that color differences among feral pigeons are probably not attributable to amino acid variation at the MC1R locus. Therefore, variation in regulatory regions of MC1R or variation in other genes may be responsible for the color polymorphism of feral pigeons.

Control of gene expression in melanocytes and RPE by conserved distal regulatory elements

Résumé Durant le développement embryonnaire, les cellules pigmentaires des mammifères se développent à partir de deux origines différentes : les melanocytes se développent à partir de la crête neurale alors que les cellules de la rétine pigmentaire (RP) ont une origine neuronale. Un grand nombre de gènes sont impliqués dans la pigmentation dont les gènes de la famille tyrosinase à savoir Tyr, Tyrp1 et Dct. Certaines études ont suggéré que les gènes de la pigmentation sont régulés de manière différentielle dans les mélanocytes et dans la RP. Dans ce travail, les gènes de la famille tyrosinase ont été étudiés comme modèle de la régulation des gènes de la pigmentation par des éléments régulateurs agissant à distance. II a été montré que le promoteur du gène Tyrp1pouvait induire l'expression d'un transgène uniquement dans la RP alors que ce gène est aussi exprimé dans les mélanocytes comme le montre le phénotype des souris mutantes pour Tyrp1. Ce résultat suggère que les éléments régulateurs du promoteur sont suffisants pour l'expression dans la RP mais pas pour l'expression dans les mélanocytes. J'ai donc cherché à identifier la séquence qui régule l'expression dans les mélanocytes. Un chromosome artificiel bactérien (CAB) contenant le gène Tyrp1 s'est avéré suffisant pour induire l'expression dans les mélanocytes, comme démontré par la correction du phénotype mutant. La séquence de ce CAB contient plusieurs régions très conservées qui pourraient représenter de nouveaux éléments régulateurs. Par la suite, j'ai focalisé mon analyse sur une séquence située à -I5 kb qui s'est révélée être un amplificateur spécifique aux mélanocytes comme démontré par des expériences de cultures cellulaire et de transgenèse. De plus, une analyse poussée de cet élément a révélé que le facteur de transcription Sox 10 représentait un transactivateur de cet amplificateur. Comme pour Tyrp1, la régulation du gène tyrosinase est contrôlée par différents éléments régulateurs dans les mélanocytes et la RP. Il a été montré que le promoteur de tyrosinase n'était pas suffisant pour une forte expression dans les mélanocytes et la RP. De plus, l'analyse de la région située en amont a révélé la présence d'un amplificateur nécessaire à l'expression dans les mélanocytes à la position -15 kb. Cet amplificateur n'est toutefois pas actif dans la RP mais agit comme un répresseur dans ces cellules. Ces résultats indiquent que certains éléments nécessaires à l'expression dans les deux types de cellules pigmentaires sont absents de ces constructions. Comme pour Tyrp1, j'ai en premier lieu démontré qu'un CAB était capable de corriger le phénotype albinique, puis ai inséré un gène reporter (lacZ) dans le CAB par recombinaison homologue et ai finalement analysé l'expression du reporter en transgenèse. Ces souris ont montré une expression forte du lacZ dans les mélanocytes et la RP, ce qui indique que le CAB contient les séquences régulatrices nécessaires à l'expression correcte de tyrosinase. Afin de localiser plus précisément les éléments régulateurs, j'ai ensuite généré des délétions dans le CAB et analysé l'expression du lacZ en transgenèse. La comparaison de séquences génomiques provenant de différentes espèces a permis par la suite d'identifier des régions représentant de nouveaux éléments régulateurs potentiels. En utilisant cette approche, j'ai identifié une région qui se comporte comme un amplificateur dans la RP et qui est nécessaire à l'expression de tyrosinase dans ce tissu. De plus, j'ai identifié les facteurs de transcription Mitf et Sox10 comme transactivateurs de l'amplificateur spécifique aux mélanocytes situé à -15 kb. L'identification et la caractérisation des ces éléments régulateurs des gènes tyrosinase et Tyrp1confirme donc que la régulation différentielle des gènes dans les mélanocytes et la RP est liée à des éléments régulateurs séparés. Summary Pigment cells of mammals originate from two different lineages: melanocytes arise from the neural crest, whereas cells of the retinal pigment epithelium (RPE) originate from the optic cup of the developing forebrain. A large set of genes are involved in pigmentation, including the members of the tyrosinase gene family, namely tyrosinase, Tyrp1 and Dct. Previous studies have suggested that pigmentation genes are differentially regulated in melanocytes and RPE. In this work, the tyrosinase gene family was used as a model for studying the involvement of distal regulatory elements in pigment cell-specific gene expression. The promoter of the Tyrp1 gene has been shown to drive detectable transgene expression only to the RPE, even though the gene is also expressed in melanocytes as evident from Tyrp1-mutant mice. This indicates that the regulatory elements responsible for Tyrp1 gene expression in the RPE are not sufficient for expression in melanocytes. I thus searched for a putative melanocyte-specific regulatory sequence and demonstrate that a bacterial artificial chromosome (BAC) containing the Tyrp1 gene and surrounding sequences is able to target transgenic expression to melanocytes and to rescue the Tyrp1 b (brown) phenotype. This BAC contains several highly conserved non-coding sequences that might represent novel regulatory elements. I further focused on a sequence located at -15 kb which I identified as amelanocyte-specific enhancer as shown by cell culture and transgenic mice. In addition, further functional analysis identified the transcription factor Sox10 as being able to bind and transactivate this enhancer. As for Tyrp1, tyrosinase gene regulation is mediated by different cis-regulatory elements in melanocytes and RPE. It was shown that the tyrosinase promoter was not sufficient to confer strong and specific expression in melanocytes and RPE. Moreover, analysis of tyrosinase upstream sequence, revealed the presence of a specific enhancer at position -15 kb which was necessary to confer strong expression in melanocytes. This enhancer element however failed to act as an enhancer in the RPE, but rather repressed expression. This indicates that some regulatory elements required for tyrosinase expression in both RPE and melanocytes are still missing from these constructs. As for Tyrp1, I first demonstrated that a BAC containing the Tyr gene is able to rescue the Tyr c (albino) phenotype in mice, then I inserted a lacZ reporter gene in the BAC by homologous recombination, and finally analysed the pattern of lacZ expression in transgenic mice. These mice showed strong lacZ expression in both RPE and melanocytes, indicating that the BAC contains the regulatory sequences required for proper tyrosinase expression. In order to localize more precisely these regulatory elements, I have then generated several deletions in the BAC and analysed lacZ expression in transgenic mice. Multi-species comparative genomic analysis then allowed identifying conserved sequences that potentially represent novel regulatory elements. Using this experimental approach, I identified a region that behaves as a RPE-specific enhancer and that is required for tyrosinase expression in the retina] pigment epithelium. In addition, I identified the transcription factors Mitf and Sox l0 as being transactivators of the melanocyte-specific enhancer located at -l5 kb. The identification and characterization of these tyrosinase and Tyrp1 distal regulatory element supports the idea that separate regulatory sequences mediate differential gene expression in melanocytes and RPE.

Detection and coat protein gene characterization of an isolate of Grapevine virus B from corky bark-affected grapevines in Southern Brazil

An isolate of Grapevine virus B (GVB), obtained by indexing Vitis labrusca and V. vinifera grapevines on the indicator LN33, was transmitted mechanically to several Nicotiana species. The virus was partially purified from N. cavicola and the coat protein estimated at 23 kDa by SDS-PAGE. In negatively stained leaf extracts of experimentally inoculated N. cavicola and N. occidentalis, flexuous particles with cross banding were observed, predominantly measuring 750-770 x 12 nm, with a modal length of 760 nm. Decoration indicated a clear, positive reaction against AS-GVB. In DAS-ELISA, GVB was detected in N. cavicola and grapevine extracts, and Western blots showed homologous and cross reaction of GVB and GVA antisera with GVB coat protein. Using specific primers for GVB, a fragment of 594 bp, comprising the coat protein gene coding for 197 amino acids, was amplified by RT-PCR with viral RNA extracted from GVB-infected N. occidentalis. The nucleotide and the deduced amino acid sequences of the coat protein gene showed high identities with Italian and Japanese isolates of GVB.

A novel polymorphism in the coding region of the vasopressin type 2 receptor gene

Nephrogenic diabetes insipidus (NDI) is a rare disease characterized by renal inability to respond properly to arginine vasopressin due to mutations in the vasopressin type 2 receptor (V2(R)) gene in affected kindreds. In most kindreds thus far reported, the mode of inheritance follows an X chromosome-linked recessive pattern although autosomal-dominant and autosomal-recessive modes of inheritance have also been described. Studies demonstrating mutations in the V2(R) gene in affected kindreds that modify the receptor structure, resulting in a dys- or nonfunctional receptor have been described, but phenotypically indistinguishable NDI patients with a structurally normal V2(R) gene have also been reported. In the present study, we analyzed exon 3 of the V2(R) gene in 20 unrelated individuals by direct sequencing. A C®T alteration in the third position of codon 331 (AGC®AGT), which did not alter the encoded amino acid, was found in nine individuals, including two unrelated patients with NDI. Taken together, these observations emphasize the molecular heterogeneity of a phenotypically homogeneous syndrome

A search for RNA insertions and NS3 gene duplication in the genome of cytopathic isolates of bovine viral diarrhea virus

Calves born persistently infected with non-cytopathic bovine viral diarrhea virus (ncpBVDV) frequently develop a fatal gastroenteric illness called mucosal disease. Both the original virus (ncpBVDV) and an antigenically identical but cytopathic virus (cpBVDV) can be isolated from animals affected by mucosal disease. Cytopathic BVDVs originate from their ncp counterparts by diverse genetic mechanisms, all leading to the expression of the non-structural polypeptide NS3 as a discrete protein. In contrast, ncpBVDVs express only the large precursor polypeptide, NS2-3, which contains the NS3 sequence within its carboxy-terminal half. We report here the investigation of the mechanism leading to NS3 expression in 41 cpBVDV isolates. An RT-PCR strategy was employed to detect RNA insertions within the NS2-3 gene and/or duplication of the NS3 gene, two common mechanisms of NS3 expression. RT-PCR amplification revealed insertions in the NS2-3 gene of three cp isolates, with the inserts being similar in size to that present in the cpBVDV NADL strain. Sequencing of one such insert revealed a 296-nucleotide sequence with a central core of 270 nucleotides coding for an amino acid sequence highly homologous (98%) to the NADL insert, a sequence corresponding to part of the cellular J-Domain gene. One cpBVDV isolate contained a duplication of the NS3 gene downstream from the original locus. In contrast, no detectable NS2-3 insertions or NS3 gene duplications were observed in the genome of 37 cp isolates. These results demonstrate that processing of NS2-3 without bulk mRNA insertions or NS3 gene duplications seems to be a frequent mechanism leading to NS3 expression and BVDV cytopathology.