Genome Segmental Duplication Research Articles

Local Genomic Instability of the SpTransformer Gene Family in the Purple Sea Urchin Inferred from BAC Insert Deletions.

The SpTransformer (SpTrf) gene family in the purple sea urchin, Strongylocentrotus purpuratus, encodes immune response proteins. The genes are clustered, surrounded by short tandem repeats, and some are present in genomic segmental duplications. The genes share regions of sequence and include repeats in the coding exon. This complex structure is consistent with putative local genomic instability. Instability of the SpTrf gene cluster was tested by 10 days of growth of Escherichia coli harboring bacterial artificial chromosome (BAC) clones of sea urchin genomic DNA with inserts containing SpTrf genes. After the growth period, the BAC DNA inserts were analyzed for size and SpTrf gene content. Clones with multiple SpTrf genes showed a variety of deletions, including loss of one, most, or all genes from the cluster. Alternatively, a BAC insert with a single SpTrf gene was stable. BAC insert instability is consistent with variations in the gene family composition among sea urchins, the types of SpTrf genes in the family, and a reduction in the gene copy number in single coelomocytes. Based on the sequence variability among SpTrf genes within and among sea urchins, local genomic instability of the family may be important for driving sequence diversity in this gene family that would be of benefit to sea urchins in their arms race with marine microbes.

Genome-wide identification of GhRLCK-VII subfamily genes in Gossypium hirsutum and investigation of their functions in resistance to Verticillium wilt

BackgroundThe receptor-like cytoplasmic kinases subfamily VII (RLCK-VII) is critical in regulating plant growth, development, and pattern-triggered immunity. However, a comprehensive exploration of these genes in the allotetraploid Gossypium hirsutum is still lacking. This study aimed to identify RLCK-VII genes in G. hirsutum and investigate their evolutionary history, structural features, expression patterns, and role in plant defense.ResultsSeventy-two RLCK-VII genes in the G. hirsutum genome were unveiled and classified into nine groups following their phylogenetic analysis with Arabidopsis thaliana. Group VII-1 was the largest, accounting for 28%, while Groups VII-2 and VII-3 had only one member each. The analysis using MCScanX revealed that these 72 genes formed 166 collinear gene pairs and were resided on 26 chromosomes of G. hirsutum, suggesting that they were derived from whole genome segmental duplication events. Their calculated Ka/Ks values were below one, implying the occurrence of purification selection during the evolution and inhibition of gene function differentiation/loss. All members of the RLCK-VII subfamily possessed two conserved domains, PKinase-Tyr and PKinase, and several conserved PBS1 kinase subdomains, individually included in one of the ten motifs identified using MEME. The RNA-Seq results showed that RLCK-VII genes exhibited different spatiotemporal expression, indicating their involvement in cotton growth, development, and defense responses to Verticillium dahliae. The transcription patterns of RLCK-VII genes found by RNA-Seq were further validated using qRT-PCR assays after inoculating “20B12” (cotton cultivar) with “V991” (V. dahliae). The virus-induced gene silencing (VIGS) assays uncovered that two RLCK-VII genes (Gohir.A13G227248 and Gohir.A10G219900) were essential to G. hirsutum resistance to Verticillium wilt.ConclusionsThese observations offer valuable insight into the attributes and roles of RLCK-VII genes in G. hirsutum, potentially enable the breeding of new cotton cultivars with enhanced resistance to Verticillium wilt.

Characterization of the MADS-Box Gene Family in Akebia trifoliata and Their Evolutionary Events in Angiosperms.

