Evidence for extensive recent intron transposition in closely related fungi

Though spliceosomal introns are a major structural component of most eukaryotic genes and intron density varies by more than three orders of magnitude among eukaryotes [1-3], the origins of introns are poorly understood, and only a few cases of unambiguous intron gain are known [4-8]. We utilized population genomic comparisons of three closely related fungi to identify crucial transitory phases of intron gain and loss. We found 74 intron positions showing intraspecific presence-absence polymorphisms (PAPs) for the entire intron. Population genetic analyses identified intron PAPs at different stages of fixation and showed that intron gain or loss was very recent. We found direct support for extensive intron transposition among unrelated genes. A substantial proportion of highly similar introns in the genome either were recently gained or showed a transient phase of intron PAP. We also identified an intron transfer among paralogous genes that created a new intron. Intron loss was due mainly to homologous recombination involving reverse-transcribed mRNA. The large number of intron positions in transient phases of either intron gain or loss shows that intron evolution is much faster than previously thought and provides an excellent model to study molecular mechanisms of intron gain.     

Genome‐wide copy number profiling using high‐density SNP array in chickens

Phenotypic diversity is a direct consequence resulting mainly from the impact of underlying genetic variation, and recent studies have shown that copy number variation (CNV) is emerging as an important contributor to both phenotypic variability and disease susceptibility. Herein, we performed a genome‐wide CNV scan in 96 chickens from 12 diversified breeds, benefiting from the high‐density Affymetrix 600 K SNP arrays. We identified a total of 231 autosomal CNV regions (CNVRs) encompassing 5.41 Mb of the chicken genome and corresponding to 0.59% of the autosomal sequence. The length of these CNVRs ranged from 2.6 to 586.2 kb with an average of 23.4 kb, including 130 gain, 93 loss and eight both gain and loss events. These CNVRs, especially deletions, had lower GC content and were located particularly in gene deserts. In particular, 102 CNVRs harbored 128 chicken genes, most of which were enriched in immune responses. We obtained 221 autosomal CNVRs after converting probe coordinates to Galgal3, and comparative analysis with previous studies illustrated that 153 of these CNVRs were regarded as novel events. Furthermore, qPCR assays were designed for 11 novel CNVRs, and eight (72.73%) were validated successfully. In this study, we demonstrated that the high‐density 600 K SNP array can capture CNVs with higher efficiency and accuracy and highlighted the necessity of integrating multiple technologies and algorithms. Our findings provide a pioneering exploration of chicken CNVs based on a high‐density SNP array, which contributes to a more comprehensive understanding of genetic variation in the chicken genome and is beneficial to unearthing potential CNVs underlying important traits of chickens.     

Ready reference to common segmental aneusomies by syndromic names, major features and chromosomal locations

Abnormal gene dosage usually results in recognizable phenotypic abnormalities, especially if it involves a series of contiguous genes. Schmickel (1986) defined contiguous gene syndromes as diseases resulting from loss or gain of a series of adjacent genes. The terms microdeletion and microduplication as well as segmental aneusomy have also been used to describe such losses or gains that may not be readily detectable by Gbanded analysis. The loss (haploinsufficiency) or gain of a series of adjoining genes may result in a direct phenotypic effect and/or cause a genetic regulatory disturbance. Such syndromic gains or losses are often detectable when in situ hybridization of fluorescent labeled DNA probes or array comparative genomic hybridization technique are used (Gersen and Keagle 2005; Stumm et al. 1999; Barch, Knutsen and Spurbeck 1997). Segmental aneusomies generally occur due to homologous pairing between non-allelic low copy repeats (LCR) followed by crossing over. The LCRs, as part of the repetitive DNA sequences range from 1-500 Kb repeats, share >97% base sequence identity and constitute up to five percent of the genomic DNA. They are distributed throughout the genome, but are more concentrated near the centromeres and telomeres. A segment of 300 bp completely identical sequence within the LCRs is adequate for mediating non-allelic homologous or paralogous pairing. This process results in generating both deletion and duplication of a defined segment.     

Quantitative comparison of cDNA-AFLP, microarrays, and GeneChip expression data in Saccharomyces cerevisiae

