OMA orthology in 2024: improved prokaryote coverage, ancestral and extant GO enrichment, a revamped synteny viewer and more in the OMA Ecosystem

In this update paper, we present the latest developments in the OMA browser knowledgebase, which aims to provide high-quality orthology inferences and facilitate the study of gene families, genomes and their evolution. First, we discuss the addition of new species in the database, particularly an expanded representation of prokaryotic species. The OMA browser now offers Ancestral Genome pages and an Ancestral Gene Order viewer, allowing users to explore the evolutionary history and gene content of ancestral genomes. We also introduce a revamped Local Synteny Viewer to compare genomic neighborhoods across both extant and ancestral genomes. Hierarchical Orthologous Groups (HOGs) are now annotated with Gene Ontology annotations, and users can easily perform extant or ancestral GO enrichments. Finally, we recap new tools in the OMA Ecosystem, including OMAmer for proteome mapping, OMArk for proteome quality assessment, OMAMO for model organism selection and Read2Tree for phylogenetic species tree construction from reads. These new features provide exciting opportunities for orthology analysis and comparative genomics. OMA is accessible at https://omabrowser.org.

The First Two Complete Mitochondrial Genomes of Neoephemeridae (Ephemeroptera): Comparative Analysis and Phylogenetic Implication for Furcatergalia

Mayflies of the family Neoephemeridae are widespread in the Holarctic and Oriental regions, and its phylogenetic position is still unstable in the group Furcatergalia (mayflies with fringed gills). In the present study, we determined the complete mitogenomes of two species, namely Potamanthellus edmundsi and Pulchephemera projecta, of this family. The lengths of two mitogenomes were 15,274 bp and 16,031 bp with an A + T content of 73.38% and 73.07%, respectively. Two neoephemerid mitogenomes had a similar gene size, base composition, and codon usage of protein-coding genes (PCGs), and the sequenced gene arrangements were consistent with the putative ancestral insect mitogenomes as understood today. The most variable gene of Furcatergalia mitogenomes was ND2, while the most conserved gene was COI. Meanwhile, the analysis of selection pressures showed that ND6 and ATP8 exhibited a relaxed purifying selection, and COI was under the strongest purifying selection. Phylogenetic trees reconstructed based on two concatenated nucleotide datasets using both maximum likelihood (ML) and Bayesian inference (BI) estimations yielded robust identical topologies. These results corroborated the monophyly of seven studied families and supported the family Leptophlebiidae as being of the basal lineage of Furcatergalia. Additionally, the sister-group relationship of Caenidae and Neoephemeridae was well supported. Methodologically, our present study provides a general reference for future phylogenetic studies of Ephemeroptera at the mitogenome level.

Sorting Signed Permutations by Inverse Tandem Duplication Random Losses

Gene order evolution of unichromosomal genomes, for example mitochondrial genomes, has been modelled mostly by four major types of genome rearrangements: inversions, transpositions, inverse transpositions, and tandem duplication random losses. Generalizing models that include all those rearrangements while admitting computational tractability are rare. In this paper, we study such a rearrangement model, namely the inverse tandem duplication random loss (iTDRL) model, where an iTDRL duplicates and inverts a continuous segment of a gene order followed by the random loss of one of the redundant copies of each gene. The iTDRL rearrangement has currently been proposed by several authors suggesting it to be a possible mechanisms of mitochondrial gene order evolution. We initiate the algorithmic study of this new model of genome rearrangement by proving that a shortest rearrangement scenario that transforms one given gene order into another given gene order can be obtained in quasilinear time. Furthermore, we show that the length of such a scenario, i.e., the minimum number of iTDRLs in the transformation, can be computed in linear time.

Complete mitogenome of endemic plum-headed parakeet Psittacula cyanocephala - characterization and phylogenetic analysis

Psittacula cyanocephala is an endemic parakeet from the Indian sub-continent that is widespread in the illegal bird trade. Previous studies on Psittacula parakeets have highlighted taxonomic ambiguities, warranting studies to resolve the issues. Since the mitochondrial genome provides useful information concerning the species evolution and phylogenetics, we sequenced the complete mitogenome of P. cyanocephala using NGS, validated 38.86% of the mitogenome using Sanger Sequencing and compared it with other available whole mitogenomes of Psittacula. The complete mitogenome of the species was 16814 bp in length with 54.08% AT composition. P. cyanocephala mitogenome comprises of 13 protein-coding genes, 2 rRNAs and 22 tRNAs. P. cyanocephala mitogenome organization was consistent with other Psittacula mitogenomes. Comparative codon usage analysis indicated the role of natural selection on Psittacula mitogenomes. Strong purifying selection pressure was observed maximum on nad1 and nad4l genes. The mitochondrial control region of all Psittacula species displayed the ancestral avian CR gene order. Phylogenetic analyses revealed the Psittacula genus as paraphyletic nature, containing at least 4 groups of species within the same genus, suggesting its taxonomic reconsideration. Our results provide useful information for developing forensic tests to control the illegal trade of the species and scientific basis for phylogenetic revision of the genus Psittacula.

OMA orthology in 2021: website overhaul, conserved isoforms, ancestral gene order and more

OMA is an established resource to elucidate evolutionary relationships among genes from currently 2326 genomes covering all domains of life. OMA provides pairwise and groupwise orthologs, functional annotations, local and global gene order conservation (synteny) information, among many other functions. This update paper describes the reorganisation of the database into gene-, group- and genome-centric pages. Other new and improved features are detailed, such as reporting of the evolutionarily best conserved isoforms of alternatively spliced genes, the inferred local order of ancestral genes, phylogenetic profiling, better cross-references, fast genome mapping, semantic data sharing via RDF, as well as a special coronavirus OMA with 119 viruses from the Nidovirales order, including SARS-CoV-2, the agent of the COVID-19 pandemic. We conclude with improvements to the documentation of the resource through primers, tutorials and short videos. OMA is accessible at https://omabrowser.org.

Close arrangement of CARK3 and PMEIL affects ABA-mediated pollen sterility in Arabidopsis thaliana

Abscisic acid (ABA) signaling is a vital plant signaling pathway for plant responses to stress conditions. ABA treatment can alter global gene expression patterns and cause significant phenotypic changes. We investigated the responses to ABA treatment during flowering in Arabidopsis thaliana. Dipping the flowers of CARK3 T-DNA mutants in ABA solution, led to less reduction of pollen fertility than in the wild type plants (Col-0). We demonstrated that PMEIL, a gene located downstream of CARK3, directly affects pollen fertility. Due to the close arrangement of CARK3 and PMEIL, CARK3 expression represses transcription of PMEIL in an ABA-dependent manner through transcriptional interference. Our study uncovers a molecular mechanism underlying ABA-mediated pollen sterility and provides an example of how transcriptional interference caused by close arrangement of genes may mediate stress responses during plant reproduction.

The complete mitochondrial genome of medicinal fungus Taiwanofungus camphoratus reveals gene rearrangements and intron dynamics of Polyporales

Taiwanofungus camphoratus is a highly valued medicinal mushroom that is endemic to Taiwan, China. In the present study, the mitogenome of T. camphoratus was assembled and compared with other published Polyporales mitogenomes. The T. camphoratus mitogenome was composed of circular DNA molecules, with a total size of 114,922 bp. Genome collinearity analysis revealed large-scale gene rearrangements between the mitogenomes of Polyporales, and T. camphoratus contained a unique gene order. The number and classes of introns were highly variable in 12 Polyporales species we examined, which proved that numerous intron loss or gain events occurred in the evolution of Polyporales. The Ka/Ks values for most core protein coding genes in Polyporales species were less than 1, indicating that these genes were subject to purifying selection. However, the rps3 gene was found under positive or relaxed selection between some Polyporales species. Phylogenetic analysis based on the combined mitochondrial gene set obtained a well-supported topology, and T. camphoratus was identified as a sister species to Laetiporus sulphureus. This study served as the first report on the mitogenome in the Taiwanofungus genus, which will provide a basis for understanding the phylogeny and evolution of this important fungus.

The SNAP hypothesis: Chromosomal rearrangements could emerge from positive Selection during Niche Adaptation

