Numerous Duplications Research Articles

Evolution of fruit development genes in flowering plants.

The genetic mechanisms regulating dry fruit development and opercular dehiscence have been identified in Arabidopsis thaliana. In the bicarpellate silique, valve elongation and differentiation is controlled by FRUITFULL (FUL) that antagonizes SHATTERPROOF1-2 (SHP1/SHP2) and INDEHISCENT (IND) at the dehiscence zone where they control normal lignification. SHP1/2 are also repressed by REPLUMLESS (RPL), responsible for replum formation. Similarly, FUL indirectly controls two other factors ALCATRAZ (ALC) and SPATULA (SPT) that function in the proper formation of the separation layer. FUL and SHP1/2 belong to the MADS-box family, IND and ALC belong to the bHLH family and RPL belongs to the homeodomain family, all of which are large transcription factor families. These families have undergone numerous duplications and losses in plants, likely accompanied by functional changes. Functional analyses of homologous genes suggest that this network is fairly conserved in Brassicaceae and less conserved in other core eudicots. Only the MADS box genes have been functionally characterized in basal eudicots and suggest partial conservation of the functions recorded for Brassicaceae. Here we do a comprehensive search of SHP, IND, ALC, SPT, and RPL homologs across core-eudicots, basal eudicots, monocots and basal angiosperms. Based on gene-tree analyses we hypothesize what parts of the network for fruit development in Brassicaceae, in particular regarding direct and indirect targets of FUL, might be conserved across angiosperms.

Extending the model of Arabidopsis telomere length and composition across Brassicaceae.

Telomeres are repetitive TG-rich DNA elements essential for maintaining the stability of genomes and replicative capacity of cells in almost all eukaryotes. Most of what is known about telomeres in plants comes from the angiosperm Arabidopsis thaliana, which has become an important comparative model for telomere biology. Arabidopsis tolerates numerous insults to its genome, many of which are catastrophic or lethal in other eukaryotic systems such as yeast and vertebrates. Despite the importance of Arabidopsis in establishing a model for the structure and regulation of plant telomeres, only a handful of studies have used this information to assay components of telomeres from across land plants, or even among the closest relatives of Arabidopsis in the plant family Brassicaceae. Here, we determined how well Arabidopsis represents Brassicaceae by comparing multiple aspects of telomere biology in species that represent major clades in the family tree. Specifically, we determined the telomeric repeat sequence, measured bulk telomere length, and analyzed variation in telomere length on syntenic chromosome arms. In addition, we used a phylogenetic approach to infer the evolutionary history of putative telomere-binding proteins, CTC1, STN1, TEN1 (CST), telomere repeat-binding factor like (TRFL), and single Myb histone (SMH). Our analyses revealed conservation of the telomeric DNA repeat sequence, but considerable variation in telomere length among the sampled species, even in comparisons of syntenic chromosome arms. We also found that the single-stranded and double-stranded telomeric DNA-binding complexes CST and TRFL, respectively, differ in their pattern of gene duplication and loss. The TRFL and SMH gene families have undergone numerous duplication events, and these duplicate copies are often retained in the genome. In contrast, CST components occur as single-copy genes in all sampled genomes, even in species that experienced recent whole genome duplication events. Taken together, our results place the Arabidopsis model in the context of other species in Brassicaceae, making the family the best characterized plant group in regard to telomere architecture.

Analysis of yeast SUMO ligase MMS21‐dependent genome maintenance pathway (925.4)