As the largest clade of modern plants, flower plants have evolved a wide variety of flowers and fruits. MADS-box genes play key roles in regulating plant morphogenesis, while basal eudicots have an evolutionarily important position of acting as an evolutionary bridge between basal angiosperms and core eudicots. Akebia trifoliata is an important member of the basal eudicot group. To study the early evolution of angiosperms, we identified and characterized the MADS-Box gene family on the whole-genome level of A. trifoliata. There were 47 MADS-box genes (13 type I and 34 type II genes) in the A. trifoliata genome; type I genes had a greater gene length and coefficient of variation and a smaller exon number than type II genes. A total of 27 (57.4%) experienced whole or segmental genome duplication and purifying selection. A transcriptome analysis suggested that three and eight genes were involved in whole fruit and seed development, respectively. The diversification and phylogenetic analysis of 1479 type II MADS-box genes of 22 angiosperm species provided some clues indicating that a γ whole genome triplication event of eudicots possibility experienced a two-step process. These results are valuable for improving A. trifoliata fruit traits and theoretically elucidating evolutionary processes of angiosperms, especially eudicots.

The Characterization of the Phloem Protein 2 Gene Family Associated with Resistance to Sclerotinia sclerotiorum in Brassica napus.

In plants, phloem is not only a vital structure that is used for nutrient transportation, but it is also the location of a response that defends against various stresses, named phloem-based defense (PBD). Phloem proteins (PP2s) are among the predominant proteins in phloem, indicating their potential functional role in PBD. Sclerotinia disease (SD), which is caused by the necrotrophic fungal pathogen S. sclerotiorum (Sclerotinia sclerotiorum), is a devastating disease that affects oil crops, especially Brassica napus (B. napus), mainly by blocking nutrition and water transportation through xylem and phloem. Presently, the role of PP2s in SD resistance is still largely estimated. Therefore, in this study, we identified 62 members of the PP2 gene family in the B. napus genome with an uneven distribution across the 19 chromosomes. A phylogenetic analysis classified the BnPP2s into four clusters (I–IV), with cluster I containing the most members (28 genes) as a consequence of its frequent genome segmental duplication. A comparison of the gene structures and conserved motifs suggested that BnPP2 genes were well conserved in clusters II to IV, but were variable in cluster I. Interestingly, the motifs in different clusters displayed unique features, such as motif 6 specifically existing in cluster III and motif 1 being excluded from cluster IV. These results indicated the possible functional specification of BnPP2s. A transcriptome data analysis showed that the genes in clusters II to IV exhibited dynamic expression alternation in tissues and the stimulation of S. sclerotiorum, suggesting that they could participate in SD resistance. A GWAS analysis of a rapeseed population comprising 324 accessions identified four BnPP2 genes that were potentially responsible for SD resistance and a transgenic study that was conducted by transiently expressing BnPP2-6 in tobacco (Nicotiana tabacum) leaves validated their positive role in regulating SD resistance in terms of reduced lesion size after inoculation with S. sclerotiorum hyphal plugs. This study provides useful information on PP2 gene functions in B. napus and could aid elaborated functional studies on the PP2 gene family.

Reconstructing the Evolutionary History of a Highly Conserved Operon Cluster in Gammaproteobacteria and Bacilli.

The evolution of gene order rearrangements within bacterial chromosomes is a fast process. Closely related species can have almost no conservation in long-range gene order. A prominent exception to this rule is a >40 kb long cluster of five core operons (secE-rpoBC-str-S10-spc-alpha) and three variable adjacent operons (cysS, tufB, and ecf) that together contain 57 genes of the transcriptional and translational machinery. Previous studies have indicated that at least part of this operon cluster might have been present in the last common ancestor of bacteria and archaea. Using 204 whole genome sequences, ∼2 Gy of evolution of the operon cluster were reconstructed back to the last common ancestors of the Gammaproteobacteria and of the Bacilli. A total of 163 independent evolutionary events were identified in which the operon cluster was altered. Further examination showed that the process of disconnecting two operons generally follows the same pattern. Initially, a small number of genes is inserted between the operons breaking the concatenation followed by a second event that fully disconnects the operons. While there is a general trend for loss of gene synteny over time, there are examples of increased alteration rates at specific branch points or within specific bacterial orders. This indicates the recurrence of relaxed selection on the gene order within bacterial chromosomes. The analysis of the alternation events indicates that segmental genome duplications and/or transposon-directed recombination play a crucial role in rearrangements of the operon cluster.

Clinical characterisation of a novel SCN5A variant associated with progressive malignant arrhythmia and dilated cardiomyopathy.

