Analysing gene function after duplication

After gene duplication, mutations cause the gene copies to diverge. The classical model predicts that these mutations will generally lead to the loss of function of one gene copy; rarely, new functions will be created and both duplicate genes are conserved. In contrast, under the subfunctionalization model both duplicates are preserved due to the partition of different functions between the duplicates. A recent study provides support for the subfunctionalization model, identifying several expressed gene duplicates common to humans and mice that contain regions conserved in one duplicate but variable in the other (and vice versa). We discuss both the methodology used in this study and also how gene phylogeny may lead to additional evidence for the importance of subfunctionalization in the evolution of new genes.     

Population-genetic models of the fates of duplicate genes

The ultimate fate of a duplicated gene is that it either silenced through inactivating mutations or both copies are maintained by selection. This later fate can occur via neofunctionalization wherein one copy acquires a new function or by subfunctionalization wherein the original function of the gene is partitioned across both copies. The relative probabilities of these three different fates involve often very subtle iterations between of population size, mutation rate, and selection. All three of these fates are critical to the expansion and diversification of gene families.     

Subfunctionalization of duplicated genes as a transition state to neofunctionalization

Background  

Gene duplication has been suggested to be an important process in the generation of evolutionary novelty. Neofunctionalization, as an adaptive process where one copy mutates into a function that was not present in the pre-duplication gene, is one mechanism that can lead to the retention of both copies. More recently, subfunctionalization, as a neutral process where the two copies partition the ancestral function, has been proposed as an alternative mechanism driving duplicate gene retention in organisms with small effective population sizes. The relative importance of these two processes is unclear.     

Coat colours in the Massese sheep breed are associated with mutations in the agouti signalling protein (ASIP) and melanocortin 1 receptor (MC1R) genes

Massese is an Italian dairy sheep breed characterized by animals with black skin and horns and black or apparent grey hairs. Owing to the presence of these two coat colour types, this breed can be considered an interesting model to evaluate the effects of coat colour gene polymorphisms on this phenotypic trait. Two main loci have been already shown to affect coat colour in sheep: Agouti and Extension coding for the agouti signalling protein (ASIP) and melanocortin 1 receptor (MC1R) genes, respectively. The Agouti locus is affected by a large duplication including the ASIP gene that may determine the Agouti white and tan allele (A(Wt)). Other disrupting or partially inactivating mutations have been identified in exon 2 (a deletion of 5 bp, D(5); and a deletion of 9 bp, D(9)) and in exon 4 (g.5172T>A, p.C126S) of the ASIP gene. Three missense mutations in the sheep MC1R gene cause the dominant black E(D) allele (p.M73K and p.D121N) and the putative recessive e allele (p.R67C). Here, we analysed these ASIP and MC1R mutations in 161 Massese sheep collected from four flocks. The presence of one duplicated copy allele including the ASIP gene was associated with grey coat colour (P = 9.4E-30). Almost all animals with a duplicated copy allele (37 out of 41) showed uniform apparent grey hair and almost all animals without a duplicated allele (117 out of 120) were completely black. Different forms of duplicated alleles were identified in Massese sheep including, in almost all cases, copies with exon 2 disrupting or partially inactivating mutations making these alleles different from the A(Wt) allele. A few exceptions were observed in the association between ASIP polymorphisms and coat colour: three grey sheep did not carry any duplicated copy allele and four black animals carried a duplicated copy allele. Of the latter four sheep, two carried the E(D) allele of the MC1R gene that may be the cause of their black coat colour. The coat colour of all other black animals may be determined by non-functional ASIP alleles (non-agouti alleles, A(a)) and in a few cases by the E(D) Extension allele. At least three frequent ASIP haplotypes ([D(5):g.5172T], [N:g.5172A] and [D(5):g.5172A]) were detected (organized into six different diplotypes). In conclusion, the results indicated that coat colours in the Massese sheep breed are mainly derived by combining ASIP and MC1R mutations.     

The evolution of the novel Sdic gene cluster in Drosophila melanogaster

