PCR-based isolation of multigene families: lessons from the avian MHC class IIB.

The amount of sequence data available today highly facilitates the access to genes from many gene families. Primers amplifying the desired genes over a range of species are readily obtained by aligning conserved gene regions, and laborious gene isolation procedures can often be replaced by quicker PCR-based approaches. However, in the case of multigene families, PCR-based approaches bear the often ignored risk of incomplete isolation of family members. This problem is most prominent in gene families with highly variable and thus unpredictable number of gene copies among species, such as in the major histocompatibility complex (MHC). In this study, we (i) report new primers for the isolation of the MHC class IIB (MHCIIB) gene family in birds and (ii) share our experience with isolating MHCIIB genes from an unprecedented number of avian species from all over the avian phylogeny. We report important and usually underappreciated problems encountered during PCR-based multigene family isolation and provide a collection of measures to help significantly improving the chance of successfully isolating complete multigene families using PCR-based approaches.

Cytogenetic Mapping of 5S and 18S rRNAs and H3 Histone Genes in 4 Ancient Proscopiidae Grasshopper Species: Contribution to Understanding the Evolutionary Dynamics of Multigene Families

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Chromosomal organization of the 18S and 5S rRNAs and histone H3 genes in Scarabaeinae coleopterans: insights into the evolutionary dynamics of multigene families and heterochromatin

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Evolution by the birth-and-death process in multigene families of the vertebrate immune system

Concerted evolution is often invoked to explain the diversity and evolution of the multigene families of major histocompatibility complex (MHC) genes and immunoglobulin (Ig) genes. However, this hypothesis has been controversial because the member genes of these families from the same species are not necessarily more closely related to one another than to the genes from different species. To resolve this controversy, we conducted phylogenetic analyses of several multigene families of the MHC and Ig systems. The results show that the evolutionary pattern of these families is quite different from that of concerted evolution but is in agreement with the birth-and-death model of evolution in which new genes are created by repeated gene duplication and some duplicate genes are maintained in the genome for a long time but others are deleted or become nonfunctional by deleterious mutations. We found little evidence that interlocus gene conversion plays an important role in the evolution of MHC and Ig multigene families.

Restriction fragment length polymorphism-coupled domain-directed differential display: a highly efficient technique for expression analysis of multigene families.

In this paper, a reverse-transcriptase PCR-based protocol suitable for efficient expression analysis of multigene families is presented. The method combines restriction fragment length polymorphism (RFLP) technology with a gene family-specific version of mRNA differential display and hence is called "RFLP-coupled domain-directed differential display. "With this method, expression of all members of a multigene family at many different developmental stages, in diverse tissues and even in different organisms, can be displayed on one gel. Moreover, bands of interest, representing gene family members, are directly accessible to sequence analysis, without the need for subcloning. The method thus enables a detailed, high-resolution expression analysis of known gene family members as well as the identification and characterization of new ones. Here the technique was used to analyze differential expression of MADS-box genes in male and female inflorescences of maize (Zea mays ssp. mays). Six different MADS-box genes could be identified, being either specifically expressed in the female sex or preferentially expressed in male or female inflorescences, respectively. Other possible applications of the method are discussed.

Heuristic algorithm for computing reversal distance with multigene families via binary integer programming

Hannenhalli and Pevzner developed the first polynomial-time algorithm for the combinatorial problem of sorting of signed genomic data. Their algorithm solves the minimum number of reversals required for rearranging a genome to another when gene duplication is nonexisting. In this paper, we show how to extend the Hannenhalli-Pevzner approach to genomes with multigene families. We propose a new heuristic algorithm to compute the reversal distance between two genomes with multigene families via the concept of binary integer programming without removing gene duplicates. The experimental results on simulated and real biological data demonstrate that the proposed algorithm is able to find the reversal distance accurately. ©2005 IEEE

Multigene Family Evolution: Perspectives from Insect Chemoreceptors.

Understanding the birth and diversification of multigene families is a fundamental evolutionary problem. I argue for the insect chemoreceptor superfamily as an outstanding model. Although these receptors are currently the preserve of neuroscientists, putative homologous genes exist in diverse animal and plant genomes, implying an ancient origin. Moreover, functional studies suggest that they act as ligand-gated ion channels in both chemosensory and non-chemosensory processes. This family permits synergism of investigations into its structural and regulatory evolution with ecological studies of the selective pressures driving these changes. In addition, sequence divergence in these receptors can be exploited through co-evolutionary and comparative genomics analyses to help to elucidate their 3D structure and signaling mechanisms, and to reveal experimentally-accessible candidate loci to explore the genetic basis of adaptation.

