Intrastrain and interstrain genetic variation within a paralogous gene family in Chlamydia pneumoniae

Background

Chlamydia pneumoniae causes human respiratory diseases and has recently been associated with atherosclerosis. Analysis of the three recently published C. pneumoniae genomes has led to the identification of a new gene family (the Cpn 1054 family) that consists of 11 predicted genes and gene fragments. Each member encodes a polypeptide with a hydrophobic domain characteristic of proteins localized to the inclusion membrane.

Results

Comparative analysis of this gene family within the published genome sequences provided evidence that multiple levels of genetic variation are evident within this single collection of paralogous genes. Frameshift mutations are found that result in both truncated gene products and pseudogenes that vary among isolates. Several genes in this family contain polycytosine (polyC) tracts either upstream or within the terminal 5' end of the predicted coding sequence. The length of the polyC stretch varies between paralogous genes and within single genes in the three genomes. Sequence analysis of genomic DNA from a collection of 12 C. pneumoniae clinical isolates was used to determine the extent of the variation in the Cpn 1054 gene family.

Conclusions

These studies demonstrate that sequence variability is present both among strains and within strains at several of the loci. In particular, changes in the length of the polyC tract associated with the different Cpn 1054 gene family members are common within each tested C. pneumoniae isolate. The variability identified within this newly described gene family may modulate either phase or antigenic variation and subsequent physiologic diversity within a C. pneumoniae population.     

Asymmetric histone modifications between the original and derived loci of human segmental duplications

Background

Sequencing and annotation of several mammalian genomes have revealed that segmental duplications are a common architectural feature of primate genomes; in fact, about 5% of the human genome is composed of large blocks of interspersed segmental duplications. These segmental duplications have been implicated in genomic copy-number variation, gene novelty, and various genomic disorders. However, the molecular processes involved in the evolution and regulation of duplicated sequences remain largely unexplored.

Results

In this study, the profile of about 20 histone modifications within human segmental duplications was characterized using high-resolution, genome-wide data derived from a ChIP-Seq study. The analysis demonstrates that derivative loci of segmental duplications often differ significantly from the original with respect to many histone methylations. Further investigation showed that genes are present three times more frequently in the original than in the derivative, whereas pseudogenes exhibit the opposite trend. These asymmetries tend to increase with the age of segmental duplications. The uneven distribution of genes and pseudogenes does not, however, fully account for the asymmetry in the profile of histone modifications.

Conclusion

The first systematic analysis of histone modifications between segmental duplications demonstrates that two seemingly 'identical' genomic copies are distinct in their epigenomic properties. Results here suggest that local chromatin environments may be implicated in the discrimination of derived copies of segmental duplications from their originals, leading to a biased pseudogenization of the new duplicates. The data also indicate that further exploration of the interactions between histone modification and sequence degeneration is necessary in order to understand the divergence of duplicated sequences.     

The ribosomal protein genes and Minute loci of Drosophila melanogaster

Background

Mutations in genes encoding ribosomal proteins (RPs) have been shown to cause an array of cellular and developmental defects in a variety of organisms. In Drosophila melanogaster, disruption of RP genes can result in the 'Minute' syndrome of dominant, haploinsufficient phenotypes, which include prolonged development, short and thin bristles, and poor fertility and viability. While more than 50 Minute loci have been defined genetically, only 15 have so far been characterized molecularly and shown to correspond to RP genes.

Results

We combined bioinformatic and genetic approaches to conduct a systematic analysis of the relationship between RP genes and Minute loci. First, we identified 88 genes encoding 79 different cytoplasmic RPs (CRPs) and 75 genes encoding distinct mitochondrial RPs (MRPs). Interestingly, nine CRP genes are present as duplicates and, while all appear to be functional, one member of each gene pair has relatively limited expression. Next, we defined 65 discrete Minute loci by genetic criteria. Of these, 64 correspond to, or very likely correspond to, CRP genes; the single non-CRP-encoding Minute gene encodes a translation initiation factor subunit. Significantly, MRP genes and more than 20 CRP genes do not correspond to Minute loci.

Conclusion

This work answers a longstanding question about the molecular nature of Minute loci and suggests that Minute phenotypes arise from suboptimal protein synthesis resulting from reduced levels of cytoribosomes. Furthermore, by identifying the majority of haplolethal and haplosterile loci at the molecular level, our data will directly benefit efforts to attain complete deletion coverage of the D. melanogaster genome.     

