Combined Genomic and Proteomic Approaches Identify Gene Clusters Involved in Anaerobic 2-Methylnaphthalene Degradation in the Sulfate-Reducing Enrichment Culture N47

The highly enriched deltaproteobacterial culture N47 anaerobically oxidizes the polycyclic aromatic hydrocarbons naphthalene and 2-methylnaphthalene, with sulfate as the electron acceptor. Combined genome sequencing and liquid chromatography-tandem mass spectrometry-based shotgun proteome analyses were performed to identify genes and proteins involved in anaerobic aromatic catabolism. Proteome analysis of 2-methylnaphthalene-grown N47 cells resulted in the identification of putative enzymes catalyzing the anaerobic conversion of 2-methylnaphthalene to 2-naphthoyl coenzyme A (2-naphthoyl-CoA), as well as the reductive ring cleavage of 2-naphthoyl-CoA, leading to the formation of acetyl-CoA and CO₂. The glycyl radical-catalyzed fumarate addition to the methyl group of 2-methylnaphthalene is catalyzed by naphthyl-2-methyl-succinate synthase (Nms), composed of α-, β-, and γ-subunits that are encoded by the genes nmsABC. Located upstream of nmsABC is nmsD, encoding the Nms-activating enzyme, which harbors the characteristic [Fe₄S₄] cluster sequence motifs of S-adenosylmethionine radical enzymes. The bns gene cluster, coding for enzymes involved in beta-oxidation reactions converting naphthyl-2-methyl-succinate to 2-naphthoyl-CoA, was found four intervening open reading frames further downstream. This cluster consists of eight genes (bnsABCDEFGH) corresponding to 8.1 kb, which are closely related to genes for enzymes involved in anaerobic toluene degradation within the denitrifiers “Aromatoleum aromaticum” EbN1, Azoarcus sp. strain T, and Thauera aromatica. Another contiguous DNA sequence harbors the gene for 2-naphthoyl-CoA reductase (ncr) and 16 additional genes that were found to be expressed in 2-methylnaphthalene-grown cells. These genes code for enzymes that were supposed to catalyze the dearomatization and ring cleavage reactions converting 2-naphthoyl-CoA to acetyl-CoA and CO₂. Comparative sequence analysis of the four encoding subunits (ncrABCD) showed the gene product to have the closest similarity to the Azoarcus type of benzoyl-CoA reductase. The present work provides the first insight into the genetic basis of anaerobic 2-methylnaphthalene metabolism and delivers implications for understanding contaminant degradation.Polycyclic aromatic hydrocarbons (PAHs) are constantly released into the environment by anthropogenic activities such as industrial use or by accidental contamination. Due to the low chemical reactivity caused by the resonance energy of the aromatic ring structure and the low bioavailability of PAHs, they are persistent in the environment (15). The understanding of microbial metabolic capabilities in terms of anaerobic PAH degradation is in its infancy. However, natural amelioration of contaminated sites relies on the degradation capacities of microorganisms, and therefore, it is an essential prerequisite to broaden knowledge about the microorganisms involved and their potentials concerning PAH breakdown.Numerous microorganisms that can degrade PAHs under aerobic conditions have already been identified, but only a small number of anaerobic cultures that degrade PAHs like naphthalene, 2-methylnaphthalene, and phenanthrene have been isolated so far (17, 20, 24, 31, 46-48, 50, 52, 66). It has been shown that these anaerobic degraders activate aromatic hydrocarbons by very unusual biochemical reactions which differ completely from those of aerobic degradation. The peripheral pathway of 2-methylnaphthalene degradation occurs in analogy to anaerobic toluene degradation by the addition of fumarate to the methyl group, catalyzed by the glycyl radical enzyme naphthyl-2-methyl-succinate synthase (Nms) (Fig. ​(Fig.1)1) (3). In subsequent reactions, naphthyl-2-methyl-succinate is activated to yield the coenzyme A (CoA) ester and oxidized to form naphthyl-2-methylene-succinyl-CoA. The following beta-oxidation of the side chain results in the formation of 2-naphthoyl-CoA and succinate (3, 53). The first three enzyme reactions of this pathway have been measured in vitro (3, 53). Recently, Musat et al. (48) identified the gene coding for the α-subunit of a putative naphthyl-2-methyl-succinate synthase (nmsA) in 2-methylnaphthalene-grown bacterial cultures. The molecular composition of the nmsA gene is analogous to that of the benzylsuccinate synthase α-subunit gene (bssA). The Bss enzyme is a well-investigated close homolog of Nms, catalyzing fumarate addition in the initial reaction of anaerobic toluene degradation (34, 40). Based on findings from comparative sequence studies, glycine radical-catalyzed fumarate addition has been shown to be a widely distributed initial reaction mechanism for anaerobic hydrocarbon degradation involving toluene and 2-methylnaphthalene, n-alkanes (12, 13, 25, 51), m-xylene (33), m- and p-cresols (9), and ethylbenzene (32).Open in a separate windowFIG. 1.Proposed pathway for anaerobic 2-methylnaphthalene degradation and reductive dearomatization of 2-naphthoyl-CoA (3, 4, 53). Genes found in the N47 genome encode the following enzymes (shown in gray boxes): NmsABC, naphthyl-2-methyl-succinate synthase; BnsEF, naphthyl-2-methyl-succinate CoA transferase; BnsG, naphthyl-2-methyl-succinyl-CoA dehydrogenase; BnsH, naphthyl-2-methylene-succinyl-CoA hydratase; BnsCD, naphthyl-2-hydroxymethyl-succinyl-CoA dehydrogenase; BnsAB, naphthyl-2-oxomethyl-succinyl-CoA thiolase; and NcrABCD, 2-naphthoyl-CoA reductase. The position of the double bond is not known for octahydro-2-naphthoyl-CoA. COSCoA, thioester of CoA and the respective carboxyl group.In a process analogous to the anaerobic benzoyl-CoA degradation pathway (7), 2-naphthoyl-CoA is subjected to aromatic ring reduction by a putative naphthoyl-CoA reductase, probably generating 5,6,7,8-tetrahydro-naphthoyl-CoA and further octahydro-2-naphthoic acid (4, 46). In the subsequent reactions, the ring system should be thiolytically cleaved and subjected to beta-oxidation, leading to the formation of acetyl-CoA and CO₂.In contrast to the first enzymatic reaction in the degradation of methylated aromatics, the first enzymatic reaction in anaerobic degradation of unsubstituted aromatic compounds such as naphthalene is still unresolved. In order to determine the initial activation reaction of anaerobic naphthalene degradation, studies based on the analysis of metabolites have been performed. Zhang and Young (66) observed the incorporation of ¹³C-labeled bicarbonate from the buffer into the carboxyl group of 2-naphthoic acid, hypothesizing that carboxylation is the initial activation reaction of anaerobic naphthalene degradation in the culture studied. Recently, Safinowski and Meckenstock (54) identified the deuterated metabolites naphthyl-2-methyl-succinate and naphthyl-2-methylene-succinate, which are exclusive intermediates of anaerobic 2-methylnaphthalene degradation, in the enrichment culture N47 when the culture was cultivated on fully deuterated naphthalene. Moreover, specific enzyme activities of the anaerobic 2-methylnaphtahlene degradation pathway have been detected in naphthalene-grown cells (54). Therefore, methylation of naphthalene to yield 2-methylnaphthalene as the initial activation reaction and subsequent degradation via the 2-methylnaphthalene pathway were proposed for this bacterial culture. The elucidation of 2-methylnaphthalene degradation may therefore reveal an important part of the naphthalene degradation pathway. However, Musat et al. (48) questioned methylation as the first reaction in naphthalene degradation for their marine naphthalene-degrading deltaproteobacterial NaphS strains.Whereas molecular components involved in anaerobic degradation of monoaromatic hydrocarbons are well known, knowledge about genes and enzymes involved in anaerobic PAH degradation is still missing (14). Here, we provide the first results of a whole-proteome- and whole-genome-based investigation of the sulfate-reducing enrichment culture N47 degrading naphthalene and 2-methylnaphtalene. We have identified some gene clusters encoding enzymes involved in 2-methylnaphthalene degradation, 2-naphthoyl-CoA dearomatization, and subsequent ring cleavage reactions in 2-methylnaphthalene-grown N47 cells.     

