RPB2 sequences reveal a close phylogenetic relationship between tetraploid Hordelymus and diploid Hordeum species in Triticeae (Poaceae)

The origin of Hordelymus genome has been debated for years, and no consensus conclusion was reached. In this study, we sequenced and analyzed the RPB2 (RNA polymerase subunit II) gene from Hordelymus europaeus (L.) Harz, and its potential diploid ancestor species those were suggested in previous studies. The focus of this study was to examine the phylogenetic relationship of Hordelymus genomes with its potential donor Hordeum, Psathyrostachys, and Taeniatherum species. Two distinguishable copies of sequences were obtained from H. europaeus. The obvious difference between the two copies of sequences is a 24 bp indel (insertion/deletion). Phylogenetic analysis showed a strong affinity between Hordeum genome and Hordelymus with 85% bootstrap support. These results suggested that one genome in tetraploid H. europaeus closely related to the genome in Hordeum species. Another genome in H. europaeus is sister to the genomes in Triticeae species examined here, which corresponds well with the recently published EF-G data. No obvious relationship was found between Hordelymus and either Ta genome donor, Taeniatherum caput-medusae or Ns genome donor, Psathyrostachys juncea. Our data does not support the presence of Ta and Ns genome in H. europaeus, and further confirms that H. europaeus is allopolyploid.     

Distinct origin of the Y and St genome in Elymus species: evidence from the analysis of a large sample of St genome species using two nuclear genes

Background

Previous cytological and single copy nuclear genes data suggested the St and Y genome in the StY-genomic Elymus species originated from different donors: the St from a diploid species in Pseudoroegneria and the Y from an unknown diploid species, which are now extinct or undiscovered. However, ITS data suggested that the Y and St genome shared the same progenitor although rather few St genome species were studied. In a recent analysis of many samples of St genome species Pseudoroegneria spicata (Pursh) À. Löve suggested that one accession of P. spicata species was the most likely donor of the Y genome. The present study tested whether intraspecific variation during sampling could affect the outcome of analyses to determining the origin of Y genome in allotetraploid StY species. We also explored the evolutionary dynamics of these species.

Methodology/Principal Findings

Two single copy nuclear genes, the second largest subunit of RNA polymerase II (RPB2) and the translation elongation factor G (EF-G) sequences from 58 accessions of Pseudoroegneria and Elymus species, together with those from Hordeum (H), Agropyron (P), Australopyrum (W), Lophopyrum (E^e), Thinopyrum (E^a), Thinopyrum (E^b), and Dasypyrum (V) were analyzed using maximum parsimony, maximum likelihood and Bayesian methods. Sequence comparisons among all these genomes revealed that the St and Y genomes are relatively dissimilar. Extensive sequence variations have been detected not only between the sequences from St and Y genome, but also among the sequences from diploid St genome species. Phylogenetic analyses separated the Y sequences from the St sequences.

Conclusions/Significance

Our results confirmed that St and Y genome in Elymus species have originated from different donors, and demonstrated that intraspecific variation does not affect the identification of genome origin in polyploids. Moreover, sequence data showed evidence to support the suggestion of the genome convergent evolution in allopolyploid StY genome species.     

The origin of the H, St, W, and Y genomes in allotetraploid species of Elymus L. and Stenostachys Turcz. (Poaceae: Triticeae)

Based on sequences from two single-copy nuclear genes (DMC1 and EF-G), four plastid genes (rbcL, rpoA, matK, and ndhF), and one mitochondrial gene (coxII), we investigate the origin of the H, St, W, and Y genomes in four allotetraploid species of Elymus and two allotetraploid species of Stenostachys. Despite significant incongruence between the two nuclear genes and between the nuclear and organelle data partitions, individual and combined analyses of the data partitions unequivocally show that the St and H genomes of the tetraploid American species of Elymus are derived from Pseudoroegneria and Hordeum, respectively, with Pseudoroegneria serving as the female parent, and that the H and W genomes of Stenostachys are derived from Hordeum and Australopyrum, respectively, with Hordeum serving as the female parent. The analyses equally clearly demonstrate that the St genome of the tetraploid Asiatic Elymus species is derived from Pseudoroegneria, with the latter serving as the female parent, but the relationship of the Y genome is less clear. Individual analyses of the nuclear genes provide conflicting results, but combined analysis of all data suggests a sister group relationship to Heteranthelium, albeit without any jackknife support.     