The relative linear order of most genes on bacterial chromosomes is not conserved over evolutionary timescales. One explanation is that selection is weak, allowing recombination to randomize gene order by genetic drift. However, most chromosomal rearrangements are deleterious to fitness. In contrast, we propose the hypothesis that rearrangements in gene order are more likely the result of selection during niche adaptation (SNAP). Partial chromosomal duplications occur very frequently by recombination between direct repeat sequences. Duplicated regions may contain tens to hundreds of genes and segregate quickly unless maintained by selection. Bacteria exposed to non-lethal selections (for example, a requirement to grow on a poor nutrient) can adapt by maintaining a duplication that includes a gene that improves relative fitness. Further improvements in fitness result from the loss or inactivation of non-selected genes within each copy of the duplication. When genes that are essential in single copy are lost from different copies of the duplication, segregation is prevented even if the original selection is lifted. Functional gene loss continues until a new genetic equilibrium is reached. The outcome is a rearranged gene order. Mathematical modelling shows that this process of positive selection to adapt to a new niche can rapidly drive rearrangements in gene order to fixation. Signature features (duplication formation and divergence) of the SNAP model were identified in natural isolates from multiple species showing that the initial two steps in the SNAP process can occur with a remarkably high frequency. Further bioinformatic and experimental analyses are required to test if and to which extend the SNAP process acts on bacterial genomes.

All life on earth has evolved from a universal common ancestor with a specific order of genes on the chromosome. This order is not maintained in modern species and the standard hypothesis is that changes reflect a lack of strong selection on gene order. Here, we propose an alternative hypothesis, SNAP. The occupation of a novel environment by bacteria is generally a trade-off situation. For example, while the bacteria may not be adapted to grow well under the new conditions, they may benefit by not having to share available resources with other microorganisms. Bacterial populations frequently acquire duplications of chromosomal segments containing genes that can help them adapt to a new environment. Other genes that are also duplicated are not required in two copies so that over time a superfluous copy can be lost. Eventually, the process of duplication and gene loss can lead to the rearrangement of the gene order in the chromosomal segment. The major benefit of this model over the standard hypothesis is that the process is driven by positive selection and can reach fixation rapidly.

Extensive intraspecific gene order and gene structural variations in upland cotton cultivars

Multiple cotton genomes (diploid and tetraploid) have been assembled. However, genomic variations between cultivars of allotetraploid upland cotton (Gossypium hirsutum L.), the most widely planted cotton species in the world, remain unexplored. Here, we use single-molecule long read and Hi-C sequencing technologies to assemble genomes of the two upland cotton cultivars TM-1 and zhongmiansuo24 (ZM24). Comparisons among TM-1 and ZM24 assemblies and the genomes of the diploid ancestors reveal a large amount of genetic variations. Among them, the top three longest structural variations are located on chromosome A08 of the tetraploid upland cotton, which account for ~30% total length of this chromosome. Haplotype analyses of the mapping population derived from these two cultivars and the germplasm panel show suppressed recombination rates in this region. This study provides additional genomic resources for the community, and the identified genetic variations, especially the reduced meiotic recombination on chromosome A08, will help future breeding.

Comparative Chloroplast Genomes of Sorghum Species: Sequence Divergence and Phylogenetic Relationships

Sorghum comprises 31 species that exhibit considerable morphological and ecological diversity. The phylogenetic relationships among Sorghum species still remain unresolved due to lower information on the traditional DNA markers, which provides a limited resolution for identifying Sorghum species. In this study, we sequenced the complete chloroplast genomes of Sorghum sudanense and S. propinquum and analyzed the published chloroplast genomes of S. bicolor and S. timorense to retrieve valuable chloroplast molecular resources for Sorghum. The chloroplast genomes ranged in length from 140,629 to 140,755 bp, and their gene contents, gene orders, and GC contents were similar to those for other Poaceae species but were slightly different in the number of SSRs. Comparative analyses among the four chloroplast genomes revealed 651 variable sites, 137 indels, and nine small inversions. Four highly divergent DNA regions (rps16-trnQ, trnG-trnM, rbcL-psaI, and rps15-ndhF), which were suitable for phylogenetic and species identification, were detected in the Sorghum chloroplast genomes. A phylogenetic analysis strongly supported that Sorghum is a monophyletic group in the tribe Andropogoneae. Overall, the genomic resources in this study could provide potential molecular markers for phylogeny and species identification in Sorghum.

Unusual dicistronic expression from closely spaced initiation codons in an umbravirus subgenomic RNA

Translation commencing at closely spaced initiation codons is common in RNA viruses with limited genome space. In the subgenomic RNA (sgRNA) of Pea enation mosaic virus 2, two closely spaced, out-of-frame start codons direct synthesis of movement/stability proteins p26 and p27. Efficient translation from AUG26/AUG27 is dependent on three 3'-proximal cap-independent translation enhancers (3'CITEs), whereas translation of the genomic (gRNA) requires only two. Contrary to strictly scanning-dependent initiation at the gRNA, sequence context of AUG26/AUG27 does not conform with Kozak requirements and insertion of efficient upstream AUGs had pronounced effects for AUG26 but only moderate effects for AUG27. Insertion of a hairpin within an extended 5' UTR did not significantly impact translation from AUG26/AUG27. Furthermore, AUG27 repressed translation from upstream AUG26 and this effect was mitigated when inter-codon spacing was reduced. Addition of a stable hairpin to the very 5' end of the sgRNA severely restricted translation, testifying that this 3'CITE-driven initiation is 5' end-dependent. Similar to gRNA, sgRNA reporter transcripts were nearly exclusively associated with light polysomes and 3'CITE-promoted long-distance interaction connecting the sgRNA ends affected the number of templates translated and not the initiation rate. We propose a non-canonical, 3'CITE-driven mechanism for efficient dicistronic expression from umbravirus sgRNAs.

The complete mitochondrial genome of the Caribbean spiny lobster Panulirus argus

Panulirus argus is a keystone species and target of the most lucrative fishery in the Caribbean region. This study reports, for the first time, the complete mitochondrial genome of Panulirus argus (average coverage depth nucleotide-1 = 70×) assembled from short Illumina 150 bp PE reads. The AT-rich mitochondrial genome of Panulirus argus was 15 739 bp in length and comprised 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, and 22 transfer RNA genes. A single 801 bp long intergenic space was assumed to be the D-loop. Most of the PCGs were encoded on the H-strand. The gene order observed in the mitochondrial genome of Panulirus argus corresponds to the presumed Pancrustacean ground pattern. KA/KS ratios calculated for all mitochondrial PCGs showed values < 1, indicating that all these PCGs are evolving under purifying selection. A maximum likelihood phylogenetic analysis (concatenated PCGs [n = 13], 154 arthropods) supported the monophyly of the Achelata and other infraorders within the Decapoda. Mitochondrial PCGs have enough phylogenetic informativeness to explore high-level genealogical relationships in the Pancrustacea. The complete mitochondrial genome of the Caribbean spiny lobster Panulirus argus will contribute to the better understanding of meta-population connectivity in this keystone overexploited species.

Sequencing of organellar genomes of Gymnomitrion concinnatum (Jungermanniales) revealed the first exception in the structure and gene order of evolutionary stable liverworts mitogenomes

BACKGROUND

Comparative analyses of chloroplast and mitochondrial genomes have shown that organelle genomes in bryophytes evolve slowly. However, in contrast to seed plants, the organellar genomes are yet poorly explored in bryophytes, especially among liverworts. Discovering another organellar genomes of liverwort species by sequencing provides new conclusions on evolution of bryophytes.

RESULTS

In this work, the organellar genomes of Gymnomitrion concinnatum liverwort were sequenced, assembled and annotated for the first time. The chloroplast genome displays, typical for most plants, quadripartite structure containing large single copy region (81,701 bp), two inverted repeat regions (8704 bp each) and small single copy region (20,179 bp). The gene order and content of chloroplast are very similar to other liverworts with minor differences observed. A total number of 739 and 222 RNA editing sites were predicted in chloroplast and mitochondrial genes of G. concinnatum. The mitochondrial genome gene content is also in accordance with liverworts except few alterations such as: intron loss in cox1 and atp1 genes. Nonetheless the analysis revealed that G. concinnatum mitogenome structure and gene order are rearranged in comparison with other mitogenomes of liverworts. The causes underlying such mitogenomic rearrangement were investigated and the probable model of recombination was proposed.