Molecular Profiling of the Lateral Habenula in a Rat Model of Depression

Objective

This study systematically investigated the effect of chronic mild stress and response to antidepressant treatment in the lateral habenula at the whole genome level.

Methods

Rat whole genome expression chips (Affymetrix) were used to detect gene expression regulations in the lateral habenula of rats subjected to chronic mild stress (mild stressors exchanged twice a day for 8 weeks). Some rats received antidepressant treatment during fifth to eights week of CMS. The lateral habenula gene expression profile was studied through the gene ontology and signal pathway analyses using bioinformatics. Real-time quantitative polymerase chain reaction (RT-PCR) was used to verify the microarray results and determine the expression of the Fcrla, Eif3k, Sec3l1, Ubr5, Abca8a, Ankrd49, Cyp2j10, Frs3, Syn2, and Znf503 genes in the lateral habenula tissue.

Results

In particular we found that stress and antidepressant treatment affected intracellular cascades like growth factor receptor signaling, G-protein-coupled receptor signaling, and Wnt signaling – processes involved in the neuroplastic changes observed during the progression of depression and antidepressant treatment.

Conclusion

The present study suggests an important role of the lateral habenula in the development of depression-like conditions and correlates to previous studies demonstrating a significant role of the lateral habenula in depressive-like conditions and antidepressant treatment.     

