Comparative population genomics in animals uncovers the determinants of genetic diversity.

Genetic diversity is the amount of variation observed between DNA sequences from distinct individuals of a given species. This pivotal concept of population genetics has implications for species health, domestication, management and conservation. Levels of genetic diversity seem to vary greatly in natural populations and species, but the determinants of this variation, and particularly the relative influences of species biology and ecology versus population history, are still largely mysterious. Here we show that the diversity of a species is predictable, and is determined in the first place by its ecological strategy. We investigated the genome-wide diversity of 76 non-model animal species by sequencing the transcriptome of two to ten individuals in each species. The distribution of genetic diversity between species revealed no detectable influence of geographic range or invasive status but was accurately predicted by key species traits related to parental investment: long-lived or low-fecundity species with brooding ability were genetically less diverse than short-lived or highly fecund ones. Our analysis demonstrates the influence of long-term life-history strategies on species response to short-term environmental perturbations, a result with immediate implications for conservation policies.

The global impact of Wolbachia on mitochondrial diversity and evolution.

The spread of maternally inherited microorganisms, such as Wolbachia bacteria, can induce indirect selective sweeps on host mitochondria, to which they are linked within the cytoplasm. The resulting reduction in effective population size might lead to smaller mitochondrial diversity and reduced efficiency of natural selection. While documented in several host species, it is currently unclear if such a scenario is common enough to globally impact the diversity and evolution of mitochondria in Wolbachia-infected lineages. Here, we address this question using a mapping of Wolbachia acquisition/extinction events on a large mitochondrial DNA tree, including over 1000 species. Our analyses indicate that on a large phylogenetic scale, other sources of variation, such as mutation rates, tend to hide the effects of Wolbachia. However, paired comparisons between closely related infected and uninfected taxa reveal that Wolbachia is associated with a twofold reduction in silent mitochondrial polymorphism, and a 13% increase in nonsynonymous substitution rates. These findings validate the conjecture that the widespread distribution of Wolbachia infections throughout arthropods impacts the effective population size of mitochondria. These effects might in part explain the disconnection between genetic diversity and demographic population size in mitochondria, and also fuel red-queen-like cytonuclear co-evolution through the fixation of deleterious mitochondrial alleles.

Regulation of dauer larva development in Caenorhabditis elegans by daf-18, a homologue of the tumour suppressor PTEN.

The tumour suppressor gene PTEN (also called MMAC1 or TEP1) is somatically mutated in a variety of cancer types [1] [2] [3] [4]. In addition, germline mutation of PTEN is responsible for two dominantly inherited, related cancer syndromes called Cowden disease and Bannayan-Ruvalcaba-Riley syndrome [4]. PTEN encodes a dual-specificity phosphatase that inhibits cell spreading and migration partly by inhibiting integrin-mediated signalling [5] [6] [7]. Furthermore, PTEN regulates the levels of phosphatidylinositol 3,4,5-trisphosphate (PIP3) by specifically dephosphorylating position 3 on the inositol ring [8]. We report here that the dauer formation gene daf-18 is the Caenorhabditis elegans homologue of PTEN. DAF-18 is a component of the insulin-like signalling pathway controlling entry into diapause and adult longevity that is regulated by the DAF-2 receptor tyrosine kinase and the AGE-1 PI 3-kinase [9]. Others have shown that mutation of daf-18 suppresses the life extension and constitutive dauer formation associated with daf-2 or age-1 mutants. Similarly, we show that inactivation of daf-18 by RNA-mediated interference mimics this suppression, and that a wild-type daf-18 transgene rescues the dauer defect. These results indicate that PTEN/daf-18 antagonizes the DAF-2-AGE-1 pathway, perhaps by catalyzing dephosphorylation of the PIP3 generated by AGE-1. These data further support the notion that mutations of PTEN contribute to the development of human neoplasia through an aberrant activation of the PI 3-kinase signalling cascade.

Highly conserved RNA sequences that are sensors of environmental stress.