CONCLUSIONS

This study provide the overview of mitochondrial and chloroplast genome structure and gene order diversity of Gymnomitrion concinnatum against the background of known organellar genomes of liverworts. The obtained results cast doubt on the idea that mitogenome structure of early land plants is highly conserved as previous studies suggested. In fact is the very first case of recombination within, evolutionary stable, mitogenomes of liverworts.

Conservation of mitochondrial genome arrangements in brittle stars (Echinodermata, Ophiuroidea)

Brittle stars are conspicuous members of benthic ecosystems, fill many ecological niches and are the most speciose of all classes of echinoderms. With high levels of biodiversity, elucidating the evolutionary history of this group is important. Understanding of higher-level relationships within Ophiuroidea has been aided by multilocus nuclear data and DNA barcoding. However, the degree of consistency between mitochondrial and nuclear data within ophiuroids remains unclear and deserves further assessment. In this study, 17 mitochondrial genomes spanning the taxonomic breadth of Ophiuroidea were utilized to explore evolutionary relationships through maximum likelihood analyses, Bayesian inference and comparative assessment of gene order. Our phylogenetic analyses, based on both nucleotide and amino acid residues, support recent findings based on multilocus nuclear data and morphology, in that the brittle star clades Ophintegrida and Euryophiurida were recovered as monophyletic with the latter comprising Euyalida, Ophiuridae and Ophiopyrgidae. Only three different arrangements of the 13 protein coding and 2 ribosomal RNA genes were observed. As expected, tRNA genes were more likely to have undergone rearrangement but the order of all 37 genes was found to be conserved in all sampled Euryalida and Ophiuridae. Both Euryalida and the clade comprised of Ophiuridae and Ophiopyrgidae, each had their own conserved rearrangement of protein coding genes and ribosomal genes, after divergence from their last common ancestor. Euryalida has a rearrangement of the two ribosomal RNA genes, rrnS and rrnL, in contrast to Ophiuridae and Ophiopyrgidae, which had an inversion of the genes nad1, nad2, and cob relative to Ophintegrida. Further, our data support the gene order found in all sampled Euryalida as the most likely ancestral order for all Ophiuroidea.

An exception among diatoms: unique organization of genes involved in isoprenoid biosynthesis in Rhizosolenia setigera CCMP 1694

The marine diatom Rhizosolenia setigera is unique among this group of microalgae given that it is only one of a handful of diatom species that can produce highly branched isoprenoid (HBI) hydrocarbons. In our efforts to determine distinguishing molecular characteristics in R. setigera CCMP 1694 that could help elucidate the underlying mechanisms for its ability to biosynthesize HBIs, we discovered the occurrence of independent genes encoding for two isopentenyl diphosphate isomerases (RsIDI1 and RsIDI2) and one squalene synthase (RsSQS), enzymes that catalyze non-consecutive steps in isoprenoid biosynthesis. These genes are peculiarly fused in all other genome-sequenced diatoms to date, making their organization in R. setigera CCMP 1694 a clear distinguishing molecular feature. Phylogenetic and sequence analysis of RsIDI1, RsIDI2, and RsSQS revealed that such an arrangement of individually transcribed genes involved in isoprenoid biosynthesis could have arisen through a secondary gene fission event. We further demonstrate that inhibition of squalene synthase (SQS) shifts the flux of exogenous isoprenoid precursors towards HBI biosynthesis suggesting the competition for isoprenoid substrates in the form of farnesyl diphosphate between the sterol and HBI biosynthetic pathways in this diatom.

Sequencing flow-sorted short arm of Haynaldia villosa chromosome 4V provides insights into its molecular structure and virtual gene order

BACKGROUND

Haynaldia villosa (H. villosa) has been recognized as a species potentially useful for wheat improvement. The availability of its genomic sequences will boost its research and application.

RESULTS

In this work, the short arm of H. villosa chromosome 4V (4VS) was sorted by flow cytometry and sequenced using Illumina platform. About 170.6 Mb assembled sequences were obtained. Further analysis showed that repetitive elements accounted for about 64.6% of 4VS, while the coding fraction, which is corresponding to 1977 annotated genes, represented 1.5% of the arm. The syntenic regions of the 4VS were searched and identified on wheat group 4 chromosomes 4AL, 4BS, 4DS, Brachypodium chromosomes 1 and 4, rice chromosomes 3 and 11, and sorghum chromosomes 1, 5 and 8. Based on genome-zipper analysis, a virtual gene order comprising 735 gene loci on 4VS genome was built by referring to the Brachypodium genome, which was relatively consistent with the scaffold order determined for Ae. tauschii chromosome 4D. The homologous alleles of several cloned genes on wheat group 4 chromosomes including Rht-1 gene were identified.

CONCLUSIONS

The sequences provided valuable information for mapping and positional-cloning genes located on 4VS, such as the wheat yellow mosaic virus resistance gene Wss1. The work on 4VS provided detailed insights into the genome of H. villosa, and may also serve as a model for sequencing the remaining parts of H. villosa genome.

Characterization of SNP and Structural Variations in the Mitochondrial Genomes of Tilletia indica and Its Closely Related Species Formed Basis for a Simple Diagnostic Assay

Tilletia indica causes the disease Karnal bunt in wheat. The disease is under international quarantine regulations. Comparative mitochondrial (mt) genome analysis of T. indica (KX394364 and DQ993184) and T. walkeri (EF536375) has found 325 to 328 SNPs, 57 to 60 short InDels (from 1 to 13 nt), two InDels (30 and 61 nt) and five (>200 nt) presence/absence variations (PAVs) between the two species. The mt genomes of both species have identical gene order. The numbers of SNPs and InDels between the mt genomes of the two species are approximately nine times of the corresponding numbers between the two T. indica isolates. There are eight SNPs between T. indica and T. walkeri that resulted in amino acid substitutions in the mt genes of cob, nad2 and nad5. In contrast, there is no amino acid substitution in the mt genes of the T. indica isolates from the SNPs found. The five PAVs present in T. indica (DQ993184) are absent in T. walkeri. Four PAVs are more than 1 kb and are not present in every T. indica isolate. Analysis of their presence and absence separates a collection of T. indica isolates into 11 subgroups. Two PAVs have ORFs for the LAGLIDAG endonuclease and two have ORFs for the GIY-YIG endonuclease family, which are representatives of homing endonuclease genes (HEGs). These intron- encoded HEGs confer intron mobility and account for their fluid distribution in T. indica isolates. The small PAV of 221 bp, present in every T. indica isolate and unique to the species, was used as the genetic fingerprint for the successful development of a rapid, highly sensitive and specific loop mediated isothermal amplification (LAMP) assay. The simple procedure of the LAMP assay and the easy detection formats will enable the assay to be automated for high throughput diagnosis.

Complete mitochondrial genome of the brown alga Sargassum vachellianum (Sargassaceae, Phaeophyceae)

Sargassum vachellianum C. Agardh is endemic to China. It inhabits in rocky intertidal zones and plays an important role in maintaining the structure and function of littoral ecosystems. In this study, we present the complete mitochondrial genome of S. vachellianum. The circular S. vachellianum mitogenome is 34,877 bp in size and contains the same set of 65 genes as the reported Sargassum mtDNAs. The overall AT content of the genome is 63.79%, and the inter-genic spacers constitute only 4.67%. The genome organization including the gene order, overlapping regions between genes, and the total length of inter-genic spacers is conserved among the known Sargassum mitogenomes. High divergence is in inter-genic spacer regions. Phylogenetic analyses indicated that S. vachellianum combined tightly with Sargassum species with strong support values (NJ/ML, 100%).

Comparative mitogenomic analysis of the superfamily Pentatomoidea (Insecta: Hemiptera: Heteroptera) and phylogenetic implications

BACKGROUND

