Genome-Wide Survey of Pseudogenes in 80 Fully Re-sequenced Arabidopsis thaliana Accessions

Pseudogenes (Ψs), including processed and non-processed Ψs, are ubiquitous genetic elements derived from originally functional genes in all studied genomes within the three kingdoms of life. However, systematic surveys of non-processed Ψs utilizing genomic information from multiple samples within a species are still rare. Here a systematic comparative analysis was conducted of Ψs within 80 fully re-sequenced Arabidopsis thaliana accessions, and 7546 genes, representing ~28% of the genomic annotated open reading frames (ORFs), were found with disruptive mutations in at least one accession. The distribution of these Ψs on chromosomes showed a significantly negative correlation between Ψs/ORFs and their local gene densities, suggesting a higher proportion of Ψs in gene desert regions, e.g. near centromeres. On the other hand, compared with the non-Ψ loci, even the intact coding sequences (CDSs) in the Ψ loci were found to have shorter CDS length, fewer exon number and lower GC content. In addition, a significant functional bias against the null hypothesis was detected in the Ψs mainly involved in responses to environmental stimuli and biotic stress as reported, suggesting that they are likely important for adaptive evolution to rapidly changing environments by pseudogenization to accumulate successive mutations.

Analysis of gene order data supports vertical inheritance of the leukotoxin operon and genome rearrangements in the 5' flanking region in genus Mannheimia

BACKGROUND: The Mannheimia subclades belong to the same bacterial genus, but have taken divergent paths toward their distinct lifestyles. For example, M. haemolytica + M. glucosida are potential pathogens of the respiratory tract in the mammalian suborder Ruminantia, whereas M. ruminalis, the supposed sister group, lives as a commensal in the ovine rumen. We have tested the hypothesis that vertical inheritance of the leukotoxin (lktCABD) operon has occurred from the last common ancestor of genus Mannheimia to any ancestor of the diverging subclades by exploring gene order data. RESULTS: We examined the gene order in the 5' flanking region of the leukotoxin operon and found that the 5' flanking gene strings, hslVU-lapB-artJ-lktC and xylAB-lktC, are peculiar to M. haemolytica + M. glucosida and M. granulomatis, respectively, whereas the gene string hslVU-lapB-lktC is present in M. ruminalis, the supposed sister group of M. haemolytica + M. glucosida, and in the most ancient subclade M. varigena. In M. granulomatis, we found remnants of the gene string hslVU-lapB-lktC in the xylB-lktC intergenic region. CONCLUSION: These observations indicate that the gene string hslVU-lapB-lktC is more ancient than the hslVU-lapB-artJ-lktC and xylAB-lktC gene strings. The presence of (remnants of) the ancient gene string hslVU-lapB-lktC among any subclades within genus Mannheimia supports that it has been vertically inherited from the last common ancestor of genus Mannheimia to any ancestor of the diverging subclades, thus reaffirming the hypothesis of vertical inheritance of the leukotoxin operon. The presence of individual 5' flanking regions in M. haemolytica + M. glucosida and M. granulomatis reflects later genome rearrangements within each subclade. The evolution of the novel 5' flanking region in M. haemolytica + M. glucosida resulted in transcriptional coupling between the divergently arranged artJ and lkt promoters. We propose that the chimeric promoter have led to high level expression of the leukotoxin operon which could explain the increased potential of certain M. haemolytica + M. glucosida strains to cause a particular type of infection.

Increased genome instability in Escherichia coli lon mutants: relation to emergence of multiple-antibiotic-resistant (Mar) mutants caused by insertion sequence elements and large tandem genomic amplifications

Thirteen spontaneous multiple-antibiotic-resistant (Mar) mutants of Escherichia coli AG100 were isolated on Luria-Bertani (LB) agar in the presence of tetracycline (4 microg/ml). The phenotype was linked to insertion sequence (IS) insertions in marR or acrR or unstable large tandem genomic amplifications which included acrAB and which were bordered by IS3 or IS5 sequences. Five different lon mutations, not related to the Mar phenotype, were also found in 12 of the 13 mutants. Under specific selective conditions, most drug-resistant mutants appearing late on the selective plates evolved from a subpopulation of AG100 with lon mutations. That the lon locus was involved in the evolution to low levels of multidrug resistance was supported by the following findings: (i) AG100 grown in LB broth had an important spontaneous subpopulation (about 3.7x10(-4)) of lon::IS186 mutants, (ii) new lon mutants appeared during the selection on antibiotic-containing agar plates, (iii) lon mutants could slowly grow in the presence of low amounts (about 2x MIC of the wild type) of chloramphenicol or tetracycline, and (iv) a lon mutation conferred a mutator phenotype which increased IS transposition and genome rearrangements. The association between lon mutations and mutations causing the Mar phenotype was dependent on the medium (LB versus MacConkey medium) and the antibiotic used for the selection. A previously reported unstable amplifiable high-level resistance observed after the prolonged growth of Mar mutants in a low concentration of tetracycline or chloramphenicol can be explained by genomic amplification.