The putative function of highly conserved regions (HCRs) within 3' untranslated regions (3'UTRs) as regulatory RNA sequences was efficiently and quantitatively assessed by using modular retroviral vectors. This strategy led to the identification of HCRs that alter gene expression in response to oxidative or mitogenic stress. Databases were screened for UTR sequences of >100 nucleotides that had retained 70% identity over more than 300 million years of evolution. The effects of 10 such HCRs on a standard reporter mRNA or protein were studied. To this end, we developed a modular retroviral vector that can allow for a direct comparison of the effects of different HCRs on gene expression independent of their gene-intrinsic 5'UTR, promoter, protein coding region, or poly(A) sequence. Five of the HCRs tested decreased mRNA steady-state levels 2- to 10-fold relative to controls, presumably by altering mRNA stability. One HCR increased translation, and one decreased translation. Elevated mitogen levels caused four HCRs to increase protein levels twofold. One HCR increased protein levels fourfold in response to hypoxia. Although nonconserved UTR sequences may also have a role, these results provide evidence that sequences that are highly conserved during evolution are good candidates for RNA motifs with posttranscriptional regulatory functions in gene expression.

Searching for regulatory elements in human noncoding sequences.

Important progress has been made in the past two years in the identification of Pol II promoters. For most other regulatory elements, however, current biological knowledge is still insufficient to allow the development of prediction tools. The phylogenetic-footprinting strategy, which is based on the comparative analysis of homologous sequences, is a very efficient approach to identify new unknown regulatory elements. The recent organization of large-scale sequencing projects for some model vertebrate organisms will be extremely valuable for the prediction of regulatory elements in the human genome.

Transposons but not retrotransposons are located preferentially in regions of high recombination rate in Caenorhabditis elegans.

We analyzed the distribution of transposable elements (TEs: transposons, LTR retrotransposons, and non-LTR retrotransposons) in the chromosomes of the nematode Caenorhabditis elegans. The density of transposons (DNA-based elements) along the chromosomes was found to be positively correlated with recombination rate, but this relationship was not observed for LTR or non-LTR retrotransposons (RNA-based elements). Gene (coding region) density is higher in regions of low recombination rate. However, the lower TE density in these regions is not due to the counterselection of TE insertions within exons since the same positive correlation between TE density and recombination rate was found in noncoding regions (both in introns and intergenic DNA). These data are not compatible with a global model of selection acting against TE insertions, for which an accumulation of elements in regions of reduced recombination is expected. We also found no evidence for a stronger selection against TE insertions on the X chromosome compared to the autosomes. The difference in distribution of the DNA and RNA-based elements along the chromosomes in relation to recombination rate can be explained by differences in the transposition processes.

Cloning of the mouse BTG3 gene and definition of a new gene family (the BTG family) involved in the negative control of the cell cycle.

It is well known that loss of tumor suppressor genes and more generally of antiproliferative genes plays a key role in the development of most tumors. We report here the cloning of the mouse BTG3 gene and show that its human counterpart maps on chromosome 21. This evolutionarily conserved gene codes for a 30 kDa protein and is expressed in most adult murine and human tissues analyzed. However, we demonstrate that its expression is cell cycle dependent and peaks at the end of the G1 phase. This gene is homologous to the human BTG1, BTG2 and TOB genes which were demonstrated to act as inhibitors of cell proliferation. Its description allowed us to define better this seven gene family (the BTG gene family) at the structural level and to speculate about its physiological role in normal and tumoral cells. This family is mainly characterized by the presence of two conserved domains (BTG boxes A and B) of as yet undetermined function which are separated by a non-conserved 20-25 amino acid sequence.

Isolation and characterization of a cDNA encoding a chicken actin-like protein.

We report the isolation and characterization of a chicken cDNA which putatively encodes an actin-like protein (chACTL). This 394-amino-acid (aa) polypeptide shares sequence homology (81, 70 and 67% identical aa, respectively) with three actin-related proteins (ARP) described for Drosophila melanogaster (ARP14D), Caenorhabditis elegans (ACTL) and Saccharomyces cerevisiae (ACT2). At least six chACTL transcripts were detected in different tissues during chick embryogenesis. Sequence analysis suggests that at least three groups of ARP have been evolutionarily conserved.

Human and nematode orthologs--lessons from the analysis of 1800 human genes and the proteome of Caenorhabditis elegans.

Recently, we have defined and analyzed over 1800 orthologous human and rodent genes. Here we extend this work to compare human and Caenorhabditis elegans coding sequences. 1880 human proteins were compared with about 20000 predicted nematode proteins presumably comprising nearly the complete proteome of C. elegans. We found that 44% of human/rodent orthologs have convincing nematode counterparts. On average, the amino acid similarity and identity between aligned human and C. elegans orthologous gene products are 69.3% and 49.1% respectively, and the nucleotide identity is 49.8%. Detailed investigation of our results suggests that some nematode gene predictions are incorrect, leading to erroneous pairing with human genes (e.g. calcineurin and polymerase II elongation factor III). Furthermore, other proteins (i.e. homologs of human ribosomal proteins S20 and L41, thymosin) are missing entirely from the nematode proteome, suggesting that it may not be complete. These results underscore the fact that metazoan gene prediction is a very challenging task and that most computer-predicted nematode genes require supporting evidence of their existence from comparative genomics and/or laboratory investigation.