The SCN5A gene is implicated in many arrhythmogenic and cardiomyopathic processes. We identified a novel SCN5A variant in a family with significant segregation in individuals affected with progressive sinus and atrioventricular nodal disease, atrial arrhythmia, dilated cardiomyopathy, and early sudden cardiac arrest. A patient pedigree was created following the clinical evaluation of three affected individuals, two monozygotic twins and a paternal half-brother, which lead to the evaluation of a paternal half-sister (four siblings with the same father and three mothers) all of whom experienced varying degrees of atrial arrhythmias, conduction disease, and dilated cardiomyopathy in addition to a paternal history of unexplained death in his 50s with similar autopsy findings. The index male underwent sequencing of 58 genes associated with cardiomyopathies. Sanger sequencing was used to provide data for bases with insufficient coverage and for bases in some known regions of genomic segmental duplications. All clinically significant and novel variants were confirmed by independent Sanger sequencing. All relatives tested were shown to have the same SCN5A variant of unknown significance (p. Asp197His) and the monozygotic twins shared a co-occurring NEXN (p. Glu575*). Segregation analysis demonstrates likely pathogenic trait for the SCN5A variant with an additional possible role for the NEXN variant in combination. There is compelling clinical evidence suggesting that the SCN5A variant p. Asp197His may be re-classified as likely pathogenic based on the segregation analysis of our family of interest. Molecular mechanism studies are pending.

Detection and analysis of ancient segmental duplications in mammalian genomes.

Although segmental duplications (SDs) represent hotbeds for genomic rearrangements and emergence of new genes, there are still no easy-to-use tools for identifying SDs. Moreover, while most previous studies focused on recently emerged SDs, detection of ancient SDs remains an open problem. We developed an SDquest algorithm for SD finding and applied it to analyzing SDs in human, gorilla, and mouse genomes. Our results demonstrate that previous studies missed many SDs in these genomes and show that SDs account for at least 6.05% of the human genome (version hg19), a 17% increase as compared to the previous estimate. Moreover, SDquest classified 6.42% of the latest GRCh38 version of the human genome as SDs, a large increase as compared to previous studies. We thus propose to re-evaluate evolution of SDs based on their accurate representation across multiple genomes. Toward this goal, we analyzed the complex mosaic structure of SDs and decomposed mosaic SDs into elementary SDs, a prerequisite for follow-up evolutionary analysis. We also introduced the concept of the breakpoint graph of mosaic SDs that revealed SD hotspots and suggested that some SDs may have originated from circular extrachromosomal DNA (ecDNA), not unlike ecDNA that contributes to accelerated evolution in cancer.

ASGART: fast and parallel genome scale segmental duplications mapping.

Segmental Duplications (SDs) are DNA fragments longer than 1 kbp, distributed within and between chromosomes and sharing more than 90% identity. Although they hold a significant role in genomic fluidity and adaptability, many key questions about their intrinsic characteristics and mutability remain unsolved due to the persistent difficulty of sequencing highly duplicated genomic regions. The recent development of long and linked-read NGS technologies will increase the need to search for SDs in genomes newly sequenced with these technics. The main limitation of SD analysis will soon be the availability of efficient detection software, to retrieve and compare SD genomic component between species or lineages. In this paper, we present the open-source ASGART, 'A Segmental duplications Gathering And Refining Tool', developed to search for segmental duplications (SDs) in any assembled sequence. We have tested and benchmarked ASGART on five models organisms. Our results demonstrate ASGART's ability to extract SDs from any genome-wide sequence, regardless of genomic size or organizational complexity and quicker than any other software available. The online version of ASGART is available at http://asgart.irit.fr. The source code of ASGART is available both on the ASGART website and at https://github.com/delehef/asgart. Supplementary data are available at Bioinformatics online.

Palindromic amplification of the ERBB2 oncogene in primary HER2-positive breast tumors