Genome instability often characterizes many cancers and is indicative of high frequency of mutations and genome rearrangements. Several pathways, including the DNA damage checkpoint, play a role in maintaining the stability of the genome by suppressing gross chromosomal rearrangements (GCRs). Additionally, eukaryotic genomes contain numerous repetitive sequences, segmental duplications and regions of copy number variation, which influence the outcomes of gross chromosomal rearrangements.Using a genetic assay in Saccharomyces cerevisiae to screen for suppressors of GCRs, we recently reported that Mms21‐dependent sumoylation specifically suppress GCRs mediated by homeologous sequences. These types of duplication‐mediated GCRs are attributed to Rad52‐dependent non‐allelic homologous recombination (NAHR). We will present findings that show unexpectedly, the DNA damage checkpoint to be required for the formation of duplication‐mediated GCRs in the absence of Mms21‐dependent sumoylation. The canonical DNA damage checkpoint genes MEC1, RAD24, MEC3, RAD9 and RAD53 are required for formation of duplication‐mediated GCRs in the mms21 mutant, while TEL1, MRC1 and CHK1 are minimally involved. Among the checkpoint genes, mutation of RAD9 strongly suppresses GCRs formed in the mms21 mutant and this suppression is independent of either SIZ1 or SIZ2, which encode two other SUMO ligases in yeast. We will also present findings on the role of DNA damage checkpoint in regulating the SUMO proteome. These results show that the outcomes of GCRs are dynamic and are governed by a fine balance of distinct genome maintenance pathways.Grant Funding Source: GM080469, T32CA009523 ,LICR

Genome Degeneration and Adaptation in a Nascent Stage of Symbiosis

Symbiotic associations between animals and microbes are ubiquitous in nature, with an estimated 15% of all insect species harboring intracellular bacterial symbionts. Most bacterial symbionts share many genomic features including small genomes, nucleotide composition bias, high coding density, and a paucity of mobile DNA, consistent with long-term host association. In this study, we focus on the early stages of genome degeneration in a recently derived insect-bacterial mutualistic intracellular association. We present the complete genome sequence and annotation of Sitophilus oryzae primary endosymbiont (SOPE). We also present the finished genome sequence and annotation of strain HS, a close free-living relative of SOPE and other insect symbionts of the Sodalis-allied clade, whose gene inventory is expected to closely resemble the putative ancestor of this group. Structural, functional, and evolutionary analyses indicate that SOPE has undergone extensive adaptation toward an insect-associated lifestyle in a very short time period. The genome of SOPE is large in size when compared with many ancient bacterial symbionts; however, almost half of the protein-coding genes in SOPE are pseudogenes. There is also evidence for relaxed selection on the remaining intact protein-coding genes. Comparative analyses of the whole-genome sequence of strain HS and SOPE highlight numerous genomic rearrangements, duplications, and deletions facilitated by a recent expansion of insertions sequence elements, some of which appear to have catalyzed adaptive changes. Functional metabolic predictions suggest that SOPE has lost the ability to synthesize several essential amino acids and vitamins. Analyses of the bacterial cell envelope and genes encoding secretion systems suggest that these structures and elements have become simplified in the transition to a mutualistic association.

Evolution of the Cation Chloride Cotransporter Family: Ancient Origins, Gene Losses, and Subfunctionalization through Duplication

The cation chloride cotransporter (CCCs) family comprises of four subfamilies-K(+)-Cl(-) cotransporters (KCCs), Na(+)-K(+)-2Cl(-) cotransporters (NKCCs), and Na(+)-Cl(-) cotransporters (NCCs)-and possibly two additional members-CCC interacting protein (CIP1) and polyamine transporters (CCC9)-as well. Altogether, CCCs can play essential physiological roles in transepithelial ion reabsorption and secretion, cell volume regulation, and inhibitory neurotransmission and so are present across all domains of life. To gain insight into the evolution of this family, we performed a comprehensive phylogenetic analysis using publically available genomic information. Our results clearly support CIP1 as being a true CCC based on shared evolutionary history. By contrast, the status of CCC9 in this regard remains equivocal. We also reveal the existence of a single ancestral CCC gene present in Archaea, from which numerous duplication events at the base of archaeans and eukaryotes lead to the divergence and subsequent neofunctionalization of the paralogous CCC subfamilies. A diversity of ensuing gene-loss events resulted in the complex distribution of CCCs present across the different taxa. Importantly, the occurrence of KCCs in "basal" metazoan taxa like sponges would allow an early formation of fast hyperpolarizing neurotransmission in metazoans. Gene duplications within the CCC subfamilies in vertebrates (in particular, KCCs, NKCCs, and NCCs) lend further evidence to the 2R hypothesis of two rounds of genome duplication at the base of the vertebrate lineage, especially in concert with our syntenic cluster analyses. This increased number of KCCs, NKCCs, and NCCs isoforms facilitates their further, important subfunctionalization in the vertebrate lineage.