The origin of new genes and of new functions for existing genes are fundamental processes in molecular evolution. Sdic is a newly evolved gene that arose recently in the D. melanogaster lineage. The gene encodes a novel sperm motility protein. It is a chimeric gene formed by duplication of two other genes followed by multiple deletions and other sequence rearrangements. The Sdic gene exists in several copies in the X chromosome, and is presumed to have undergone several duplications to form a tandemly arrayed gene cluster. Given the very recent origin of the gene and the gene cluster, the analysis of the composition of this gene cluster represents an excellent opportunity to study the origin and evolution of new gene functions and the fate of gene duplications. We have analyzed the nucleotide sequence of this region and reconstructed the evolutionary history of this gene cluster. We found that the cluster is composed by four tandem copies of Sdic; these duplicates are very similar but can be distinguished by the unique pattern of insertions, deletions, and point mutations in each copy. The oldest gene copy in the array has a 3' exon that has undergone accelerated diversification, and also shows divergent regulatory sequences. Moreover, there is evidence that this might be the only gene copy in the tandem array that is transcribed at a significant level, expressing a novel sperm-specific protein. There is also a retrotransposon located at the 3' end of each Sdic gene copy. We argue that this gene cluster was formed in the last two million years by at least three tandem duplications and one retrotransposition event.     

The probability of preservation of a newly arisen gene duplicate.

Newly emerging data from genome sequencing projects suggest that gene duplication, often accompanied by genetic map changes, is a common and ongoing feature of all genomes. This raises the possibility that differential expansion/contraction of various genomic sequences may be just as important a mechanism of phenotypic evolution as changes at the nucleotide level. However, the population-genetic mechanisms responsible for the success vs. failure of newly arisen gene duplicates are poorly understood. We examine the influence of various aspects of gene structure, mutation rates, degree of linkage, and population size (N) on the joint fate of a newly arisen duplicate gene and its ancestral locus. Unless there is active selection against duplicate genes, the probability of permanent establishment of such genes is usually no less than 1/(4N) (half of the neutral expectation), and it can be orders of magnitude greater if neofunctionalizing mutations are common. The probability of a map change (reassignment of a key function of an ancestral locus to a new chromosomal location) induced by a newly arisen duplicate is also generally >1/(4N) for unlinked duplicates, suggesting that recurrent gene duplication and alternative silencing may be a common mechanism for generating microchromosomal rearrangements responsible for postreproductive isolating barriers among species. Relative to subfunctionalization, neofunctionalization is expected to become a progressively more important mechanism of duplicate-gene preservation in populations with increasing size. However, even in large populations, the probability of neofunctionalization scales only with the square of the selective advantage. Tight linkage also influences the probability of duplicate-gene preservation, increasing the probability of subfunctionalization but decreasing the probability of neofunctionalization.     

Molecular evolution of the ligands of the insulin-signaling pathway: dilp genes in the genus Drosophila

Drosophila melanogaster, unlike mammals, has seven insulin-like peptides (DILPS). In Drosophila, all seven genes (dilp1-7) are single copy in the 12 species studied, except for D. grimshawi with two tandem copies of dilp2. Our comparative analysis revealed that genes dilp1-dilp7 exhibit differential functional constraint, which is indicative of some functional divergence. Species of the subgenera Sophophora and Drosophila differ in some traits likely affected by the insulin-signaling pathway, such as adult body size. It is in the branch connecting the two subgenera that we found the footprint left by positive selection driving nonsynonymous changes at some dilp1 codons to fixation. Finally, the similar rate at which the two dilp2 copies of D. grimshawi have evolved since their duplication and the presence of a putative regulatory region highly conserved between the two paralogs would suggest that both copies were preserved either because of subfunctionalization or dose dependency rather than by the neofunctionalization of one of the two copies.     

Retention of a Duplicate Gene Through Changes in Subcellular Targeting: An Electron Transport Protein Homologue Localizes to the Golgi