cDNA-AFLP is a genome-wide expression analysis technology that does not require any prior knowledge of gene sequences. This PCR-based technique combines a high sensitivity with a high specificity, allowing detection of rarely expressed genes and distinguishing between homologous genes. In this report, we validated quantitative expression data of 110 cDNA-AFLP fragments in yeast with DNA microarrays and GeneChip data. The best correlation was found between cDNA-AFLP and GeneChip data. The cDNA-AFLP data revealed a low number of inconsistent profiles that could be explained by gel artifact, overexposure, or mismatch amplification. In addition, 18 cDNA-AFLP fragments displayed homology to genomic yeast DNA, but could not be linked unambiguously to any known ORF. These fragments were most probably derived from 5' or 3' noncoding sequences or might represent previously unidentified ORFs. Genes liable to cross hybridization showed identical results in cDNA-AFLP and GeneChip analysis. Three genes, which were readily detected with cDNA-AFLP, showed no significant expression in GeneChip experiments. We show that cDNA-AFLP is a very good alternative to microarrays and since no preexisting biological or sequence information is required, it is applicable to any species.     

In vitro and in vivo complementation of the Helicobacter pylori arginase mutant using an intergenic chromosomal site

BACKGROUND: Gene complementation strategies are important in validating the roles of genes in specific phenotypes. Complementation systems in Helicobacter pylori include shuttle vectors, which transform H. pylori at relatively low frequencies, and chromosomally based approaches. Chromosomal complementation strategies are susceptible to polar effects and disruption of other H. pylori genes, leading to unwanted pleiotropic effects. MATERIALS AND METHODS: A new complementation strategy was developed for H. pylori by utilizing a suicide plasmid vector that contains fragments of an H. pylori intergenic region (hp0203-hp0204), a chloramphenicol acetyltransferase cassette (cat), and a multiple-cloning site. Genes of interest could be cloned into the intergenic plasmid and the genes integrated into H. pylori by homologous recombination into the intergenic chromosomal region without disrupting any annotated H. pylori gene. The complementation system was validated using the gene encoding arginase (rocF). RESULTS: A rocF mutant unable to hydrolyze or consume l-arginine regained these functions by complementation with the wild-type rocF gene. Complemented strains also had restored arginase protein as determined by Western blot analysis. The complementation system could be successfully applied to multiple H. pylori strains. The intergenic region varied in length and sequence across 17 H. pylori strains, but the flanking-3' ends of the hp0203 and hp0204 coding regions were highly conserved. Inserting a cat cassette and wild-type rocF into the intergenic region did not alter the ability of strain SS1 to colonize mice. CONCLUSIONS: This complementation strategy should greatly facilitate genetic experiments in H. pylori.     

Genome of Helicobacter pylori strain 908

Helicobacter pylori is a genetically diverse and coevolved pathogen inhabiting human gastric niches and leading to a spectrum of gastric diseases in susceptible populations. We describe the genome sequence of H. pylori 908, which was originally isolated from an African patient living in France who suffered with recrudescent duodenal ulcer disease. The strain was found to be phylogenetically related to H. pylori J99, and its comparative analysis revealed several specific genome features and novel insertion-deletion and substitution events. The genome sequence revealed several strain-specific deletions and/or gain of genes exclusively present in HP908 compared with different sequenced genomes already available in the public domain. Comparative and functional genomics of HP908 and its subclones will be important in understanding genomic plasticity and the capacity to colonize and persist in a changing host environment.     

Complete genomic sequence of bacteriophage B3, a Mu-like phage of Pseudomonas aeruginosa

Bacteriophage B3 is a transposable phage of Pseudomonas aeruginosa. In this report, we present the complete DNA sequence and annotation of the B3 genome. DNA sequence analysis revealed that the B3 genome is 38,439 bp long with a G+C content of 63.3%. The genome contains 59 proposed open reading frames (ORFs) organized into at least three operons. Of these ORFs, the predicted proteins from 41 ORFs (68%) display significant similarity to other phage or bacterial proteins. Many of the predicted B3 proteins are homologous to those encoded by the early genes and head genes of Mu and Mu-like prophages found in sequenced bacterial genomes. Only two of the predicted B3 tail proteins are homologous to other well-characterized phage tail proteins; however, several Mu-like prophages and transposable phage D3112 encode approximately 10 highly similar proteins in their predicted tail gene regions. Comparison of the B3 genomic organization with that of Mu revealed evidence of multiple genetic rearrangements, the most notable being the inversion of the proposed B3 immunity/early gene region, the loss of Mu-like tail genes, and an extreme leftward shift of the B3 DNA modification gene cluster. These differences illustrate and support the widely held view that tailed phages are genetic mosaics arising by the exchange of functional modules within a diverse genetic pool.     

Comparison of eukaryotic phylogenetic profiling approaches using species tree aware methods

Background  

Phylogenetic profiling encompasses an important set of methodologies for in silico high throughput inference of functional relationships between genes. The simplest profiles represent the distribution of gene presence-absence in a set of species as a sequence of 0's and 1's, and it is assumed that functionally related genes will have more similar profiles. The methodology has been successfully used in numerous studies of prokaryotic genomes, although its application in eukaryotes appears problematic, with reported low accuracy due to the complex genomic organization within this domain of life. Recently some groups have proposed an alternative approach based on the correlation of homologous gene group sizes, taking into account all potentially informative genetic events leading to a change in group size, regardless of whether they result in a de novo group gain or total gene group loss.     