Insect mitochondrial genomes (mitogenomes) are the most extensively used genetic marker for evolutionary and population genetics studies of insects. The Pentatomoidea superfamily is economically important and the largest superfamily within Pentatomomorpha with over 7,000 species. To better understand the diversity and evolution of pentatomoid species, we sequenced and annotated the mitogenomes of Eurydema gebleri and Rubiconia intermedia, and present the first comparative analysis of the 11 pentatomoid mitogenomes that have been sequenced to date.

RESULTS

We obtained the complete mitogenome of Eurydema gebleri (16,005 bp) and a nearly complete mitogenome of Rubiconia intermedia (14,967 bp). Our results show that gene content, gene arrangement, base composition, codon usage, and mitochondrial transcription termination factor sequences are highly conserved in pentatomoid species, especially for species in the same family. Evolutionary rate analyses of protein-coding genes reveal that the highest and lowest rates are found in atp8 and cox1 and distinctive evolutionary patterns are significantly correlated with the G + C content of genes. We inferred the secondary structures for two rRNA genes for eleven pentatomoid species, and identify some conserved motifs of RNA structures in Pentatomidea. All tRNA genes in pentatomoid mitogenomes have a canonical cloverleaf secondary structure, except for two tRNAs (trnS1 and trnV) which appear to lack the dihydrouridine arm. Regions that are A + T-rich have several distinct characteristics (e.g. size variation and abundant tandem repeats), and have potential as species or population level molecular markers. Phylogenetic analyses based on mitogenomic data strongly support the monophyly of Pentatomoidea, and the estimated phylogenetic relationships are: (Urostylididae + (Plataspidae + (Pentatomidae + (Cydnidae + (Dinidoridae + Tessaratomidae))))).

CONCLUSIONS

This comparative mitogenomic analysis sheds light on the architecture and evolution of mitogenomes in the superfamily Pentatomoidea. Mitogenomes can be effectively used to resolve phylogenetic relationships of pentatomomorphan insects at various taxonomic levels. Sequencing more mitogenomes at various taxonomic levels, particularly from closely related species, will improve the annotation accuracy of mitochondrial genes, as well as greatly enhance our understanding of mitogenomic evolution and phylogenetic relationships in pentatomoids.

Whole genome analyses of marine fish pathogenic isolate, Mycobacterium sp. 012931

Mycobacterium is a genus within the order Actinomycetales that comprises of a large number of well-characterized species, several of which includes pathogens known to cause serious disease in human and animal. Here, we report the whole genome sequence of Mycobacterium sp. strain 012931 isolated from the marine fish, yellowtail (Seriola quinqueradiata). Mycobacterium sp. 012931 is a fish pathogen causing serious damage to aquaculture farms in Japan. DNA dot plot analysis showed that Mycobacterium sp. 012931 was more closely related to Mycobacterium marinum when compared across several Mycobacterium species. However, little conservation of the gene order was observed between Mycobacterium sp. 012931 and M. marinum genome. The annotated 5,464 genes of Mycobacterium sp. 012931 was classified into 26 subsystems. The insertion/deletion gene analysis shows Mycobacterium sp. 012931 had 643 unique genes that were not found in the M. marinum strains. In the virulence, disease, and defense subsystem, both insertion and deletion genes of Mycobacterium sp. 012931 were associated with the PPE gene cluster of Mycobacteria. Of seven plcB genes in Mycobacterium sp. 012931, plcB_2 and plcB_3 showed low identities with those of M. marinum strains. Therefore, Mycobacterium sp. 012931 has differences on genetic and virulence from M. marinum and may induce different interaction mechanisms between host and pathogen.

Genome-wide analysis of trans-splicing in the nematode Pristionchus pacificus unravels conserved gene functions for germline and dauer development in divergent operons

Discovery of trans-splicing in multiple metazoan lineages led to the identification of operon-like gene organization in diverse organisms, including trypanosomes, tunicates, and nematodes, but the functional significance of such operons is not completely understood. To see whether the content or organization of operons serves similar roles across species, we experimentally defined operons in the nematode model Pristionchus pacificus. We performed affinity capture experiments on mRNA pools to specifically enrich for transcripts that are trans-spliced to either the SL1- or SL2-spliced leader, using spliced leader-specific probes. We obtained distinct trans-splicing patterns from the analysis of three mRNA pools (total mRNA, SL1 and SL2 fraction) by RNA-seq. This information was combined with a genome-wide analysis of gene orientation and spacing. We could confirm 2219 operons by RNA-seq data out of 6709 candidate operons, which were predicted by sequence information alone. Our gene order comparison of the Caenorhabditis elegans and P. pacificus genomes shows major changes in operon organization in the two species. Notably, only 128 out of 1288 operons in C. elegans are conserved in P. pacificus. However, analysis of gene-expression profiles identified conserved functions such as an enrichment of germline-expressed genes and higher expression levels of operonic genes during recovery from dauer arrest in both species. These results provide support for the model that a necessity for increased transcriptional efficiency in the context of certain developmental processes could be a selective constraint for operon evolution in metazoans. Our method is generally applicable to other metazoans to see if similar functional constraints regulate gene organization into operons.

Up-regulating CYP3A4 expression in C3A cells by transfection with a novel chimeric regulator of hPXR-p53-AD

Most hepatoma cell lines lack proper expression and induction of CYP3A4 enzyme, which limits their use for predicting drug metabolism and toxicity. Nuclear receptor pregnane X receptor (PXR) has been well recognized for its critical role in regulating expression of CYP3A4 gene. However, its physiological activity of binding to the particular site of promoter is significantly weakened in hepatic cell lines. To address this problem, we created “chimeric PXR” constructs by appending a strong activation domain (AD) from p53 subunit to either N- or C- termini of the human PXR (hPXR), that is, hPXR-p53 and p53-hPXR. C3A, a hepatoma cell line, was used as the cell model to test the regulation effect of chimeric hPXR over wild type (WT) hPXR on CYP3A4 expression at gene, protein, and metabolism levels, respectively. Compared with C3A cells transiently transfected with WT hPXR, the activity of CYP3A4.XREM.luc reporter gene in C3A cells transfected with hPXR-p53 or p53-hPXR increased 5- and 9-fold respectively, and the levels of CYP3A4 mRNA expression increased 3.5- and 2.6-fold, respectively. C3A cells stably transfected with hPXR-p53-AD exhibited an improved expression of CYP3A4 at both gene (2-fold) and protein (1.5-fold) levels compared to WT C3A cells. Testosterone, a CYP3A4-specific substrate, was used for detecting the metabolism activity of CYP3A4. No testosterone metabolite could be detected in microsomes from WT C3A cells and WT C3A cells-based array, while the formation of 6β-hydroxytestosterone metabolite in the transfected cells was 714 and 55 pmol/mg protein/min, respectively. In addition, all the above expression levels in the transfected cell models could be further induced with additional treatment of Rifampicin, a specific inducer for CYP3A4. In conclusion, our study established a proof-of-principle example that genetic modification with chimeric hPXR-p53-AD could improve CYP3A4 metabolism ability in hepatic cell line.

Complete mitochondrial genome of the anadromous river pufferfish, Takifugu obscurus (Tetraodontiformes: Tetraodontidae)

The complete mitochondrial genome was sequenced from the anadromous river pufferfish, Takifugu obscurus. The genome sequence was 16,446 bp in size, and the gene order and contents were identical with those of congeneric species in the genus Takifugu. Of 13 protein-coding genes (PGCs), 2 genes (CO2 and ND4) had incomplete stop codons as shown in T. obscurus. Furthermore, the stop codon of CO1 and ND6 genes was AGG. The base composition of T. obscurus mitogenome showed anti-G bias (13.2% and 6.4%) on the second and third positions of protein-coding genes (PCGs), respectively.

Complete mitochondrial genome of Southern Smiler, Opistognathus jacksoniensis (Perciformes: Opistognathidae)

The first complete mitochondrial genome (mitogenome) of Opistognathidae (Teleostei: Perciformes) has been determined using the Southern Smiler, Opistognathus jacksoniensis, as a representative species of the family. This mitogenome was 16,650 bp in length, including 13 typical vertebrate protein-coding genes, 22 transfer RNA (tRNA) genes, 2 ribosomal RNA genes and 1 control region. All genes were encoded on the heavy strain except for ND6 and eight tRNA genes. The overall base composition of the heavy strain was 27.4% for A, 30.2% for C, 25.7% for T and 16.7% for G. The mitogenome data of O. jacksoniensis should contribute to clarify phylogenetic position of Opistognathidae under mitogenome-based phylogenetic reconstruction of perciform fishes.