Oncogene amplification confers a growth advantage to tumor cells for clonal expansion. There are several, recurrently amplified oncogenes throughout the human genome. However, it remains unclear whether this recurrent amplification is solely a manifestation of increased fitness resulting from random amplification mechanisms, or if a genomic locus-specific amplification mechanism plays a role. Here we show that the ERBB2 oncogene at 17q12 is susceptible to palindromic gene amplification, a mechanism characterized by the inverted (palindromic) duplication of genomic segments, in HER2-positive breast tumors. We applied two genomic approaches to investigate amplification mechanisms: sequencing of DNA libraries enriched with tumor-derived palindromic DNA (Genome-wide Analysis of Palindrome Formation) and whole genome sequencing (WGS). We observed significant enrichment of palindromic DNA within amplified ERBB2 genomic segments. Palindromic DNA was particularly enriched at amplification peaks and at boundaries between amplified and normal copy-number regions. Thus, palindromic gene amplification shaped the amplified ERBB2 locus. The enrichment of palindromic DNA throughout the amplified segments leads us to propose that the ERBB2 locus is amplified through the mechanism that repeatedly generates palindromic DNA, such as Breakage-Fusion-Bridge cycles. The genomic architecture surrounding ERBB2 in the normal genome, such as segmental duplications, could promote the locus-specific mechanism.

Genome-wide investigation and expression analyses of the pentatricopeptide repeat protein gene family in foxtail millet.

BackgroundPentatricopeptide repeat (PPR) proteins are encoded by a large gene family of approximately 450 members in Arabidopsis and 477 in rice, which characterized by tandem repetitions of a degenerate 35 amino acid characteristic sequence motifs. A large majority of the PPR genes in the higher plants are localized in organelles. Their functions remain as yet largely unknown. The majority of characterized PPR proteins have been found to function in modulating the expression plastid and mitochondrial genes in plants.ResultsHere, a genome-wide identification and comparison of the PPR genes from 5 organisms was performed, including the moss Physcomitrella patens, the lycophyte Selaginella moellendorffii, the eudicot Arabidopsis, and the monocots rice and foxtail millet. It appears that the expansion of this gene family prior to the divergence of the euphyllophytes and the lycophytes in land plants. The duplication and divergence rates of the foxtail millet PPR genes (SiPPRs) showed that the expansion period of this gene family around 400 Mya, and indicated that genome segmental duplication was very likely the primary mechanism underlying the expansion of the PPR gene family in vascular plants. An analysis of a complete set of SiPPR genes/proteins that included classification, chromosomal location, orthologous relationships, duplication analysis, and auxiliary motifs is presented. Expression analysis of the SiPPR genes under stress conditions revealed that the expression of 24 SiPPR genes was responsive to abiotic stress. Subcellular localization analysis of 11 PPR proteins indicated that 5 proteins were localized to chloroplasts, that 4 were localized to mitochondria, and that 2 were localized to the cytoplasm.ConclusionsOur results contribute to a more comprehensive understanding the roles of PPR proteins and will be useful in the prioritization of particular PPR proteins for subsequent functional validation studies in foxtail millet.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-3184-2) contains supplementary material, which is available to authorized users.

De novo deletions and duplications of 17q25.3 cause susceptibility to cardiovascular malformations.

BackgroundGenomic disorders resulting from deletion or duplication of genomic segments are known to be an important cause of cardiovascular malformations (CVMs). In our previous study, we identified a unique individual with a de novo 17q25.3 deletion from a study of 714 individuals with CVM.MethodsTo understand the contribution of this locus to cardiac malformations, we reviewed the data on 60,000 samples submitted for array comparative genomic hybridization (CGH) studies to Medical Genetics Laboratories at Baylor College of Medicine, and ascertained seven individuals with segmental aneusomy of 17q25. We validated our findings by studying another individual with a de novo submicroscopic deletion of this region from Cytogenetics Laboratory at Cincinnati Children’s Hospital. Using bioinformatic analyses including protein-protein interaction network, human tissue expression patterns, haploinsufficiency scores, and other annotation systems, including a training set of 251 genes known to be linked to human cardiac disease, we constructed a pathogenicity score for cardiac phenotype for each of the 57 genes within the terminal 2.0 Mb of 17q25.3.ResultsWe found relatively high penetrance of cardiovascular defects (~60 %) with five deletions and three duplications, observed in eight unrelated individuals. Distinct cardiac phenotypes were present in four of these subjects with non-recurrent de novo deletions (range 0.08 Mb–1.4 Mb) in the subtelomeric region of 17q25.3. These included coarctation of the aorta (CoA), total anomalous pulmonary venous return (TAPVR), ventricular septal defect (VSD) and atrial septal defect (ASD). Amongst the three individuals with variable size duplications of this region, one had patent ductus arteriosus (PDA) at 8 months of age.ConclusionThe distinct cardiac lesions observed in the affected patients and the bioinformatics analyses suggest that multiple genes may be plausible drivers of the cardiac phenotype within this gene-rich critical interval of 17q25.3.Electronic supplementary materialThe online version of this article (doi:10.1186/s13023-015-0291-0) contains supplementary material, which is available to authorized users.