Anaerobic degradation of naphthalene and 2-methylnaphthalene by strains of marine sulfate-reducing bacteria

The anaerobic biodegradation of naphthalene, an aromatic hydrocarbon in tar and petroleum, has been repeatedly observed in environments but scarcely in pure cultures. To further explore the relationships and physiology of anaerobic naphthalene-degrading microorganisms, sulfate-reducing bacteria (SRB) were enriched from a Mediterranean sediment with added naphthalene. Two strains (NaphS3, NaphS6) with oval cells were isolated which showed naphthalene-dependent sulfate reduction. According to 16S rRNA gene sequences, both strains were Deltaproteobacteria and closely related to each other and to a previously described naphthalene-degrading sulfate-reducing strain (NaphS2) from a North Sea habitat. Other close relatives were SRB able to degrade alkylbenzenes, and phylotypes enriched anaerobically with benzene. If in adaptation experiments the three naphthalene-grown strains were exposed to 2-methylnaphthalene, this compound was utilized after a pronounced lag phase, indicating that naphthalene did not induce the capacity for 2-methylnaphthalene degradation. Comparative denaturing gel electrophoresis of cells grown with naphthalene or 2-methylnaphthalene revealed a striking protein band which was only present upon growth with the latter substrate. Peptide sequences from this band perfectly matched those of a protein predicted from genomic libraries of the strains. Sequence similarity (50% identity) of the predicted protein to the large subunit of the toluene-activating enzyme (benzylsuccinate synthase) from other anaerobic bacteria indicated that the detected protein is part of an analogous 2-methylnaphthalene-activating enzyme. The absence of this protein in naphthalene-grown cells together with the adaptation experiments as well as isotopic metabolite differentiation upon growth with a mixture of d(8)-naphthalene and unlabelled 2-methylnaphthalene suggest that the marine strains do not metabolize naphthalene by initial methylation via 2-methylnaphthalene, a previously suggested mechanism. The inability to utilize 1-naphthol or 2-naphthol also excludes these compounds as free intermediates. Results leave open the possibility of naphthalene carboxylation, another previously suggested activation mechanism.     

Nutrient reserves may allow for genome size increase: evidence from comparison of geophytes and their sister non-geophytic relatives

Background and Aims