Phylogenetic relationships and reticulation among Asian Elymus (Poaceae) allotetraploids: Analyses of three nuclear gene trees

This phylogenetic study focuses on a subset of the species in Elymus—specifically, the endemic Asian tetraploids presumed to combine the St genome from Pseudoroegneria with the Y genome from an unknown donor. The primary goals were to (1) determine whether the St and Y genomes are derived from phylogenetically distinct donors; (2) identify the closest relative, and potentially the likely donor, of the Y genome; and (3) interpret variation among StStYY species in terms of multiple origins and/or introgression. The goals were addressed using phylogenetic analyses of sequences from three low-copy nuclear genes: phosphoenolpyruvate carboxylase, β-amylase, and granule-bound starch synthase I. Data sets include 16 StStYY individuals representing nine species, along with a broad sample of representatives from most of the monogenomic (i.e., non-allopolyploid) genera in the tribe. To briefly summarize the results: (1) the data clearly support an allopolyploid origin for the Asian tetraploids, involving two distinct donors; (2) the Y genome was contributed by a single donor, or multiple closely-related donors; (3) the phylogenetic position of the Elymus Y genome varies among the three trees and its position is not strongly supported, so the identity of the donor remains a mystery; and (4) conflicts among the gene trees with regard to the St-genome sequences suggest introgression involving both Elymus and Pseudoroegneria.     

Phylogenetic analysis in some Hordeum species (Triticeae; Poaceae) based on two single-copy nuclear genes encoding acetyl-CoA carboxylase

We investigate the phylogenetic relationship, and evolutionary history of 18 diploid and polyploid Hordeum species including 22 taxa based on two single-copy nuclear ACC1 and ACC2 genes using maximum parsimony, maximum likelihood, and Bayesian inference. The results of molecular phylogenetic analysis demonstrated genetic relationships among taxa and origin of polyploids. Our phylogenetic analyses revealed a clear alloploid origin of Hordeum capense, with Eurasian Hordeum marinum subsp. gussoneanum as the Xa genome donor and diploid Asian Hordeum roshevitzii as the H genome donor. The formation of hexaploid Hordeum lechleri likely involves hybridization between tetraploid Hordeum brachyantherum subsp. brachyantherum and a diploid possessing the I genome. The Acc1 and Acc2 gene data analyses suggested that Siberian Hordeum bogdanii might have be the common ancestor of the diploid New World Hordeum species. Perennial diploid South American species, Hordeum comosum was the first-diverging group within the clade of diploid American species in the analyses.     

Molecular Evidence for the Hybrid Origin of Ilex dabieshanensis (Aquifoliaceae)

Ilex, the largest genus of dioecious woody plants, is a good study system to assess the role of hybridization in speciation and evolution. Ilex dabieshanensis, a tree endemic to Dabieshan Mountains region, was initially described as a new species. Based on morphological intermediacy and sympatric distribution with its putative parental species, I. cornuta and I. latifolia, we proposed it as a natural hybrid between them. In this study, we sequenced one chloroplast intergenic spacer (trnH-psbA) and two nuclear genes (gapC and nepGS) in I. dabieshanensis and its putative parental species to test the hybrid origin hypothesis. Our results showed that there were one to two differentially fixed sequence differences between I. cornuta and I. latifolia at the two nuclear genes. Twelve of the 14 individuals of I. dabieshanensis exhibited additivity in chromatograms on these differentially fixed sites at both nuclear genes, and the remaining two exhibited additivity in chromatograms on the fixed site at only the nepGS gene. Except one haplotype of I. cornuta at the nepGS gene, all of the haplotypes of I. cornuta at the two nuclear genes were well separated from those of I. latifolia, and most haplotypes of I. dabieshanensis were shared with those of I. cornuta and I. latifolia. Phylogenetic analysis of these haplotypes was largely consistent with haplotype network analysis. I. cornuta and I. latifolia differed by two nucleotide substitutions in the chloroplast intergenic spacer, and 12 individuals of I. dabieshanensis had the same sequences as I. latifolia, while the remaining two were identical with I. cornuta. The molecular data provide convincing evidence for the hybrid origin of I. dabieshanensis and asymmetrical direction of hybridization. One haplotype of I. cornuta at the nepGS gene was nested with those of I. latifolia, indicating introgression to I. cornuta.     

