Multigene families and the evolution of complexity

Summary Higher organisms are complex, and their developmental processes are controlled by the sequential expression of genes that often form multigene families. Facts are surveyed on how functional diversity of genes is related to duplication of genes or segments of genes, by emphasizing that diversity is often enhanced by alternate splicing and proteolytic cleavage involving duplicated genes or gene segments. Analyses of a population genetics model for the origin of gene families suggest that positive Darwinian selection is needed for acquiring gene families with desirable functions. Based on these considerations, examples that show acceleration of amino acid substitution relative to synonymous change during evolutionary processes are surveyed. Some of such examples strongly suggest that positive selection has worked. In other cases it is difficult to judge whether or not acceleration is caused by positive Darwinian selection. As a general pattern, acceleration of amino acid substitution is often found to be related to gene duplication. It is thought that complexity and diversity of gene function have been advantageous in the long evolutionary course of higher organisms.     

Critical review of NGS analyses for de novo genotyping multigene families

The genotyping of highly polymorphic multigene families across many individuals used to be a particularly challenging task because of methodological limitations associated with traditional approaches. Next‐generation sequencing (NGS) can overcome most of these limitations, and it is increasingly being applied in population genetic studies of multigene families. Here, we critically review NGS bioinformatic approaches that have been used to genotype the major histocompatibility complex (MHC) immune genes, and we discuss how the significant advances made in this field are applicable to population genetic studies of gene families. Increasingly, approaches are introduced that apply thresholds of sequencing depth and sequence similarity to separate alleles from methodological artefacts. We explain why these approaches are particularly sensitive to methodological biases by violating fundamental genotyping assumptions. An alternative strategy that utilizes ultra‐deep sequencing (hundreds to thousands of sequences per amplicon) to reconstruct genotypes and applies statistical methods on the sequencing depth to separate alleles from artefacts appears to be more robust. Importantly, the ‘degree of change’ (DOC) method avoids using arbitrary cut‐off thresholds by looking for statistical boundaries between the sequencing depth for alleles and artefacts, and hence, it is entirely repeatable across studies. Although the advances made in generating NGS data are still far ahead of our ability to perform reliable processing, analysis and interpretation, the community is developing statistically rigorous protocols that will allow us to address novel questions in evolution, ecology and genetics of multigene families. Future developments in third‐generation single molecule sequencing may potentially help overcome problems that still persist in de novo multigene amplicon genotyping when using current second‐generation sequencing approaches.     

Clade classification of monolignol biosynthesis gene family members reveals target genes to decrease lignin in Lolium perenne

In monocots, lignin content has a strong impact on the digestibility of the cell wall fraction. Engineering lignin biosynthesis requires a profound knowledge of the role of paralogues in the multigene families that constitute the monolignol biosynthesis pathway. We applied a bioinformatics approach for genome‐wide identification of candidate genes in Lolium perenne that are likely to be involved in the biosynthesis of monolignols. More specifically, we performed functional subtyping of phylogenetic clades in four multigene families: 4CL, COMT, CAD and CCR. Essential residues were considered for functional clade delineation within these families. This classification was complemented with previously published experimental evidence on gene expression, gene function and enzymatic activity in closely related crops and model species. This allowed us to assign functions to novel identified L. perenne genes, and to assess functional redundancy among paralogues. We found that two 4CL paralogues, two COMT paralogues, three CCR paralogues and one CAD gene are prime targets for genetic studies to engineer developmentally regulated lignin in this species. Based on the delineation of sequence conservation between paralogues and a first analysis of allelic diversity, we discuss possibilities to further study the roles of these paralogues in lignin biosynthesis, including expression analysis, reverse genetics and forward genetics, such as association mapping. We propose criteria to prioritise paralogues within multigene families and certain SNPs within these genes for developing genotyping assays or increasing power in association mapping studies. Although L. perenne was the target of the analyses presented here, this functional subtyping of phylogenetic clades represents a valuable tool for studies investigating monolignol biosynthesis genes in other monocot species.     

An Extension of a Model for the Evolution of Multigene Families by Unequal Crossing over

Evolution of a multigene family is studied from the standpoint of population genetics. It is assumed that the multigene family is undergoing continuous interchromosomal unequal crossing over, mutation and random frequency drift. The equilibrium properties of the probability of gene identity (clonality) are investigated, using two measures: identity probability within and between chromosomes. The measures represent homogeneity of genes within a family in one chromosome and similarity of gene families between two homologous chromosomes. The means, the variances and the covariance of these two measures of identity probability are obtained by using the diffusion equation method. It is shown that the means and the variances are generally smaller than those predicted in the previous model assuming intrachromosomal (sister chromatid) unequal crossing over (Ohta 1978a,b).     