RNaseIII and T4 Polynucleotide Kinase sequence biases and solutions during RNA-seq library construction

BackgroundRNA-seq is a next generation sequencing method with a wide range of applications including single nucleotide polymorphism (SNP) detection, splice junction identification, and gene expression level measurement. However, the RNA-seq sequence data can be biased during library constructions resulting in incorrect data for SNP, splice junction, and gene expression studies. Here, we developed new library preparation methods to limit such biases.ResultsA whole transcriptome library prepared for the SOLiD system displayed numerous read duplications (pile-ups) and gaps in known exons. The pile-ups and gaps of the whole transcriptome library caused a loss of SNP and splice junction information and reduced the quality of gene expression results. Further, we found clear sequence biases for both 5' and 3' end reads in the whole transcriptome library. To remove this bias, RNaseIII fragmentation was replaced with heat fragmentation. For adaptor ligation, T4 Polynucleotide Kinase (T4PNK) was used following heat fragmentation. However, its kinase and phosphatase activities introduced additional sequence biases. To minimize them, we used OptiKinase before T4PNK. Our study further revealed the specific target sequences of RNaseIII and T4PNK.ConclusionsOur results suggest that the heat fragmentation removed the RNaseIII sequence bias and significantly reduced the pile-ups and gaps. OptiKinase minimized the T4PNK sequence biases and removed most of the remaining pile-ups and gaps, thus maximizing the quality of RNA-seq data.ReviewersThis article was reviewed by Dr. A. Kolodziejczyk (nominated by Dr. Sarah Teichmann), Dr. Eugene Koonin, and Dr. Christoph Adami. For the full reviews, see the Reviewers' Comments section.

The numerical duplication of a numerical semigroup

In this paper we present and study the numerical duplication of a numerical semigroup, a construction that, starting with a numerical semigroup S and a semigroup ideal E⊆S, produces a new numerical semigroup, denoted by S⋈bE (where b is any odd integer belonging to S), such that S=(S⋈bE)/2. In particular, we characterize the ideals E such that S⋈bE is almost symmetric and we determine its type.

Evolutionary Analysis of the MIXTA Gene Family Highlights Potential Targets for the Study of Cellular Differentiation

Differentiated epidermal cells such as trichomes and conical cells perform numerous essential functions in plant biology and are important for our understanding of developmental patterning and cell shape regulation. Many are also commercially significant, such as cotton fibers and trichomes that secrete pharmaceutically useful or herbivore-deterring compounds. Here, we focus on the phylogeny and evolution of the subgroup 9 R2R3 MYB gene transcription factors, which include the MIXTA gene, and that are important for the specification and regulation of plant cellular differentiation. We have sequenced 49 subgroup 9 R2R3 MYB genes from key experimental taxa and combined these sequences with those identified by an exhaustive bioinformatic search, to compile a data set of 223 subgroup 9 R2R3 MYB genes. Our phylogenetic analyses demonstrate, for the first time, the complex evolutionary history of the subgroup 9 R2R3 MYB genes. A duplication event is inferred before the origin of seed plants giving rise to two major gene lineages, here termed SBG9-A and SBG9-B. The evolutionary conservation of the SBG9-B gene lineage has not been previously recognized and its role in cellular differentiation is unknown, thus an entire clade of potential candidate genes for epidermal cell regulation remains to be explored. Using a heterologous transformation bioassay, we provide functional data that implicate members of the SBG9-B lineage in the specification of epidermal projections. Furthermore, we reveal numerous putative duplication events in both SBG9-A and SBG9-B lineages, resolving uncertainty about orthology and paralogy among the subgroup 9 R2R3 MYB genes. Finally, we provide a robust framework over which to interpret existing functional data and to direct ongoing comparative genetic research into the evolution of plant cellular diversity.

Early events in the evolution of spider silk genes.