Complete Mitochondrial Genome of the Medicinal Mushroom Ganoderma lucidum

Ganoderma lucidum is one of the well-known medicinal basidiomycetes worldwide. The mitochondrion, referred to as the second genome, is an organelle found in most eukaryotic cells and participates in critical cellular functions. Elucidating the structure and function of this genome is important to understand completely the genetic contents of G. lucidum. In this study, we assembled the mitochondrial genome of G. lucidum and analyzed the differential expressions of its encoded genes across three developmental stages. The mitochondrial genome is a typical circular DNA molecule of 60,630 bp with a GC content of 26.67%. Genome annotation identified genes that encode 15 conserved proteins, 27 tRNAs, small and large rRNAs, four homing endonucleases, and two hypothetical proteins. Except for genes encoding trnW and two hypothetical proteins, all genes were located on the positive strand. For the repeat structure analysis, eight forward, two inverted, and three tandem repeats were detected. A pair of fragments with a total length around 5.5 kb was found in both the nuclear and mitochondrial genomes, which suggests the possible transfer of DNA sequences between two genomes. RNA-Seq data for samples derived from three stages, namely, mycelia, primordia, and fruiting bodies, were mapped to the mitochondrial genome and qualified. The protein-coding genes were expressed higher in mycelia or primordial stages compared with those in the fruiting bodies. The rRNA abundances were significantly higher in all three stages. Two regions were transcribed but did not contain any identified protein or tRNA genes. Furthermore, three RNA-editing sites were detected. Genome synteny analysis showed that significant genome rearrangements occurred in the mitochondrial genomes. This study provides valuable information on the gene contents of the mitochondrial genome and their differential expressions at various developmental stages of G. lucidum. The results contribute to the understanding of the functions and evolution of fungal mitochondrial DNA.     

Evolutionary History of Triticum petropavlovskyi Udacz. et Migusch. Inferred from the Sequences of the 3-Phosphoglycerate Kinase Gene

Single- and low-copy genes are less likely to be subject to concerted evolution. Thus, they are appropriate tools to study the origin and evolution of polyploidy plant taxa. The plastid 3-phosphoglycerate kinase gene (Pgk-1) sequences from 44 accessions of Triticum and Aegilops, representing diploid, tetraploid, and hexaploid wheats, were used to estimate the origin of Triticum petropavlovskyi. Our phylogenetic analysis was carried out on exon+intron, exon and intron sequences, using maximum likelihood, Bayesian inference and haplotype networking. We found the D genome sequences of Pgk-1 genes from T. petropavlovskyi are similar to the D genome orthologs in T. aestivum, while their relationship with Ae. tauschii is more distant. The A genome sequences of T. petropavlovskyi group with those of T. polonicum, but its Pgk-1 B genome sequences to some extent diverge from those of other species of Triticum. Our data do not support for the origin of T. petropavlovskyi either as an independent allopolyploidization event between Ae. tauschii and T. polonicum, or as a monomendelian mutation in T. aestivum. We suggest that T. petropavlovskyi originated via spontaneous introgression from T. polonicum into T. aestivum. The dating of this introgression indicates an age of 0.78 million years; a further mutation event concerning the B genome occurred 0.69 million years ago.     

Origin and Reticulate Evolutionary Process of Wheatgrass Elymus trachycaulus (Triticeae: Poaceae)