The Cohesive Population Genetics of Molecular Drive

The long-term population genetics of multigene families is influenced by several biased and unbiased mechanisms of nonreciprocal exchanges (gene conversion, unequal exchanges, transposition) between member genes, often distributed on several chromosomes. These mechanisms cause fluctuations in the copy number of variant genes in an individual and lead to a gradual replacement of an original family of n genes (A) in N number of individuals by a variant gene (a). The process for spreading a variant gene through a family and through a population is called molecular drive. Consideration of the known slow rates of nonreciprocal exchanges predicts that the population variance in the copy number of gene a per individual is small at any given generation during molecular drive. Genotypes at a given generation are expected only to range over a small section of all possible genotypes from one extreme (n number of A) to the other (n number of a). A theory is developed for estimating the size of the population variance by using the concept of identity coefficients. In particular, the variance in the course of spreading of a single mutant gene of a multigene family was investigated in detail, and the theory of identity coefficients at the state of steady decay of genetic variability proved to be useful. Monte Carlo simulations and numerical analysis based on realistic rates of exchange in families of known size reveal the correctness of the theoretical prediction and also assess the effect of bias in turnover. The population dynamics of molecular drive in gradually increasing the mean copy number of a variant gene without the generation of a large variance (population cohesion) is of significance regarding potential interactions between natural selection and molecular drive.     

The phylogenetic origins of natural killer receptors and recognition: relationships,possibilities, and realities

Natural killer (NK) cells affect a form of innate immunity that recognizes and eliminates cells that are infected with certain viruses or have undergone malignant transformation. In mammals, this recognition can be mediated through immunoglobulin- (Ig) and/or lectin-type NK receptors (NKRs). NKR genes in mammals range from minimally polymorphic single-copy genes to complex multigene families that exhibit high levels of haplotypic complexity and exhibit significant interspecific variation. Certain single-copy NKR genes that are present in one mammal are present as expanded multigene families in other mammals. These observations highlight NKRs as one of the most rapidly evolving eukaryotic gene families and likely reflect the influence of pathogens, especially viruses, on their evolution. Although well characterized in human and mice, cytotoxic cells that are functionally similar to NK cells have been identified in species ranging from birds to reptiles, amphibians and fish. Although numerous receptors have been identified in non-mammalian vertebrates that share structural relationships with mammalian NKRs, functionally defining these lower vertebrate molecules as NKRs is confounded by methodological and interpretive complexities. Nevertheless, several lines of evidence suggest that NK-type function or its equivalent has sustained a long evolutionary history throughout vertebrate species.     

Theoretical population genetics of repeated genes forming a multigene family

The evolution of repeated genes forming a multigene family in a finite population is studied with special reference to the probability of gene identity, i.e., the identity probability of two gene units chosen from the gene family. This quantity is called clonality and is defined as the sum of squares of the frequencies of gene lineages in the family. The multigene family is undergoing continuous unequal somatic crossing over, ordinary interchromosomal crossing over, mutation and random frequency drift. Two measures of clonality are used: clonality within one chromosome and that between two different chromosomes. The equilibrium properties of the means, the variances and the covariance of the two measures of clonality are investigated by using the diffusion equation method under the assumption of constant number of gene units in the multigene family. Some models of natural selection based on clonality are considered. The possible significance of the variance and covariance of clonality among the chromosomes on the adaptive differentiation of gene families such as those producing antibodies is discussed.     

Genome-wide Comparative Analysis of Annexin Superfamily in Plants

Most annexins are calcium-dependent, phospholipid-binding proteins with suggested functions in response to environmental stresses and signaling during plant growth and development. They have previously been identified and characterized in Arabidopsis and rice, and constitute a multigene family in plants. In this study, we performed a comparative analysis of annexin gene families in the sequenced genomes of Viridiplantae ranging from unicellular green algae to multicellular plants, and identified 149 genes. Phylogenetic studies of these deduced annexins classified them into nine different arbitrary groups. The occurrence and distribution of bona fide type II calcium binding sites within the four annexin domains were found to be different in each of these groups. Analysis of chromosomal distribution of annexin genes in rice, Arabidopsis and poplar revealed their localization on various chromosomes with some members also found on duplicated chromosomal segments leading to gene family expansion. Analysis of gene structure suggests sequential or differential loss of introns during the evolution of land plant annexin genes. Intron positions and phases are well conserved in annexin genes from representative genomes ranging from Physcomitrella to higher plants. The occurrence of alternative motifs such as K/R/HGD was found to be overlapping or at the mutated regions of the type II calcium binding sites indicating potential functional divergence in certain plant annexins. This study provides a basis for further functional analysis and characterization of annexin multigene families in the plant lineage.     