Characterization of lipooligosaccharide-biosynthetic loci of Campylobacter jejuni reveals new lipooligosaccharide classes: evidence of mosaic organizations

The lipooligosaccharide (LOS) biosynthesis region is one of the more variable genomic regions between strains of Campylobacter jejuni. Indeed, eight classes of LOS biosynthesis loci have been established previously based on gene content and organization. In this study, we characterize additional classes of LOS biosynthesis loci and analyze various mechanisms that result in changes to LOS structures. To gain further insights into the genomic diversity of C. jejuni LOS biosynthesis region, we sequenced the LOS biosynthesis loci of 15 strains that possessed gene content that was distinct from the eight classes. This analysis identified 11 new classes of LOS loci that exhibited examples of deletions and insertions of genes and cassettes of genes found in other LOS classes or capsular biosynthesis loci leading to mosaic LOS loci. The sequence analysis also revealed both missense mutations leading to "allelic" glycosyltransferases and phase-variable and non-phase-variable gene inactivation by the deletion or insertion of bases. Specifically, we demonstrated that gene inactivation is an important mechanism for altering the LOS structures of strains possessing the same class of LOS biosynthesis locus. Together, these observations suggest that LOS biosynthesis region is a hotspot for genetic exchange and variability, often leading to changes in the LOS produced.     

Microarrays: how many do you need?

We estimate the number of microarrays that is required in order to gain reliable results from a common type of study: the pairwise comparison of different classes of samples. We show that current knowledge allows for the construction of models that look realistic with respect to searches for individual differentially expressed genes and derive prototypical parameters from real data sets. Such models allow investigation of the dependence of the required number of samples on the relevant parameters: the biological variability of the samples within each class, the fold changes in expression that are desired to be detected, the detection sensitivity of the microarrays, and the acceptable error rates of the results. We supply experimentalists with general conclusions as well as a freely accessible Java applet at www.scai.fhg.de/special/bio/howmanyarrays/ for fine tuning simulations to their particular settings.     

Phylogenomic analysis reveals dynamic evolutionary history of the Drosophila heterochromatin protein 1 (HP1) gene family

Heterochromatin is the gene-poor, satellite-rich eukaryotic genome compartment that supports many essential cellular processes. The functional diversity of proteins that bind and often epigenetically define heterochromatic DNA sequence reflects the diverse functions supported by this enigmatic genome compartment. Moreover, heterogeneous signatures of selection at chromosomal proteins often mirror the heterogeneity of evolutionary forces that act on heterochromatic DNA. To identify new such surrogates for dissecting heterochromatin function and evolution, we conducted a comprehensive phylogenomic analysis of the Heterochromatin Protein 1 gene family across 40 million years of Drosophila evolution. Our study expands this gene family from 5 genes to at least 26 genes, including several uncharacterized genes in Drosophila melanogaster. The 21 newly defined HP1s introduce unprecedented structural diversity, lineage-restriction, and germline-biased expression patterns into the HP1 family. We find little evidence of positive selection at these HP1 genes in both population genetic and molecular evolution analyses. Instead, we find that dynamic evolution occurs via prolific gene gains and losses. Despite this dynamic gene turnover, the number of HP1 genes is relatively constant across species. We propose that karyotype evolution drives at least some HP1 gene turnover. For example, the loss of the male germline-restricted HP1E in the obscura group coincides with one episode of dramatic karyotypic evolution, including the gain of a neo-Y in this lineage. This expanded compendium of ovary- and testis-restricted HP1 genes revealed by our study, together with correlated gain/loss dynamics and chromosome fission/fusion events, will guide functional analyses of novel roles supported by germline chromatin.     

Homology Requirements for Unequal Crossing over in Humans