To study origin and evolutionary dynamics of tetraploid Elymus trachycaulus that has been cytologically defined as containing StH genomes, thirteen accessions of E. trachycaulus were analyzed using two low-copy nuclear gene Pepc (phosphoenolpyruvate carboxylase) and Rpb2 (the second largest subunit of RNA polymerase II), and one chloroplast region trnL–trnF (spacer between the tRNA Leu (UAA) gene and the tRNA-Phe (GAA) gene). Our chloroplast data indicated that Pseudoroegneria (St genome) was the maternal donor of E. trachycaulus. Rpb2 data indicated that the St genome in E. trachycaulus was originated from either P. strigosa, P. stipifolia, P. spicata or P. geniculate. The Hordeum (H genome)-like sequences of E. trachycaulus are polyphyletic in the Pepc tree, suggesting that the H genome in E. trachycaulus was contributed by multiple sources, whether due to multiple origins or introgression resulting from subsequent hybridization. Failure to recovering St copy of Pepc sequence in most accessions of E. trachycaulus might be caused by genome convergent evolution in allopolyploids. Multiple copies of H-like Pepc sequence from each accession with relative large deletions and insertions might be caused by either instability of Pepc sequence in H- genome or incomplete concerted evolution. Our results highlighted complex evolutionary history of E. trachycaulus.     

RECONSTRUCTING THE COMPLEX EVOLUTIONARY ORIGIN OF WILD ALLOPOLYPLOID TOBACCOS (NICOTIANA SECTION SUAVEOLENTES)

Nicotiana (Solanaceae) provides an ideal system for understanding polyploidization, a pervasive and powerful evolutionary force in plants, as this genus contains several groups of allotetraploids that formed at different times from different diploid progenitors. However, the parental lineages of the largest group of allotetraploids, Nicotiana section Suaveolentes, have been problematic to identify. Using data from four regions of three low‐copy nuclear genes, nuclear ribosomal DNA, and regions of the plastid genome, we have reconstructed the evolutionary origin of sect. Suaveolentes and identified the most likely diploid progenitors by using a combination of gene trees and network approaches to uncover the most strongly supported evidence of species relationships. Our analyses best support a scenario where a member of the sect. Sylvestres lineage acted as the paternal progenitor and a member of either sect. Petunioides or sect. Noctiflorae that also contained introgressed DNA from the other, or a hypothetical hybrid species between these two sections, was the maternal progenitor. Nicotiana exemplifies many of the factors that can complicate the reconstruction of polyploid evolutionary history and highlights how reticulate evolution at the diploid level can add even greater complexity to allopolyploid genomes.     

Phylogenetic analysis of PISTILLATA sequences in Neillia (Rosaceae)

Putative PISTILLATA genes were generated in nine species of Neillia, to examine the phylogenetic relationships among the species, and to test the hypothesis of hybrid origin within the genus. The PI genes determined in Neillia have two introns in the I-box region, which is consistent with PI genes in other Rosaceae. Phylogenetic analyses of the I-box region, including the introns, indicated that the species formerly classified in Stephanandra were nested within Neillia, supporting the taxonomic merger of the two genera. The PI data do not have a sufficiently strong signal to reject the hypothesis that Stephanandra is a hybrid in origin. The PI data, in conjunction with nuclear LEAFY, ribosomal DNA, and chloroplast DNA data, suggest that N. affinis might have been derived from hybridization between N. thibetica and N. gracilis. The phylogenetic position of N. affinis in the N. thibetica clade is supported by the PI and rDNA data, whereas N. affinis is also supported as a sister to N. gracilis in the LEAFY and cpDNA data. The pattern of phylogenetic placements of N. affinis in two different clades in two different sets of data suggests that the genome of the species might be comprised of a combination of the putative parental species.     

MOLECULAR EVIDENCE FOR HOMOPLOID RETICULATE EVOLUTION AMONG AUSTRALIAN SPECIES OF GOSSYPIUM