The genome size of an organism is determined by its capacity to tolerate genome expansion, given the species'' life strategy and the limits of a particular environment, and the ability for retrotransposon suppression and/or removal. In some giant-genomed bulb geophytes, this tolerance is explained by their ability to pre-divide cells in the dormant stages or by the selective advantage of larger cells in the rapid growth of their fleshy body. In this study, a test shows that the tendency for genome size expansion is a more universal feature of geophytes, and is a subject in need of more general consideration.

Methods

Differences in monoploid genome sizes were compared using standardized phylogenetically independent contrasts in 47 sister pairs of geophytic and non-geophytic taxa sampled across all the angiosperms. The genome sizes of 96 species were adopted from the literature and 53 species were newly measured using flow cytometry with propidium iodide staining.

Key Results

The geophytes showed increased genome sizes compared with their non-geophytic relatives, regardless of the storage organ type and regardless of whether or not vernal geophytes, polyploids or annuals were included in the analyses.

Conclusions

The universal tendency of geophytes to possess a higher genome size suggests the presence of a universal mechanism allowing for genome expansion. It is assumed that this is primarily due to the nutrient and energetic independence of geophytes perhaps allowing continuous synthesis of DNA, which is known to proceed in the extreme cases of vernal geophytes even in dormant stages. This independence may also be assumed as a reason for allowing large genomes in some parasitic plants, as well as the nutrient limitation of small genomes of carnivorous plants.     

Horizontal transfer of a non-autonomous Helitron among insect and viral genomes

Background

The movement of mobile elements among species by horizontal transposon transfer (HTT) influences the evolution of genomes through the modification of structure and function. Helitrons are a relatively new lineage of DNA-based (class II) transposable elements (TEs) that propagate by rolling-circle replication, and are capable of acquiring host DNA. The rapid spread of Helitrons among animal lineages by HTT is facilitated by shuttling in viral particles or by unknown mechanisms mediated by close organism associations (e.g. between hosts and parasites).

Results

A non-autonomous Helitron independently annotated as BmHel-2 from Bombyx mori and the MITE01 element from Ostrinia nubilalis was predicted in the genomes of 24 species in the insect Order Lepidoptera. Integrated Helitrons retained ≥ 65% sequence identity over a 250 bp consensus, and were predicted to retain secondary structures inclusive of a 3′-hairpin and a 5′-subterminal inverted repeat. Highly similar Hel-2 copies were predicted in the genomes of insects and associated viruses, which along with a previous documented case of real-time virus-insect cell line transposition suggests that this Helitron has likely propagated by HTT.

Conclusions

These findings provide evidence that insect virus may mediate the HTT of Helitron-like TEs. This movement may facilitate the shuttling of DNA elements among insect genomes. Further sampling is required to determine the putative role of HTT in insect genome evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1318-6) contains supplementary material, which is available to authorized users.     

Modeling the fitness consequences of a cyanophage-encoded photosynthesis gene

Background

Phages infecting marine picocyanobacteria often carry a psbA gene, which encodes a homolog to the photosynthetic reaction center protein, D1. Host encoded D1 decays during phage infection in the light. Phage encoded D1 may help to maintain photosynthesis during the lytic cycle, which in turn could bolster the production of deoxynucleoside triphosphates (dNTPs) for phage genome replication.

Methodology / Principal Findings

To explore the consequences to a phage of encoding and expressing psbA, we derive a simple model of infection for a cyanophage/host pair — cyanophage P-SSP7 and Prochlorococcus MED4— for which pertinent laboratory data are available. We first use the model to describe phage genome replication and the kinetics of psbA expression by host and phage. We then examine the contribution of phage psbA expression to phage genome replication under constant low irradiance (25 µE m⁻² s⁻¹). We predict that while phage psbA expression could lead to an increase in the number of phage genomes produced during a lytic cycle of between 2.5 and 4.5% (depending on parameter values), this advantage can be nearly negated by the cost of psbA in elongating the phage genome. Under higher irradiance conditions that promote D1 degradation, however, phage psbA confers a greater advantage to phage genome replication.

Conclusions / Significance

These analyses illustrate how psbA may benefit phage in the dynamic ocean surface mixed layer.     