Complete mitochondrial genome of the freshwater gudgeon, Pseudopungtungia tenuicorpa (Cypriniformes, Gobioninae)

The complete mitochondrial genome was sequenced from the freshwater fish, Pseudopungtungia tenuicorpa. The genome sequence was 16,590 bp in size, and the gene order and contents were identical with P. nigra in the same genus Pseudopungtungia without any peculiar features. Of 13 protein-coding genes (PGCs), 3 genes (CO2, CO3 and Cytb) had incomplete stop codon as shown in P. nigra. The base composition of P. tenuicorpa showed anti-G bias (13.94% and 12.73%) on the second and third positions of PCGs, indicating that P. tenuicorpa showed slightly lower anti-G bias than P. nigra (13.71% and 9.48%).

Complete mitochondrial genome of the Bloch's gizzard shad Nematalosa nasus (Clupeiformes: Clupeidae)

The complete mitochondrial genome sequence of Nematalosa nasus was determined in this paper. The entire sequence was 16,674 bp in length and the overall base compositions were 26.5% of A, 28.4% of C, 25.4% of T and 19.7% of G. The mitogenome contained 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes and two non-coding regions (the control region and the light strand replication origin). Gene arrangement was similar to that observed in most other vertebrates. Except for the ND6 genes and eight tRNA genes, all other genes were encoded on the heavy strand. Within the control region, the extended termination associated sequence was identified, as well as four central conserved sequence blocks (CSB-F, CSB-E, CSB-D and CSB-C) and three conserved sequence blocks (CSB-1, CSB-2 and CSB-3).

Complete mitochondrial genome of squaretail coralgrouper Plectropomus areolatus (Perciformes, Epinephelidae)

The complete mitochondrial genome (mitogenome) sequence of squaretail coralgrouper Plectropomus areolatus was determined with long PCR approach. The genome was 16,770 bp in length and contained 13 protein-coding genes, 22 transfer RNAs (tRNAs), 2 ribosomal RNAs (rRNAs), and 2 non-coding regions, with the gene content and order being similar to most other teleosts. The overall base composition of the mitogenome was estimated to be 28.81% A, 27.63% T, 27.24% C, and 16.33% G, with a high AT content of 56.44%, indicating an obvious anti-guanine bias commonly observed in fishes. It shared 90.7%, 78.8%, and 78.7% mitogenome sequence with Plectropomus leopardus, Epinephelus moara, and Epinephelus coioides, respectively.

Mitochondrial genome of the Luciogobius platycephalus (Perciformes, Gobioidei)

In this paper, the complete mitochondrial genome of Luciogobius platycephalus was first determined. The genome is 16,478 bp in length and consists of 13 protein-coding genes, 22 tRNA genes, 2 ribosomal RNA genes, and 2 main non-coding regions (the control region and the origin of the light-strand replication). The overall base composition of L. platycephalus is 29.1% for T, 26.7% for C, 29.2% for A, and 15.1% for G, with a slight A+T bias of 58.3%. It has the typical vertebrate mitochondrial gene arrangement.

The complete mitochondrial genome of the hybrid of Ctenopharyngodon idella (♀) × Squaliobarbus curriculus (♂)

The complete mitochondrial genome sequence of the hybrid of Ctenopharyngodon idella (♀) × Squaliobarbus curriculus (♂) was determined using PCR-based method. The total length of the mitogenome is 16,609 bp. It contains the typical structure as that of most other vertebrates, including 2 ribosomal RNA genes, 13 protein-coding genes, 22 transfer RNA genes, and 1 non-coding control region (D-loop region). The overall composition of the mitogenome was estimated to be 31.86% for A, 26.08% for T, 26.38% for C, and 15.68% for G, respectively, indicating that an A+T (57.94%)-rich feature occurs in the hybrid mitogenome. Both the termination-associated sequence and three conserved sequence blocks (CSB1, CSB2, and CSB3) were also detected in the D-loop region.

Complete mitochondrial genome of Chinese sleeper, Perccottus glenii

Chinese sleeper (Perccottus glenii) belongs to the largest vertebrate order, Perciformes. In this study, the complete mitochondrial genome of P. glenii was determined to be 16,487 bp in length, including 13 protein-coding genes, 22 transfer RNA genes and 2 ribosomal RNA genes. We also analysed the sequence structure of non-coding control region. Comparing the mitochondrial genome of P. glenii with its congener Rhyacichthys aspro showed sequence similarity and the identical gene arrangement. The complete mitochondrial genome of Chinese sleeper provides the basis for the studies in Perciformes evolution and conservation.

Pro-inflammatory flagellin proteins of prevalent motile commensal bacteria are variably abundant in the intestinal microbiome of elderly humans

Some Eubacterium and Roseburia species are among the most prevalent motile bacteria present in the intestinal microbiota of healthy adults. These flagellate species contribute “cell motility” category genes to the intestinal microbiome and flagellin proteins to the intestinal proteome. We reviewed and revised the annotation of motility genes in the genomes of six Eubacterium and Roseburia species that occur in the human intestinal microbiota and examined their respective locus organization by comparative genomics. Motility gene order was generally conserved across these loci. Five of these species harbored multiple genes for predicted flagellins. Flagellin proteins were isolated from R. inulinivorans strain A2-194 and from E. rectale strains A1-86 and M104/1. The amino-termini sequences of the R. inulinivorans and E. rectale A1-86 proteins were almost identical. These protein preparations stimulated secretion of interleukin-8 (IL-8) from human intestinal epithelial cell lines, suggesting that these flagellins were pro-inflammatory. Flagellins from the other four species were predicted to be pro-inflammatory on the basis of alignment to the consensus sequence of pro-inflammatory flagellins from the β- and γ- proteobacteria. Many fliC genes were deduced to be under the control of σ²⁸. The relative abundance of the target Eubacterium and Roseburia species varied across shotgun metagenomes from 27 elderly individuals. Genes involved in the flagellum biogenesis pathways of these species were variably abundant in these metagenomes, suggesting that the current depth of coverage used for metagenomic sequencing (3.13–4.79 Gb total sequence in our study) insufficiently captures the functional diversity of genomes present at low (≤1%) relative abundance. E. rectale and R. inulinivorans thus appear to synthesize complex flagella composed of flagellin proteins that stimulate IL-8 production. A greater depth of sequencing, improved evenness of sequencing and improved metagenome assembly from short reads will be required to facilitate in silico analyses of complete complex biochemical pathways for low-abundance target species from shotgun metagenomes.

The complete mitochondrial genome of the land snail Cornu aspersum (Helicidae: Mollusca): intra-specific divergence of protein-coding genes and phylogenetic considerations within Euthyneura

The complete sequences of three mitochondrial genomes from the land snail Cornu aspersum were determined. The mitogenome has a length of 14050 bp, and it encodes 13 protein-coding genes, 22 transfer RNA genes and two ribosomal RNA genes. It also includes nine small intergene spacers, and a large AT-rich intergenic spacer. The intra-specific divergence analysis revealed that COX1 has the lower genetic differentiation, while the most divergent genes were NADH1, NADH3 and NADH4. With the exception of Euhadra herklotsi, the structural comparisons showed the same gene order within the family Helicidae, and nearly identical gene organization to that found in order Pulmonata. Phylogenetic reconstruction recovered Basommatophora as polyphyletic group, whereas Eupulmonata and Pulmonata as paraphyletic groups. Bayesian and Maximum Likelihood analyses showed that C. aspersum is a close relative of Cepaea nemoralis, and with the other Helicidae species form a sister group of Albinaria caerulea, supporting the monophyly of the Stylommatophora clade.

The "fossilized" mitochondrial genome of Liriodendron tulipifera: ancestral gene content and order, ancestral editing sites, and extraordinarily low mutation rate

BACKGROUND