Characterization of the gene conversions between the multigene family members of the yeast genome

Stanley Sawyer's gene conversion detection method, implemented in his GENECONV computer program, was used to detect and characterize the gene conversions between the multigene family members of the yeast genome. This method gave different gene conversion frequencies and size distribution for gene families with two members and multigene families with more than two members. The 69 gene conversions detected in multigene families with more than two members occur at a frequency of 7.8% gene conversion/pair of genes compared and have an average size of 173+/-220 nucleotides. Larger gene conversions are found only between more similar genes, the genes involved in gene conversions are distributed almost randomly among the 16 yeast chromosomes, and the frequency of gene conversions increases as the distance between repeated genes decreases. In contrast to previous studies, no relationship was observed between the level of expression of a gene and its involvement in gene conversions. These analyses also suggest that gene conversions might occur by different mechanisms in closely linked genes and unlinked genes. The excess of converted regions at the 3? end of unlinked genes suggests that recombination with incomplete cDNA molecules is the main mechanism responsible for gene conversions between such genes.     

Molecular Evolution of Drosophila Cuticular Protein Genes

Several multigene families have been described that together encode scores of structural cuticular proteins in Drosophila, although the functional significance of this diversity remains to be explored. Here I investigate the evolutionary histories of several multigene families (CPR, Tweedle, CPLCG, and CPF/CPFL) that vary in age, size, and sequence complexity, using sequenced Drosophila genomes and mosquito outgroups. My objective is to describe the rates and mechanisms of ‘cuticle-ome’ divergence, in order to identify conserved and rapidly evolving elements. I also investigate potential examples of interlocus gene conversion and concerted evolution within these families during Drosophila evolution. The absolute rate of change in gene number (per million years) is an order of magnitude lower for cuticular protein families within Drosophila than it is among Drosophila and the two mosquito taxa, implying that major transitions in the cuticle proteome have occurred at higher taxonomic levels. Several hotspots of intergenic conversion and/or gene turnover were identified, e.g. some gene pairs have independently undergone intergenic conversion within different lineages. Some gene conversion hotspots were characterized by conversion tracts initiating near nucleotide repeats within coding regions, and similar repeats were found within concertedly evolving cuticular protein genes in Anopheles gambiae. Rates of amino-acid substitution were generally severalfold higher along the branch connecting the Sophophora and Drosophila species groups, and 13 genes have Ka/Ks significantly greater than one along this branch, indicating adaptive divergence. Insect cuticular proteins appear to be a source of adaptive evolution within genera and, at higher taxonomic levels, subject to periods of gene-family expansion and contraction followed by quiescence. However, this relative stasis is belied by hotspots of molecular evolution, particularly concerted evolution, during the diversification of Drosophila. The prominent association between interlocus gene conversion and repeats within the coding sequence of interacting genes suggests that the latter promote strand exchange.     

Role of gene duplication in evolution

It is now known that many multigene and supergene families exist in eukaryote genomes: multigene families with uniform copy members like genes for ribosomal RNA, those with variable members like immunoglobulin genes, and supergene families such as those for various growth factor and hormone receptors. Many such examples indicate that gene duplication and subsequent differentiation are extremely important for organismal evolution. In particular, gene duplication could well have been the primary mechanism for the evolution of complexity in higher organisms. Population genetic models for the origin of gene families with diverse functions are presented, in which natural selection favors those genomes with more useful mutants in duplicated genes. Since any gene has a certain probability of degenerating by mutation, success versus failure in acquiring a new gene by duplication may be expressed as the ratio of probabilities of spreading of useful versus detrimental mutations in redundant gene copies. Also examined are the effects of gene duplication on evolution by compensatory advantageous mutations. Results of the analyses show that both natural selection and random drift are important for the origin of gene families. In addition, interaction between molecular mechanisms such as unequal crossing-over and gene conversion, and selection or drift is found to have a large effect on evolution by gene duplication.     

Multigene organization of structural genes as a phenomenon of eukaryotic organisms

The paper contains a short review of the data on multigene families of protein genes. Structural features and organization of a number of multigene systems coding for proteins are described as models of eukaryotic gene loci arrangements. Based on the extensive information of multigene families, the principles of eukaryotic genome organization are discussed.     

Two chromosomes with multigene families

Watterson's formulae for the distribution, mean, and variance of the number of alleles in common on two chromosomes with multigene families are derived as simpler forms, and extended to chromosomes with an infinite number of genes, each evolving as in an infinitely many alleles model.     

