Evidence for a genetic sex determination in Cnidaria, the Mediterranean red coral (Corallium rubrum).

Sexual reproduction is widespread among eukaryotes, and the sex-determining processes vary greatly among species. While genetic sex determination (GSD) has been intensively described in bilaterian species, no example has yet been recorded among non-bilaterians. However, the quasi-ubiquitous repartition of GSD among multicellular species suggests that similar evolutionary forces can promote this system, and that these forces could occur also in non-bilaterians. Studying sex determination across the range of Metazoan diversity is indeed important to understand better the evolution of this mechanism and its lability. We tested the existence of sex-linked genes in the gonochoric red coral (Corallium rubrum, Cnidaria) using restriction site-associated DNA sequencing. We analysed 27 461 single nucleotide polymorphisms (SNPs) in 354 individuals from 12 populations including 53 that were morphologically sexed. We found a strong association between the allele frequencies of 472 SNPs and the sex of individuals, suggesting an XX/XY sex-determination system. This result was confirmed by the identification of 435 male-specific loci. An independent test confirmed that the amplification of these loci enabled us to identify males with absolute certainty. This is the first demonstration of a GSD system among non-bilaterian species and a new example of its convergence in multicellular eukaryotes.

Genomic resources and their influence on the detection of the signal of positive selection in genome scans.

Genome scans represent powerful approaches to investigate the action of natural selection on the genetic variation of natural populations and to better understand local adaptation. This is very useful, for example, in the field of conservation biology and evolutionary biology. Thanks to Next Generation Sequencing, genomic resources are growing exponentially, improving genome scan analyses in non-model species. Thousands of SNPs called using Reduced Representation Sequencing are increasingly used in genome scans. Besides, genome sequences are also becoming increasingly available, allowing better processing of short-read data, offering physical localization of variants, and improving haplotype reconstruction and data imputation. Ultimately, genome sequences are also becoming the raw material for selection inferences. Here, we discuss how the increasing availability of such genomic resources, notably genome sequences, influences the detection of signals of selection. Mainly, increasing data density and having the information of physical linkage data expand genome scans by (i) improving the overall quality of the data, (ii) helping the reconstruction of demographic history for the population studied to decrease false-positive rates and (iii) improving the statistical power of methods to detect the signal of selection. Of particular importance, the availability of a high-quality reference genome can improve the detection of the signal of selection by (i) allowing matching the potential candidate loci to linked coding regions under selection, (ii) rapidly moving the investigation to the gene and function and (iii) ensuring that the highly variable regions of the genomes that include functional genes are also investigated. For all those reasons, using reference genomes in genome scan analyses is highly recommended.

The red coral (Corallium rubrum) transcriptome: a new resource for population genetics and local adaptation studies.

The question of species survival and evolution in heterogeneous environments has long been a subject for study. Indeed, it is often difficult to identify the molecular basis of adaptation to contrasted environments, and nongenetic effects increase the difficulty to disentangle fixed effects, such as genetic adaptation, from variable effects, such as individual phenotypic plasticity, in adaptation. Nevertheless, this question is also of great importance for understanding the evolution of species in a context of climate change. The red coral (Corallium rubrum) lives in the Mediterranean Sea, where at depths ranging from 5 to 600 m, it meets very contrasted thermal conditions. The shallowest populations of this species suffered from mortality events linked with thermal anomalies that have highlighted thermotolerance differences between individuals. We provide here a new transcriptomic resource, as well as candidate markers for the study of local adaptation. We sequenced the transcriptome of six individuals from 5 m and six individuals from 40 m depth at the same site of the Marseilles bay, after a period of common garden acclimatization. We found differential expression maintained between the two depths even after common garden acclimatization, and we analysed the polymorphism pattern of these samples. We highlighted contigs potentially implicated in the response to thermal stress, which could be good candidates for the study of thermal adaptation for the red coral. Some of these genes are also involved in the response to thermal stress in other corals. Our method enables the identification of candidate loci of local adaptation useful for other nonmodel organisms.

Compound Poisson approximation and testing for gene clusters with multigene families.

We present in this article a compound Poisson approximation for computing probabilities involved in significance tests for conserved genomic regions between different species. We consider the case when the conserved genomic regions are found by the reference region approach. An important aspect of our computations is the fact that we are taking into account the existence of multigene families. We obtain convergence results for the error of our approximation by using the Stein-Chen method for compound Poisson approximation.