PR gene families of citrus: their organ specific-biotic and abiotic inducible expression profiles based on ESTs approach

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Molecular evolution and phylogeny of elapid snake venom three-finger toxins

Animal venom components are of considerable interest to researchers across a wide variety of disciplines, including molecular biology, biochemistry, medicine, and evolutionary genetics. The three-finger family of snake venom peptides is a particularly interesting and biochemically complex group of venom peptides, because they are encoded by a large multigene family and display a diverse array of functional activities. In addition, understanding how this complex and highly varied multigene family evolved is an interesting question to researchers investigating the biochemical diversity of these peptides and their impact on human health. Therefore, the purpose of our study was to investigate the long-term evolutionary patterns exhibited by these snake venom toxins to understand the mechanisms by which they diversified into a large, biochemically diverse, multigene family. Our results show a much greater diversity of family members than was previously known, including a number of subfamilies that did not fall within any previously identified groups with characterized activities. In addition, we found that the long-term evolutionary processes that gave rise to the diversity of three-finger toxins are consistent with the birth-and-death model of multigene family evolution. It is anticipated that this three-finger toxin toolkit will prove to be useful in providing a clearer picture of the diversity of investigational ligands or potential therapeutics available within this important family.

Isolation of a novel G-protein gamma-subunit from Arabidopsis thaliana and its interaction with G beta

There is increasing evidence that heterotrimeric G-proteins (G-proteins) are involved in many plant processes including phytohormone response, pathogen defence and stomatal control. In animal systems, each of the three G-protein subunits belong to large multigene families; however, few subunits have been isolated from plants. Here we report the cloning of a second plant G-protein γ-subunit (AGG2) from Arabidopsis thaliana. The predicted AGG2 protein sequence shows 48% identity to the first identified Arabidopsis Gγ-subunit, AGG1. Furthermore, AGG2 contains all of the conserved characteristics of γ-subunits including a small size (100 amino acids, 11.1 kDa), C-terminal CAAX box and a N-terminal α-helix region capable of forming a coiled-coil interaction with the β-subunit. A strong interaction between AGG2 and both the tobacco (TGB1) and Arabidopsis (AGB1) β-subunits was observed in vivo using the yeast two-hybrid system. The strong association between AGG2 and AGB1 was confirmed in vitro. Southern and Northern analyses showed that AGG2 is a single copy gene in Arabidopsis producing two transcripts that are present in all tissues tested. The isolation of a second γ-subunit from A. thaliana indicates that plant G-proteins, like their mammalian counterparts, may form different heterotrimer combinations that presumably regulate multiple signal transduction pathways.

New insights into the long-term evolutionary history of the avian MHC