The mitochondrial genomes of flowering plants vary greatly in size, gene content, gene order, mutation rate and level of RNA editing. However, the narrow phylogenetic breadth of available genomic data has limited our ability to reconstruct these traits in the ancestral flowering plant and, therefore, to infer subsequent patterns of evolution across angiosperms.

RESULTS

We sequenced the mitochondrial genome of Liriodendron tulipifera, the first from outside the monocots or eudicots. This 553,721 bp mitochondrial genome has evolved remarkably slowly in virtually all respects, with an extraordinarily low genome-wide silent substitution rate, retention of genes frequently lost in other angiosperm lineages, and conservation of ancestral gene clusters. The mitochondrial protein genes in Liriodendron are the most heavily edited of any angiosperm characterized to date. Most of these sites are also edited in various other lineages, which allowed us to polarize losses of editing sites in other parts of the angiosperm phylogeny. Finally, we added comprehensive gene sequence data for two other magnoliids, Magnolia stellata and the more distantly related Calycanthus floridus, to measure rates of sequence evolution in Liriodendron with greater accuracy. The Magnolia genome has evolved at an even lower rate, revealing a roughly 5,000-fold range of synonymous-site divergence among angiosperms whose mitochondrial gene space has been comprehensively sequenced.

CONCLUSIONS

Using Liriodendron as a guide, we estimate that the ancestral flowering plant mitochondrial genome contained 41 protein genes, 14 tRNA genes of mitochondrial origin, as many as 7 tRNA genes of chloroplast origin, >700 sites of RNA editing, and some 14 colinear gene clusters. Many of these gene clusters, genes and RNA editing sites have been variously lost in different lineages over the course of the ensuing ∽200 million years of angiosperm evolution.

Complete mitochondrial genome of the free-living earwig, Challia fletcheri (Dermaptera: Pygidicranidae) and phylogeny of Polyneoptera

The insect order Dermaptera, belonging to Polyneoptera, includes ∼2,000 extant species, but no dermapteran mitochondrial genome has been sequenced. We sequenced the complete mitochondrial genome of the free-living earwig, Challia fletcheri, compared its genomic features to other available mitochondrial sequences from polyneopterous insects. In addition, the Dermaptera, together with the other known polyneopteran mitochondrial genome sequences (protein coding, ribosomal RNA, and transfer RNA genes), were employed to understand the phylogeny of Polyneoptera, one of the least resolved insect phylogenies, with emphasis on the placement of Dermaptera. The complete mitochondrial genome of C. fletcheri presents the following several unusual features: the longest size in insects is 20,456 bp; it harbors the largest tandem repeat units (TRU) among insects; it displays T- and G-skewness on the major strand and A- and C-skewness on the minor strand, which is a reversal of the general pattern found in most insect mitochondrial genomes, and it possesses a unique gene arrangement characterized by a series of gene translocations and/or inversions. The reversal pattern of skewness is explained in terms of inversion of replication origin. All phylogenetic analyses consistently placed Dermaptera as the sister to Plecoptera, leaving them as the most basal lineage of Polyneoptera or sister to Ephemeroptera, and placed Odonata consistently as the most basal lineage of the Pterygota.

High confidence prediction of essential genes in Burkholderia cenocepacia

Background

Essential genes are absolutely required for the survival of an organism. The identification of essential genes, besides being one of the most fundamental questions in biology, is also of interest for the emerging science of synthetic biology and for the development of novel antimicrobials. New antimicrobial therapies are desperately needed to treat multidrug-resistant pathogens, such as members of the Burkholderia cepacia complex.

Methodology/Principal Findings

We hypothesize that essential genes may be highly conserved within a group of evolutionary closely related organisms. Using a bioinformatics approach we determined that the core genome of the order Burkholderiales consists of 649 genes. All but two of these identified genes were located on chromosome 1 of Burkholderia cenocepacia. Although many of the 649 core genes of Burkholderiales have been shown to be essential in other bacteria, we were also able to identify a number of novel essential genes present mainly, or exclusively, within this order. The essentiality of some of the core genes, including the known essential genes infB, gyrB, ubiB, and valS, as well as the so far uncharacterized genes BCAL1882, BCAL2769, BCAL3142 and BCAL3369 has been confirmed experimentally in B. cenocepacia.

Conclusions/Significance

We report on the identification of essential genes using a novel bioinformatics strategy and provide bioinformatics and experimental evidence that the large majority of the identified genes are indeed essential. The essential genes identified here may represent valuable targets for the development of novel antimicrobials and their detailed study may shed new light on the functions required to support life.

Conservation of gene order and content in the circular chromosomes of 'Candidatus Liberibacter asiaticus' and other Rhizobiales

'Ca. Liberibacter asiaticus,' an insect-vectored, obligate intracellular bacterium associated with citrus-greening disease, also called "HLB," is a member of the Rhizobiales along with nitrogen-fixing microsymbionts Sinorhizobium meliloti and Bradyrhizobium japonicum, plant pathogen Agrobacterium tumefaciens and facultative intracellular mammalian pathogen Bartonella henselae. Comparative analyses of their circular chromosomes identified 514 orthologous genes shared among all five species. Shared among all five species are 50 identical blocks of microsyntenous orthologous genes (MOGs), containing a total of 283 genes. While retaining highly conserved genomic blocks of microsynteny, divergent evolution, horizontal gene transfer and niche specialization have disrupted macrosynteny among the five circular chromosomes compared. Highly conserved microsyntenous gene clusters help define the Rhizobiales, an order previously defined by 16S RNA gene similarity and herein represented by the three families: Bartonellaceae, Bradyrhizobiaceae and Rhizobiaceae. Genes without orthologs in the other four species help define individual species. The circular chromosomes of each of the five Rhizobiales species examined had genes lacking orthologs in the other four species. For example, 63 proteins are encoded by genes of 'Ca. Liberibacter asiaticus' not shared with other members of the Rhizobiales. Of these 63 proteins, 17 have predicted functions related to DNA replication or RNA transcription, and some of these may have roles related to low genomic GC content. An additional 17 proteins have predicted functions relevant to cellular processes, particularly modifications of the cell surface. Seventeen unshared proteins have specific metabolic functions including a pathway to synthesize cholesterol encoded by a seven-gene operon. The remaining 12 proteins encoded by 'Ca. Liberibacter asiaticus' genes not shared with other Rhizobiales are of bacteriophage origin. 'Ca. Liberibacter asiaticus' shares 11 genes with only Sinorhizobium meliloti and 12 genes are shared with only Bartonella henselae.

Evaluation of the role of functional constraints on the integrity of an ultraconserved region in the genus Drosophila

Why gene order is conserved over long evolutionary timespans remains elusive. A common interpretation is that gene order conservation might reflect the existence of functional constraints that are important for organismal performance. Alteration of the integrity of genomic regions, and therefore of those constraints, would result in detrimental effects. This notion seems especially plausible in those genomes that can easily accommodate gene reshuffling via chromosomal inversions since genomic regions free of constraints are likely to have been disrupted in one or more lineages. Nevertheless, no empirical test has been performed to this notion. Here, we disrupt one of the largest conserved genomic regions of the Drosophila genome by chromosome engineering and examine the phenotypic consequences derived from such disruption. The targeted region exhibits multiple patterns of functional enrichment suggestive of the presence of constraints. The carriers of the disrupted collinear block show no defects in their viability, fertility, and parameters of general homeostasis, although their odorant perception is altered. This change in odorant perception does not correlate with modifications of the level of expression and sex bias of the genes within the genomic region disrupted. Our results indicate that even in highly rearranged genomes, like those of Diptera, unusually high levels of gene order conservation cannot be systematically attributed to functional constraints, which raises the possibility that other mechanisms can be in place and therefore the underpinnings of the maintenance of gene organization might be more diverse than previously thought.

The complete mitochondrial genome and novel gene arrangement of the unique-headed bug Stenopirates sp. (Hemiptera: Enicocephalidae)