Comparative genomics of mutualistic viruses of Glyptapanteles parasitic wasps

Background

Polydnaviruses, double-stranded DNA viruses with segmented genomes, have evolved as obligate endosymbionts of parasitoid wasps. Virus particles are replication deficient and produced by female wasps from proviral sequences integrated into the wasp genome. These particles are co-injected with eggs into caterpillar hosts, where viral gene expression facilitates parasitoid survival and, thereby, survival of proviral DNA. Here we characterize and compare the encapsidated viral genome sequences of bracoviruses in the family Polydnaviridae associated with Glyptapanteles gypsy moth parasitoids, along with near complete proviral sequences from which both viral genomes are derived.

Results

The encapsidated Glyptapanteles indiensis and Glyptapanteles flavicoxis bracoviral genomes, each composed of 29 different size segments, total approximately 517 and 594 kbp, respectively. They are generated from a minimum of seven distinct loci in the wasp genome. Annotation of these sequences revealed numerous novel features for polydnaviruses, including insect-like sugar transporter genes and transposable elements. Evolutionary analyses suggest that positive selection is widespread among bracoviral genes.

Conclusions

The structure and organization of G. indiensis and G. flavicoxis bracovirus proviral segments as multiple loci containing one to many viral segments, flanked and separated by wasp gene-encoding DNA, is confirmed. Rapid evolution of bracovirus genes supports the hypothesis of bracovirus genes in an 'arms race' between bracovirus and caterpillar. Phylogenetic analyses of the bracoviral genes encoding sugar transporters provides the first robust evidence of a wasp origin for some polydnavirus genes. We hypothesize transposable elements, such as those described here, could facilitate transfer of genes between proviral segments and host DNA.     

The Dictyostelium genome encodes numerous RasGEFs with multiple biological roles

Background

Dictyostelium discoideum is a eukaryote with a simple lifestyle and a relatively small genome whose sequence has been fully determined. It is widely used for studies on cell signaling, movement and multicellular development. Ras guanine-nucleotide exchange factors (RasGEFs) are the proteins that activate Ras and thus lie near the top of many signaling pathways. They are particularly important for signaling in development and chemotaxis in many organisms, including Dictyostelium.

Results

We have searched the genome for sequences encoding RasGEFs. Despite its relative simplicity, we find that the Dictyostelium genome encodes at least 25 RasGEFs, with a few other genes encoding only parts of the RasGEF consensus domains. All appear to be expressed at some point in development. The 25 genes include a wide variety of domain structures, most of which have not been seen in other organisms. The LisH domain, which is associated with microtubule binding, is seen particularly frequently; other domains that confer interactions with the cytoskeleton are also common. Disruption of a sample of the novel genes reveals that many have clear phenotypes, including altered morphology and defects in chemotaxis, slug phototaxis and thermotaxis.

Conclusion

These results suggest that the unexpectedly large number of RasGEF genes reflects an evolutionary expansion of the range of Ras signaling rather than functional redundancy or the presence of multiple pseudogenes.     

The human HOX gene family.

We report the identification of 10 new human homeobox sequences. Altogether, we have isolated and sequenced 30 human homeoboxes clustered in 4 chromosomal regions called HOX loci. HOX1 includes 8 homeoboxes in 90 kb of DNA on chromosome 7. HOX2 includes 9 homeoboxes in 180 kb on chromosome 17. HOX3 contains at least 7 homeoboxes in 160 kb on chromosome 12. Finally, HOX4 includes 6 homeoboxes in 70 kb on chromosome 2. Homeodomains obtained from the conceptual translation of the isolated homeoboxes can be attributed to 13 homology groups on the basis of their primary peptide sequence. Moreover, it is possible to align the 4 HOX loci so that corresponding homeodomains in all loci share the maximal sequence identity. The complex of these observations supports and extends an evolutionary hypothesis concerning the origin of mammalian and fly homeobox gene complexes. We also determined the coding region present in 3 HOX2 cDNA clones corresponding to HOX2G, HOX2H and HOX2I.     

Evolution and diversity of Rickettsia bacteria

Background