Interspecific hybridization and introgression are important evolutionary processes in plants, but their full significance with respect to speciation at the diploid level remains unresolved. In this study, molecular markers from the plastid and nuclear genomes were used to document an unusual evolutionary history of Gossypium bickii Prokh. (Malvaceae). This species is one of three morphologically similar Australian cottons (along with G. austrate F. Muell. and G. nelsonii Fryx.) in section Hibiscoidea. In contrast to expectations based on previous morphological data, cladistic analysis of maternally inherited cpDNA restriction site mutations unites G. bickii with G. sturtianum J. H. Willis, a morphologically distant species in a different taxonomic section (Sturtia). Few restriction site mutations distinguish the plastomes of G. bickii and G. sturtianum, but these two cpDNAs are differentiated from those of G. australe and G. nelsonii by a minimum of 33 mutations (out of 640 sites scored). These two highly distinct clades are not supported by phylogenetic analyses of allozyme markers (from 58 populations) and restriction site mutations in nuclear ribosomal DNAs. Rather, phylogenies based on 83 nuclear markers indicate that G. bickii shares a more recent common ancestor with G. australe and G. nelsonii than it does with G. sturtianum. We suggest that the striking discrepancy between independent molecular phylogenies from two different genomes indicates a biphyletic ancestry of G. bickii. Our preferred hypothesis involves an ancient hybridization, in which G. sturtianum, or a similar species, served as the maternal parent with a paternal donor from the lineage leading to G. australe and G. nelsonii. Because we detected no G. sturtianum nuclear genes in G. bickii, we suggest that the nuclear genomic contribution of the maternal parent was subsequently eliminated from the hybrid or its descendent maternal lineage. Several possible mechanisms of cytoplasm transfer are suggested, including repeated backcrossing of the hybrid, as female, into the paternal donor lineage, selection against recombinant nuclear genomes and a form of apomixis known as semigamy. This example, and several others in Gossypium as well as other genera, attest to the evolutionary possibility of interspecific cytoplasmic transfer, and perhaps the origin of diploid species via reticulate speciation. In addition, this study offers an example of natural cytoplasmic introgression without long-term survival of nuclear genes from the maternal progenitor.     

Next generation sequencing from Hepatozoon canis (Apicomplexa: Coccidia: Adeleorina): Complete apicoplast genome and multiple mitochondrion-associated sequences

Extrachromosomal genomes of the adeleorinid parasite Hepatozoon canis infecting an Israeli dog were investigated using next-generation and standard sequencing technologies. A complete apicoplast genome and several mitochondrion-associated sequences were generated. The apicoplast genome (31,869?bp) possessed two copies of both large subunit (23S) and small subunit (16S) ribosomal RNA genes (rDNA) within an inverted repeat region, as well as 22 protein-coding sequences, 25 transfer RNA genes (tDNA) and seven open reading frames of unknown function. Although circular-mapping, the apicoplast genome was physically linear according to next-generation data. Unlike other apicoplast genomes, genes encoding ribosomal protein S19 and tDNAs for alanine, aspartic acid, histidine, threonine and valine were not identified. No complete mitochondrial genome was recovered using next-generation data or directed PCR amplifications. Eight mitochondrion-associated (215–3523?bp) contigs assembled from next-generation data encoded a complete cytochrome c oxidase subunit I coding sequence, a complete cytochrome c oxidase subunit III coding sequence, two complete cytochrome B coding sequences, a non-coding, pseudogene for cytochrome B and multiple fragmented mitochondrial rDNA genes (SSUA, SSUB, SSUD, LSUC, LSUG, RNA6, RNA10, RNA14, RNA18). The paucity of NGS reads generating each of the mitochondrion-like sequences suggested that a complete mitochondrial genome at typically high copy number was absent in H. canis. In contrast, the complete nuclear rDNA unit sequence of H. canis (18S rDNA to 28S rDNA, 6977?bp) had >1000-fold next-generation coverage. Multiple divergent (from 93.6% to 99.9% pairwise identities) nuclear 18S rDNA contigs were generated (three types with 10 subtypes total). To our knowledge this is the first apicoplast genome sequenced from any adeleorinid coccidium and the first mitochondrion-associated sequences from this serious pathogen of wild and domestic canids. These newly generated sequences may provide useful genetic loci for high-resolution species-level genotyping that is currently impossible using existing nuclear rDNA targets.     

The phylogeny of the Anderson's White‐bellied Rat (Niviventer andersoni) based on complete mitochondrial genomes