Nuclear and plastid markers reveal the persistence of genetic identity: A new perspective on the evolutionary history of Petunia exserta

Recently divergent species that can hybridize are ideal models for investigating the genetic exchanges that can occur while preserving the species boundaries. Petunia exserta is an endemic species from a very limited and specific area that grows exclusively in rocky shelters. These shaded spots are an inhospitable habitat for all other Petunia species, including the closely related and widely distributed species P. axillaris. Individuals with intermediate morphologic characteristics have been found near the rocky shelters and were believed to be putative hybrids between P. exserta and P. axillaris, suggesting a situation where Petunia exserta is losing its genetic identity. In the current study, we analyzed the plastid intergenic spacers trnS/trnG and trnH/psbA and six nuclear CAPS markers in a large sampling design of both species to understand the evolutionary process occurring in this biological system. Bayesian clustering methods, cpDNA haplotype networks, genetic diversity statistics, and coalescence-based analyses support a scenario where hybridization occurs while two genetic clusters corresponding to two species are maintained. Our results reinforce the importance of coupling differentially inherited markers with an extensive geographic sample to assess the evolutionary dynamics of recently diverged species that can hybridize. (C) 2013 Elsevier Inc. All rights reserved.

Molecular phylogenetic analyses identify Alpine differentiation and dysploid chromosome number changes as major forces for the evolution of the European endemic Phyteuma (Campanulaceae)

Phyteuma is a chromosomally and ecologically diverse vascular plant genus and constitutes an excellent system for studying both the role of chromosomal change for species diversification and the evolution of high-mountain biota. This kind of research is, however, hampered by the lack of a sound phylogenetic framework exacerbated by the notoriously low predictive power of traditional taxonomy with respect to phylogenetic relationships in Campanulaceae. Based on a comprehensive taxon sampling and analyses of nuclear and plastid sequence and AFLP fingerprint data, Phyteuma is confirmed as a monophyletic group sister to the monotypic Physoplexis, which is in line with their peculiar flower morphologies. Within Phyteuma two clades, largely corresponding to previously recognized sections, are consistently found. The traditional circumscription of taxonomic series is largely rejected. Whereas distinctness of the currently recognized species is mostly corroborated, some interspecific relationships remain ambiguous due to incongruences between nuclear and plastid data. Major forces for diversification and evolution of Phyteuma are descending dysploidy (i.e., a decrease in chromosome base number) as well as allopatric and ecological differentiation within the Alps, the genus' center of species diversity. (C) 2013 Elsevier Inc. All rights reserved.

The nuclear mRNA export receptor Mex67-Mtr2 of Trypanosoma brucei contains a unique and essential zinc finger motif.

Trypanosoma brucei is the causative agent of Human African Trypanosomiasis. Trypanosomes are early diverged protozoan parasites and show significant differences in their gene expression compared with higher eukaryotes. Due to a lack of individual gene promoters, large polycistronic transcripts are produced and individual mRNAs mature by trans-splicing and polyadenylation. In the absence of transcriptional control, regulation of gene expression occurs post-transcriptionally mainly by control of transcript stability and translation. Regulation of mRNA export from the nucleus to the cytoplasm might be an additional post-transcriptional event involved in gene regulation. However, our knowledge about mRNA export in trypanosomes is very limited. Although export factors of higher eukaryotes are reported to be conserved, only a few orthologues can be readily identified in the genome of T. brucei. Hence, biochemical approaches are needed to identify the export machinery of trypanosomes. Here, we report the functional characterization of the essential mRNA export factor TbMex67. TbMex67 contains a unique and essential N-terminal zinc finger motif. Furthermore, we could identify two interacting export factors namely TbMtr2 and the karyopherin TbIMP1. Our data show that the general heterodimeric export receptor Mex67-Mtr2 is conserved throughout the eukaryotic kingdom albeit exhibiting parasite-specific features.

Nuclear organization in the nematode C. elegans.

With its invariant cell lineage, easy genetics and small genome, the nematode Caenorhabditis elegans has emerged as one of the prime models in developmental biology over the last 50 years. Surprisingly however, until a decade ago very little was known about nuclear organization in worms, even though it is an ideal model system to explore the link between nuclear organization and cell fate determination. Here, we review the latest findings that exploit the repertoire of genetic tools developed in worms, leading to the identification of important sequences and signals governing the changes in chromatin tridimensional architecture. We also highlight parallels and differences to other model systems.