Facioscapulohumeral muscular dystrophy (FSHD) is an autosomal dominant neuromuscular disorder associated with the partial deletion of integral numbers of 3.3 kb D4Z4 DNA repeats within the subtelomere of chromosome 4q. A number of candidate FSHD genes, adenine nucleotide translocator 1 gene (ANT1), FSHD-related gene 1 (FRG1), FRG2 and DUX4c, upstream of the D4Z4 array (FSHD locus), and double homeobox chromosome 4 (DUX4) within the repeat itself, are upregulated in some patients, thus suggesting an underlying perturbation of the chromatin structure. Furthermore, a mouse model overexpressing FRG1 has been generated, displaying skeletal muscle defects.

Results

In the context of myogenic differentiation, we compared the chromatin structure and tridimensional interaction of the D4Z4 array and FRG1 gene promoter, and FRG1 expression, in control and FSHD cells. The FRG1 gene was prematurely expressed during FSHD myoblast differentiation, thus suggesting that the number of D4Z4 repeats in the array may affect the correct timing of FRG1 expression. Using chromosome conformation capture (3C) technology, we revealed that the FRG1 promoter and D4Z4 array physically interacted. Furthermore, this chromatin structure underwent dynamic changes during myogenic differentiation that led to the loosening of the FRG1/4q-D4Z4 array loop in myotubes. The FRG1 promoter in both normal and FSHD myoblasts was characterized by H3K27 trimethylation and Polycomb repressor complex binding, but these repression signs were replaced by H3K4 trimethylation during differentiation. The D4Z4 sequences behaved similarly, with H3K27 trimethylation and Polycomb binding being lost upon myogenic differentiation.

Conclusion

We propose a model in which the D4Z4 array may play a critical chromatin function as an orchestrator of in cis chromatin loops, thus suggesting that this repeat may play a role in coordinating gene expression.     

Gene organization of murine homeobox-containing gene clusters

A chromosomal walk which links a previously described and a new homeobox to the Hox-2 murine homeobox gene cluster is described, and the nucleotide sequence of the new homeobox is presented. With these new data the Hox-2 gene cluster contains seven loci on an approximately 100-kb-long physical map. Homology comparisons reveal that a significant number of vertebrate homeoboxes are in fact analogous. We also find that the linear order of homologous homeoboxes is similar in the two murine gene complexes, Hox-1 and Hox-2, and among the human homeobox loci on chromosome 17. Conservation of the homeo-domain and the linear gene order of homeobox-containing genes in vertebrates is discussed in light of the interactions and the anteroposterior linear order of homeotic loci in insects.     

Comparison of Francisella tularensis genomes reveals evolutionary events associated with the emergence of human pathogenic strains

Background

Francisella tularensis subspecies tularensis and holarctica are pathogenic to humans, whereas the two other subspecies, novicida and mediasiatica, rarely cause disease. To uncover the factors that allow subspecies tularensis and holarctica to be pathogenic to humans, we compared their genome sequences with the genome sequence of Francisella tularensis subspecies novicida U112, which is nonpathogenic to humans.

Results

Comparison of the genomes of human pathogenic Francisella strains with the genome of U112 identifies genes specific to the human pathogenic strains and reveals pseudogenes that previously were unidentified. In addition, this analysis provides a coarse chronology of the evolutionary events that took place during the emergence of the human pathogenic strains. Genomic rearrangements at the level of insertion sequences (IS elements), point mutations, and small indels took place in the human pathogenic strains during and after differentiation from the nonpathogenic strain, resulting in gene inactivation.

Conclusion

The chronology of events suggests a substantial role for genetic drift in the formation of pseudogenes in Francisella genomes. Mutations that occurred early in the evolution, however, might have been fixed in the population either because of evolutionary bottlenecks or because they were pathoadaptive (beneficial in the context of infection). Because the structure of Francisella genomes is similar to that of the genomes of other emerging or highly pathogenic bacteria, this evolutionary scenario may be shared by pathogens from other species.     

Complex Evolution of a Y-Chromosomal Double Homeobox 4 (DUX4)-Related Gene Family in Hominoids