The phylogenetic structure of the genus Niviventer has been studied based on several individual mitochondrial and nuclear genes, but the results seem to be inconsistent. In order to clarify the phylogeny of Niviventer, we sequenced the complete mitochondrial genome of white‐bellied rat (Niviventer andersoni of the family Muridae) by next‐generation sequencing. The 16,291 bp mitochondrial genome consists of 22 transfer RNA genes, 13 protein‐coding genes (PCGs), two ribosomal RNA genes, and one noncoding control region (D‐Loop). Phylogenetic analyses of the nucleotide sequences of all 13 PCGs, PCGs minus ND6, and the entire mitogenome sequence except for the D‐loop revealed well‐resolved topologies supporting that N. andersoni was clustered with N. excelsior forming a sister division with N. confucianus, which statistically rejected the hypothesis based on the tree of cytochrome b (cytb) gene that N. confucianus is sister to N. fulvescens. Our research provides the first annotated complete mitochondrial genome of N. andersoni, extending the understanding about taxonomy and mitogenomic evolution of the genus Niviventer.     

Phylogeny and evolutionary history of Leymus (Triticeae; Poaceae) based on a single-copy nuclear gene encoding plastid acetyl-CoA carboxylase

Background  

Single- and low- copy genes are less likely subject to concerted evolution, thus making themselves ideal tools for studying the origin and evolution of polyploid taxa. Leymus is a polyploid genus with a diverse array of morphology, ecology and distribution in Triticeae. The genomic constitution of Leymus was assigned as NsXm, where Ns was presumed to be originated from Psathyrostachys, while Xm represented a genome of unknown origin. In addition, little is known about the evolutionary history of Leymus. Here, we investigate the phylogenetic relationship, genome donor, and evolutionary history of Leymus based on a single-copy nuclear Acc1 gene.     

A unique tRNA gene family and a novel,highly expressed ORF in the mitochondrial genome of the silver-lip pearl oyster, Pinctada maxima (Bivalvia: Pteriidae)

Characteristics of mitochondrial (mt) DNA such as gene content and arrangement, as well as mt tRNA secondary structure, are frequently used in comparative genomic analyses because they provide valuable phylogenetic information. However, most analyses do not characterize the relationship of tRNA genes from the same mt genome and, in some cases, analyses overlook possible novel open reading frames (ORFs) when the 13 expected protein-coding genes are already annotated. In this study, we describe the sequence and characterization of the complete mt genome of the silver-lip pearl oyster, Pinctada maxima. The 16,994-bp mt genome contains the same 13 protein-coding genes (PCGs) and two ribosomal RNA genes typical of metazoans. The gene arrangement, however, is completely distinct from that of all other available bivalve mt genomes, and a unique tRNA gene family is observed in this genome. The unique tRNA gene family includes two trnS^− AGY and trnQ genes, a trnM isomerism, but it lacks trnS^− CUN. We also report the first clear evidence of alloacceptor tRNA gene recruitment (trnP → trnS^− AGY) in mollusks. In addition, a novel ORF (orfUR1) expressed at high levels is present in the mt genome of this pearl oyster. This gene contains a conserved domain, “Oxidored_q1_N”, which is a member of Complex I and thus may play an important role in key biological functions. Because orfUR1 has a very similar nucleotide composition and codon bias to that of other genes in this genome, we hypothesize that this gene may have been moved to the mt genome via gene transfer from the nuclear genome at an early stage of speciation of P. maxima, or it may have evolved as a result of gene duplication, followed by rapid sequence divergence. Lastly, a 319-bp region was identified as the possible control region (CR) even though it does not correspond to the longest non-coding region in the genome. Unlike other studies of mt genomes, this study compares the evolutionary patterns of all available bivalve mt tRNA and atp8 genes.     

Chromosomal-level genome assembly of the orchid tree Bauhinia variegata (Leguminosae; Cercidoideae) supports the allotetraploid origin hypothesis of Bauhinia

Cercidoideae, one of the six subfamilies of Leguminosae, contains one genus Cercis with its chromosome number 2n = 14 and all other genera with 2n = 28. An allotetraploid origin hypothesis for the common ancestor of non-Cercis genera in this subfamily has been proposed; however, no chromosome-level genomes from Cercidoideae have been available to test this hypothesis. Here, we conducted a chromosome-level genome assembly of Bauhinia variegata to test this hypothesis. The assembled genome is 326.4 Mb with the scaffold N50 of 22.1 Mb and contains 37,996 protein-coding genes. The Ks distribution between gene pairs in the syntenic regions indicates two whole-genome duplications (WGDs): one is B. variegata-specific, and the other is shared among core eudicots. Although Ks between gene pairs generated by the recent WGD in Bauhinia is greater than that between Bauhinia and Cercis, the WGD was not detected in Cercis, which can be explained by an accelerated evolutionary rate in Bauhinia after divergence from Cercis. Ks distribution and phylogenetic analysis for gene pairs generated by the recent WGD in Bauhinia and their corresponding orthologs in Cercis support the allopolyploidy origin hypothesis of Bauhinia. The genome of B. variegata also provides a genomic resource for dissecting genetic basis of its ornamental traits.     