Many of true bugs are important insect pests to cultivated crops and some are important vectors of human diseases, but few cladistic analyses have addressed relationships among the seven infraorders of Heteroptera. The Enicocephalomorpha and Nepomorpha are consider the basal groups of Heteroptera, but the basal-most lineage remains unresolved. Here we report the mitochondrial genome of the unique-headed bug Stenopirates sp., the first mitochondrial genome sequenced from Enicocephalomorpha. The Stenopirates sp. mitochondrial genome is a typical circular DNA molecule of 15, 384 bp in length, and contains 37 genes and a large non-coding fragment. The gene order differs substantially from other known insect mitochondrial genomes, with rearrangements of both tRNA genes and protein-coding genes. The overall AT content (82.5%) of Stenopirates sp. is the highest among all the known heteropteran mitochondrial genomes. The strand bias is consistent with other true bugs with negative GC-skew and positive AT-skew for the J-strand. The heteropteran mitochondrial atp8 exhibits the highest evolutionary rate, whereas cox1 appears to have the lowest rate. Furthermore, a negative correlation was observed between the variation of nucleotide substitutions and the GC content of each protein-coding gene. A microsatellite was identified in the putative control region. Finally, phylogenetic reconstruction suggests that Enicocephalomorpha is the sister group to all the remaining Heteroptera.

Discovering patterns in gene order

Various genetic events during the process of natural evolution shape the landscape of the genomes. In this chapter, we explore an approach to investigating multiple genomes in order to unravel their complex relationships that go beyond their placement on a phylogeny. To this end, we treat genes as the smallest syntactic unit on the genome and explore their relative organization across multiple genomes. In the first half of the chapter, we discuss mathematical models to capture the combinatorial structures of this relative organization and statistical models to study their distributions. In the second half of the chapter, we apply these models to analyze the relationship between three closely related plant genomes.

Intron gains and losses in the evolution of Fusarium and Cryptococcus fungi

The presence of spliceosomal introns in eukaryotic genes poses a major puzzle for the study of genome evolution. Intron densities vary enormously among distant lineages. However, the mechanisms driving intron gains are poorly understood and very few intron gains and losses have been documented over short evolutionary time spans. Fungi emerged recently as excellent models to study intron evolution and "reverse splicing" was found to be a major driver of recent intron gains in a clade of ascomycete fungi. We screened a total of 38 genomes from two fungal clades important in medicine and agriculture to identify intron gains and losses both within and between species. We detected 86 and 198 variable intron positions in the Cryptococcus and Fusarium clades, respectively. Some genes underwent extensive changes in their exon-intron structure, with up to six variable intron positions per gene. We identified a very recently gained intron in a group of tomato-infecting strains belonging to the F. oxysporum species complex. In the human pathogen C. gattii, we found recent intron losses in subtypes of the species. The two studied fungal clades provided evidence for extensive changes in their exon-intron structure within and among closely related species. We show that both intronization of previously coding DNA and insertion of exogenous DNA are the major drivers of intron gains.

Analysis of gene order evolution beyond single-copy genes

The purpose of this chapter is to provide a comprehensive review of the field of genome rearrangement, i.e., comparative genomics, based on the representation of genomes as ordered sequences of signed genes. We specifically focus on the "hard part" of genome rearrangement, how to handle duplicated genes. The main questions are: how have present-day genomes evolved from a common ancestor? What are the most realistic evolutionary scenarios explaining the observed gene orders? What was the content and structure of ancestral genomes? We aim to provide a concise but complete overview of the field, starting with the practical problem of finding an appropriate representation of a genome as a sequence of ordered genes or blocks, namely the problems of orthology, paralogy, and synteny block identification. We then consider three levels of gene organization: the gene family level (evolution by duplication, loss, and speciation), the cluster level (evolution by tandem duplications), and the genome level (all types of rearrangement events, including whole genome duplication).

Identification of Anaplasma marginale type IV secretion system effector proteins

BACKGROUND

Anaplasma marginale, an obligate intracellular alphaproteobacterium in the order Rickettsiales, is a tick-borne pathogen and the leading cause of anaplasmosis in cattle worldwide. Complete genome sequencing of A. marginale revealed that it has a type IV secretion system (T4SS). The T4SS is one of seven known types of secretion systems utilized by bacteria, with the type III and IV secretion systems particularly prevalent among pathogenic Gram-negative bacteria. The T4SS is predicted to play an important role in the invasion and pathogenesis of A. marginale by translocating effector proteins across its membrane into eukaryotic target cells. However, T4SS effector proteins have not been identified and tested in the laboratory until now.

RESULTS

By combining computational methods with phylogenetic analysis and sequence identity searches, we identified a subset of potential T4SS effectors in A. marginale strain St. Maries and chose six for laboratory testing. Four (AM185, AM470, AM705 [AnkA], and AM1141) of these six proteins were translocated in a T4SS-dependent manner using Legionella pneumophila as a reporter system.

CONCLUSIONS

The algorithm employed to find T4SS effector proteins in A. marginale identified four such proteins that were verified by laboratory testing. L. pneumophila was shown to work as a model system for A. marginale and thus can be used as a screening tool for A. marginale effector proteins. The first T4SS effector proteins for A. marginale have been identified in this work.

Genomic sequence analysis of granulovirus isolated from the tobacco cutworm, Spodoptera litura

BACKGROUND

Spodoptera litura is a noctuid moth that is considered an agricultural pest. The larvae feed on a wide range of plants and have been recorded on plants from 40 plant families (mostly dicotyledons). It is a major pest of many crops. To better understand Spodoptera litura granulovirus (SpliGV), the nucleotide sequence of the SpliGV DNA genome was determined and analyzed.

METHODOLOGY/PRINCIPAL FINDINGS

The genome of the SpliGV was completely sequenced. The nucleotide sequence of the SpliGV genome was 124,121 bp long with 61.2% A+T content and contained 133 putative open reading frames (ORFs) of 150 or more nucleotides. The 133 putative ORFs covered 86.3% of the genome. Among these, 31 ORFs were conserved in most completely sequenced baculovirus genomes, 38 were granulovirus (GV)-specific, and 64 were present in some nucleopolyhedroviruses (NPVs) and/or GVs. We proved that 9 of the ORFs were SpliGV specific.

CONCLUSIONS/SIGNIFICANCE

The genome of SpliGV is 124,121 bp in size. One hundred thirty-three ORFs that putatively encode proteins of 50 or more amino acid residues with minimal overlap were determined. No chitinase or cathepsin genes, which are involved in the liquefaction of the infected host, were found in the SpliGV genome, explaining why SpliGV-infected insects do not degrade in a typical manner. The DNA photolyase gene was first found in the genus Granulovirus. When phylogenic relationships were analyzed, the SpliGV was most closely related to Trichoplusia ni granulovirus (TnGV) and Xestia c-nigrum granulovirus (XecnGV), which belong to the Type I-granuloviruses (Type I-GV).

First survey of the wheat chromosome 5A composition through a next generation sequencing approach

Wheat is one of the world's most important crops and is characterized by a large polyploid genome. One way to reduce genome complexity is to isolate single chromosomes using flow cytometry. Low coverage DNA sequencing can provide a snapshot of individual chromosomes, allowing a fast characterization of their main features and comparison with other genomes. We used massively parallel 454 pyrosequencing to obtain a 2x coverage of wheat chromosome 5A. The resulting sequence assembly was used to identify TEs, genes and miRNAs, as well as to infer a virtual gene order based on the synteny with other grass genomes. Repetitive elements account for more than 75% of the genome. Gene content was estimated considering non-redundant reads showing at least one match to ESTs or proteins. The results indicate that the coding fraction represents 1.08% and 1.3% of the short and long arm respectively, projecting the number of genes of the whole chromosome to approximately 5,000. 195 candidate miRNA precursors belonging to 16 miRNA families were identified. The 5A genes were used to search for syntenic relationships between grass genomes. The short arm is closely related to Brachypodium chromosome 4, sorghum chromosome 8 and rice chromosome 12; the long arm to regions of Brachypodium chromosomes 4 and 1, sorghum chromosomes 1 and 2 and rice chromosomes 9 and 3. From these similarities it was possible to infer the virtual gene order of 392 (5AS) and 1,480 (5AL) genes of chromosome 5A, which was compared to, and found to be largely congruent with the available physical map of this chromosome.

The mitochondrial genomes of the early land plants Treubia lacunosa and Anomodon rugelii: dynamic and conservative evolution