SummaryGene duplication and neofunctidnalization are important processes in the evolution of phenotypic complexity. They account for important evolutionary novelties that confer ecological adaptation, such as the major histocompatibility complex (MHC), a multigene family with a central role in vertebrates' adaptive immune system. Multigene families, which evolved in large part through duplication, represent promising systems to study the still strongly depbated relative roles of neutral and adaptive processes in the evolution of phenotypic complexity. Detailed knowledge on ecological function and a well-characterized evolutionary history place the mammals' MHC amongst ideal study systems. However mammalian MHCs usually encompass several million base pairs and hold a large number of functional and non-functional duplicate genes, which makes their study complex. Avian MHCs on the other hand are usually way more compact, but the reconstruction of. their evolutionary history has proven notoriously difficult. However, no focused attempt has been undertaken so far to study the avian MHC evolutionary history in a broad phylogenetic context and using adequate gene regions.In the present PhD, we were able to make important contributions to the understanding of the long-term evolution of the avian MHC class II Β (MHCI1B). First, we isolated and characterized MHCIIB genes in barn owl (Tyto alba?, Strigiformes, Tytonidae), a species from an avian lineage in which MHC has not been studied so far. Our results revealed that with only two functional MHCIIB genes the MHC organization of barn owl may be similar to the 'minimal essential' MHC of chicken (Gallus gallus), indicating that simple MHC organization may be ancestral to birds. Taking advantage of the sequence information from barn owl, we studied the evolution of MHCIIB genes in 13 additional species of 'typical' owls (Strigiformes, Strigidae). Phylogenetic analyses revealed that according to their function, in owls the peptide-binding region (PBR) encoding exon 2 and the non-PBR encoding exon 3 evolve by different patterns. Exon 2 exhibited an evolutionary history of positive selection and recombination, while exon 3 traced duplication history and revealed two paralogs evolving divergently from each other in owls, and in a shorebird, the great snipe {Gallinago media). The results from exon 3 were the first ever from birds to demonstrate gene orthology in species that diverged tens of millions of years ago, and strongly questioned whether the taxa studied before provided an adequate picture of avian MHC evolution. In a follow-up study, we aimed at explaining a striking pattern revealed by phylogenetic trees analyzing the owl sequences along with MHCIIB sequences from other birds: One owl paralog (termed DAB1) grouped with sequences of passerines and falcons, while the other (DAB2) grouped with wildfowl, penguins and birds of prey. This could be explained by either a duplication event preceding the evolution of these bird orders, or by convergent evolution of similar sequences in a number of orders. With extensive phylogenetic analyses we were able to show, that indeed a duplication event preceeded the major avian radiation -100 my ago, and that following this duplication, the paralogs evolved under positive selection. Furthermore, we showed that the divergently evolving amino acid residues in the MHCIIB-encoded β-chain potentially interact with the MHCI I α-chain, and that molecular coevolution of the interacting residues may have been involved in the divergent evolution of the MHCIIB paralogs.The findings of this PhD are of particular interest to the understanding of the evolutionary history of the avian MHC and, by providing essential information on long-term gene history in the avian MHC, open promising perspectives for advances in the understanding of the evolution of multigene families in general, and for avian MHC organization in particular. Amongst others I discuss the importance of including protein structure in the phylogenetic study of multigene families, and the roles of ecological versus molecular selection pressures. I conclude by providing a population genomic perspective on avian MHC, which may serve as a basis for future research to investigate the relative roles of neutral processes involving effective population size effects and of adaptation in the evolution of avian MHC diversity and organization.RésuméLa duplication de gènes et leur néo-fonctionnalisation sont des processus importants dans l'évolution de la complexité phénotypique. Ils sont impliqués dans l'apparition d'importantes nouveautés évolutives favorisant l'adaptation écologique, comme c'est le cas pour le complexe majeur d'histocompatibilité

Identification of a new cancer/testis gene family, CT47, among expressed multicopy genes on the human X chromosome.

Cancer/testis (CT) genes are normally expressed in germ cells only, yet are reactivated and expressed in some tumors. Of the approximately 40 CT genes or gene families identified to date, 20 are on the X chromosome and are present as multigene families, many with highly conserved members. This indicates that novel CT gene families may be identified by detecting duplicated expressed genes on chromosome X. By searching for transcript clusters that map to multiple locations on the chromosome, followed by in silico analysis of their gene expression profiles, we identified five novel gene families with testis-specific expression and >98% sequence identity among family members. The expression of these genes in normal tissues and various tumor cell lines and specimens was evaluated by qualitative and quantitative RT-PCR, and a novel CT gene family with at least 13 copies was identified on Xq24, designated as CT47. mRNA expression of CT47 was found mainly in the testes, with weak expression in the placenta. Brain tissue was the only positive somatic tissue tested, with an estimated CT47 transcript level 0.09% of that found in testis. Among the tumor specimens tested, CT47 expression was found in approximately 15% of lung cancer and esophageal cancer specimens, but not in colorectal cancer or breast cancer. The putative CT47 protein consists of 288 amino acid residues, with a C-terminus rich in alanine and glutamic acid. The only species other than human in which a gene homologous to CT47 has been detected is the chimpanzee, with the predicted protein showing approximately 80% identity in its carboxy terminal region.

Comparative genomic analysis of the odorant-binding protein family in 12 Drosophila genomes: purifying selection and birth-and-death evolution