Building silent compartments at the nuclear periphery: a recurrent theme.

In eukaryotes, the genetic material is stored in the nucleus, which is enclosed in a double lipid bilayer, the nuclear envelope (NE). It protects the genome from physical stress and separates it from the rest of the cell. On top of this physical function, growing evidence shows that the nuclear periphery contributes to the 3D organization of the genome. In turn, tridimensional organization of chromatin in the nuclear space influences genome expression. Here we review recent findings on the function of this physical barrier in gene repression and latest models on how silent subnuclear compartments at the NE are built in yeast as well as in the nematode C. elegans and mammalian cells; trying to draw parallels between the three systems.

The evolution of host associations in the parasitic wasp genus Ichneumon (Hymenoptera: Ichneumonidae): convergent adaptations to host pupation sites

Background: The diversification of organisms with a parasitic lifestyle is often tightly linked to the evolution of their host associations. If a tight host association exists, closely related species tend to attack closely related hosts; host associations are less stable if associations are determined by more plastic traits like parasitoid searching and oviposition behaviour. The pupal-parasitoids of the genus Ichneumon attack a variety of macrolepidopteran hosts.They are either monophagous or polyphagous, and therefore offer a promissing system to investigate the evolution of host associations. Ichneumon was previously divided into two groups based on general body shape; however, a stout shape has been suggested as an adaptation to buried host pupation sites, and might thus not represent a reliable phylogenetic character. Results: We here reconstruct the first molecular phylogeny of the genus Ichneumon using two mitochondrial (CO1 and NADH1) and one nuclear marker (28S). The resulting phylogeny only supports monophyly of Ichneumon when Ichneumon lugens Gravenhorst, 1829 (formerly in Chasmias, stat. rev.) and Ichneumon deliratorius Linnaeus, 1758 (formerly Coelichneumon) are included. Neither parasitoid species that attack hosts belonging to one family nor those attacking butterflies (Rhopalocera) form monophyletic clades. Ancestral state reconstructions suggest multiple transitions between searching for hosts above versus below ground and between a stout versus elongated body shape. A model assuming correlated evolution between the two characters was preferred over independent evolution of host-searching niche and body shape. Conclusions: Host relations, both in terms of phylogeny and ecology, evolved at a high pace in the genus Ichneumon. Numerous switches between hosts of different lepidopteran families have occurred, a pattern that seems to be the rule among idiobiont parasitoids. A stout body and antennal shape in the parasitoid female is confirmed as an ecological adaptation to host pupation sites below ground and has evolved convergently several times. Morphological characters that might be involved in adaptation to hosts should be avoided as diagnostic characters for phylogeny and classification, as they can be expected to show high levels of homoplasy.

The Oxytricha trifallax macronuclear genome: a complex eukaryotic genome with 16,000 tiny chromosomes.

The macronuclear genome of the ciliate Oxytricha trifallax displays an extreme and unique eukaryotic genome architecture with extensive genomic variation. During sexual genome development, the expressed, somatic macronuclear genome is whittled down to the genic portion of a small fraction (∼5%) of its precursor "silent" germline micronuclear genome by a process of "unscrambling" and fragmentation. The tiny macronuclear "nanochromosomes" typically encode single, protein-coding genes (a small portion, 10%, encode 2-8 genes), have minimal noncoding regions, and are differentially amplified to an average of ∼2,000 copies. We report the high-quality genome assembly of ∼16,000 complete nanochromosomes (∼50 Mb haploid genome size) that vary from 469 bp to 66 kb long (mean ∼3.2 kb) and encode ∼18,500 genes. Alternative DNA fragmentation processes ∼10% of the nanochromosomes into multiple isoforms that usually encode complete genes. Nucleotide diversity in the macronucleus is very high (SNP heterozygosity is ∼4.0%), suggesting that Oxytricha trifallax may have one of the largest known effective population sizes of eukaryotes. Comparison to other ciliates with nonscrambled genomes and long macronuclear chromosomes (on the order of 100 kb) suggests several candidate proteins that could be involved in genome rearrangement, including domesticated MULE and IS1595-like DDE transposases. The assembly of the highly fragmented Oxytricha macronuclear genome is the first completed genome with such an unusual architecture. This genome sequence provides tantalizing glimpses into novel molecular biology and evolution. For example, Oxytricha maintains tens of millions of telomeres per cell and has also evolved an intriguing expansion of telomere end-binding proteins. In conjunction with the micronuclear genome in progress, the O. trifallax macronuclear genome will provide an invaluable resource for investigating programmed genome rearrangements, complementing studies of rearrangements arising during evolution and disease.

Modern techniques for the analysis of chromatin and nuclear organization in C. elegans.

