Accelerated evolution of the mitochondrial genome in an alloplasmic line of durum wheat

Background

Wheat is an excellent plant species for nuclear mitochondrial interaction studies due to availability of large collection of alloplasmic lines. These lines exhibit different vegetative and physiological properties than their parents. To investigate the level of sequence changes introduced into the mitochondrial genome under the alloplasmic condition, three mitochondrial genomes of the Triticum-Aegilops species were sequenced: 1) durum alloplasmic line with the Ae. longissima cytoplasm that carries the T. turgidum nucleus designated as (lo) durum, 2) the cytoplasmic donor line, and 3) the nuclear donor line.

Results

The mitochondrial genome of the T. turgidum was 451,678 bp in length with high structural and nucleotide identity to the previously characterized T. aestivum genome. The assembled mitochondrial genome of the (lo) durum and the Ae. longissima were 431,959 bp and 399,005 bp in size, respectively. The high sequence coverage for all three genomes allowed analysis of heteroplasmy within each genome. The mitochondrial genome structure in the alloplasmic line was genetically distant from both maternal and paternal genomes. The alloplasmic durum and the Ae. longissima carry the same versions of atp6, nad6, rps19-p, cob and cox2 exon 2 which are different from the T. turgidum parent. Evidence of paternal leakage was also observed by analyzing nad9 and orf359 among all three lines. Nucleotide search identified a number of open reading frames, of which 27 were specific to the (lo) durum line.

Conclusions

Several heteroplasmic regions were observed within genes and intergenic regions of the mitochondrial genomes of all three lines. The number of rearrangements and nucleotide changes in the mitochondrial genome of the alloplasmic line that have occurred in less than half a century was significant considering the high sequence conservation between the T. turgidum and the T. aestivum that diverged from each other 10,000 years ago. We showed that the changes in genes were not limited to paternal leakage but were sufficiently significant to suggest that other mechanisms, such as recombination and mutation, were responsible. The newly formed ORFs, differences in gene sequences and copy numbers, heteroplasmy, and substoichiometric changes show the potential of the alloplasmic condition to accelerate evolution towards forming new mitochondrial genomes.

Electronic supplementary material

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

Deletion Genotypes Reduce Occlusion Body Potency but Increase Occlusion Body Production in a Colombian Spodoptera frugiperda Nucleopolyhedrovirus Population

A Colombian field isolate (SfCOL-wt) of Spodoptera frugiperda multiple nucleopolyhedrovirus (SfMNPV) is a mixture of different genotypes. To evaluate the insecticidal properties of the different genotypic variants, 83 plaque purified virus were characterized. Ten distinct genotypes were identified (named A through J). SfCOL-A was the most prevalent (71±2%; mean ± SE) showing a PstI restriction profile indistinguishable to that of SfCOL-wt. The remaining nine genotypes presented genomic deletions of 3.8 - 21.8 Kb located mainly between nucleotides 11,436 and 33,883 in the reference genome SfMNPV-B, affecting the region between open reading frames (ORFs) sf20 and sf33. The insecticidal activity of each genotype from SfCOL-wt and several mixtures of genotypes was compared to that of SfCOL-wt. The potency of SfCOL-A occlusion bodies (OBs) was 4.4-fold higher than SfCOL-wt OBs, whereas the speed of kill of SfCOL-A was similar to that of SfCOL-wt. Deletion genotype OBs were similarly or less potent than SfCOL-wt but six deletion genotypes were faster killing than SfCOL-wt. The potency of genotype mixtures co-occluded within OBs were consistently reduced in two-genotype mixtures involving equal proportions of SfCOL-A and one of three deletion genotypes (SfCOL-C, -D or -F). Speed of kill and OB production were improved only when the certain genotype mixtures were co-occluded, although OB production was higher in the SfCOL-wt isolate than in any of the component genotypes, or mixtures thereof. Deleted genotypes reduced OB potency but increased OB production of the SfCOL-wt population, which is structured to maximize the production of OBs in each infected host.     

Genomic Sequence of Heliothis virescens Ascovirus 3g Isolated from Spodoptera exigua