Early land plant mitochondrial genomes captured important changes of mitochondrial genome evolution when plants colonized land. The chondromes of seed plants show several derived characteristics, e.g., large genome size variation, rapid intra-genomic rearrangement, abundant introns, and highly variable levels of RNA editing. On the other hand, the chondromes of charophytic algae are still largely ancestral in these aspects, resembling those of early eukaryotes. When the transition happened has been a long-standing question in studies of mitochondrial genome evolution. Here we report complete mitochondrial genome sequences from an early-diverging liverwort, Treubia lacunosa, and a late-evolving moss, Anomodon rugelii. The two genomes, 151,983 and 104,239 base pairs in size respectively, contain standard sets of protein coding genes for respiration and protein synthesis, as well as nearly full sets of rRNA and tRNA genes found in the chondromes of the liverworts Marchantia polymorpha and Pleurozia purpurea and the moss Physcomitrella patens. The gene orders of these two chondromes are identical to those of the other liverworts and moss. Their intron contents, with all cis-spliced group I or group II introns, are also similar to those in the previously sequenced liverwort and moss chondromes. These five chondromes plus the two from the hornworts Phaeoceros laevis and Megaceros aenigmaticus for the first time allowed comprehensive comparative analyses of structure and organization of mitochondrial genomes both within and across the three major lineages of bryophytes. These analyses led to the conclusion that the mitochondrial genome experienced dynamic evolution in genome size, gene content, intron acquisition, gene order, and RNA editing during the origins of land plants and their major clades. However, evolution of this organellar genome has remained rather conservative since the origin and initial radiation of early land plants, except within vascular plants.

Genomic sequence around butterfly wing development genes: annotation and comparative analysis

BACKGROUND

Analysis of genomic sequence allows characterization of genome content and organization, and access beyond gene-coding regions for identification of functional elements. BAC libraries, where relatively large genomic regions are made readily available, are especially useful for species without a fully sequenced genome and can increase genomic coverage of phylogenetic and biological diversity. For example, no butterfly genome is yet available despite the unique genetic and biological properties of this group, such as diversified wing color patterns. The evolution and development of these patterns is being studied in a few target species, including Bicyclus anynana, where a whole-genome BAC library allows targeted access to large genomic regions.

METHODOLOGY/PRINCIPAL FINDINGS

We characterize ∼1.3 Mb of genomic sequence around 11 selected genes expressed in B. anynana developing wings. Extensive manual curation of in silico predictions, also making use of a large dataset of expressed genes for this species, identified repetitive elements and protein coding sequence, and highlighted an expansion of Alcohol dehydrogenase genes. Comparative analysis with orthologous regions of the lepidopteran reference genome allowed assessment of conservation of fine-scale synteny (with detection of new inversions and translocations) and of DNA sequence (with detection of high levels of conservation of non-coding regions around some, but not all, developmental genes).

CONCLUSIONS

The general properties and organization of the available B. anynana genomic sequence are similar to the lepidopteran reference, despite the more than 140 MY divergence. Our results lay the groundwork for further studies of new interesting findings in relation to both coding and non-coding sequence: 1) the Alcohol dehydrogenase expansion with higher similarity between the five tandemly-repeated B. anynana paralogs than with the corresponding B. mori orthologs, and 2) the high conservation of non-coding sequence around the genes wingless and Ecdysone receptor, both involved in multiple developmental processes including wing pattern formation.

The mitochondrial genome of Sinentomon erythranum (Arthropoda: Hexapoda: Protura): an example of highly divergent evolution

BACKGROUND

The phylogenetic position of the Protura, traditionally considered the most basal hexapod group, is disputed because it has many unique morphological characters compared with other hexapods. Although mitochondrial genome information has been used extensively in phylogenetic studies, such information is not available for the Protura. This has impeded phylogenetic studies on this taxon, as well as the evolution of the arthropod mitochondrial genome.

RESULTS

In this study, the mitochondrial genome of Sinentomon erythranum was sequenced, as the first proturan species to be reported. The genome contains a number of special features that differ from those of other hexapods and arthropods. As a very small arthropod mitochondrial genome, its 14,491 nucleotides encode 37 typical mitochondrial genes. Compared with other metazoan mtDNA, it has the most biased nucleotide composition with T = 52.4%, an extreme and reversed AT-skew of -0.351 and a GC-skew of 0.350. Two tandemly repeated regions occur in the A+T-rich region, and both could form stable stem-loop structures. Eighteen of the 22 tRNAs are greatly reduced in size with truncated secondary structures. The gene order is novel among available arthropod mitochondrial genomes. Rearrangements have involved in not only small tRNA genes, but also PCGs (protein-coding genes) and ribosome RNA genes. A large block of genes has experienced inversion and another nearby block has been reshuffled, which can be explained by the tandem duplication and random loss model. The most remarkable finding is that trnL2(UUR) is not located between cox1 and cox2 as observed in most hexapod and crustacean groups, but is between rrnL and nad1 as in the ancestral arthropod ground pattern. The "cox1-cox2" pattern was further confirmed in three more representative proturan species. The phylogenetic analyses based on the amino acid sequences of 13 mitochondrial PCGs suggest S. erythranum failed to group with other hexapod groups.

CONCLUSIONS

The mitochondrial genome of S. erythranum shows many different features from other hexapod and arthropod mitochondrial genomes. It underwent highly divergent evolution. The "cox1-cox2" pattern probably represents the ancestral state for all proturan mitogenomes, and suggests a long evolutionary history for the Protura.

The complete mitochondrial genome of the grand jackknife clam, Solen grandis (Bivalvia: Solenidae): a novel gene order and unusual non-coding region

Molluscs in general, and bivalves in particular, exhibit an extraordinary degree of mitochondrial gene order variation when compared with other metazoans. The complete mitochondrial genome of Solen grandis (Bivalvia: Solenidae) was determined using long-PCR and genome walking techniques. The entire mitochondrial genome sequence of S. grandis is 16,784 bp in length, and contains 36 genes including 12 protein-coding genes (atp8 is absent), 2 ribosomal RNAs, and 22 tRNAs. All genes are encoded on the same strand. Compared with other species, it bears a novel gene order. Besides these, we find a peculiar non-coding region of 435 bp with a microsatellite-like (TA)(12) element, poly-structures and many hairpin structures. In contrast to the available heterodont mitochondrial genomes from GenBank, the complete mtDNA of S. grandis has the shortest cox3 gene, and the longest atp6, nad4, nad5 genes.

Sex-linked pheromone receptor genes of the European corn borer, Ostrinia nubilalis, are in tandem arrays

Background

Tuning of the olfactory system of male moths to conspecific female sex pheromones is crucial for correct species recognition; however, little is known about the genetic changes that drive speciation in this system. Moths of the genus Ostrinia are good models to elucidate this question, since significant differences in pheromone blends are observed within and among species. Odorant receptors (ORs) play a critical role in recognition of female sex pheromones; eight types of OR genes expressed in male antennae were previously reported in Ostrinia moths.

Methodology/Principal Findings

We screened an O. nubilalis bacterial artificial chromosome (BAC) library by PCR, and constructed three contigs from isolated clones containing the reported OR genes. Fluorescence in situ hybridization (FISH) analysis using these clones as probes demonstrated that the largest contig, which contained eight OR genes, was located on the Z chromosome; two others harboring two and one OR genes were found on two autosomes. Sequence determination of BAC clones revealed the Z-linked OR genes were closely related and tandemly arrayed; moreover, four of them shared 181-bp direct repeats spanning exon 7 and intron 7.

Conclusions/Significance

This is the first report of tandemly arrayed sex pheromone receptor genes in Lepidoptera. The localization of an OR gene cluster on the Z chromosome agrees with previous findings for a Z-linked locus responsible for O. nubilalis male behavioral response to sex pheromone. The 181-bp direct repeats might enhance gene duplications by unequal crossovers. An autosomal locus responsible for male response to sex pheromone in Heliothis virescens and H. subflexa was recently reported to contain at least four OR genes. Taken together, these findings support the hypothesis that generation of additional copies of OR genes can increase the potential for male moths to acquire altered specificity for pheromone components, and accordingly, facilitate differentiation of sex pheromones.