In recent years, Caenorhabditis elegans has emerged as a new model to investigate the relationships between nuclear architecture, cellular differentiation, and organismal development. On one hand, C. elegans with its fixed lineage and transparent body is a great model organism to observe gene functions in vivo in specific cell types using microscopy. On the other hand, two different techniques have been applied in nematodes to identify binding sites for chromatin-associated proteins genome-wide: chromatin immunoprecipitation (ChIP), and Dam-mediated identification (DamID). We summarize here all three techniques together as they are complementary. We also highlight strengths and differences of the individual approaches.

Trypanosoma brucei RRM1 is a nuclear RNA-binding protein and modulator of chromatin structure

TbRRM1 of Trypanosoma brucei is a nucleoprotein that was previously identified in a search for splicing factors in T. brucei. We show that TbRRM1 associates with mRNAs and with the auxiliary splicing factor polypyrimidine tract-binding protein 2, but not with components of the core spliceosome. TbRRM1 also interacts with several retrotransposon hot spot (RHS) proteins and histones. RNA immunoprecipitation of a tagged form of TbRRM1 from procyclic (insect) form trypanosomes identified ca. 1,500 transcripts that were enriched and 3,000 transcripts that were underrepresented compared to cellular mRNA. Enriched transcripts encoded RNA-binding proteins, including TbRRM1 itself, several RHS transcripts, mRNAs with long coding regions, and a high proportion of stage-regulated mRNAs that are more highly expressed in bloodstream forms. Transcripts encoding ribosomal proteins, other factors involved in translation, and procyclic-specific transcripts were underrepresented. Knockdown of TbRRM1 by RNA interference caused widespread changes in mRNA abundance, but these changes did not correlate with the binding of the protein to transcripts, and most splice sites were unchanged, negating a general role for TbRRM1 in splice site selection. When changes in mRNA abundance were mapped across the genome, regions with many downregulated mRNAs were identified. Two regions were analyzed by chromatin immunoprecipitation, both of which exhibited increases in nucleosome occupancy upon TbRRM1 depletion. In addition, subjecting cells to heat shock resulted in translocation of TbRRM1 to the cytoplasm and compaction of chromatin, consistent with a second role for TbRRM1 in modulating chromatin structure. IMPORTANCE: Trypanosoma brucei, the parasite that causes human sleeping sickness, is transmitted by tsetse flies. The parasite progresses through different life cycle stages in its two hosts, altering its pattern of gene expression in the process. In trypanosomes, protein-coding genes are organized as polycistronic units that are processed into monocistronic mRNAs. Since genes in the same unit can be regulated independently of each other, it is believed that gene regulation is essentially posttranscriptional. In this study, we investigated the role of a nuclear RNA-binding protein, TbRRM1, in the insect stage of the parasite. We found that TbRRM1 binds nuclear mRNAs and also affects chromatin status. Reduction of nuclear TbRRM1 by RNA interference or heat shock resulted in chromatin compaction. We propose that TbRRM1 regulates RNA polymerase II-driven gene expression both cotranscriptionally, by facilitating transcription and efficient splicing, and posttranscriptionally, via its interaction with nuclear mRNAs.

Genome and transcriptome sequencing identifies breeding targets in the orphan crop tef (Eragrostis tef)

Background: Tef (Eragrostis tef), an indigenous cereal critical to food security in the Horn of Africa, is rich in minerals and protein, resistant to many biotic and abiotic stresses and safe for diabetics as well as sufferers of immune reactions to wheat gluten. We present the genome of tef, the first species in the grass subfamily Chloridoideae and the first allotetraploid assembled de novo. We sequenced the tef genome for marker-assisted breeding, to shed light on the molecular mechanisms conferring tef's desirable nutritional and agronomic properties, and to make its genome publicly available as a community resource. Results: The draft genome contains 672 Mbp representing 87% of the genome size estimated from flow cytometry. We also sequenced two transcriptomes, one from a normalized RNA library and another from unnormalized RNASeq data. The normalized RNA library revealed around 38000 transcripts that were then annotated by the SwissProt group. The CoGe comparative genomics platform was used to compare the tef genome to other genomes, notably sorghum. Scaffolds comprising approximately half of the genome size were ordered by syntenic alignment to sorghum producing tef pseudo-chromosomes, which were sorted into A and B genomes as well as compared to the genetic map of tef. The draft genome was used to identify novel SSR markers, investigate target genes for abiotic stress resistance studies, and understand the evolution of the prolamin family of proteins that are responsible for the immune response to gluten. Conclusions: It is highly plausible that breeding targets previously identified in other cereal crops will also be valuable breeding targets in tef. The draft genome and transcriptome will be of great use for identifying these targets for genetic improvement of this orphan crop that is vital for feeding 50 million people in the Horn of Africa.