Heliothis virescens ascovirus 3a (HvAV-3a), a member of the family Ascoviridae, has the highest diversity among ascovirus species that have been reported in Australia, Indonesia, China, and the United States. To understand the diversity and origin of this important ascovirus, the complete genome of the HvAV Indonesia strain (HvAV-3g), isolated from Spodoptera exigua, was determined to be 199,721 bp, with a G+C content of 45.9%. Therefore, HvAV-3g has the largest genome among the reported ascovirus genomes to date. There are 194 predicted open reading frames (ORFs) encoding proteins of 50 or more amino acid residues. In comparison to HvAV-3e reported from Australia, HvAV-3g has all the ORFs in HvAV-3e with 6 additional ORFs unique to HvAV-3g, including 1 peptidase C26 gene with the highest identity to Drosophila spp. and 2 gas vesicle protein U (GvpU) genes with identities to Bacillus megaterium. The five unique homologous regions (hrs) and 25 baculovirus repeat ORFs (bro) of HvAV-3g are highly variable.     

Analysis of three leafminers' complete mitochondrial genomes

Liriomyza trifolii (Burgess), Liriomyza huidobrensis (Blanchard), and Liriomyza bryoniae (Kaltenbach), are three closely related and economically important leafminer pests in the world. This study examined the complete mitochondrial genomes of L. trifolii, L. huidobrensis and L. bryoniae, which were 16141 bp, 16236 bp and 16183 bp in length, respectively. All of them displayed 37 typical animal mitochondrial genes and an A + T-rich region. The genomes were highly compact with only 60–68 bp of non-coding intergenic spacer. However, considerable differences in the A + T-rich region were detected among the three species. Results of this study also showed the two ribosomal RNA genes of the three species had very limited variable sites and thus should not provide much information in the study of population genetics of these species. Data generated from three leafminers' complete mitochondrial genomes should provide valuable information in studying phylogeny of Diptera, and developing genetic markers for species identification in leafminers.     

Analysis of a naturally-occurring deletion mutant of Spodoptera frugiperda multiple nucleopolyhedrovirus reveals sf58 as a new per os infectivity factor of lepidopteran-infecting baculoviruses

The Nicaraguan population of Spodoptera frugiperda multiple nucleopolyhedrovirus, SfMNPV-NIC, is structured as a mixture of nine genotypes (A-I). Occlusion bodies (OBs) of SfMNPV-C, -D and -G pure genotypes are incapable of oral transmission; a phenotype which in SfMNPV-C and -D is due to the absence of pif1 and pif2 genes. The complete sequence of the SfMNPV-G genome was determined to identify possible factors involved in this phenotype. Deletions of 4860 bp (22,366-27,225) and 60 bp (119,759-119,818) were observed in SfMNPV-G genome compared with that of the predominant complete genotype SfMNPV-B (132,954 bp). However no genes homologous to previously described per os infectivity factors were located within the deleted sequences. Significant differences were detected in the nucleotide sequence in sf58 gene (unknown function) that produced changes in the amino acid sequence and the predicted secondary structure of the corresponding protein. This gene is conserved only in lepidopteran baculoviruses (alpha- and betabaculoviruses). To determine the role of sf58 in peroral infectivity a deletion mutant was constructed using bacmid technology. OBs of the deletion mutant (Sf58null) were not orally infectious for S. frugiperda larvae, whereas Sf58null rescue virus OBs recovered oral infectivity. Sf58null DNA and occlusion derived virions (ODVs) were as infective as SfMNPV bacmid DNA and ODVs in intrahemocelically infected larvae or cell culture, indicating that defects in ODV or OB morphogenesis were not involved in the loss of peroral infectivity. Addition of optical brightener or the presence of the orally infectious SfMNPV-B OBs in mixtures with SfMNPV-G OBs did not recover Sf58null OB infectivity. According to these results sf58 is a new per os infectivity factor present only in lepidopteran baculoviruses.     

Complete mitochondrial genome sequences of three bats species and whole genome mitochondrial analyses reveal patterns of codon bias and lend support to a basal split in Chiroptera