To gain insight into mechanisms of unequal homologous recombination in vivo, genes generated by homologous unequal crossovers in the human beta-globin gene cluster were examined by nucleotide sequencing and hybridization experiments. The naturally occurring genes studied included one delta-beta Lepore-Baltimore fusion gene, one delta-beta Lepore-Hollandia fusion gene, 12 delta-beta Lepore-Boston genes, one A gamma-beta fusion Kenya gene, one A gamma-G gamma fusion (the central gene of a triplication) and one G gamma-A gamma fusion. A comparison of the nucleotide sequences of three Lepore-Boston genes indicates that they were derived from at least two independent homologous but unequal crossover events, although the crossovers occurred within the same 58-bp region. Nine additional Lepore-Boston genes from individuals of various ethnic origins were shown, by hybridization to specific oligonucleotide probes, to have been generated by a crossover in the same region as the sequenced genes. Evidence for gene conversion accompanying a homologous unequal crossover event was found in only one case (although some of the single nucleotide differences observed in other genes in this study may be related to the crossover events in ways that we do not presently understand). Thus, as judged by this limited sample, concurrent gene conversions are not commonly associated with homologous but unequal exchange in humans in vivo. Classification of the recombinant chromosomes by their polymorphic restriction sites in the beta-globin gene cluster indicated that the Lepore-Boston genes are found in at least six different haplotype backgrounds. Therefore the total number of independent examples in this study is at least 6, and at most 12. We have shown that in at least six cases of genes that have arisen by homologous but unequal crossing over in vivo, each event occurred in a relatively extensive region of uninterrupted identity between the parental genes. This preference cannot be explained by a mechanism whereby crossovers occur at random within misaligned related but not identical genes. In general, crossovers occur in regions that are among the largest available stretches of identity for a particular pair of mismatched genes. Our data are in agreement with those of other types of studies of homologous recombination, and support the idea that sequence identity, rather than general homology, is a critical factor in homologous recombination.     

Gain and loss of multiple genes during the evolution of Helicobacter pylori

Sequence diversity and gene content distinguish most isolates of Helicobacter pylori. Even greater sequence differences differentiate distinct populations of H. pylori from different continents, but it was not clear whether these populations also differ in gene content. To address this question, we tested 56 globally representative strains of H. pylori and four strains of Helicobacter acinonychis with whole genome microarrays. Of the weighted average of 1,531 genes present in the two sequenced genomes, 25% are absent in at least one strain of H. pylori and 21% were absent or variable in H. acinonychis. We extrapolate that the core genome present in all isolates of H. pylori contains 1,111 genes. Variable genes tend to be small and possess unusual GC content; many of them have probably been imported by horizontal gene transfer. Phylogenetic trees based on the microarray data differ from those based on sequences of seven genes from the core genome. These discrepancies are due to homoplasies resulting from independent gene loss by deletion or recombination in multiple strains, which distort phylogenetic patterns. The patterns of these discrepancies versus population structure allow a reconstruction of the timing of the acquisition of variable genes within this species. Variable genes that are located within the cag pathogenicity island were apparently first acquired en bloc after speciation. In contrast, most other variable genes are of unknown function or encode restriction/modification enzymes, transposases, or outer membrane proteins. These seem to have been acquired prior to speciation of H. pylori and were subsequently lost by convergent evolution within individual strains. Thus, the use of microarrays can reveal patterns of gene gain or loss when examined within a phylogenetic context that is based on sequences of core genes.     

Helicobacter pylori evolution and phenotypic diversification in a changing host

Helicobacter pylori colonizes the stomachs of more than 50% of the world's population, making it one of the most successful of all human pathogens. One striking characteristic of H. pylori biology is its remarkable allelic diversity and genetic variability. Not only does almost every infected person harbour their own individual H. pylori strain, but strains can undergo genetic alteration in vivo, driven by an elevated mutation rate and frequent intraspecific recombination. This genetic variability, which affects both housekeeping and virulence genes, has long been thought to contribute to host adaptation, and several recently published studies support this concept. We review the available knowledge relating to the genetic variation of H. pylori, with special emphasis on the changes that occur during chronic colonization, and argue that H. pylori uses mutation and recombination processes to adapt to its individual host by modifying molecules that interact with the host. Finally, we put forward the hypothesis that the lack of opportunity for intraspecies recombination as a result of the decreasing prevalence of H. pylori could accelerate its disappearance from Western populations.     

Suppression of homologous and homeologous recombination by the bacterial MutS2 protein

In addition to their role in DNA repair, recombination events are associated with processes aimed at providing the genetic variability needed for adaptation and evolution of a population. In bacteria, recombination is involved in the appearance of new variants by allowing the incorporation of exogenous DNA or the reshuffling of endogenous sequences. Here we show that HpMutS2, a protein belonging to the MutS2 family in Helicobacter pylori, is not involved in mismatch repair but inhibits homologous and homeologous recombination. Disruption of HpmutS2 leads to an increased efficiency of exogenous DNA incorporation. HpMutS2 has a selective affinity for DNA structures mimicking recombination intermediates with no specificity for homoduplex DNA or mismatches. The purified protein has an ATPase activity stimulated by the same DNA structures. Finally, we show that HpMutS2 inhibits DNA strand exchange reactions in vitro. Thus, MutS2 proteins are candidates for controlling recombination and therefore genetic diversity in bacteria.