The Mitochondrial Genome of an Aquatic Plant,Spirodela polyrhiza

Background

Spirodela polyrhiza is a species of the order Alismatales, which represent the basal lineage of monocots with more ancestral features than the Poales. Its complete sequence of the mitochondrial (mt) genome could provide clues for the understanding of the evolution of mt genomes in plant.

Methods

Spirodela polyrhiza mt genome was sequenced from total genomic DNA without physical separation of chloroplast and nuclear DNA using the SOLiD platform. Using a genome copy number sensitive assembly algorithm, the mt genome was successfully assembled. Gap closure and accuracy was determined with PCR products sequenced with the dideoxy method.

Conclusions

This is the most compact monocot mitochondrial genome with 228,493 bp. A total of 57 genes encode 35 known proteins, 3 ribosomal RNAs, and 19 tRNAs that recognize 15 amino acids. There are about 600 RNA editing sites predicted and three lineage specific protein-coding-gene losses. The mitochondrial genes, pseudogenes, and other hypothetical genes (ORFs) cover 71,783 bp (31.0%) of the genome. Imported plastid DNA accounts for an additional 9,295 bp (4.1%) of the mitochondrial DNA. Absence of transposable element sequences suggests that very few nuclear sequences have migrated into Spirodela mtDNA. Phylogenetic analysis of conserved protein-coding genes suggests that Spirodela shares the common ancestor with other monocots, but there is no obvious synteny between Spirodela and rice mtDNAs. After eliminating genes, introns, ORFs, and plastid-derived DNA, nearly four-fifths of the Spirodela mitochondrial genome is of unknown origin and function. Although it contains a similar chloroplast DNA content and range of RNA editing as other monocots, it is void of nuclear insertions, active gene loss, and comprises large regions of sequences of unknown origin in non-coding regions. Moreover, the lack of synteny with known mitochondrial genomic sequences shed new light on the early evolution of monocot mitochondrial genomes.     

A novel mitochondrial genome architecture in thrips (Insecta: Thysanoptera): extreme size asymmetry among chromosomes and possible recent control region duplication

Background

Multipartite mitochondrial genomes are very rare in animals but have been found previously in two insect orders with highly rearranged genomes, the Phthiraptera (parasitic lice), and the Psocoptera (booklice/barklice).

Results

We provide the first report of a multipartite mitochondrial genome architecture in a third order with highly rearranged genomes: Thysanoptera (thrips). We sequenced the complete mitochondrial genomes of two divergent members of the Scirtothrips dorsalis cryptic species complex. The East Asia 1 species has the single circular chromosome common to animals while the South Asia 1 species has a genome consisting of two circular chromosomes. The fragmented South Asia 1 genome exhibits extreme chromosome size asymmetry with the majority of genes on the large, 14.28 kb, chromosome and only nad6 and trnC on the 0.92 kb mini-circle chromosome. This genome also features paralogous control regions with high similarity suggesting a very recent origin of the nad6 mini-circle chromosome in the South Asia 1 cryptic species.

Conclusions

Thysanoptera, along with the other minor paraenopteran insect orders should be considered models for rapid mitochondrial genome evolution, including fragmentation. Continued use of these models will facilitate a greater understanding of recombination and other mitochondrial genome evolutionary processes across eukaryotes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1672-4) contains supplementary material, which is available to authorized users.     

High correlation between genotypes and phenotypes of environmental bacteria Comamonas testosteroni strains

Background

Members of Comamonas testosteroni are environmental microorganisms that are usually found in polluted environment samples. They utilize steroids and aromatic compounds but rarely sugars, and show resistance to multiple heavy metals and multiple drugs. However, comprehensive genomic analysis among the C. testosteroni strains is lacked.

Results