Order Chiroptera is a unique group of mammals whose members have attained self-powered flight as their main mode of locomotion. Much speculation persists regarding bat evolution; however, lack of sufficient molecular data hampers evolutionary and conservation studies. Of ~ 1200 species, complete mitochondrial genome sequences are available for only eleven. Additional sequences should be generated if we are to resolve many questions concerning these fascinating mammals. Herein, we describe the complete mitochondrial genomes of three bats: Corynorhinus rafinesquii, Lasiurus borealis and Artibeus lituratus. We also compare the currently available mitochondrial genomes and analyze codon usage in Chiroptera. C. rafinesquii, L. borealis and A. lituratus mitochondrial genomes are 16438 bp, 17048 bp and 16709 bp, respectively. Genome organization and gene arrangements are similar to other bats. Phylogenetic analyses using complete mitochondrial genome sequences support previously established phylogenetic relationships and suggest utility in future studies focusing on the evolutionary aspects of these species. Comprehensive analyses of available bat mitochondrial genomes reveal distinct nucleotide patterns and synonymous codon preferences corresponding to different chiropteran families. These patterns suggest that mutational and selection forces are acting to different extents within Chiroptera and shape their mitochondrial genomes.     

Sequence analysis of the genome of the Neodiprion sertifer nucleopolyhedrovirus

The genome of the Neodiprion sertifer nucleopolyhedrovirus (NeseNPV), which infects the European pine sawfly, N. sertifer (Hymenoptera: Diprionidae), was sequenced and analyzed. The genome was 86,462 bp in size. The C+G content of 34% was lower than that of the majority of baculoviruses. A total of 90 methionine-initiated open reading frames (ORFs) with more than 50 amino acids and minimal overlapping were found. From those, 43 ORFs were homologous to other baculovirus ORFs, and 29 of these were from the 30 conserved core genes among all baculoviruses. A NeseNPV homolog to the ld130 gene, which is present in all other baculovirus genomes sequenced to date, could not be identified. Six NeseNPV ORFs were similar to non-baculovirus-related genes, one of which was a trypsin-like gene. Only one iap gene, containing a single BIR motif and a RING finger, was found in NeseNPV. Two NeseNPV ORFs (nese18 and nese19) were duplicates transcribed in opposite orientations from each other. NeseNPV did not have an AcMNPV ORF 2 homolog characterized as the baculovirus repeat ORF (bro). Six homologous regions (hrs) were located within the NeseNPV genome, each containing small palindromes embedded within direct repeats. A phylogenetic analysis was done to root the tree based upon the sequences of DNA polymerase genes of NeseNPV, 23 other baculoviruses, and other phyla. Baculovirus phylogeny was then constructed with 29 conserved genes from 24 baculovirus genomes. Culex nigripalpus nucleopolyhedrovirus (CuniNPV) was the most distantly related baculovirus, branching to the hymenopteran NeseNPV and the lepidopteran nucleopolyhedroviruses and granuloviruses.     

Otx2 expression in anterior neuroectoderm and forebrain/midbrain is directed by more than six enhancers

Otx2 plays essential roles in each site at each step of head development. We previously identified the AN1 enhancer at 91 kb 5' upstream for the Otx2 expressions in anterior neuroectoderm (AN) at neural plate stage before E8.5, and the FM1 enhancer at 75 kb 5' upstream and the FM2 enhancer at 122 kb 3' downstream for the expression in forebrain/midbrain (FM) at brain vesicle stage after E8.5. The present study identified a second AN enhancer (AN2) at 88 kb 5' upstream; the AN2 enhancer also recapitulates the endogenous Otx2 expression in choroid plexus, cortical hem and choroidal roof. However, the enhancer mutants indicated the presence of another AN enhancer. The study also identified a third FM enhancer (FM3) at 153 kb 5' upstream. Thus, the Otx2 expressions in anterior neuroectoderm and forebrain/midbrain are regulated by more than six enhancers located far from the coding region. The enhancers identified are differentially conserved among vertebrates; none of the AN enhancers has activities in caudal forebrain and midbrain at brain vesicle stage after E8.5, nor do any of the FM enhancers in anterior neuroectoderm at neural plate stage before E8.5.     

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.     