Functional roles of evolutionary conserved motifs and residues in vertebrate chemokine receptors.

Chemokine receptors regulate cell migration and homing. They belong to the rhodopsin-like family of GPCRs. Their ancestor genes emerged in the early stages of vertebrate evolution. Since then, the family has been greatly expanded through whole and segmental genome duplication events. During evolution, many amino acid changes have been introduced in individual chemokine receptors, but certain motifs and residues are highly conserved. Previously, we proposed a nomenclature system of the vertebrate chemokine receptors based on their evolutionary history and phylogenetic analyses. With the use of this classification system, we are now able to confidently assign the species orthologs of vertebrate chemokine receptors. Here, we systematically analyze conserved motifs and residues of each group of orthologous chemokine receptors that may play important roles in their signaling and biologic functions. Our present analysis may provide useful information on how individual chemokine receptors are activated upon ligand binding.

Characterization of small HSPs from Anemonia viridis reveals insights into molecular evolution of alpha crystallin genes among cnidarians.

Gene family encoding small Heat-Shock Proteins (sHSPs containing α-crystallin domain) are found both in prokaryotic and eukaryotic organisms; however, there is limited knowledge of their evolution. In this study, two small HSP genes termed AvHSP28.6 and AvHSP27, both organized in one intron and two exons, were characterised in the Mediterranean snakelocks anemone Anemonia viridis. The release of the genome sequence of Hydra magnipapillata and Nematostella vectensis enabled a comprehensive study of the molecular evolution of α-crystallin gene family among cnidarians. Most of the H. magnipapillata sHSP genes share the same gene organization described for AvHSP28.6 and AvHSP27, differing from the sHSP genes of N. vectensis which mainly show an intronless architecture. The different genomic organization of sHSPs, the phylogenetic analyses based on protein sequences, and the relationships among Cnidarians, suggest that the A.viridis sHSPs represent the common ancestor from which H. magnipapillata genes directly evolved through segmental genome duplication. Additionally retroposition events may be considered responsible for the divergence of sHSP genes of N. vectensis from A. viridis. Analyses of transcriptional expression profile showed that AvHSP28.6 was constitutively expressed among different tissues from both ectodermal and endodermal layers of the adult sea anemones, under normal physiological conditions and also under different stress condition. Specifically, we profiled the transcriptional activation of AvHSP28.6 after challenges with different abiotic/biotic stresses showing induction by extreme temperatures, heavy metals exposure and immune stimulation. Conversely, no AvHSP27 transcript was detected in such dissected tissues, in adult whole body cDNA library or under stress conditions. Hence, the involvement of AvHSP28.6 gene in the sea anemone defensome is strongly suggested.

Characterization and Evolution of microRNA Genes Derived from Repetitive Elements and Duplication Events in Plants