Cytochrome c oxidase (COX), the terminal enzyme complex of the electron transport chain, contains 13 subunits, 3 encoded by mitochondrial DNA and 10 by nuclear. Several of the nuclear subunits, including subunit VIIa, are known to have two tissue- and development-specific isoforms in mammals. A recently identified third member of the gene family, COX7AR, encodes a protein previously thought to function in mitochondria. However, observation of fluorescent pCOX7AR C-terminal fusion proteins in HeLa cells showed that pCOX7AR is localized to the Golgi apparatus. Sequence analyses indicate that the duplication of COX7AR occurred prior to the origin of the Euteleostomi (bony vertebrates) and that pCOX7AR is more highly conserved than the two other isoforms. These results indicate that, after gene duplication and modification of the mitochondrial targeting signal, pCOX7AR was evolutionarily altered to a new and apparently important function in the Golgi. These results also suggest that predictions of function from homology can be misleading and show that specialization and modification of subcellular localization are similar to cis-element subfunctionalization. In cis-element subfunctionalization, complementary null mutations occur to the cis-elements of the descendents of a gene duplication, causing both descendent genes to be obligate. In the process described in this paper, which could be termed subcellular subfunctionalization, complementary null mutations can occur to the subcellular localization signals of the descendants of a gene duplication, causing both descendent genes to be similarly obligate. Noncomplementary null mutations could also uncover an alternate localization, which is the more likely case for pCOX7AR.     

The limits of subfunctionalization

Background  

The duplication-degeneration-complementation (DDC) model has been proposed as an explanation for the unexpectedly high retention of duplicate genes. The hypothesis proposes that, following gene duplication, the two gene copies degenerate to perform complementary functions that jointly match that of the single ancestral gene, a process also known as subfunctionalization. We distinguish between subfunctionalization at the regulatory level and at the product level (e.g within temporal or spatial expression domains).     

Dual duplication of neurohypophysial hormones in an Australian marsupial: mesotocin, oxytocin, lysine vasopressin and arginine vasopressin in a single gland of the northern bandicoot (Isoodon macrourus)

Neurohypophysial hormones of an Australian marsupial, the Northern bandicoot (Isoodon macrourus), have been identified by their retention times in high-pressure reverse-phase liquid chromatography using two solvent systems and by their molar pressor or uterotonic activities. Two pressor peptides, arginine vasopressin and lysipressin, and two uterotonic peptides, mesotocin and oxytocin, have been characterized. Because mesotocin and arginine vasopressin have been identified in three other Australian marsupial families, it is assumed that a duplication of each ancestral gene occurred in Peramelidae and subsequent mutations in one copy led to the additional oxytocin and lysipressin. A similar dual duplication of neurohypophysial hormones has previously been discovered in the North-American opossum (Didelphis virginiana) so that the duplication propensity seems peculiar to marsupials in contrast to placental mammals.     

Rapid evolution through gene duplication and subfunctionalization of the testes-specific alpha4 proteasome subunits in Drosophila

Gene duplication is an important mechanism for acquiring new genes and creating genetic novelty in organisms. Evidence suggests that duplicated genes are retained at a much higher rate than originally thought and that functional divergence of gene copies is a major factor promoting their retention in the genome. We find that two Drosophila testes-specific alpha4 proteasome subunit genes (alpha4-t1 and alpha4-t2) have a higher polymorphism within species and are significantly more diverged between species than the somatic alpha4 gene. Our data suggest that following gene duplication, the alpha4-t1 gene experienced relaxed selective constraints, whereas the alpha4-t2 gene experienced positive selection acting on several codons. We report significant heterogeneity in evolutionary rates among all three paralogs at homologous codons, indicating that functional divergence has coincided with genic divergence. Reproductive subfunctionalization may allow for a more rapid evolution of reproductive traits and a greater specialization of testes function. Our data add to the increasing evidence that duplicated genes experience lower selective constraints and in some cases positive selection following duplication. Newly duplicated genes that are freer from selective constraints may provide a mechanism for developing new interactions and a pathway for the evolution of new genes.     

Questioning the Ubiquity of Neofunctionalization

Gene duplication provides much of the raw material from which functional diversity evolves. Two evolutionary mechanisms have been proposed that generate functional diversity: neofunctionalization, the de novo acquisition of function by one duplicate, and subfunctionalization, the partitioning of ancestral functions between gene duplicates. With protein interactions as a surrogate for protein functions, evidence of prodigious neofunctionalization and subfunctionalization has been identified in analyses of empirical protein interactions and evolutionary models of protein interactions. However, we have identified three phenomena that have contributed to neofunctionalization being erroneously identified as a significant factor in protein interaction network evolution. First, self-interacting proteins are underreported in interaction data due to biological artifacts and design limitations in the two most common high-throughput protein interaction assays. Second, evolutionary inferences have been drawn from paralog analysis without consideration for concurrent and subsequent duplication events. Third, the theoretical model of prodigious neofunctionalization is unable to reproduce empirical network clustering and relies on untenable parameter requirements. In light of these findings, we believe that protein interaction evolution is more persuasively characterized by subfunctionalization and self-interactions.     