The human Y chromosome carries four human Y-chromosomal euchromatin/heterochromatin transition regions, all of which are characterized by the presence of interchromosomal segmental duplications. The Yq11.1/Yq11.21 transition region harbours a peculiar segment composed of an imperfectly organized tandem-repeat structure encoding four members of the double homeobox (DUX) gene family. By comparative fluorescence in situ hybridization (FISH) analysis we have documented the primary appearance of Y-chromosomal DUX genes (DUXY) on the gibbon Y chromosome. The major amplification and dispersal of DUXY paralogs occurred after the gibbon and hominid lineages had diverged. Orthologous DUXY loci of human and chimpanzee show a highly similar structural organization. Sequence alignment survey, phylogenetic reconstruction and recombination detection analyses of human and chimpanzee DUXY genes revealed the existence of all copies in a common ancestor. Comparative analysis of the circumjacent beta-satellites indicated that DUXY genes and beta-satellites evolved in concert. However, evolutionary forces acting on DUXY genes may have induced amino acid sequence differences in the orthologous chimpanzee and human DUXY open reading frames (ORFs). The acquisition of complete ORFs in human copies might relate to evolutionary advantageous functions indicating neo-functionalization. We propose an evolutionary scenario in which an ancestral tandem array DUX gene cassette transposed to the hominoid Y chromosome followed by lineage-specific chromosomal rearrangements paved the way for a species-specific evolution of the Y-chromosomal members of a large highly diverged homeobox gene family.     

Quantitative trait loci of stripe rust resistance in wheat

Key message

Over 140 QTLs for resistance to stripe rust in wheat have been published and through mapping flanking markers on consensus maps, 49 chromosomal regions are identified.

Abstract

Over thirty publications during the last 10 years have identified more than 140 QTLs for stripe rust resistance in wheat. It is likely that many of these QTLs are identical genes that have been spread through plant breeding into diverse backgrounds through phenotypic selection under stripe rust epidemics. Allelism testing can be used to differentiate genes in similar locations but in different genetic backgrounds; however, this is problematic for QTL studies where multiple loci segregate from any one parent. This review utilizes consensus maps to illustrate important genomic regions that have had effects against stripe rust in wheat, and although this methodology cannot distinguish alleles from closely linked genes, it does highlight the extent of genetic diversity for this trait and identifies the most valuable loci and the parents possessing them for utilization in breeding programs. With the advent of cheaper, high throughput genotyping technologies, it is envisioned that there will be many more publications in the near future describing ever more QTLs. This review sets the scene for the coming influx of data and will quickly enable researchers to identify new loci in their given populations.     

Genome-wide analyses implicate 33 loci in heritable dog osteosarcoma, including regulatory variants near CDKN2A/B

Background

Canine osteosarcoma is clinically nearly identical to the human disease, but is common and highly heritable, making genetic dissection feasible.

Results

Through genome-wide association analyses in three breeds (greyhounds, Rottweilers, and Irish wolfhounds), we identify 33 inherited risk loci explaining 55% to 85% of phenotype variance in each breed. The greyhound locus exhibiting the strongest association, located 150 kilobases upstream of the genes CDKN2A/B, is also the most rearranged locus in canine osteosarcoma tumors. The top germline candidate variant is found at a >90% frequency in Rottweilers and Irish wolfhounds, and alters an evolutionarily constrained element that we show has strong enhancer activity in human osteosarcoma cells. In all three breeds, osteosarcoma-associated loci and regions of reduced heterozygosity are enriched for genes in pathways connected to bone differentiation and growth. Several pathways, including one of genes regulated by miR124, are also enriched for somatic copy-number changes in tumors.

Conclusions

Mapping a complex cancer in multiple dog breeds reveals a polygenic spectrum of germline risk factors pointing to specific pathways as drivers of disease.     

Complete reannotation of the Arabidopsis genome: methods, tools, protocols and the final release

Background

Since the initial publication of its complete genome sequence, Arabidopsis thaliana has become more important than ever as a model for plant research. However, the initial genome annotation was submitted by multiple centers using inconsistent methods, making the data difficult to use for many applications.

Results

Over the course of three years, TIGR has completed its effort to standardize the structural and functional annotation of the Arabidopsis genome. Using both manual and automated methods, Arabidopsis gene structures were refined and gene products were renamed and assigned to Gene Ontology categories. We present an overview of the methods employed, tools developed, and protocols followed, summarizing the contents of each data release with special emphasis on our final annotation release (version 5).

Conclusion

Over the entire period, several thousand new genes and pseudogenes were added to the annotation. Approximately one third of the originally annotated gene models were significantly refined yielding improved gene structure annotations, and every protein-coding gene was manually inspected and classified using Gene Ontology terms.