To understand the genome bases of the features of C. testosteroni, we sequenced 10 strains of this species and analyzed them together with other related published genome sequences. The results revealed that: 1) the strains of C. testosteroni have genome sizes ranging from 5.1 to 6.0 Mb and G + C contents ranging from 61.1% to 61.8%. The pan-genome contained 10,165 gene families and the core genome contained 3,599 gene families. Heap’s law analysis indicated that the pan-genome of C. testosteroni may be open (α = 0.639); 2) by analyzing 31 phenotypes of 11 available C. testosteroni strains, 99.4% of the genotypes (putative genes) were found to be correlated to the phenotypes, indicating a high correlation between phenotypes and genotypes; 3) gene clusters for nitrate reduction, steroids degradation and metal and multi-drug resistance were found and were highly conserved among all the genomes of this species; 4) the genome similarity of C. testosteroni may be related to the geographical distances.

Conclusions

This work provided an overview on the genomes of C. testosteroni and new genome resources that would accelerate the further investigations of this species. Importantly, this work focused on the analysis of potential genetic determinants for the typical characters and found high correlation between the phenotypes and their corresponding genotypes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1314-x) contains supplementary material, which is available to authorized users.     

Characterization of Transposable Elements in the Ectomycorrhizal Fungus Laccaria bicolor

Background

The publicly available Laccaria bicolor genome sequence has provided a considerable genomic resource allowing systematic identification of transposable elements (TEs) in this symbiotic ectomycorrhizal fungus. Using a TE-specific annotation pipeline we have characterized and analyzed TEs in the L. bicolor S238N-H82 genome.

Methodology/Principal Findings

TEs occupy 24% of the 60 Mb L. bicolor genome and represent 25,787 full-length and partial copy elements distributed within 171 families. The most abundant elements were the Copia-like. TEs are not randomly distributed across the genome, but are tightly nested or clustered. The majority of TEs exhibits signs of ancient transposition except some intact copies of terminal inverted repeats (TIRS), long terminal repeats (LTRs) and a large retrotransposon derivative (LARD) element. There were three main periods of TE expansion in L. bicolor: the first from 57 to 10 Mya, the second from 5 to 1 Mya and the most recent from 0.5 Mya ago until now. LTR retrotransposons are closely related to retrotransposons found in another basidiomycete, Coprinopsis cinerea.

Conclusions

This analysis 1) represents an initial characterization of TEs in the L. bicolor genome, 2) contributes to improve genome annotation and a greater understanding of the role TEs played in genome organization and evolution and 3) provides a valuable resource for future research on the genome evolution within the Laccaria genus.     

Mating system shifts and transposable element evolution in the plant genus Capsella

Background

Despite having predominately deleterious fitness effects, transposable elements (TEs) are major constituents of eukaryote genomes in general and of plant genomes in particular. Although the proportion of the genome made up of TEs varies at least four-fold across plants, the relative importance of the evolutionary forces shaping variation in TE abundance and distributions across taxa remains unclear. Under several theoretical models, mating system plays an important role in governing the evolutionary dynamics of TEs. Here, we use the recently sequenced Capsella rubella reference genome and short-read whole genome sequencing of multiple individuals to quantify abundance, genome distributions, and population frequencies of TEs in three recently diverged species of differing mating system, two self-compatible species (C. rubella and C. orientalis) and their self-incompatible outcrossing relative, C. grandiflora.

Results

We detect different dynamics of TE evolution in our two self-compatible species; C. rubella shows a small increase in transposon copy number, while C. orientalis shows a substantial decrease relative to C. grandiflora. The direction of this change in copy number is genome wide and consistent across transposon classes. For insertions near genes, however, we detect the highest abundances in C. grandiflora. Finally, we also find differences in the population frequency distributions across the three species.

Conclusion

Overall, our results suggest that the evolution of selfing may have different effects on TE evolution on a short and on a long timescale. Moreover, cross-species comparisons of transposon abundance are sensitive to reference genome bias, and efforts to control for this bias are key when making comparisons across species.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-602) contains supplementary material, which is available to authorized users.     

Evolutionary mode and functional divergence of vertebrate NMDA receptor subunit 2 genes

Background

Ionotropic glutamate receptors in the central nervous system play a major role in numerous brain functions including learning and memory in many vertebrate species. NR2 subunits have been regarded as rate-limiting molecules in controlling the optimal N-methyl-D-aspartate (NMDA) receptor''s coincidence-detection property and subsequent learning and memory function across multi-species. However, its evolutionary mode among vertebrate species remains unclear.