Silk spinning is essential to spider ecology and has had a key role in the expansive diversification of spiders. Silk is composed primarily of proteins called spidroins, which are encoded by a multi-gene family. Spidroins have been studied extensively in the derived clade, Orbiculariae (orb-weavers), from the suborder Araneomorphae (‘true spiders’). Orbicularians produce a suite of different silks, and underlying this repertoire is a history of duplication and spidroin gene divergence. A second class of silk proteins, Egg Case Proteins (ECPs), is known only from the orbicularian species, Lactrodectus hesperus (Western black widow). In L. hesperus, ECPs bond with tubuliform spidroins to form egg case silk fibers. Because most of the phylogenetic diversity of spiders has not been sampled for their silk genes, there is limited understanding of spidroin gene family history and the prevalence of ECPs. Silk genes have not been reported from the suborder Mesothelae (segmented spiders), which diverged from all other spiders >380 million years ago, and sampling from Mygalomorphae (tarantulas, trapdoor spiders) and basal araneomorph lineages is sparse. In comparison to orbicularians, mesotheles and mygalomorphs have a simpler silk biology and thus are hypothesized to have less diversity of silk genes. Here, we present cDNAs synthesized from the silk glands of six mygalomorph species, a mesothele, and a non-orbicularian araneomorph, and uncover a surprisingly rich silk gene diversity. In particular, we find ECP homologs in the mesothele, suggesting that ECPs were present in the common ancestor of extant spiders, and originally were not specialized to complex with tubuliform spidroins. Furthermore, gene-tree/species-tree reconciliation analysis reveals that numerous spidroin gene duplications occurred after the split between Mesothelae and Opisthothelae (Mygalomorphae plus Araneomorphae). We use the spidroin gene tree to reconstruct the evolution of amino acid compositions of spidroins that perform different ecological functions.

Complex evolutionary history and diverse domain organization of SET proteins suggest divergent regulatory interactions

• Plants and animals possess very different developmental processes, yet share conserved epigenetic regulatory mechanisms, such as histone modifications. One of the most important forms of histone modification is methylation on lysine residues of the tails, carried out by members of the SET protein family, which are widespread in eukaryotes. • We analyzed molecular evolution by comparative genomics and phylogenetics of the SET genes from plant and animal genomes, grouping SET genes into several subfamilies and uncovering numerous gene duplications, particularly in the Suv, Ash, Trx and E(z) subfamilies. • Domain organizations differ between different subfamilies and between plant and animal SET proteins in some subfamilies, and support the grouping of SET genes into seven main subfamilies, suggesting that SET proteins have acquired distinctive regulatory interactions during evolution. We detected evidence for independent evolution of domain organization in different lineages, including recruitment of new domains following some duplications. • More recent duplications in both vertebrates and land plants are probably the result of whole-genome or segmental duplications. The evolution of the SET gene family shows that gene duplications caused by segmental duplications and other mechanisms have probably contributed to the complexity of epigenetic regulation, providing insights into the evolution of the regulation of chromatin structure.

Additional genomic duplications in AZFc underlie the b2/b3 deletion-associated risk of spermatogenic impairment in Han Chinese population

The azoospermia factor c (AZFc) region on the Y chromosome is a genetically dynamic locus in the human genome. Numerous genomic rearrangements, including deletion, duplication and inversion, have been identified in AZFc. The complete deletion of AZFc can cause spermatogenic impairment. However, the roles of partial AZFc deletions (e.g. b2/b3 deletion) in spermatogenesis are controversial and variable among human populations. Secondary duplication has been hypothesized to be a compensatory factor for partial AZFc deletions. To further study genomic duplications in AZFc as a potential genetic modifier underlying the phenotypic variations of partial AZFc deletions in spermatogenesis, we conducted comprehensive molecular analyses in 711 idiopathic infertile men and 390 healthy controls. Unexpectedly, we found that additional AZFc duplications accompanying the b2/b3 deletion, instead of the b2/b3 deletion alone, led to the b2/b3 deletion-associated risk of spermatogenic impairment previously reported in Han Chinese population. In addition, partial AZFc duplication also rendered a risk factor in the non-deletion patients. DAZ is a multi-copy AZFc gene (DAZ1-DAZ4) implicated in spermatogenesis. Genetic variations do exist between DAZ copies. Intriguingly, we found that the DAZ1/2 cluster was the main duplicated copies in the partial AZFc duplications associated with spermatogenic impairment, suggesting a potential different role of spermatogenesis between DAZ copies. Our findings demonstrated that additional AZFc duplications did not compensate but convey the susceptibility of the b2/b3 deletion to spermatogenic impairment in the tested population. Notably, genomic duplications and deletions in AZFc deserve comprehensive investigations to uncover spermatogenic roles of the AZFc region.