MicroRNAs (miRNAs) are a major class of small non-coding RNAs that act as negative regulators at the post-transcriptional level in animals and plants. In this study, all known miRNAs in four plant species (Arabidopsis thaliana, Populus trichocarpa, Oryza sativa and Sorghum bicolor) have been analyzed, using a combination of computational and comparative genomic approaches, to systematically identify and characterize the miRNAs that were derived from repetitive elements and duplication events. The study provides a complete mapping, at the genome scale, of all the miRNAs found on repetitive elements in the four test plant species. Significant differences between repetitive element-related miRNAs and non-repeat-derived miRNAs were observed for many characteristics, including their location in protein-coding and intergenic regions in genomes, their conservation in plant species, sequence length of their hairpin precursors, base composition of their hairpin precursors and the minimum free energy of their hairpin structures. Further analysis showed that a considerable number of miRNA families in the four test plant species arose from either tandem duplication events, segmental duplication events or a combination of the two. However, comparative analysis suggested that the contribution made by these two duplication events differed greatly between the perennial tree species tested and the other three annual species. The expansion of miRNA families in A. thaliana, O. sativa and S. bicolor are more likely to occur as a result of tandem duplication events than from segmental duplications. In contrast, genomic segmental duplications contributed significantly more to the expansion of miRNA families in P. trichocarpa than did tandem duplication events. Taken together, this study has successfully characterized miRNAs derived from repetitive elements and duplication events at the genome scale and provides comprehensive knowledge and deeper insight into the origins and evolution of miRNAs in plants.

Advances in plant genome sequencing

The study of plant biology in the 21st century is, and will continue to be, vastly different from that in the 20th century. One driver for this has been the use of genomics methods to reveal the genetic blueprints for not one but dozens of plant species, as well as resolving genome differences in thousands of individuals at the population level. Genomics technology has advanced substantially since publication of the first plant genome sequence, that of Arabidopsis thaliana, in 2000. Plant genomics researchers have readily embraced new algorithms, technologies and approaches to generate genome, transcriptome and epigenome datasets for model and crop species that have permitted deep inferences into plant biology. Challenges in sequencing any genome include ploidy, heterozygosity and paralogy, all which are amplified in plant genomes compared to animal genomes due to the large genome sizes, high repetitive sequence content, and rampant whole- or segmental genome duplication. The ability to generate de novo transcriptome assemblies provides an alternative approach to bypass these complex genomes and access the gene space of these recalcitrant species. The field of genomics is driven by technological improvements in sequencing platforms; however, software and algorithm development has lagged behind reductions in sequencing costs, improved throughput, and quality improvements. It is anticipated that sequencing platforms will continue to improve the length and quality of output, and that the complementary algorithms and bioinformatic software needed to handle large, repetitive genomes will improve. The future is bright for an exponential improvement in our understanding of plant biology.

Genome change in wheat observed through the structure and expression of α/β-gliadin genes

To better understand genome structure and the expression of α/β-gliadin multigenes in hexaploid wheat, bacterial artificial chromosome (BAC) clones containing α/β-gliadin genes from the three loci, Gli-A2, Gli-B2, and Gli-D2, were screened. Based on their restriction fragment patterns, we selected five BAC clones, namely, two clones for Gli-A2, two clones for Gli-B2, and one clone for Gli-D2, to fully sequence. Approximately 200 kb was sequenced for each locus. In total, twelve α/β-gliadin intact genes and four pseudogenes were found, and retrotransposons or other transposons existed in each BAC clone. Dot-plot analysis revealed the pattern of genome segmental duplication within each BAC. We calculated time since duplication of each set of α/β-gliadin genes and insertion of retrotransposons. Duplication of all adjacent genes within the same BAC clone took place before or after allotetrapolyploidization, but duplication of certain genes occurred before diploid differentiation of wheat species. Retrotransposons were also inserted before and after the segmental duplication events. Furthermore, translocation of α/β-gliadin genes from chromosomes 1 to 6 apparently occurred before the diversification of various wheat genomes. Duplication of genome segments containing α/β-gliadin genes and retrotransposons were brought about through unequal crossing-over or saltatory replication and α/β-gliadin genes per se were duplicated without any recombination events. Out of twelve intact α/β-gliadin genes detected from their sequences, nine were expressed, although their patterns of expression were distinct. Since they have similar cis-elements and promoter structures, the mechanisms underlying their distinct gene expression and possible applications are discussed.

Efficient algorithms for analyzing segmental duplications with deletions and inversions in genomes