Results

With extensive analysis of phylogeny, exon structure, protein domain, paralogon and synteny, we demonstrated that two-round genome duplication generated quartet GRIN2 genes and the third-round fish-specific genome duplication generated extra copies of fish GRIN2 genes. In addition, in-depth investigation has enabled the identification of three novel genes, GRIN2C_Gg, GRIN2D-1_Ol and GRIN2D-2_Tr in the chicken, medaka and fugu genome, respectively. Furthermore, we showed functional divergence of NR2 genes mostly occurred at the first-round duplication, amino acid residues located at the N-terminal Lig_chan domain were responsible for type I functional divergence between these GRIN2 subfamilies and purifying selection has been the prominent natural pressure operating on these diversified GRIN2 genes.

Conclusion and Significance

These findings provide intriguing subjects for testing the 2R and 3R hypothesis and we expect it could provide new insights into the underlying evolution mechanisms of cognition in vertebrate.     

Genomic correlates of recombination rate and its variability across eight recombination maps in the western honey bee (Apis mellifera L.)

Background

Meiotic recombination has traditionally been explained based on the structural requirement to stabilize homologous chromosome pairs to ensure their proper meiotic segregation. Competing hypotheses seek to explain the emerging findings of significant heterogeneity in recombination rates within and between genomes, but intraspecific comparisons of genome-wide recombination patterns are rare. The honey bee (Apis mellifera) exhibits the highest rate of genomic recombination among multicellular animals with about five cross-over events per chromatid.

Results

Here, we present a comparative analysis of recombination rates across eight genetic linkage maps of the honey bee genome to investigate which genomic sequence features are correlated with recombination rate and with its variation across the eight data sets, ranging in average marker spacing ranging from 1 Mbp to 120 kbp. Overall, we found that GC content explained best the variation in local recombination rate along chromosomes at the analyzed 100 kbp scale. In contrast, variation among the different maps was correlated to the abundance of microsatellites and several specific tri- and tetra-nucleotides.

Conclusions

The combined evidence from eight medium-scale recombination maps of the honey bee genome suggests that recombination rate variation in this highly recombining genome might be due to the DNA configuration instead of distinct sequence motifs. However, more fine-scale analyses are needed. The empirical basis of eight differing genetic maps allowed for robust conclusions about the correlates of the local recombination rates and enabled the study of the relation between DNA features and variability in local recombination rates, which is particularly relevant in the honey bee genome with its exceptionally high recombination rate.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1281-2) contains supplementary material, which is available to authorized users.     

Differential Gene Repertoire in Mycobacterium ulcerans Identifies Candidate Genes for Patho-Adaptation

Background

Based on large genomic sequence polymorphisms, several haplotypes belonging to two major lineages of the human pathogen Mycobacterium ulcerans could be distinguished among patient isolates from various geographic origins. However, the biological relevance of insertional/deletional diversity is not understood.

Methodology

Using comparative genomics, we have investigated the genes located in regions of difference recently identified by DNA microarray based hybridisation analysis. The analysed regions of difference comprise ∼7% of the entire M. ulcerans genome.

Principal Findings

Several different mechanisms leading to loss of functional genes were identified, ranging from pseudogenization, caused by frame shift mutations or mobile genetic element interspersing, to large sequence polymorphisms. Four hot spot regions for genetic instability were unveiled. Altogether, 229 coding sequences were found to be differentially inactivated, constituting a repertoire of coding sequence variation in the rather monomorphic M. ulcerans.

Conclusions/Significance

The differential gene inactivation patterns associated with the M. ulcerans haplotypes identified candidate genes that may confer enhanced adaptation upon ablation of expression. A number of gene conversions confined to the classical lineage may contribute to particular virulence of this group comprising isolates from Africa and Australia. Identification of this spectrum of anti-virulence gene candidates expands our understanding of the pathogenicity and ecology of the emerging infectious disease Buruli ulcer.