A rigorous method for multigenic families' functional annotation: the peptidyl arginine deiminase (PADs) proteins family example.

BACKGROUND: large scale and reliable proteins' functional annotation is a major challenge in modern biology. Phylogenetic analyses have been shown to be important for such tasks. However, up to now, phylogenetic annotation did not take into account expression data (i.e. ESTs, Microarrays, SAGE, ...). Therefore, integrating such data, like ESTs in phylogenetic annotation could be a major advance in post genomic analyses. We developed an approach enabling the combination of expression data and phylogenetic analysis. To illustrate our method, we used an example protein family, the peptidyl arginine deiminases (PADs), probably implied in Rheumatoid Arthritis. RESULTS: the analysis was performed as follows: we built a phylogeny of PAD proteins from the NCBI's NR protein database. We completed the phylogenetic reconstruction of PADs using an enlarged sequence database containing translations of ESTs contigs. We then extracted all corresponding expression data contained in EST database This analysis allowed us 1/To extend the spectrum of homologs-containing species and to improve the reconstruction of genes' evolutionary history. 2/To deduce an accurate gene expression pattern for each member of this protein family. 3/To show a correlation between paralogous sequences' evolution rate and pattern of tissular expression. CONCLUSION: coupling phylogenetic reconstruction and expression data is a promising way of analysis that could be applied to all multigenic families to investigate the relationship between molecular and transcriptional evolution and to improve functional annotation.

DOK4 and DOK5: new Dok-related genes expressed in human T cells.

Dok proteins are adapter proteins involved in signal transduction. Several intracellular proteins expressed in lymphocytes meet the criteria of membrane-associated adapter proteins such as members of the Dok family. To understand the role and the formation of multiprotein networks involving Dok proteins in T lymphocytes, we search for potential additional members of this family. Here, we describe the two new human dok-related genes DOK4 and DOK5 and present data showing the expression of DOK4 and DOK5 genes in T cells. These genes are the orthologues of mouse Dok4 and Dok5 genes. Based on analysis of phylogenetic trees and exon/intron structure of Dok family members, DOK4 and DOK5 define a subfamily within dok genes distinct from DOK1, DOK2 and DOK3. So, Dok-4 and Dok-5 molecules constitute a new group of adapter proteins in T cells, requiring further functional analysis.

Identification and characterization of a beta proteasome subunit cluster in the Japanese pufferfish (Fugu rubripes).

The low molecular mass polypeptide (LMP2, LMP7, and MECL-1) genes code for beta-type subunits of the proteasome, a multimeric complex that degrades proteins into peptides as part of the MHC class I-mediated Ag-presenting pathway. These gene products are up-regulated in response to infection by IFN-gamma and replace the corresponding constitutively expressed subunits (X, Y, and Z) during the immune response. In humans, the LMP2 and LMP7 genes both reside within the class II region of the MHC (6p21.3), while MECL-1 is located at 16q22.1. In the present study, we have identified all three IFN-gamma-regulated beta-type proteasome subunits in Fugu, which are present as a cluster within the Fugu MHC class I region. We show that in this species, LMP7, LMP2, and MECL-1 are linked. Also within this cluster is an LMP2-like subunit (which seems specific to all teleosts tested to date) and a closely linked LMP7 pseudogene, indicating that within Fugu and potentially other teleosts, there has been an additional regional duplication involving these genes.

"Paleogenomics": looking in the past to the future.

The complete sequence of the human and other vertebrate and nonvertebrate genomes provide a wealth of information on the organization, relationships and evolution of the metazoans. Soon the fine structure of our innermost biological identity will be unveiled and what has so far remained deep and secret will shine like an unearthed treasure and shape and fuel our future quests. A key treasure, for many molecular scientists interested in molecular evolution and development would be the knowledge of the genome of the ancestral precursor of all metazoans. In the absence of fossil DNA, this knowledge will forever remain a yearning for dreamy molecular biologists. And yet, will not the power of deduction and reconstitution of information gained through man's sophisticated technologies one day recreate a "virtual" metazoan ancestor?

let-756, a C. elegans fgf essential for worm development.