Going green in Cryptococcus neoformans: The recycling of a selectable drug marker

Cryptococcus neoformans is an opportunistic fungal pathogen that primarily affects immunocompromised individuals. Reverse genetics is commonly used to identify and characterize genes involved in a variety of cellular processes. In C. neoformans there is a limited set of positive selectable markers available to make gene deletions or other genetic manipulations. This has hampered the application of reverse genetics in this organism. We have adapted the Bacteriophage P1 Cre-loxP system for use in C. neoformans and successfully excised and reused the same drug marker, G418, to make two sequential gene deletions, lac1Δ and cap59Δ, in the same strain. This tool will allow investigators to make multiple sequential gene deletions in the same strain, which should facilitate the analysis of multigene families.     

Coevolution of immunoglobulin heavy- and light-chain variable-region gene families

The gene families encoding the immunoglobulin variable regions of heavy (VH) and light (VL) chains in vertebrates are composed of many genes. However, the gene number and the extent of diversity among VH and VL gene copies vary with species. To examine the causes of this variation and the evolutionary forces for these multigene families, we conducted a phylogenetic analysis of VH and VL genes from the species of amniotes. The results of our analysis showed that for each species, VH and VL genes have the same pattern of clustering in the trees, and, according to this clustering pattern, the species can be divided into two groups. In the first group of species (humans and mice), VH and VL genes were extensively intermingled with genes from other organisms; in the second group of species (chickens, rabbits, cattle, sheep, swine, and horses), the genes tended to form clusters within the same group of organisms. These results suggest that the VH and VL multigene families have evolved in the same fashion: they have undergone coordinated contraction and expansion of gene repertoires such that each group of organisms is characterized by a certain level of diversity of VH and VL genes. The extent of diversity among copies of VH and VL genes in each species is related to the mechanism of generation of antibody variety. In humans and mice, DNA rearrangement of immunoglobulin variable, diversity, and joining-segment genes is a main source of antibody diversity, whereas in chickens, rabbits, cattle, sheep, swine, and horses, somatic hypermutation and somatic gene conversion play important roles. The evolutionary pattern of VH and VL multigene families is consistent with the birth-and-death model of evolution, yet different levels of diversifying selection seem to operate in the VH and VL genes of these two groups of species.      

Chromosomal Mapping of Repetitive DNAs in the Grasshopper Abracris flavolineata Reveal Possible Ancestry of the B Chromosome and H3 Histone Spreading

Supernumerary chromosomes (B chromosomes) occur in approximately 15% of eukaryote species. Although these chromosomes have been extensively studied, knowledge concerning their specific molecular composition is lacking in most cases. The accumulation of repetitive DNAs is one remarkable characteristic of B chromosomes, and the occurrence of distinct types of multigene families, satellite DNAs and some transposable elements have been reported. Here, we describe the organization of repetitive DNAs in the A complement and B chromosome system in the grasshopper species Abracris flavolineata using classical cytogenetic techniques and FISH analysis using probes for five multigene families, telomeric repeats and repetitive C₀t-1 DNA fractions. The 18S rRNA and H3 histone multigene families are highly variable and well distributed in A. flavolineata chromosomes, which contrasts with the conservation of U snRNA genes and less variable distribution of 5S rDNA sequences. The H3 histone gene was an extensively distributed with clusters occurring in all chromosomes. Repetitive DNAs were concentrated in C-positive regions, including the pericentromeric region and small chromosomal arms, with some occurrence in C-negative regions, but abundance was low in the B chromosome. Finally, the first demonstration of the U2 snRNA gene in B chromosomes in A. flavolineata may shed light on its possible origin. These results provide new information regarding chromosomal variability for repetitive DNAs in grasshoppers and the specific molecular composition of B chromosomes.     

Molecular evolution and population genetics of duplicated accessory gland protein genes in Drosophila

To investigate the potential importance of gene duplication in D. melanogaster accessory gland protein (Acp) gene evolution we carried out a computational analysis comparing annotated D. melanogaster Acp genes to the entire D. melanogaster genome. We found that two known Acp genes are actually members of small multigene families. Polymorphism and divergence data from these duplicated genes suggest that in at least four cases, protein divergence between D. melanogaster and D. simulans is a result of directional selection. One putative Acp revealed by our computational analysis shows evidence of a recent selective sweep in a non-African population (but not in an African population). These data support the idea that selection on reproduction-related genes may drive divergence of populations within species, and strengthen the conclusion that Acps may often be under directional selection in Drosophila.