BackgroundSegmental duplications, or low-copy repeats, are common in mammalian genomes. In the human genome, most segmental duplications are mosaics comprised of multiple duplicated fragments. This complex genomic organization complicates analysis of the evolutionary history of these sequences. One model proposed to explain this mosaic patterns is a model of repeated aggregation and subsequent duplication of genomic sequences.ResultsWe describe a polynomial-time exact algorithm to compute duplication distance, a genomic distance defined as the most parsimonious way to build a target string by repeatedly copying substrings of a fixed source string. This distance models the process of repeated aggregation and duplication. We also describe extensions of this distance to include certain types of substring deletions and inversions. Finally, we provide a description of a sequence of duplication events as a context-free grammar (CFG).ConclusionThese new genomic distances will permit more biologically realistic analyses of segmental duplications in genomes.

Ubiquitous internal gene duplication and intron creation in eukaryotes

Duplication of genomic segments provides a primary resource for the origin of evolutionary novelties. However, most previous studies have focused on duplications of complete protein-coding genes, whereas little is known about the significance of duplication segments that are entirely internal to genes. Our examination of six fully sequenced genomes reveals that internal duplications of gene segments occur at a high frequency (0.001-0.013 duplications/gene per million years), similar to that of complete gene duplications, such that 8-17% of the genes in a genome carry duplicated intronic and/or exonic regions. At least 7-30% of such genes have acquired novel introns, either because a prior intron in the same gene has been duplicated, or more commonly, because a spatial change has activated a latent splice site. These results strongly suggest a major evolutionary role for internal gene duplications in the origin of genomic novelties, particularly as a mechanism for intron gain.

Evidence and evolutionary analysis of ancient whole-genome duplication in barley predating the divergence from rice

BackgroundWell preserved genomic colinearity among agronomically important grass species such as rice, maize, Sorghum, wheat and barley provides access to whole-genome structure information even in species lacking a reference genome sequence. We investigated footprints of whole-genome duplication (WGD) in barley that shaped the cereal ancestor genome by analyzing shared synteny with rice using a ~2000 gene-based barley genetic map and the rice genome reference sequence.ResultsBased on a recent annotation of the rice genome, we reviewed the WGD in rice and identified 24 pairs of duplicated genomic segments involving 70% of the rice genome. Using 968 putative orthologous gene pairs, synteny covered 89% of the barley genetic map and 63% of the rice genome. We found strong evidence for seven shared segmental genome duplications, corresponding to more than 50% of the segmental genome duplications previously determined in rice. Analysis of synonymous substitution rates (Ks) suggested that shared duplications originated before the divergence of these two species. While major genome rearrangements affected the ancestral genome of both species, small paracentric inversions were found to be species specific.ConclusionWe provide a thorough analysis of comparative genome evolution between barley and rice. A barley genetic map of approximately 2000 non-redundant EST sequences provided sufficient density to allow a detailed view of shared synteny with the rice genome. Using an indirect approach that included the localization of WGD-derived duplicated genome segments in the rice genome, we determined the current extent of shared WGD-derived genome duplications that occurred prior to species divergence.

PseudoGeneQuest – Service for identification of different pseudogene types in the human genome

BackgroundPseudogenes, nonfunctional copies of genes, evolve fast due the lack of evolutionary pressures and thus appear in several different forms. PseudoGeneQuest is an online tool to search the human genome for a given query sequence and to identify different types of pseudogenes as well as novel genes and gene fragments.DescriptionThe service can detect pseudogenes, that have arisen either by retrotransposition or segmental genome duplication, many of which are not listed in the public pseudogene databases. The service has a user-friendly web interface and uses a powerful computer cluster in order to perform parallel searches and provide relatively fast runtimes despite exhaustive database searches and analyses.ConclusionPseudoGeneQuest is a versatile tool for detecting novel pseudogene candidates from the human genome. The service searches human genome sequences for five types of pseudogenes and provides an output that allows easy further analysis of observations. In addition to the result file the system provides visualization of the results linked to Ensembl Genome Browser. PseudoGeneQuest service is freely available.