Sequence analysis reveals that the BTG1 anti-proliferative gene is conserved throughout evolution in its coding and 3' non-coding regions.

The human BTG1 gene (expressing an anti-proliferative function) is an evolutionarily conserved gene homologous to the murine PC3/TIS21 genes. Here, we report the cloning and sequencing of the murine BTG1 coding region and chicken BTG1 cDNA. The putative human and mouse BTG1 proteins are 100% identical; the chicken BTG1 cDNA contains an open reading frame of 170 amino acids with a 91% identity to its human and murine counterparts. The 3'-untranslated region of BTG1 is also highly conserved (82% homology between human and chicken), suggesting that it plays a key role in the regulation of BTG1 expression. These data confirm that BTG1 is phylogenetically highly conserved and that BTG1 and PC3/TIS21 may constitute the first members of a new family of functionally related genes.

Identification and molecular analysis of BANP.

BTG3 belongs to a family of structurally related genes whose biochemical functions remain elusive. In order to investigate the mechanism underlying BTG3-mediated functions, we tried to identify BTG3 potential partners. The use of the yeast 'two-hybrid system', with BTG3 as bait, enabled us to isolate BANP (BTG3 Associated Nuclear Protein). Other commonly used protein-binding assays did not confirm this yeast interaction. However, BANP had never been described before, and this prompted us to further characterise this gene. In this paper, we present data on its molecular organization in mouse, then we speculate on the nature of this nuclear protein, and finally we localise BANP on the human chromosome 16q24 subregion; we discuss the fact that frequent loss of heterozygosity within this region has been observed in different tumours.

Cloning and characterization of a gene encoding a novel immunodominant antigen of Trypanosoma cruzi.

A Trypanosoma cruzi genomic expression library was screened with a pool of sera obtained from chronic chagasic patients. The recombinant antigen (Tc40) isolated from this library reacted with a large number of serum samples of chronic chagasic patients, suggesting that the presence of anti-Tc40 antibodies may be specifically associated to Chagas' disease. The full-length sequence of the Tc40 gene was determined after isolation of genomic and cDNA clones. The Tc40 cDNA includes a large open reading frame (2745 bp-long) that encodes a polypeptide of 100 kDa without any homology with previously described T. cruzi sequences. In contrast with other T. cruzi antigens whose immunodominant B-cell epitopes are composed by amino acid repetitive motifs, Tc40 does not show any amino acid repetition. Antibodies against the Tc40 recombinant protein reacted with three native polypeptides of 100, 41 and 38 kDa which are tightly associated with membranes or cytoskeleton and expressed in all developmental stages of the parasite life cycle. A transcript of 3.9-kb was detected in Northern blot analysis which is large enough to encode a 100 kDa polypeptide. Tc40 genes were mapped on a chromosomal band of 1.1 Mbp and in a few copies per haploid genome in the G strain.

Chromosomal distribution and coding capacity of the human endogenous retrovirus HERV-W family.

Some genomic elements of the multicopy HERV-W endogenous retroviral family have been previously identified in databases. One of them, located on chromosome 7, contains a single complete open reading frame (ORF) putatively encoding an envelope protein. We have experimentally investigated the genomic complexity and coding capacity of the HERV-W family. The human haploid genome contains at least 70, 100, and 30 HERV-W-related gag, pro, and env regions, respectively, widely and heterogeneously dispersed among chromosomes. Using in vitro transcription-translation procedures, three putative HERV-W gag, pro, and env ORFs were detected on chromosomes 3, 6, and 7, respectively, and their sequences analyzed. A 363 amino acid gag ORF containing matrix and carboxy-terminal truncated capsid domains encoded a putative 45-kDa protein. No gag-pro ORF was found, but a pro sequence containing a DTG active site was detected. Finally, the previously described 538 amino acid HERV-W env ORF, located on chromosome 7, was shown to be unique and encoded a putative 80-kDa glycosylated protein. Proteins of molecular mass identical to the one obtained by an in vitro transcription-translation procedure were detected in human placenta, using anti HERV-W Gag- and Env-specific antibodies. The absence of an HERV-W replication-competent provirus versus the existence of HERV-W-related Gag and Env proteins in healthy human placenta is discussed with respect to particle formation, physiology, and pathology.