Prevalence of laterally transferred Wolbachia genes in Japanese pine sawyer,Monochamus alternatus (Coleoptera: Cerambycidae)

Lateral gene transfer is any process in which an organism incorporates genetic material from a phylogenetically distant and reproductively isolated organism. A previous survey demonstrated that two Japanese populations of Monochamus alternatus Hope carried many Wolbachia genes that had been acquired by lateral gene transfer in their own genomes, but were not infected with Wolbachia. To understand the prevalence of Wolbachia infection and laterally transferred Wolbachia genes in this beetle, we performed a broader survey of natural populations covering two subspecies of M. alternatus from Japan (M. alternatus endai) and Taiwan (M. alternatus alternatus). We detected laterally transferred Wolbachia genes in all Japanese and Taiwanese populations of M. alternatus, but no Wolbachia infection in any population. In addition, we confirmed the absence of Wolbachia infection and of transferred Wolbachia genes in two Japanese populations of Monochamus saltuarius Gebler, which is a congeneric relative of M. alternatus. Our findings suggest that the Wolbachia endosymbiont as gene donor has disappeared from M. alternatus, and that the transfer of part of its genome to M. alternatus occurred in the ancestor of M. alternatus before the subspeciation event.     

Comparative analysis of mitochondrial genomes between the hau cytoplasmic male sterility (CMS) line and its iso-nuclear maintainer line in Brassica juncea to reveal the origin of the CMS-associated gene orf288

Background

Cytoplasmic male sterility (CMS) is not only important for exploiting heterosis in crop plants, but also as a model for investigating nuclear-cytoplasmic interaction. CMS may be caused by mutations, rearrangement or recombination in the mitochondrial genome. Understanding the mitochondrial genome is often the first and key step in unraveling the molecular and genetic basis of CMS in plants. Comparative analysis of the mitochondrial genome of the hau CMS line and its maintainer line in B. juneca (Brassica juncea) may help show the origin of the CMS-associated gene orf288.

Results

Through next-generation sequencing, the B. juncea hau CMS mitochondrial genome was assembled into a single, circular-mapping molecule that is 247,903 bp in size and 45.08% in GC content. In addition to the CMS associated gene orf288, the genome contains 35 protein-encoding genes, 3 rRNAs, 25 tRNA genes and 29 ORFs of unknown function. The mitochondrial genome sizes of the maintainer line and another normal type line “J163-4” are both 219,863 bp and with GC content at 45.23%. The maintainer line has 36 genes with protein products, 3 rRNAs, 22 tRNA genes and 31 unidentified ORFs. Comparative analysis the mitochondrial genomes of the hau CMS line and its maintainer line allowed us to develop specific markers to separate the two lines at the seedling stage. We also confirmed that different mitotypes coexist substoichiometrically in hau CMS lines and its maintainer lines in B. juncea. The number of repeats larger than 100 bp in the hau CMS line (16 repeats) are nearly twice of those found in the maintainer line (9 repeats). Phylogenetic analysis of the CMS-associated gene orf288 and four other homologous sequences in Brassicaceae show that orf288 was clearly different from orf263 in Brassica tournefortii despite of strong similarity.

Conclusion

The hau CMS mitochondrial genome was highly rearranged when compared with its iso-nuclear maintainer line mitochondrial genome. This study may be useful for studying the mechanism of natural CMS in B. juncea, performing comparative analysis on sequenced mitochondrial genomes in Brassicas, and uncovering the origin of the hau CMS mitotype and structural and evolutionary differences between different mitotypes.

Electronic supplementary material

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