A New Highly Conserved Antibiotic Sensing/Resistance Pathway in Firmicutes Involves an ABC Transporter Interplaying with a Signal Transduction System

Signal transduction systems and ABC transporters often contribute jointly to adaptive bacterial responses to environmental changes. In Bacillus subtilis, three such pairs are involved in responses to antibiotics: BceRSAB, YvcPQRS and YxdJKLM. They are characterized by a histidine kinase belonging to the intramembrane sensing kinase family and by a translocator possessing an unusually large extracytoplasmic loop. It was established here using a phylogenomic approach that systems of this kind are specific but widespread in Firmicutes, where they originated. The present phylogenetic analyses brought to light a highly dynamic evolutionary history involving numerous horizontal gene transfers, duplications and lost events, leading to a great variety of Bce-like repertories in members of this bacterial phylum. Based on these phylogenetic analyses, it was proposed to subdivide the Bce-like modules into six well-defined subfamilies. Functional studies were performed on members of subfamily IV comprising BceRSAB from B. subtilis, the expression of which was found to require the signal transduction system as well as the ABC transporter itself. The present results suggest, for the members of this subfamily, the occurrence of interactions between one component of each partner, the kinase and the corresponding translocator. At functional and/or structural levels, bacitracin dependent expression of bceAB and bacitracin resistance processes require the presence of the BceB translocator loop. Some other members of subfamily IV were also found to participate in bacitracin resistance processes. Taken together our study suggests that this regulatory mechanism might constitute an important common antibiotic resistance mechanism in Firmicutes. [Supplemental material is available online at http://www.genome.org.]

Genome duplication events have led to a diversification in the CPT I gene family in fish

The enzyme carnitine palmitoyltransferase (CPT) I is a major regulator of mitochondrial fatty acid oxidation in vertebrates. Numerous genome duplication events throughout evolution have given rise to three (in mammals) or multiple (in fish) genetically and functionally different isoforms of this enzyme. In particular, these isoforms represent a diversification of kinetic and regulatory properties stemming from mutations at the genomic and proteomic levels. Phylogenetic reconstructions reveal a comprehensive view of the CPT I family in vertebrates and genomic modifications leading to structural changes in proteins and functional differences between tissues and taxa. In a model fish species (rainbow trout), the presence of five CPT I isoforms suggests repeated duplication events in bony fishes and salmonids. Subsequently, an array of nucleotide and amino acid substitutions in the isoforms may contribute to a tissue-specific and a previously observed species-specific difference in the IC(50) for malonyl-CoA. Moreover, all five isoforms are expressed in trout at the mRNA level in skeletal muscle, heart, liver, kidney, and intestine. In general, transcript levels of the beta-isoforms were higher in muscle tissues, while levels of the alpha-isoforms were higher in other tissues. Rainbow trout also exhibit developmental plasticity in relative mRNA expression of CPT I isoforms from fry to juvenile to adult stage. Thus the evolution of CPT I has resulted in a very diverse family of isoforms. These differences represent a degree of specificity in the ability of species to regulate function at the protein and tissue levels, which, in turn, may allow for precise control of lipid oxidation in individual tissues during physiological perturbations.

Phylogenetic detection of numerous gene duplications shared by animals, fungi and plants

A phylogentic analysis reveals that many gene duplications occurred prior to the split in animals, fungi and plants.

Genome-wide analysis of cyclins in maize (Zea mays)

Cyclins are primary regulators of the activity of cyclin-dependent kinases and play crucial roles in cell cycle progression in eukaryotes. Although extensive studies have revealed the roles of some cyclins and underlying mechanisms in plants, relatively few cyclins have been functionally analyzed in maize. We identified 59 cyclins in the maize genome, distributed on 10 chromosomes; these were grouped into six types by phylogenetic analysis. The cyclin genes in the maize genome went through numerous tandem gene duplications on five chromosomes. However, no segmental duplications, which occur in rice, were found on maize chromosomes. This information allows us to assess the position of plant cyclin genes in terms of evolution and classification, which will be useful for functional studies of maize cyclins.

The CRF domain defines Cytokinin Response Factor proteins in plants

BackgroundCytokinin Response Factors (CRFs) are a small subset of AP2/ERF transcription factor genes shown in Arabidopsis to regulate leaf development as part of the cytokinin signal transduction pathway. This study examines the phylogenetic distribution of CRF genes in other plant species, and attempts to identify the extent of sequence conservation and potential gene function among all CRF genes.ResultsWe identified CRF genes in representatives of all major land plant lineages, including numerous flowering plant taxa in addition to the model systems in which ERF genes have been catalogued. Comparative analysis across this broader sampling has identified strongly conserved amino acid motifs other than the AP2/ERF domain for all CRF proteins as well as signature sequences unique to specific clades of CRF genes. One of these motifs, here designated as the CRF domain, is conserved in and unique to CRF proteins distinguishing them from related genes. We show that this novel domain of approximately 65 amino acids is found in CRF proteins from all groups of land plants and only in CRF genes. Phylogenetic analyses suggest that the evolution of CRF genes has included numerous duplication events. In this phylogenetic context we examine protein evolution including the gain and loss of accessory domains, correlate these molecular evolutionary events with experimental data on cytokinin regulation and speculate on the function and evolution of the CRF domain within AP2/ERF transcription factor proteins. We also tested a prediction drawn from the phylogenetic analyses that four CRF domain containing genes from Tomato, previously unexamined for cytokinin response, are transcriptionally inducible by cytokinin, supporting the link between CRF genes, CRF-specific domains and cytokinin regulation.ConclusionCRF genes can be identified in all lineages of land plants, as a distinct subset of AP2/ERF proteins containing a specific and unique CRF domain. The CRF domain can be used to identify previously unclassified predicted genes or genes identified only as members of the AP2/ERF protein family. CRF domain presence and phylogenetic relatedness to known Arabidopsis CRF genes predicts gene function to some extent.

Phylogenomic analysis of the cystatin superfamily in eukaryotes and prokaryotes

BackgroundThe cystatin superfamily comprises cysteine protease inhibitors that play key regulatory roles in protein degradation processes. Although they have been the subject of many studies, little is known about their genesis, evolution and functional diversification. Our aim has been to obtain a comprehensive insight into their origin, distribution, diversity, evolution and classification in Eukaryota, Bacteria and Archaea.ResultsWe have identified in silico the full complement of the cystatin superfamily in more than 2100 prokaryotic and eukaryotic genomes. The analysis of numerous eukaryotic genomes has provided strong evidence for the emergence of this superfamily in the ancestor of eukaryotes. The progenitor of this superfamily was most probably intracellular and lacked a signal peptide and disulfide bridges, much like the extant Giardia cystatin. A primordial gene duplication produced two ancestral eukaryotic lineages, cystatins and stefins. While stefins remain encoded by a single or a small number of genes throughout the eukaryotes, the cystatins have undergone a more complex and dynamic evolution through numerous gene and domain duplications. In the cystatin superfamily we discovered twenty vertebrate-specific and three angiosperm-specific orthologous families, indicating that functional diversification has occurred only in multicellular eukaryotes. In vertebrate orthologous families, the prevailing trends were loss of the ancestral inhibitory activity and acquisition of novel functions in innate immunity. Bacterial cystatins and stefins may be emergency inhibitors that enable survival of bacteria in the host, defending them from the host's proteolytic activity.ConclusionThis study challenges the current view on the classification, origin and evolution of the cystatin superfamily and provides valuable insights into their functional diversification. The findings of this comprehensive study provide guides for future structural and evolutionary studies of the cystatin superfamily as well as of other protease inhibitors and proteases.

Unique genes in plants: specificities and conserved features throughout evolution

BackgroundPlant genomes contain a high proportion of duplicated genes as a result of numerous whole, segmental and local duplications. These duplications lead up to the formation of gene families, which are the usual material for many evolutionary studies. However, all characterized genomes include single-copy (unique) genes that have not received much attention. Unlike gene duplication, gene loss is not an unspecific mechanism but is rather influenced by a functional selection. In this context, we have established and used stringent criteria in order to identify suitable sets of unique genes present in plant proteomes. Comparisons of unique genes in the green phylum were used to characterize the gene and protein features exhibited by both conserved and species-specific unique genes.ResultsWe identified the unique genes within both A. thaliana and O. sativa genomes and classified them according to the number of homologs in the alternative species: none (U{1:0}), one (U{1:1}) or several (U{1:m}). Regardless of the species, all the genes in these groups present some conserved characteristics, such as small average protein size and abnormal intron number. In order to understand the origin and function of unique genes, we further characterized the U{1:1} gene pairs. The possible involvement of sequence convergence in the creation of U{1:1} pairs was discarded due to the frequent conservation of intron positions. Furthermore, an orthology relationship between the two members of each U{1:1} pair was strongly supported by a high conservation in the protein sizes and transcription levels. Within the promoter of the unique conserved genes, we found a number of TATA and TELO boxes that specifically differed from their mean number in the whole genome. Many unique genes have been conserved as unique through evolution from the green alga Ostreococcus lucimarinus to higher plants. Plant unique genes may also have homologs in bacteria and we showed a link between the targeting towards plastids of proteins encoded by plant nuclear unique genes and their homology with a bacterial protein.ConclusionMany of the A. thaliana and O. sativa unique genes are conserved in plants for which the ancestor diverged at least 725 million years ago (MYA). Half of these genes are also present in other eukaryotic and/or prokaryotic species. Thus, our results indicate that (i) a strong negative selection pressure has conserved a number of genes as unique in genomes throughout evolution, (ii) most unique genes are subjected to a low divergence rate, (iii) they have some features observed in housekeeping genes but for most of them there is no functional annotation and (iv) they may have an ancient origin involving a possible gene transfer from ancestral chloroplasts or bacteria to the plant nucleus.

Analysis of segmental duplications via duplication distance

Segmental duplications are common in mammalian genomes, but their evolutionary origins remain mysterious. A major difficulty in analyzing segmental duplications is that many duplications are complex mosaics of fragments of numerous other segmental duplications. We introduce a novel measure called duplication distance that describes the minimum number of duplications necessary to create a target string by repeated insertions of fragments of a source string. We derive an efficient algorithm to compute duplication distance, and we use the algorithm to analyze segmental duplications in the human genome. Our analysis reveals possible ancestral relationships between segmental duplications including numerous examples of duplications that contain multiple, nested insertions of fragments from one or more other duplications. Using duplication distance, we also identify a small number of segmental duplications that appear to have seeded many other duplications in the genome, lending support to a two-step model of segmental duplication in the genome. Software for computing duplication distance is available upon request.

Extensive Copy-Number Variation of the Human Olfactory Receptor Gene Family

As much as a quarter of the human genome has been reported to vary in copy number between individuals, including regions containing about half of the members of the olfactory receptor (OR) gene family. We have undertaken a detailed study of copy-number variation of ORs to elucidate the selective and mechanistic forces acting on this gene family and the true impact of copy-number variation on human OR repertoires. We argue that the properties of copy-number variants (CNVs) and other sets of large genomic regions violate the assumptions of statistical methods that are commonly used in the assessment of gene enrichment. Using more appropriate methods, we provide evidence that OR enrichment in CNVs is not due to positive selection but is because of OR preponderance in segmentally duplicated regions, which are known to be frequently copy-number variable, and because purifying selection against CNVs is lower in OR-containing regions than in regions containing essential genes. We also combine multiplex ligation-dependent probe amplification (MLPA) and PCR to assay the copy numbers of 37 candidate CNV ORs in a panel of approximately 50 human individuals. We confirm copy-number variation of 18 ORs but find no variation in this human-diversity panel for 16 other ORs, highlighting the caveat that reported intervals often overrepresent true CNVs. The copy-number variation we describe is likely to underpin significant variation in olfactory abilities among human individuals. Finally, we show that both homology-based and homology-independent processes have played a recent role in remodeling the OR family.