Comparative Genome Sequence Analysis of Choristoneura occidentalis Freeman and C. rosaceana Harris (Lepidoptera: Tortricidae) Alphabaculoviruses

The complete genome sequences of Choristoneura occidentalis and C. rosaceana nucleopolyhedroviruses (ChocNPV and ChroNPV, respectively) (Baculoviridae: Alphabaculovirus) were determined and compared with each other and with those of other baculoviruses, including the genome of the closely related C. fumiferana NPV (CfMNPV). The ChocNPV genome was 128,446 bp in length (1147 bp smaller than that of CfMNPV), had a G+C content of 50.1%, and contained 148 open reading frames (ORFs). In comparison, the ChroNPV genome was 129,052 bp in length, had a G+C content of 48.6% and contained 149 ORFs. ChocNPV and ChroNPV shared 144 ORFs in common, and had a 77% sequence identity with each other and 96.5% and 77.8% sequence identity, respectively, with CfMNPV. Five homologous regions (hrs), with sequence similarities to those of CfMNPV, were identified in ChocNPV, whereas the ChroNPV genome contained three hrs featuring up to 14 repeats. Both genomes encoded three inhibitors of apoptosis (IAP-1, IAP-2, and IAP-3), as reported for CfMNPV, and the ChocNPV IAP-3 gene represented the most divergent functional region of this genome relative to CfMNPV. Two ORFs were unique to ChocNPV, and four were unique to ChroNPV. ChroNPV ORF chronpv38 is a eukaryotic initiation factor 5 (eIF-5) homolog that has also been identified in the C. occidentalis granulovirus (ChocGV) and is believed to be the product of horizontal gene transfer from the host. Based on levels of sequence identity and phylogenetic analysis, both ChocNPV and ChroNPV fall within group I alphabaculoviruses, where ChocNPV appears to be more closely related to CfMNPV than does ChroNPV. Our analyses suggest that it may be appropriate to consider ChocNPV and CfMNPV as variants of the same virus species.     

Isolation, Sequence Analysis, and Comparison of Two Plasmids (28 and 29 Kilobases) from the Biomining Bacterium Leptospirillum ferrooxidans ATCC 49879

Two plasmids, of 28,878 bp and 28,012 bp, were isolated from Leptospirillum ferrooxidans ATCC 49879. Altogether, a total of 67 open reading frames (ORFs) were identified on both plasmids, of which 32 had predicted products with high homology to proteins of known function, while 11 ORFs had predicted products with homology to previously identified proteins of unknown function. Twenty-four ORFs had products with no homologues in the GenBank/NCBI database. An analysis of the ORFs and other features of the two plasmids, the first to be isolated from a bacterium of the genus Leptospirillum, is presented.     

Fast assembly of the mitochondrial genome of a plant parasitic nematode (Meloidogyne graminicola) using next generation sequencing

Little is known about the variations of nematode mitogenomes (mtDNA). Sequencing a complete mtDNA using a PCR approach remains a challenge due to frequent genome reorganizations and low sequence similarities between divergent nematode lineages. Here, a genome skimming approach based on HiSeq sequencing (shotgun) was used to assemble de novo the first complete mtDNA sequence of a root-knot nematode (Meloidogyne graminicola). An AT-rich genome (84.3%) of 20,030 bp was obtained with a mean sequencing depth superior to 300. Thirty-six genes were identified with a semi-automated approach. A comparison with a gene map of the M. javanica mitochondrial genome indicates that the gene order is conserved within this nematode lineage. However, deep genome rearrangements were observed when comparing with other species of the superfamily Hoplolaimoidea. Repeat elements of 111 bp and 94 bp were found in a long non-coding region of 7.5 kb, as similarly reported in M. javanica and M. hapla. This study points out the power of next generation sequencing to produce complete mitochondrial genomes, even without a reference sequence, and possibly opening new avenues for species/race identification, phylogenetics and population genetics of nematodes.     

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.