En route to a genome-based classification of Archaea and Bacteria?

Given the considerable promise whole-genome sequencing offers for phylogeny and classification, it is surprising that microbial systematics and genomics have not yet been reconciled. This might be due to the intrinsic difficulties in inferring reasonable phylogenies from genomic sequences, particularly in the light of the significant amount of lateral gene transfer in prokaryotic genomes. However, recent studies indicate that the species tree and the hierarchical classification based on it are still meaningful concepts, and that state-of-the-art phylogenetic inference methods are able to provide reliable estimates of the species tree to the benefit of taxonomy. Conversely, we suspect that the current lack of completely sequenced genomes for many of the major lineages of prokaryotes and for most type strains is a major obstacle in progress towards a genome-based classification of microorganisms. We conclude that phylogeny-driven microbial genome sequencing projects such as the Genomic Encyclopaedia of Archaea and Bacteria (GEBA) project are likely to rectify this situation.     

Standard operating procedure for calculating genome-to-genome distances based on high-scoring segment pairs

DNA-DNA hybridization (DDH) is a widely applied wet-lab technique to obtain an estimate of the overall similarity between the genomes of two organisms. To base the species concept for prokaryotes ultimately on DDH was chosen by microbiologists as a pragmatic approach for deciding about the recognition of novel species, but also allowed a relatively high degree of standardization compared to other areas of taxonomy. However, DDH is tedious and error-prone and first and foremost cannot be used to incrementally establish a comparative database. Recent studies have shown that in-silico methods for the comparison of genome sequences can be used to replace DDH. Considering the ongoing rapid technological progress of sequencing methods, genome-based prokaryote taxonomy is coming into reach. However, calculating distances between genomes is dependent on multiple choices for software and program settings. We here provide an overview over the modifications that can be applied to distance methods based in high-scoring segment pairs (HSPs) or maximally unique matches (MUMs) and that need to be documented. General recommendations on determining HSPs using BLAST or other algorithms are also provided. As a reference implementation, we introduce the GGDC web server (http://ggdc.gbdp.org).     

Taxonomic evaluation of the Streptomyces hygroscopicus clade using multilocus sequence analysis and DNA-DNA hybridization, validating the MLSA scheme for systematics of the whole genus

Streptomyces hygroscopicus and related species are the most well known candidate producers of antibiotics and many other industrially and agronomically important secondary metabolites in the genus Streptomyces. Multilocus sequence analysis (MLSA) has shown to be a powerful and pragmatic molecular method for unraveling streptomycete diversities. In this investigation, a multilocus phylogeny of 58 representatives of the S. hygroscopicus 16S rRNA gene clade including S. violaceusniger and related species was examined. The result demonstrated that the MLSA data were helpful in defining members of the S. hygroscopicus clade, providing further evidence that the MLSA scheme of five housekeeping genes (atpD, gyrB, recA, rpoB and trpB) is a valuable alternative for creating and maintaining operational protocols for the Streptomyces species assignment. DNA-DNA hybridization (DDH) between strains with representative MLSA evolutionary distances, combined with previous data from S. griseus and S. albidoflavus clades, revealed a high correlation between MLSA and DDH, and sustains that the five-gene nucleotide sequence distance of 0.007 could be considered as the species cut-off for the whole genus. This significant correlation thus makes the MLSA scheme applicable to construction of a theory-based taxonomy for both ecology and bioprospecting of streptomycetes. Based on the MLSA and DDH data, as well as phenotypic characteristics, 10 species and three subspecies of the S. hygroscopicus clade are considered to be later heterotypic synonyms of eight genomic species, and Streptomyces glebosus sp. nov., comb. nov. (type strain CGMCC 4.1873(T)=LMG 19950(T)=DSM 40823(T)) and Streptomyces ossamyceticus sp. nov., comb. nov. (type strain CGMCC 4.1866(T)=LMG 19951(T)=DSM 40824(T)) are also proposed.     

Mitochondrial phylogenomics and genetic relationships of closely related pine moth (Lasiocampidae: Dendrolimus) species in China,using whole mitochondrial genomes

Background

Pine moths (Lepidoptera; Bombycoidea; Lasiocampidae: Dendrolimus spp.) are among the most serious insect pests of forests, especially in southern China. Although COI barcodes (a standardized portion of the mitochondrial cytochrome c oxidase subunit I gene) can distinguish some members of this genus, the evolutionary relationships of the three morphospecies Dendrolimus punctatus, D. tabulaeformis and D. spectabilis have remained largely unresolved. We sequenced whole mitochondrial genomes of eight specimens, including D. punctatuswenshanensis. This is an unambiguous subspecies of D. punctatus, and was used as a reference for inferring the relationships of the other two morphospecies of the D. punctatus complex. We constructed phylogenetic trees from this data, including twelve published mitochondrial genomes of other Bombycoidea species, and examined the relationships of the Dendrolimus taxa using these trees and the genomic features of the mitochondrial genome.

Results

The eight fully sequenced mitochondrial genomes from the three morphospecies displayed similar genome structures as other Bombycoidea species in terms of gene content, base composition, level of overall AT-bias and codon usage. However, the Dendrolimus genomes possess a unique feature in the large ribosomal 16S RNA subunits (rrnL), which are more than 60 bp longer than other members of the superfamily and have a higher AC proportion. The eight mitochondrial genomes of Dendrolimus were highly conservative in many aspects, for example with identical stop codons and overlapping regions. But there were many differences in start codons, intergenic spacers, and numbers of mismatched base pairs of tRNA (transfer RNA genes).Our results, based on phylogenetic trees, genetic distances, species delimitation and genomic features (such as intergenic spacers) of the mitochondrial genome, indicated that D. tabulaeformis is as close to D. punctatus as is D. punctatus wenshanensis, whereas D. spectabilis evolved independently from D. tabulaeformis and D. punctatus. Whole mitochondrial DNA phylogenies showed that D. spectabilis formed a well-supported monophyletic clade, with a clear species boundary separating it from the other congeners examined here. However, D. tabulaeformis often clustered with D. punctatus and with the subspecies D. punctatus wenshanensis. Genetic distance analyses showed that the distance between D. tabulaeformis and D. punctatus is generally less than the intraspecific distance of D. punctatus and its subspecies D. punctatus wenshanensis. In the species delimitation analysis of Poisson Tree Processes (PTP), D. tabulaeformis, D. punctatus and D. punctatus wenshanensis clustered into a putative species separated from D. spectabilis. In comparison with D. spectabilis, D. tabulaeformis and D. punctatus also exhibit a similar structure in intergenic spacer characterization. These different types of evidence suggest that D. tabulaeformis is very close to D. punctatus and its subspecies D. punctatus wenshanensis, and is likely to be another subspecies of D. punctatus.

Conclusions

Whole mitochondrial genomes possess relatively rich genetic information compared with the traditional use of single or multiple genes for phylogenetic purposes. They can be used to better infer phylogenetic relationships and degrees of relatedness of taxonomic groups, at least from the aspect of maternal lineage: caution should be taken due to the maternal-only inheritance of this genome. Our results indicate that D. spectabilis is an independent lineage, while D. tabulaeformis shows an extremely close relationship to D. punctatus.

Electronic supplementary material

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

Genome-wide comparative analysis of simple sequence coding repeats among 25 insect species

We present a detailed genome-scale comparative analysis of simple sequence repeats within protein coding regions among 25 insect genomes. The repetitive sequences in the coding regions primarily represented single codon repeats and codon pair repeats. The CAG triplet is highly repetitive in the coding regions of insect genomes. It is frequently paired with the synonymous codon CAA to code for polyglutamine repeats. The codon pairs that are least repetitive code for polyalanine repeats. The frequency of hexanucleotide and dinucleotide motifs of codon pair repeats is significantly (p<0.001) different in the Drosophila species compared to the non-Drosophila species. However, the frequency of synonymous and non-synonymous codon pair repeats varies in a correlated manner (r(2)=0.79) among all the species. Results further show that perfect and imperfect repeats have significant association with the trinucleotide and hexanucleotide coding repeats in most of these insects. However, only select species show significant association between the numbers of perfect/imperfect hexamers and repeat coding for single amino acid/amino acid pair runs. Our data further suggests that genes containing simple sequence coding repeats may be under negative selection as they tend to be poorly conserved across species. The sequences of coding repeats of orthologous genes vary according to the known phylogeny among the species. In conclusion, the study shows that simple sequence coding repeats are important features of genome diversity among insects.     

Trypanosomatid comparative genomics: Contributions to the study of parasite biology and different parasitic diseases

In 2005, draft sequences of the genomes of Trypanosoma brucei, Trypanosoma cruzi and Leishmania major, also known as the Tri-Tryp genomes, were published. These protozoan parasites are the causative agents of three distinct insect-borne diseases, namely sleeping sickness, Chagas disease and leishmaniasis, all with a worldwide distribution. Despite the large estimated evolutionary distance among them, a conserved core of ~6,200 trypanosomatid genes was found among the Tri-Tryp genomes. Extensive analysis of these genomic sequences has greatly increased our understanding of the biology of these parasites and their host-parasite interactions. In this article, we review the recent advances in the comparative genomics of these three species. This analysis also includes data on additional sequences derived from other trypanosmatid species, as well as recent data on gene expression and functional genomics. In addition to facilitating the identification of key parasite molecules that may provide a better understanding of these complex diseases, genome studies offer a rich source of new information that can be used to define potential new drug targets and vaccine candidates for controlling these parasitic infections.     

Non-referenced genome assembly from epigenomic short-read data

Current computational methods used to analyze changes in DNA methylation and chromatin modification rely on sequenced genomes. Here we describe a pipeline for the detection of these changes from short-read sequence data that does not require a reference genome. Open source software packages were used for sequence assembly, alignment, and measurement of differential enrichment. The method was evaluated by comparing results with reference-based results showing a strong correlation between chromatin modification and gene expression. We then used our de novo sequence assembly to build the DNA methylation profile for the non-referenced Psammomys obesus genome. The pipeline described uses open source software for fast annotation and visualization of unreferenced genomic regions from short-read data.     

Molecular Approaches to Taenia asiatica

Taenia solium, T. saginata, and T. asiatica are taeniid tapeworms that cause taeniasis in humans and cysticercosis in intermediate host animals. Taeniases remain an important public health concerns in the world. Molecular diagnostic methods using PCR assays have been developed for rapid and accurate detection of human infecting taeniid tapeworms, including the use of sequence-specific DNA probes, PCR-RFLP, and multiplex PCR. More recently, DNA diagnosis using PCR based on histopathological specimens such as 10% formalin-fixed paraffin-embedded and stained sections mounted on slides has been applied to cestode infections. The mitochondrial gene sequence is believed to be a very useful molecular marker for not only studying evolutionary relationships among distantly related taxa, but also for investigating the phylo-biogeography of closely related species. The complete sequence of the human Taenia tapeworms mitochondrial genomes were determined, and its organization and structure were compared to other human-tropic Taenia tapeworms for which complete mitochondrial sequence data were available. The multiplex PCR assay with the Ta4978F, Ts5058F, Tso7421F, and Rev7915 primers will be useful for differential diagnosis, molecular characterization, and epidemiological surveys of human Taenia tapeworms.     

Identification and experimental validation of splicing regulatory elements in Drosophila melanogaster reveals functionally conserved splicing enhancers in metazoans

RNA sequence elements involved in the regulation of pre-mRNA splicing have previously been identified in vertebrate genomes by computational methods. Here, we apply such approaches to predict splicing regulatory elements in Drosophila melanogaster and compare them with elements previously found in the human, mouse, and pufferfish genomes. We identified 99 putative exonic splicing enhancers (ESEs) and 231 putative intronic splicing enhancers (ISEs) enriched near weak 5' and 3' splice sites of constitutively spliced introns, distinguishing between those found near short and long introns. We found that a significant proportion (58%) of fly enhancer sequences were previously reported in at least one of the vertebrates. Furthermore, 20% of putative fly ESEs were previously identified as ESEs in human, mouse, and pufferfish; while only two fly ISEs, CTCTCT and TTATAA, were identified as ISEs in all three vertebrate species. Several putative enhancer sequences are similar to characterized binding-site motifs for Drosophila and mammalian splicing regulators. To provide additional evidence for the function of putative ISEs, we separately identified 298 intronic hexamers significantly enriched within sequences phylogenetically conserved among 15 insect species. We found that 73 putative ISEs were among those enriched in conserved regions of the D. melanogaster genome. The functions of nine enhancer sequences were verified in a heterologous splicing reporter, demonstrating that these sequences are sufficient to enhance splicing in vivo. Taken together, these data identify a set of predicted positive-acting splicing regulatory motifs in the Drosophila genome and reveal regulatory sequences that are present in distant metazoan genomes.     

高榕榕果内Eupristina属两种榕小蜂的遗传进化关系

对生长在福州地区的高榕进行长期追踪观察,发现高榕榕果内仅生活着Eupristina altissima和Eupristina sp.榕小蜂,前者为高榕的传粉小蜂,后者无传粉行为,两者雌蜂之间在体色、触角、花粉袋和花粉刷等部位存在细微的差异,而两者雄蜂之间无形态差异。通过克隆福建地区5个样地的高榕榕果内收集到的E. altissima和Eupristina sp.榕小蜂,以及细叶榕的传粉小蜂Eupristina verticillata(外群)的Cytb及COI基因,并进行碱基组成及遗传距离分析,用邻接法构建系统发育树,分析两榕小蜂群体之间的遗传进化关系,结果显示:(1)榕小蜂COI及Cytb序列碱基组成中A+T的含量(Cytb序列中A+T=75.3%,COI序列中A+T=75.5%)显著高于G+C,符合膜翅目昆虫线粒体基因碱基组成特征。(2)对两群体小蜂进行遗传距离分析显示,Cytb序列中E. altissima和Eupristina sp. 群体内各样本之间的平均遗传距离分别为0.0092和0.0030,而E. altissima与Eupristina sp. 群体间的平均距离为0.1588;COI序列中E. altissima 与Eupristina sp. 群体内各样本之间的平均遗传距离分别为0.0065和0.0205,而二者群体间的平均遗传距离为0.1043,表明两者群体间的遗传距离明显大于各自群体内各样本间的遗传距离。统计GenBank中下载的6个属34种榕小蜂Cytb序列的种间遗传距离为0.0811-0.1723,6个属28种榕小蜂COI序列的种间遗传距离为0.0939-0.1986。由此认为E. altissima和Eupristina sp.之间的遗传距离差异已经达到了种间水平,即E. altissima与Eupristina sp.为两个不同的种。(3)在形态上,两种小蜂的雌蜂之间有微小差异,而二者雄蜂之间无差异,但Cytb与COI序列分析结果一致表明:E. altissima与Eupristina sp.雄蜂之间,以及二者雌蜂之间的遗传距离均差异显著,表明形态变异滞后于基因变异。雌蜂在表型上进化快于雄蜂,可能是由于雌蜂羽化后从榕果出飞,受到外界环境因素的影响较大,且两种雌蜂在传粉功能上存在差异,故二者之间的形态差异较大,而雄蜂寿命短,又终生生活在黑暗封闭、环境变化相对恒定的榕果内,两种雄蜂在行为上不存在差异,故二者表型变异较为缓慢。E. altissima和Eupristina sp.小蜂对宿主的专一性不强,在榕-蜂协同进化过程中,可能发生过宿主转移事件。     

Comparative genome analysis of Solanum lycopersicum and Solanum tuberosum

Solanum lycopersicum and Solanum tuberosum are agriculturally important crop species as they are rich sources of starch, protein, antioxidants, lycopene, beta-carotene, vitamin C, and fiber. The genomes of S. lycopersicum and S. tuberosum are currently available. However the linear strings of nucleotides that together comprise a genome sequence are of limited significance by themselves. Computational and bioinformatics approaches can be used to exploit the genomes for fundamental research for improving their varieties. The comparative genome analysis, Pfam analysis of predicted reviewed paralogous proteins was performed. It was found that S. lycopersicum proteins belong to more families, domains and clans in comparison with S. tuberosum. It was also found that mostly intergenic regions are conserved in two genomes followed by exons, intron and UTR. This can be exploited to predict regions between genomes that are similar to each other and to study the evolutionary relationship between two genomes, leading towards the development of disease resistance, stress tolerance and improved varieties of tomato.     

An assembly and alignment-free method of phylogeny reconstruction from next-generation sequencing data

Background

Next-generation sequencing technologies are rapidly generating whole-genome datasets for an increasing number of organisms. However, phylogenetic reconstruction of genomic data remains difficult because de novo assembly for non-model genomes and multi-genome alignment are challenging.

Results

To greatly simplify the analysis, we present an Assembly and Alignment-Free (AAF) method (https://sourceforge.net/projects/aaf-phylogeny) that constructs phylogenies directly from unassembled genome sequence data, bypassing both genome assembly and alignment. Using mathematical calculations, models of sequence evolution, and simulated sequencing of published genomes, we address both evolutionary and sampling issues caused by direct reconstruction, including homoplasy, sequencing errors, and incomplete sequencing coverage. From these results, we calculate the statistical properties of the pairwise distances between genomes, allowing us to optimize parameter selection and perform bootstrapping. As a test case with real data, we successfully reconstructed the phylogeny of 12 mammals using raw sequencing reads. We also applied AAF to 21 tropical tree genome datasets with low coverage to demonstrate its effectiveness on non-model organisms.

Conclusion

Our AAF method opens up phylogenomics for species without an appropriate reference genome or high sequence coverage, and rapidly creates a phylogenetic framework for further analysis of genome structure and diversity among non-model organisms.

Electronic supplementary material

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

GC content dependency of open reading frame prediction via stop codon frequencies

A frequently used approach for detecting potential coding regions is to search for stop codons. In the standard genetic code 3 out of 64 trinucleotides are stop codons. Hence, in random or non-coding DNA one can expect every 21st trinucleotide to have the same sequence as a stop codon. In contrast, the open reading frames (ORFs) of most protein-coding genes are considerably longer. Thus, the stop codon frequency in coding sequences deviates from the background frequency of the corresponding trinucleotides. This has been utilized for gene prediction, in particular, in detecting protein-coding ORFs. Traditional methods based on stop codon frequency are based on the assumption that the GC content is about 50%. However, many genomes show significant deviations from that value. With the presented method we can describe the effects of GC content on the selection of appropriate length thresholds of potentially coding ORFs. Conversely, for a given length threshold, we can calculate the probability of observing it in a random sequence. Thus, we can derive the maximum GC content for which ORF length is practicable as a feature for gene prediction methods and the resulting false positive rates. A rough estimate for an upper limit is a GC content of 80%. This estimate can be made more precise by including further parameters and by taking into account start codons as well. We demonstrate the feasibility of this method by applying it to the genomes of the bacteria Rickettsia prowazekii, Escherichia coli and Caulobacter crescentus, exemplifying the effect of GC content variations according to our predictions. We have adapted the method for predicting coding ORFs by stop codon frequency to the case of GC contents different from 50%. Usually, several methods for gene finding need to be combined. Thus, our results concern a specific part within a package of methods. Interestingly, for genomes with low GC content such as that of R. prowazekii, the presented method provides remarkably good results even when applied alone.     

Diversity array technology markers: genetic diversity analyses and linkage map construction in rapeseed (Brassica napus L.)

We developed Diversity Array Technology (DArT) markers for application in genetic studies of Brassica napus and other Brassica species with A or C genomes. Genomic representation from 107 diverse genotypes of B. napus L. var. oleifera (rapeseed, AACC genomes) and B. rapa (AA genome) was used to develop a DArT array comprising 11 520 clones generated using PstI/BanII and PstI/BstN1 complexity reduction methods. In total, 1547 polymorphic DArT markers of high technical quality were identified and used to assess molecular diversity among 89 accessions of B. napus, B. rapa, B. juncea, and B. carinata collected from different parts of the world. Hierarchical cluster and principal component analyses based on genetic distance matrices identified distinct populations clustering mainly according to their origin/pedigrees. DArT markers were also mapped in a new doubled haploid population comprising 131 lines from a cross between spring rapeseed lines 'Lynx-037DH' and 'Monty-028DH'. Linkage groups were assigned on the basis of previously mapped simple sequence repeat (SSRs), intron polymorphism (IP), and gene-based markers. The map consisted of 437 DArT, 135 SSR, 6 IP, and 6 gene-based markers and spanned 2288 cM. Our results demonstrate that DArT markers are suitable for genetic diversity analysis and linkage map construction in rapeseed.     

Analysis of genetic population structure in Acacia caven (Leguminosae,Mimosoideae), comparing one exploratory and two Bayesian-model-based methods

Bayesian clustering as implemented in STRUCTURE or GENELAND software is widely used to form genetic groups of populations or individuals. On the other hand, in order to satisfy the need for less computer-intensive approaches, multivariate analyses are specifically devoted to extracting information from large datasets. In this paper, we report the use of a dataset of AFLP markers belonging to 15 sampling sites of Acacia caven for studying the genetic structure and comparing the consistency of three methods: STRUCTURE, GENELAND and DAPC. Of these methods, DAPC was the fastest one and showed accuracy in inferring the K number of populations (K = 12 using the find.clusters option and K = 15 with a priori information of populations). GENELAND in turn, provides information on the area of membership probabilities for individuals or populations in the space, when coordinates are specified (K = 12). STRUCTURE also inferred the number of K populations and the membership probabilities of individuals based on ancestry, presenting the result K = 11 without prior information of populations and K = 15 using the LOCPRIOR option. Finally, in this work all three methods showed high consistency in estimating the population structure, inferring similar numbers of populations and the membership probabilities of individuals to each group, with a high correlation between each other.     

Non-contiguous finished genome sequence of the opportunistic oral pathogen Prevotella multisaccharivorax type strain (PPPA20)

Prevotella multisaccharivorax Sakamoto et al. 2005 is a species of the large genus Prevotella, which belongs to the family Prevotellaceae. The species is of medical interest because its members are able to cause diseases in the human oral cavity such as periodontitis, root caries and others. Although 77 Prevotella genomes have already been sequenced or are targeted for sequencing, this is only the second completed genome sequence of a type strain of a species within the genus Prevotella to be published. The 3,388,644 bp long genome is assembled in three non-contiguous contigs, harbors 2,876 protein-coding and 75 RNA genes and is a part of the Genomic Encyclopedia of Bacteria and Archaea project.     

Resistance of Auto- and Allotetraploid Triticeae Species and Accessions to Meloidogyne chitwoodi based on Genome Composition

The Columbia root-knot nematode Meloidogyne chitwoodi parasitizes several plant species, including grasses that have been developed for semiarid environments, and substantially reduces the productivity of cereals and the longevity of perennial grasses growing under semiarid conditions throughout the intermountain region. Thirty-two auto- and allotetraploid (2n = 28) taxa in the perennial Triticeae were evaluated as possible sources of resistance to M. chitwoodi. Low levels of root galling were observed on roots of all accessions; root-gall indices ranged from 0 (no galls) to 1.95 in the grasses compared to 4.67 for the susceptible ''Ranger'' alfalfa check on a scale of 1 to 6. Even though the gall ratings were low, significant (P < 0.01) differences among accessions of the same species, among species, and among genera with different genomes were observed. Within the reproductive indices, which ranged from 0.01 to 1.20 in the grasses compared to 65.38 for the alfalfa check, there was no difference among genera with different genomes and accessions within the same species and genome; however, there was a significant (P < 0.05) difference among species with the same genomes. This variation can be traced to Thinopyrum nodosum (Jaaska-19), which was the only accession with a reproductive factor greater than 1.00. Based on the data, all auto- and allotetraploids are considered resistant to M. chitwoodi.     

The genus Allochromatium (Chromatiales Chromatiaceae) revisited: a study on its intragenic structure based on multilocus sequence analysis (MLSA) and DNA-DNA hybridization (DDH)

In this study, the taxonomic status of anoxygenic photosynthetic bacteria belonging to the genus Allochromatium is revisited. The inter- and intraspecies relationship of seven Allochromatium strains, including a set of well described type strains, were examined by DNA-DNA hybridization (DDH) and multilocus sequence analysis (MLSA) using segments of seven protein-coding genes. The re-sequencing of the 16S rRNA, the internal transcriber spacer (ITS), multi-gene analysis and DDH comparison indicated that both type strains Allochromatium vinosum DSM 180^T and Allochromatium minutissimum DSM 1376^T are closely related to each other forming an independent cluster together with the strains A. vinosum DSM 183 and DSM 1686. The internal comparison of members of this A. vinosum phylogroup showed values of DDH relatedness above 80% and concatenated sequence similarities (4744 bp) above 98%. In contrast, the MLSA scheme has identified A. vinosum strain BH-2 as a separate lineage. Strain BH-2 was first classified as a member of the species A. vinosum based on DDH comparison. However, this strain showed the lowest similarity values of the 16S rRNA gene and concatenated sequences, as well as amino acid identity (AAI) when compared to other Allochromatium strains, suggesting that strain BH-2 may represent a new species.     

Disentangling the effects of mating systems and mutation rates on cytoplamic diversity in gynodioecious Silene nutans and dioecious Silene otites

Many flowering plant species exhibit a variety of distinct sexual morphs, the two most common cases being the co-occurrence of females and males (dioecy) or the co-occurrence of hermaphrodites and females (gynodioecy). In this study, we compared DNA sequence variability of the three genomes (nuclear, mitochondrial and chloroplastic) of a gynodioecious species, Silene nutans, with that of a closely related dioecious species, Silene otites. In the light of theoretical models, we expect cytoplasmic diversity to differ between the two species due to the selective dynamics that acts on cytoplasmic genomes in gynodioecious species: under an epidemic scenario, the gynodioecious species is expected to exhibit lower cytoplasmic diversity than the dioecious species, while the opposite is expected in the case of balancing selection maintaining sterility cytoplasms in the gynodioecious species. We found no difference between the species for nuclear gene diversity, but, for the cytoplasmic loci, the gynodioecious S. nutans had more haplotypes, and higher nucleotide diversity, than the dioecious relative, S. otites, even though the latter has a relatively high rate of mitochondrial synonymous substitutions, and therefore presumably a higher mutation rate. Therefore, as the mitochondrial mutation rate cannot account for the higher cytoplasmic diversity found in S. nutans, our findings support the hypothesis that gynodioecy in S. nutans has been maintained by balancing selection rather than by epidemic-like dynamics.