Copy number of DAZ genes in Slovenian and Bosnian general population

Deletions of two of four DAZ (Deleted in AZoospermia) gene copies located on the Y chromosome were associated with spermatogenic failure, but the information on DAZ copy number is still very scarce. The aim of this study was to determine the frequency of partial DAZ gene deletions and to analyze the existence of duplications in general Slovenian and Bosnian population. To answer these questions, we used real time PCR. We analyzed 100 male samples from Slovenian and Bosnian general population. The incidence of two DAZ gene copies was 6% (3/50) in Slovenian population. The incidence of more than four DAZ genes was 2% (1/50) in Slovenian population and 8% (4/50) in Bosnian population. Observed differences have not reached statistical significance. In conclusion we demonstrate that DAZ genes are not only prone to deletions but also to duplication events. Further studies are needed to estimate the prevalence of these mutations and its' relevance to male infertility.     

Polyploidy and DNA methylation: new tools available

Most plant species are recent or ancient polyploids (displaying at least one round of genome duplication in their history). Cultivated species (e.g. wheat, cotton, canola, sugarcane, coffee) and invasive species are often relatively recent polyploids, and frequently of hybrid origin (i.e. allopolyploids). Despite the genetic bottleneck occurring during the allopolyploid speciation process, the formation of such species from two divergent lineages leads to fixed heterozygosity decisive to their success. New phenotypes and new niche occupation are usually associated with this mode of speciation, as a result of both genomic rearrangements and gene expression changes of different magnitudes depending on the different polyploid species investigated. These gene expression changes affecting newly formed polyploid species may result from various, interconnected mechanisms, including (i) functional interactions between the homoeologous copies and between their products, that are reunited in the same nucleus and cell; (ii) the fate of duplicated copies, selective pressure on one of the parental copy being released which could lead to gene loss, pseudogenization, or alternatively, to subfunctionalization or neofunctionalization; and (iii) epigenetic landscape changes that in turn affect gene expression. As one of the interrelated processes leading to epigenetic regulation of gene expression, the DNA methylation status of newly formed species appears to be consistently affected following both hybridization and genome doubling. In this issue, Verhoeven et al. have investigated the fate of DNA methylation patterns that could affect naturally occurring new asexual triploid lineages of dandelions. As a result of such a ploidy level change, the authors demonstrate stably transmitted DNA methylation changes leading to unique DNA methylation patterns in each newly formed lineage. Most studies published to date on plant DNA methylation polymorphism were performed using restriction enzymes sensitive to methylation. Recently, new high‐throughput methods were made available, thanks to the development of ‘next‐generation sequencing’ techniques. The combination of these methods offers powerful and promising tools to investigate epigenetic variation in both model and non‐model systems.     

Two Sequence-Ready Contigs Spanning the Two Copies of a 200-kb Duplication on Human 21q: Partial Sequence and Polymorphisms

Physical mapping across a duplication can be a tour de force if the region is larger than the size of a bacterial clone. This was the case of the 170- to 275-kb duplication present on the long arm of chromosome 21 in normal human at 21q11.1 (proximal region) and at 21q22.1 (distal region), which we described previously. We have constructed sequence-ready contigs of the two copies of the duplication of which all the clones are genuine representatives of one copy or the other. This required the identification of four duplicon polymorphisms that are copy-specific and nonallelic variations in the sequence of the STSs. Thirteen STSs were mapped inside the duplicated region and 5 outside but close to the boundaries. Among these STSs 10 were end clones from YACs, PACs, or cosmids, and the average interval between two markers in the duplicated region was 16 kb. Eight PACs and cosmids showing minimal overlaps were selected in both copies of the duplication. Comparative sequence analysis along the duplication showed three single-basepair changes between the two copies over 659 bp sequenced (4 STSs), suggesting that the duplication is recent (less than 4 mya). Two CpG islands were located in the duplication, but no genes were identified after a 36-kb cosmid from the proximal copy of the duplication was sequenced. The homology of this chromosome 21 duplicated region with the pericentromeric regions of chromosomes 13, 2, and 18 suggests that the mechanism involved is probably similar to pericentromeric-directed mechanisms described in interchromosomal duplications.