In vertebrates, Fibroblast Growth Factors (FGFs) and their receptors are involved in various developmental and pathological processes, including neoplasia. The number of FGFs and their large range of activities have made the understanding of their precise functions difficult. Investigating their biology in other species might be enlightening. A sequence encoding a putative protein presenting 30-40% identity with the conserved core of vertebrate FGFs has been identified by the C. elegans sequencing consortium. We show here that this gene is transcribed and encodes a putative protein of 425 amino acids (aa). The gene is expressed at all stages of development beyond late embryogenesis, peaking at the larval stages. Loss-of-function mutants of the let-756 gene are rescued by the wild type fgf gene in germline transformation experiments. Two partial loss-of-function alleles, s2613 and s2809, have a mutation that replaces aa 317 by a stop. The truncated protein retains the FGF core but lacks a C-termins portion. These worms are small and develop slowly into clear and scrawny, yet viable and fertile adults. A third allele, s2887, is inactivated by an inversion that disrupts the first exon. It causes a developmental arrest early in the larval stages. Thus, in contrast to the other nematode fgf gene egl-17, let-756/fgf is essential for worm development.

The family of Caenorhabditis elegans tyrosine kinase receptors: similarities and differences with mammalian receptors.

Transmembrane receptors with tyrosine kinase activity (RTK) constitute a superfamily of proteins present in all metazoans that is associated with the control and regulation of cellular processes. They have been the focus of numerous studies and are a good subject for comparative analyses of multigene families in different species aimed at understanding metazoan evolution. The sequence of the genome of the nematode worm Caenorhabditis elegans is available. This offers a good opportunity to study the superfamily of nematode RTKs in its entirety and to compare it with its mammalian counterpart. We show that the C. elegans RTKs constitute various groups with different phylogenetic relationships with mammalian RTKs. A group of four RTKs show structural similarity with the three mammalian receptors for the vascular endothelial growth factors. Another group comprises RTKs with a short extracellular region, a feature not known in mammals; the genes encoding these RTKs are clustered on chromosome II with other gene families, including genes encoding chitinase-like proteins. Most of the C. elegans RTKs have no direct orthologous relationship with any mammalian RTK, providing an illustration of the importance of the separate evolution of the different phyla.

A susceptibility region for myasthenia gravis extending into the HLA-class I sector telomeric to HLA-C.

We have analyzed a series of HLA region markers in 207 UK Caucasoids with early-onset myasthenia gravis (EOMG, onset before age 40), where there is a strong female bias. The well known associations with HLA-DR3 and -B8 have now proved to be significantly stronger in the 165 females than in the 42 males. In patients (of either sex) lacking -DR3, there was also a significant increase in HLA-DR2. Although the muscle weakness in EOMG is clearly mediated by autoantibodies, the associations are consistently stronger with HLA-B8 (in class I) than with HLADR3 (in class II), as confirmed here. We therefore typed 87-137 cases for polymorphisms at four loci in the intervening class III region, and also at three in the adjacent stretch of class I. At each locus, one allele tended to co-occur with HLA-B8 and showed strong and highly significant associations in the patients. There appeared to be a region of maximal susceptibility extending from HSP70 (in class III) past HLA-B and HLA-C at least 600 kb telomerically into the class I region, which is now being mapped in detail. Any candidate genes here that act shortly after puberty may allow more precise localization of susceptibility.

Structure and evolution of the extended B7 family.

Here, Joëlle Henry and colleagues explore structural and evolutionary relationships between the B7 costimulator molecules and a growing number of molecules encoded within the major histocompatibility complex. They propose that B7 and MHC genes are derived from a common ancestor, with several members of this large gene family possibly having pivotal influences on T-cell activation.

The MHC big bang.

The human Major Histocompatibility Complex (MHC) shares similarities with three other chromosome regions in human. This could be the vestige of ancestral large scale duplications. We discuss here the possibility i) that these duplications occurred during two rounds of tetraploidization supposed to have taken place during chordate evolution before the jawed vertebrate radiation, and ii) that one of the quadruplicate regions, relaxed of functional constraints, gave rise to the vertebrate MHC by a quick round of gene cis-duplication and cis-exon shuffling. These different rounds of cis-duplications and exon shufflings allowed the emergence of new genes participating in novel biological functions i.e. adaptive immune responses. Cis-duplications and cis-exon shufflings are ongoing processes in the evolution of some of these genes in this region as they have occurred and were fixed at different times and in different lineages during vertebrate evolution. In contrast, other genes within the MHC have remained stable since the emergence of jawed vertebrates.