The covariation between TpA deficiency, CpG deficiency, and G+C content of human isochores is due to a mathematical artifact.

CpG and TpA dinucleotides are underrepresented in the human genome. The CpG deficiency is due to the high mutation rate from C to T in methylated CpG's. The TpA suppression was thought to reflect a counterselection against TpA's destabilizing effect in RNA. Unexpectedly, the TpA and CpG deficiencies vary according to the G+C contents of sequences. It has been proposed that the variation in CpG suppression was correlated with a particular chromatin organization in G+C-rich isochores. Here, we present an improved model of dinucleotide evolution accounting for the overlap between successive dinucleotides. We show that an increased mutation rate from CpG to TpG or CpA induces both an apparent TpA deficiency and a correlation between CpG and TpA deficiencies and G+C content. Moreover, this model shows that the ratio of observed over expected CpG frequency underestimates the real CpG deficiency in G+C-rich sequences. The predictions of our model fit well with observed frequencies in human genomic data. This study suggests that previously published selectionist interpretations of patterns of dinucleotide frequencies should be taken with caution. Moreover, we propose new criteria to identify unmethylated CpG islands taking into account this bias in the measure of CpG depletion.

Synonymous codon usage, accuracy of translation, and gene length in Caenorhabditis elegans.

In many unicellular organisms, invertebrates, and plants, synonymous codon usage biases result from a coadaptation between codon usage and tRNAs abundance to optimize the efficiency of protein synthesis. However, it remains unclear whether natural selection acts at the level of the speed or the accuracy of mRNAs translation. Here we show that codon usage can improve the fidelity of protein synthesis in multicellular species. As predicted by the model of selection for translational accuracy, we find that the frequency of codons optimal for translation is significantly higher at codons encoding for conserved amino acids than at codons encoding for nonconserved amino acids in 548 genes compared between Caenorhabditis elegans and Homo sapiens. Although this model predicts that codon bias correlates positively with gene length, a negative correlation between codon bias and gene length has been observed in eukaryotes. This suggests that selection for fidelity of protein synthesis is not the main factor responsible for codon biases. The relationship between codon bias and gene length remains unexplained. Exploring the differences in gene expression process in eukaryotes and prokaryotes should provide new insights to understand this key question of codon usage.

Determinants of substitution rates in mammalian genes: expression pattern affects selection intensity but not mutation rate.

To determine whether gene expression patterns affect mutation rates and/or selection intensity in mammalian genes, we studied the relationships between substitution rates and tissue distribution of gene expression. For this purpose, we analyzed 2,400 human/rodent and 834 mouse/rat orthologous genes, and we measured (using expressed sequence tag data) their expression patterns in 19 tissues from three development states. We show that substitution rates at nonsynonymous sites are strongly negatively correlated with tissue distribution breadth: almost threefold lower in ubiquitous than in tissue-specific genes. Nonsynonymous substitution rates also vary considerably according to the tissues: the average rate is twofold lower in brain-, muscle-, retina- and neuron-specific genes than in lymphocyte-, lung-, and liver-specific genes. Interestingly, 5' and 3' untranslated regions (UTRs) show exactly the same trend. These results demonstrate that the expression pattern is an essential factor in determining the selective pressure on functional sites in both coding and noncoding regions. Conversely, silent substitution rates do not vary with expression pattern, even in ubiquitously expressed genes. This latter result thus suggests that synonymous codon usage is not constrained by selection in mammals. Furthermore, this result also indicates that there is no reduction of mutation rates in genes expressed in the germ line, contrary to what had been hypothesized based on the fact that transcribed DNA is more efficiently repaired than nontranscribed DNA.

Expression pattern and, surprisingly, gene length shape codon usage in Caenorhabditis, Drosophila, and Arabidopsis.

We measured the expression pattern and analyzed codon usage in 8,133, 1,550, and 2,917 genes, respectively, from Caenorhabditis elegans, Drosophila melanogaster, and Arabidopsis thaliana. In those three species, we observed a clear correlation between codon usage and gene expression levels and showed that this correlation is not due to a mutational bias. This provides direct evidence for selection on silent sites in those three distantly related multicellular eukaryotes. Surprisingly, there is a strong negative correlation between codon usage and protein length. This effect is not due to a smaller size of highly expressed proteins. Thus, for a same-expression pattern, the selective pressure on codon usage appears to be lower in genes encoding long rather than short proteins. This puzzling observation is not predicted by any of the current models of selection on codon usage and thus raises the question of how translation efficiency affects fitness in multicellular organisms.