Background: Chemoreception is a widespread mechanism that is involved in critical biologic processes, including individual and social behavior. The insect peripheral olfactory system comprises three major multigene families: the olfactory receptor (Or), the gustatory receptor (Gr), and the odorant-binding protein (OBP) families. Members of the latter family establish the first contact with the odorants, and thus constitute the first step in the chemosensory transduction pathway.Results: Comparative analysis of the OBP family in 12 Drosophila genomes allowed the identification of 595 genes that encode putative functional and nonfunctional members in extant species, with 43 gene gains and 28 gene losses (15 deletions and 13 pseudogenization events). The evolution of this family shows tandem gene duplication events, progressive divergence in DNA and amino acid sequence, and prevalence of pseudogenization events in external branches of the phylogenetic tree. We observed that the OBP arrangement in clusters is maintained across the Drosophila species and that purifying selection governs the evolution of the family; nevertheless, OBP genes differ in their functional constraints levels. Finally, we detect that the OBP repertoire evolves more rapidly in the specialist lineages of the Drosophila melanogaster group (D. sechellia and D. erecta) than in their closest generalists.Conclusion: Overall, the evolution of the OBP multigene family is consistent with the birth-and-death model. We also found that members of this family exhibit different functional constraints, which is indicative of some functional divergence, and that they might be involved in some of the specialization processes that occurred through the diversification of the Drosophila genus.

Functional analysis of Plasmodium vivax VIR proteins reveals different subcellular localizations and cytoadherence to the ICAM-1 endothelial receptor.

The subcellular localization and function of variant subtelomeric multigene families in Plasmodium vivax remain vastly unknown. Among them, the vir superfamily is putatively involved in antigenic variation and in mediating adherence to endothelial receptors. In the absence of a continuous in vitro culture system for P. vivax, we have generated P. falciparum transgenic lines expressing VIR proteins to infer location and function. We chose three proteins pertaining to subfamilies A (VIR17), C (VIR14) and D (VIR10), with domains and secondary structures that predictably traffic these proteins to different subcellular compartments. Here, we showed that VIR17 remained inside the parasite and around merozoites, whereas VIR14 and VIR10 were exported to the membrane of infected red blood cells (iRBCs) in an apparent independent pathway of Maurer's clefts. Remarkably, VIR14 was exposed at the surface of iRBCs and mediated adherence to different endothelial receptors expressed in CHO cells under static conditions. Under physiological flow conditions, however, cytoadherence was only observed to ICAM-1, which was the only receptor whose adherence was specifically and significantly inhibited by antibodies against conserved motifs of VIR proteins. Immunofluorescence studies using these antibodies also showed different subcellular localizations of VIR proteins in P. vivax-infected reticulocytes from natural infections. These data suggest that VIR proteins are trafficked to different cellular compartments and functionally demonstrates that VIR proteins can specifically mediate cytoadherence to the ICAM-1 endothelial receptor.

Functional and evolutionary analysis of DXL1, a non-essential gene encoding a 1-deoxy-D-xylulose 5-phosphate synthase like protein in Arabidopsis thaliana

The synthesis of 1-deoxy-D-xylulose 5-phosphate (DXP), catalyzed by the enzyme DXP synthase (DXS), represents a key regulatory step of the 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway for isoprenoid biosynthesis. In plants DXS is encoded by small multigene families that can be classified into, at least, three specialized subfamilies. Arabidopsis thaliana contains three genes encoding proteins with similarity to DXS, including the well-known DXS1/CLA1 gene, which clusters within subfamily I. The remaining proteins, initially named DXS2 and DXS3, have not yet been characterized. Here we report the expression and functional analysis of A. thaliana DXS2. Unexpectedly, the expression of DXS2 failed to rescue Escherichia coli and A. thaliana mutants defective in DXS activity. Coherently, we found that DXS activity was negligible in vitro, being renamed as DXL1 following recent nomenclature recommendation. DXL1 is targeted to plastids as DXS1, but shows a distinct expression pattern. The phenotypic analysis of a DXL1 defective mutant revealed that the function of the encoded protein is not essential for growth and development. Evolutionary analyses indicated that DXL1 emerged from DXS1 through a recent duplication apparently specific of the Brassicaceae lineage. Divergent selective constraints would have affected a significant fraction of sites after diversification of the paralogues. Furthermore, amino acids subjected to divergent selection and likely critical for functional divergence through the acquisition of a novel, although not yet known, biochemical function, were identified. Our results provide with the first evidences of functional specialization at both the regulatory and biochemical level within the plant DXS family.