B30.2-like domain proteins: update and new insights into a rapidly expanding family of proteins.

The B30.2 domain is a conserved region of around 170 amino acids associated with several different protein domains, including the immunoglobulin folds of butyrophilin and the RING finger domain of ret finger protein. We recently reported several novel members of this family as well as previously undescribed protein families possessing the B30.2 domain. Many proteins have subsequently been found to possess this domain, including pyrin/marenostrin and the midline 1 (MID1) protein. Mutations in the B30.2 domain of pyrin/marenostrin are implicated in familial Mediterranean fever, and partial loss of the B30.2 domain of MID1 is responsible for Opitz G/BBB syndrome, characterized by developmental midline defects. In this study, we scrutinized the available sequence data bases for the identification of novel B30.2 domain proteins using highly sensitive database-searching tools. In addition, we discuss the chromosomal localization of genes in the B30.2 family, since the encoded proteins are likely to be involved in other forms of periodic fever, autoimmune, and genetic diseases.

Ancient large-scale genome duplications: phylogenetic and linkage analyses shed light on chordate genome evolution.

Paralogous genes from several families were found in four human chromosome regions (4p16, 5q33-35, 8p12-21, and 10q24-26), suggesting that their common ancestral region underwent several rounds of large-scale duplication. Searches in the EMBL databases, followed by phylogenetic analyses, showed that cognates (orthologs) of human duplicated genes can be found in other vertebrates, including bony fishes. In contrast, within each family, only one gene showing the same high degree of similarity with all the duplicated mammalian genes was found in nonvertebrates (echinoderms, insects, nematodes). This indicates that large-scale duplications occurred after the echinoderms/chordates split and before the bony vertebrate radiation. It has been suggested that two rounds of gene duplication occurred in the vertebrate lineage after the separation of Amphioxus and craniate (vertebrates + Myxini) ancestors. Before these duplications, the genes that have led to the families of paralogous genes in vertebrates must have been physically linked in the craniate ancestor. Linkage of some of these genes can be found in the Drosophila melanogaster and Caenorhabditis elegans genomes, suggesting that they were linked in the triploblast Metazoa ancestor.

A first map of the porcine major histocompatibility complex class I region.

A map of the SLA complex, or swine major histocompatibility complex (MHC), class I region was constructed by alignment of yeast artificial chromosomes (YACs) harboring MHC class I genes as well as anchor genes already mapped within the human MHC complex (HLA). Five YACs containing 9 anchor genes built a contig of about 1.0-1.2 Mb between the SLA class III BAT1 locus and the olfactory receptor-like genes OLF42. Ten different SLA class I sequences, including putative allelic forms of published classical and non-classical SLA class I genes, were assigned to the 400-kb enclosing centromeric part of the contig. Three additional YACs comprising the OLF89 genes and two YACs containing the butyrophilin gene were located telomeric to the contig. Comparison between the human and porcine MHC complexes showed a perfect conserved order of anchor genes, whereas no orthologous relationships were found for the class I loci.

Cloning, structural analysis, and mapping of the B30 and B7 multigenic families to the major histocompatibility complex (MHC) and other chromosomal regions.

We present the cloning, structural analysis, and mapping of new members belonging to two multigenic families, the B30-RING finger family and the B7.1-B7.2 family, as well as two genes derived by exon shuffling from members of these families. Eight new members were found and three of them map to the human major histocompatibilitiy complex (MHC) region. Phylogenic and physical mapping analysis allowed us to decipher the evolutionary story of these two multigenic families and to shed light on the evolution of the MHC region. We also show that a deductive analysis can be used to predict the existence of a given gene.

B30.2-like domain proteins: a growing family.

The B30.2 domain is a conserved domain of around 170 amino acids. It is found associated with different protein domains: immunoglobulin domain in the case of butyrophilin and Ring Finger domain in the case of Ret Finger Protein. B30.2 should therefore be considered a migratory domain. We here report new members of these families as well as new protein families having the B30.2 domain, and we tentatively propose a general function for this domain.