Base Composition Bias Research Articles

Evolutionary Relationships of Omani Macrotermes subhyalinus, Macrotermitinae.

A study was conducted to investigate the evolutionary relationships of Macrotermes subhyalinus from Oman, in the southeastern part of the Arabian Peninsula. Sequences of the mitochondrial COI and the nuclear large-subunit ribosomal RNA (LSU rRNA, 28S) genes were used to investigate the populations of M. subhyalinus across their distribution in Oman to determine their relationships with other Macrotermes species. Our findings indicate that M. subhyalinus in Oman is a member of an East African clade, distinct from those in West Africa. Analyses of the COI showed that there is base composition bias among the taxa (non-stationarity) that has not been considered in earlier studies. We provide the first report of pseudogene copies of 28S in M. subhyalinus that are differentially amplified.

Enriching barcoding markers in environmental samples utilizing a phylogenetic probe design: Insights from mock communities

AbstractHybridization capture is an emerging method making use of short oligonucleotide baits to enrich DNA libraries for genomic fragments of specific organisms thus enabling detection of their presence in environmental samples. Although it offers a primer‐independent alternative to metabarcoding, little empirical work has been dedicated to characterizing the underlying biases and coupled implications for biological interpretation. Moreover, few published bioinformatic pipelines are available for designing polynucleotide capture baits from a reference sequence collection. We designed RNA‐baits specifically targeting two chloroplast barcoding genes matK and rbcL to reveal the plant taxonomic diversity present in a given environmental sample. Our approach leverages the sensitivity of hybridization capture and the capacity of high‐throughput DNA sequencing instruments. It builds on a new and universal method based on ancestral sequence reconstruction, ultimately limiting the number of bait‐probes required and reducing experimental costs, while accessing high levels of taxonomic diversity. Our bait‐set selectively targets four main plant orders (Fagales, Pinales, Asterales, and Poales), representing ~18% of all described vascular plants. This is achieved through the use of only 4084 baits, each 80 nucleotides in length (80‐mer), capturing ~1.0–1.6 k nucleotide sequences from each taxon. Tests on mock communities revealed important factors influencing capture efficiency and relative abundance estimates, including GC‐content, the overall target length per taxa, and the bait density and mean number of mismatches to the bait sequence. Our results show that hybridization capture, like metabarcoding, requires caution when interpreting results quantitatively within (paleo)‐ecological studies. Biases detected in this work have the potential to be mitigated with bait designs that avoid extreme base compositional biases and balancing bait targets across taxa. However, we strongly recommend the use of mock communities and read simulations to quantify the accuracy of taxonomic representation when using new bait designs.

The complete mitochondrial genome of a species of Cirrhipathes de Blainville, 1830 from Kauaʻi, Hawaiʻi (Hexacorallia: Antipatharia)

This study reports the first mitogenome from the antipatharian (black coral) genus Cirrhipathes (GenBank accession number ON653414). The 20,452 bp mitochondrial genome of Cirrhipathes cf. anguina LS-2022 consists of 13 protein-coding genes, two rRNA genes, and two tRNA genes (trnM and trnW). The mitogenome is typical of other antipatharian families, including an A + T biased (64.1%) base composition and cytochrome c oxidase subunit I (COX1) intron with embedded homing endonuclease gene (HEG). A phylogenetic tree based on complete mitogenome sequences of currently available antipatharians indicates Cirrhipathes cf. anguina LS-2022 is sister and closely related to Stichopathes sp. SCBUCN-8849. However, it seems unlikely that intergeneric taxa share 99.97% similarity across their complete mitogenomes, raising questions about the current taxonomy of this group. This study highlights the need for additional vouchered antipatharian species to be sequenced so phylogenetic relationships can be compared with accepted taxonomy.

GCparagon: evaluating and correcting GC biases in cell-free DNA at the fragment level.

Analyses of cell-free DNA (cfDNA) are increasingly being employed for various diagnostic and research applications. Many technologies aim to increase resolution, e.g. for detecting early-stage cancer or minimal residual disease. However, these efforts may be confounded by inherent base composition biases of cfDNA, specifically the over - and underrepresentation of guanine (G) and cytosine (C) sequences. Currently, there is no universally applicable tool to correct these effects on sequencing read-level data. Here, we present GCparagon, a two-stage algorithm for computing and correcting GC biases in cfDNA samples. In the initial step, length and GC base count parameters are determined. Here, our algorithm minimizes the inclusion of known problematic genomic regions, such as low-mappability regions, in its calculations. In the second step, GCparagon computes weights counterbalancing the distortion of cfDNA attributes (correction matrix). These fragment weights are added to a binary alignment map (BAM) file as alignment tags for individual reads. The GC correction matrix or the tagged BAM file can be used for downstream analyses. Parallel computing allows for a GC bias estimation below 1 min. We demonstrate that GCparagon vastly improves the analysis of regulatory regions, which frequently show specific GC composition patterns and will contribute to standardized cfDNA applications.

Base composition, codon usage and patterns of gene sequence evolution in butterflies.

Coding sequence evolution is influenced by both natural selection and neutral evolutionary forces. In many species, the effects of mutation bias, codon usage and GC-biased gene conversion (gBGC) on gene sequence evolution have not been detailed. Quantification of how these forces shape substitution patterns is therefore necessary to understand the strength and direction of natural selection. Here, we used comparative genomics to investigate the association between base composition and codon usage bias on gene sequence evolution in butterflies and moths (Lepidoptera), including an in-depth analysis of underlying patterns and processes in one species, Leptidea sinapis. The data revealed significant G/C to A/T substitution bias at third codon position with some variation in the strength among different butterfly lineages. However, the substitution bias was lower than expected from previously estimated mutation rate ratios, partly due to the influence of gBGC. We found that A/T-ending codons were overrepresented in most species, but there was a positive association between the magnitude of codon usage bias and GC-content in third codon positions. In addition, the tRNA-gene population in L. sinapis showed higher GC-content at third codon positions compared to coding sequences in general and less overrepresentation of A/T-ending codons. There was an inverse relationship between synonymous substitutions and codon usage bias indicating selection on synonymous sites. We conclude that the evolutionary rate in Lepidoptera is affected by a complex interaction between underlying G/C -> A/T mutation bias and partly counteracting fixation biases, predominantly conferred by overall purifying selection, gBGC and selection on codon usage.

Mitochondrial phylogenomics provides conclusive evidence that the family Ancyrocephalidae is deeply paraphyletic

BackgroundUnresolved taxonomic classification and paraphyly pervade the flatworm class Monogenea: the class itself may be paraphyletic and split into Polyopisthocotylea and Monopisthocotylea; there are some indications that the monopisthocotylean order Dactylogyridea may also be paraphyletic; single-gene markers and some morphological traits indicate that the family Ancyrocephalidae is paraphyletic and intertwined with the family Dactylogyridae.MethodsTo attempt to study the relationships of Ancyrocephalidae and Monopisthocotylea using a phylogenetic marker with high resolution, we sequenced mitochondrial genomes of two fish ectoparasites from the family Dactylogyridae: Dactylogyrus simplex and Dactylogyrus tuba. We conducted phylogenetic analyses using three datasets and three methods. Datasets were ITS1 (nuclear) and nucleotide and amino acid sequences of almost complete mitogenomes of almost all available Monopisthocotylea mitogenomes. Methods were maximum likelihood (IQ-TREE), Bayesian inference (MrBayes) and CAT-GTR (PhyloBayes).ResultsBoth mitogenomes exhibited the ancestral gene order for Neodermata, and both were compact, with few and small intergenic regions and many and large overlaps. Gene sequences were remarkably divergent for nominally congeneric species, with only trnI exhibiting an identity value > 80%. Both mitogenomes had exceptionally low A + T base content and AT skews. We found evidence of pervasive compositional heterogeneity in the dataset and indications that base composition biases cause phylogenetic artefacts. All six mitogenomic analyses produced unique topologies, but all nine analyses produced topologies that rendered Ancyrocephalidae deeply paraphyletic. Mitogenomic data consistently resolved the order Capsalidea as nested within the Dactylogyridea.ConclusionsThe analyses indicate that taxonomic revisions are needed for multiple Polyopisthocotylea lineages, from genera to orders. In combination with previous findings, these results offer conclusive evidence that Ancyrocephalidae is a paraphyletic taxon. The most parsimonious solution to resolve this is to create a catch-all Dactylogyridae sensu lato clade comprising the current Ancyrocephalidae, Ancylodiscoididae, Pseudodactylogyridae and Dactylogyridae families, but the revision needs to be confirmed by another marker with a sufficient resolution.Graphical

Detecting and Removing Sample Contamination in Phylogenomic Data: An Example and its Implications for Cicadidae Phylogeny (Insecta: Hemiptera).

Contamination of a genetic sample with DNA from one or more nontarget species is a continuing concern of molecular phylogenetic studies, both Sanger sequencing studies and next-generation sequencing studies. We developed an automated pipeline for identifying and excluding likely cross-contaminated loci based on the detection of bimodal distributions of patristic distances across gene trees. When contamination occurs between samples within a data set, a comparison between a contaminated sample and its contaminant taxon will yield bimodal distributions with one peak close to zero patristic distance. This new method does not rely on a priori knowledge of taxon relatedness nor does it determine the causes(s) of the contamination. Exclusion of putatively contaminated loci from a data set generated for the insect family Cicadidae showed that these sequences were affecting some topological patterns and branch supports, although the effects were sometimes subtle, with some contamination-influenced relationships exhibiting strong bootstrap support. Long tip branches and outlier values for one anchored phylogenomic pipeline statistic (AvgNHomologs) were correlated with the presence of contamination. While the anchored hybrid enrichment markers used here, which target hemipteroid taxa, proved effective in resolving deep and shallow level Cicadidae relationships in aggregate, individual markers contained inadequate phylogenetic signal, in part probably due to short length. The cleaned data set, consisting of 429 loci, from 90 genera representing 44 of 56 current Cicadidae tribes, supported three of the four sampled Cicadidae subfamilies in concatenated-matrix maximum likelihood (ML) and multispecies coalescent-based species tree analyses, with the fourth subfamily weakly supported in the ML trees. No well-supported patterns from previous family-level Sanger sequencing studies of Cicadidae phylogeny were contradicted. One taxon (Aragualna plenalinea) did not fall with its current subfamily in the genetic tree, and this genus and its tribe Aragualnini is reclassified to Tibicininae following morphological re-examination. Only subtle differences were observed in trees after the removal of loci for which divergent base frequencies were detected. Greater success may be achieved by increased taxon sampling and developing a probe set targeting a more recent common ancestor and longer loci. Searches for contamination are an essential step in phylogenomic analyses of all kinds and our pipeline is an effective solution. [Auchenorrhyncha; base-composition bias; Cicadidae; Cicadoidea; Hemiptera; phylogenetic conflict.].

Inverted base composition skews and discontinuous mitochondrial genome architecture evolution in the Enoplea (Nematoda)

BackgroundWithin the class Enoplea, the earliest-branching lineages in the phylum Nematoda, the relatively highly conserved ancestral mitochondrial architecture of Trichinellida is in stark contrast to the rapidly evolving architecture of Dorylaimida and Mermithida. To better understand the evolution of mitogenomic architecture in this lineage, we sequenced the mitogenome of a fish parasite Pseudocapillaria tomentosa (Trichinellida: Capillariidae) and compared it to all available enoplean mitogenomes.ResultsP. tomentosa exhibited highly reduced noncoding regions (the largest was 98 bp), and a unique base composition among the Enoplea. We attributed the latter to the inverted GC skew (0.08) in comparison to the ancestral skew in Trichinellidae (-0.43 to -0.37). Capillariidae, Trichuridae and Longidoridae (Dorylaimida) generally exhibited low negative or low positive skews (-0.1 to 0.1), whereas Mermithidae exhibited fully inverted low skews (0 to 0.05). This is indicative of inversions in the strand replication order or otherwise disrupted replication mechanism in the lineages with reduced/inverted skews. Among the Trichinellida, Trichinellidae and Trichuridae have almost perfectly conserved architecture, whereas Capillariidae exhibit multiple rearrangements of tRNA genes. In contrast, Mermithidae (Mermithida) and Longidoridae (Dorylaimida) exhibit almost no similarity to the ancestral architecture.ConclusionsLongidoridae exhibited more rearranged mitogenomic architecture than the hypervariable Mermithidae. Similar to the Chromadorea, the evolution of mitochondrial architecture in enoplean nematodes exhibits a strong discontinuity: lineages possessing a mostly conserved architecture over tens of millions of years are interspersed with lineages exhibiting architectural hypervariability. As Longidoridae also have some of the smallest metazoan mitochondrial genomes, they contradict the prediction that compact mitogenomes should be structurally stable. Lineages exhibiting inverted skews appear to represent the intermediate phase between the Trichinellidae (ancestral) and fully derived skews in Chromadorean mitogenomes (GC skews = 0.18 to 0.64). Multiple lines of evidence (CAT-GTR analysis in our study, a majority of previous mitogenomic results, and skew disruption scenarios) support the Dorylaimia split into two sister-clades: Dorylaimida + Mermithida and Trichinellida. However, skew inversions produce strong base composition biases, which can hamper phylogenetic and other evolutionary studies, so enoplean mitogenomes have to be used with utmost care in evolutionary studies.

Chloroplast phylogenomics of unicellular and colonial Volvocales provides perspectives on the evolution of morphological characters

AbstractVolvocales forms a species‐rich clade with wide morphological variety and is regarded as an ideal model for tracing the evolutionary transitions in multicellularity. The phylogenetic relationships among the colonial volvocine algae and its relatives are important for investigating the origin of multicellularity in the clade Reinhardtinia. Therefore, a robust phylogenetic framework of the unicellular and colonial volvocine algae with broad taxon and gene sampling is essential for illuminating the evolution of multicellularity. Recent chloroplast phylogenomic studies have uncovered five major orders in the Chlorophyceae, but the family‐level relationships within Sphaeropleales and Volvocales remain elusive due to the uncertain positions of some incertae sedis taxa. In this study, we contributed six newly sequenced chloroplast genomes in the Volvocales and analyzed a dataset with 91 chlorophycean taxa and 58 protein‐coding genes. Conflicting phylogenetic signals were detected among chloroplast genes that resulted in discordant tree topologies among different analyses. We compared the phylogenetic trees inferred from original nucleotide, RY‐coding, codon‐degenerate, and amino acid datasets, and improved the robustness of phylogenetic inference in the Chlorophyceae by reducing base compositional bias. Our analyses indicate that the unicellular Chlamydomonas and Vitreochlamys are close to or nested within the colonial taxa, and all the incertae sedis taxa are nested within the monophyletic Sphaeropleales s.l. We propose that the colonial taxa in the Reinhardtinia are paraphyletic and multicellularity evolved once in the volvocine green algae and might be lost in Chlamydomonas and Vitreochlamys.

Triplex and other DNA motifs show motif-specific associations with mitochondrial DNA deletions and species lifespan

The “theory of resistant biomolecules” posits that long-lived species show resistance to molecular damage at the level of their biomolecules. Here, we test this hypothesis in the context of mitochondrial DNA (mtDNA) as it implies that predicted mutagenic DNA motifs should be inversely correlated with species maximum lifespan (MLS).First, we confirmed that guanine-quadruplex and direct repeat (DR) motifs are mutagenic, as they associate with mtDNA deletions in the human major arc of mtDNA, while also adding mirror repeat (MR) and intramolecular triplex motifs to a growing list of potentially mutagenic features. What is more, triplex motifs showed disease-specific associations with deletions and an apparent interaction with guanine-quadruplex motifs.Surprisingly, even though DR, MR and guanine-quadruplex motifs were associated with mtDNA deletions, their correlation with MLS was explained by the biased base composition of mtDNA. Only triplex motifs negatively correlated with MLS even after adjusting for body mass, phylogeny, mtDNA base composition and effective number of codons.Taken together, our work highlights the importance of base composition for the comparative biogerontology of mtDNA and suggests that future research on mitochondrial triplex motifs is warranted.

Sequencing and comparative analysis of the chloroplast genome of Ribes odoratum provide insights for marker development and phylogenetics in Ribes.

The online version of this article (10.1007/s12298-021-00932-4) contains supplementary material, which is available to authorized users.

Molecular phylogeny of the orb-weaving spider genus Leucauge and the intergeneric relationships of Leucauginae (Araneae, Tetragnathidae)

The tetragnathid genus Leucauge includes some of the most common orb-weaving spiders in the tropics. Although some species in this genus have attained relevance as model systems for several aspects of spider biology, our understanding of the generic diversity and evolutionary relationships among the species is poor. In this study we present the first attempt to determine the phylogenetic structure within Leucauge and the relationship of this genus with other genera of Leucauginae. This is based on DNA sequences from the five loci commonly used and Histone H4, used for the first time in spider phylogenetics. We also assess the informativeness of the standard markers and test for base composition biases in the dataset. Our results suggest that Leucauge is not monophyletic since species of the genera Opas, Opadometa, Mecynometa and Alcimosphenus are included within the current circumscription of the genus. Based on a phylogenetic re-circumscription of the genus to fulfil the requirement for monophyly of taxa, Leucauge White, 1841 is deemed to be a senior synonym of the genera Opas Pickard-Cambridge, 1896 revalidated synonymy, Mecynometa Simon, 1894 revalidated synonymy, Opadometa Archer, 1951 new synonymy and Alcimosphenus Simon, 1895 new synonymy. We identify groups of taxa critical for resolving relationships within Leucauginae and describe the limitations of the standard loci for accomplishing these resolutions.

Information theoretic perspective on genome clustering.

Shannon’s information theoretic perspective of communication helps one to understand the storage and processing of information in one-dimensional sequences. An information theoretic analysis of 937 available completely sequenced prokaryotic genomes and 238 eukaryotic chromosomes is presented. Information content (Id) values were used to cluster these chromosomes. Chargaff’s second parity rule i.e compositional self-complementarity, an empirical fact is observed in all the genomes, except for the proteobacteria Candidatus Hodgkinia cicadicola. High information content, arising out of biased base composition in all the 14 chromosomes of Plasmodium falciparum is found among two other genomes of prokaryotes viz. Buchnera aphidicola str. Cc (Cinara cedri) and Candidatus Carsonella ruddii PV. Despite size and compositional variations, both prokaryotic and eukaryotic genomes do not deviate significantly from an equiprobable and random situation. Eukaryotic chromosomes of an organism tend to have similar informational restraints as seen when a simple distance based method is used to cluster them. In eukaryotes, in certain cases, Id values are also similar for the two arms (p and q arm) of the chromosomes. The results of this current study confirm that the information content can provide insights into the clustering of genomes and the evolution of messaging strategies of the genomes. An efficient and robust Perl CGI standalone tool is created based on this information theory algorithm for the analysis of the whole genomes and is made available at https://github.com/AlagurajVeluchamy/InformationTheory.

Phylogenomics of Parasitic and Nonparasitic Lice (Insecta: Psocodea): Combining Sequence Data and Exploring Compositional Bias Solutions in Next Generation Data Sets.

The insect order Psocodea is a diverse lineage comprising both parasitic (Phthiraptera) and nonparasitic members (Psocoptera). The extreme age and ecological diversity of the group may be associated with major genomic changes, such as base compositional biases expected to affect phylogenetic inference. Divergent morphology between parasitic and nonparasitic members has also obscured the origins of parasitism within the order. We conducted a phylogenomic analysis on the order Psocodea utilizing both transcriptome and genome sequencing to obtain a data set of 2370 orthologous genes. All phylogenomic analyses, including both concatenated and coalescent methods suggest a single origin of parasitism within the order Psocodea, resolving conflicting results from previous studies. This phylogeny allows us to propose a stable ordinal level classification scheme that retains significant taxonomic names present in historical scientific literature and reflects the evolution of the group as a whole. A dating analysis, with internal nodes calibrated by fossil evidence, suggests an origin of parasitism that predates the K-Pg boundary. Nucleotide compositional biases are detected in third and first codon positions and result in the anomalous placement of the Amphientometae as sister to Psocomorpha when all nucleotide sites are analyzed. Likelihood-mapping and quartet sampling methods demonstrate that base compositional biases can also have an effect on quartet-based methods.[Illumina; Phthiraptera; Psocoptera; quartet sampling; recoding methods.].

Utility of targeted sequence capture for phylogenomics in rapid, recent angiosperm radiations: Neotropical Burmeistera bellflowers as a case study

Targeted sequence capture is a promising approach for large-scale phylogenomics. However, rapid evolutionary radiations pose significant challenges for phylogenetic inference (e.g. incomplete lineages sorting (ILS), phylogenetic noise), and the ability of targeted nuclear loci to resolve species trees despite such issues remains poorly studied. We test the utility of targeted sequence capture for inferring phylogenetic relationships in rapid, recent angiosperm radiations, focusing on Burmeistera bellflowers (Campanulaceae), which diversified into ~130 species over less than 3 million years. We compared phylogenies estimated from supercontig (exons plus flanking sequences), exon-only, and flanking-only datasets with 506–546 loci (~4.7 million bases) for 46 Burmeistera species/lineages and 10 outgroup taxa. Nuclear loci resolved backbone nodes and many congruent internal relationships with high support in concatenation and coalescent-based species tree analyses, and inferences were largely robust to effects of missing taxa and base composition biases. Nevertheless, species trees were incongruent between datasets, and gene trees exhibited remarkably high levels of conflict (~4–60% congruence, ~40–99% conflict) not simply driven by poor gene tree resolution. Higher gene tree heterogeneity at shorter branches suggests an important role of ILS, as expected for rapid radiations. Phylogenetic informativeness analyses also suggest this incongruence has resulted from low resolving power at short internal branches, consistent with ILS, and homoplasy at deeper nodes, with exons exhibiting much greater risk of incorrect topologies due to homoplasy than other datasets. Our findings suggest that targeted sequence capture is feasible for resolving rapid, recent angiosperm radiations, and that results based on supercontig alignments containing nuclear exons and flanking sequences have higher phylogenetic utility and accuracy than either alone. We use our results to make practical recommendations for future target capture-based studies of Burmeistera and other rapid angiosperm radiations, including that such studies should analyze supercontigs to maximize the phylogenetic information content of loci.

Structural and phylogenetic implications of the complete mitochondrial genome of Ledra auditura

We sequenced and annotated the first complete mitochondrial genome (mitogenome) of Ledra auditura (Hemiptera: Cicadellidae: Ledrinae) and reconstructed phylogenetic relationships among 47 species (including 2 outgroup species) on the basis of 3 datasets using maximum likelihood (ML) and Bayesian inference (BI) analyses. The complete L. auditura mitogenome (length, 16,094 bp) comprises 37 genes [13 protein-coding genes (PCGs), 22 tRNAs, and 2 rRNAs], 1 control region, and 2 long non-coding regions. The first long non-coding region (length, 211 bp) is located between tRNA-I and tRNA-Q and the second region (length, 994 bp) between tRNA-S2 and ND1. All PCGs show ATN (Met/Ile) as their start codon and TAR as their stop codon. Except tRNA-S1 (AGN), which lacks the dihydrouridine arm, all tRNAs can fold into the typical cloverleaf secondary structure. The complete L. auditura mitogenome shows a base composition bias of 76.3% A + T (A = 29.9%, T = 46.4%, G = 13.3%, and C = 10.5%), negative AT skew of −0.22, and positive GC skew of 0.12. In ML and BI analyses, L. auditura was clustered with Evacanthus heimianus (Hemiptera: Cicadellidae: Evacanthinae) with strong branch support.

Comprehensive analysis of synonymous codon usage patterns in orf3 gene of porcine epidemic diarrhea virus in China

The ORF3 protein of porcine epidemic diarrhea virus (PEDV) is found to function as an ion channel which influences virus virulence and production. Taking consideration of the importance of PEDV orf3 gene, we have performed comprehensive analysis to investigate its synonymous codon usage patterns. In this study, the results of base composition analysis showed A/T rich and G/C poor in PEDV orf3 genes, and the most abundant base was nucleotide T. The relative synonymous codon usage value in each codon revealed that codon usage bias existed. The mean ENC value of each gene was 48.75, indicating a low codon usage bias, as well as a relatively instable change in PEDV orf3 genes. The general correlation analysis between base composition and codon usage bias indicated that mutational bias has an impact on the PEDV codon usage bias. Neutral analysis suggested that natural selection pressure takes a more important influence than mutational bias in shaping codon usage bias. Moreover, other factors including hydrophobicity and aromaticity have been also found to influence the codon usage variation among the PEDV orf3 genes. This study not only represents the most systematic analysis of codon usage patterns in PEDV orf3 genes, but also provides a basic shaping mechanism of the codon usage bias.

Genetic Diversity of atp6 and cox3 Gene in Wild Drosophila melanogaster

Mitochondria are the energy factory in cells, and more than 95% of the energy in cells is produced by the mitochondrial oxidative phosphorylation (OXPHOS). Therefore, mitochondria likely plays an important role in the process of high-altitude adaptation. Wild Drosophila melanogaster population structure provides a good precondition to detect positive selection on mitochondrial genes related to environmental or climatic variation, specifically hypoxia and high altitudes. In this paper, we sequenced atp6 (ATPase synthase 6) from 27 individuals and cox3 (cytochrome c oxidase III) from 26 individuals of wild Drosophila melanogaster. The percent G+C content of atp6 and cox3 genes has a value of 24 and 29.3%, respectively, exhibiting an extreme bias in base composition. 5 single nucleotide polymorphisms (SNPs) were detected and 2 nonsynonymous substitutions (m.555T > A and m.577A > G) were detected in atp6, only one SNP (m.126C > T) was found from highland Drosophila melanogaster. 6 SNPs were detected in cox3, we found 2 nonsynonymous substitutions (m.745G > A and m.427G > A) of only one highland strain from Yun Nan province. Substitution rates and dN/dS ratios were relatively high for cox3 (ω = 0.34851) compared to atp6 (ω = 0.07329), both of them inferring is purifying selection by using PAML. Only one amino acid site 255I (p G and m.764T > C). From these analyses we conclude that selection may not have played a role in shaping Drosophila melanogaster regional mtDNA variation, the changes in selection is very likely due to some environmental stressors other than hypoxia and high altitudes, and more broadly, our results add to an emerging body of research in Drosophila.

Accelerated Laboratory Evolution Reveals the Influence of Replication on the GC Skew in Escherichia coli.

Most bacterial genomes display contrasting strand asymmetry in a variety of features, such as nucleotide composition and gene orientation, of the two replichores separated by the replication origin and terminus. The cause for the polarization is often attributed to mutations arising from the asymmetric replication machinery. Notably, a base compositional bias known as a GC skew is focused on as a footprint of the bacterial genome evolution driven by DNA replication. Previously, although a replication driven mutation pattern responsible for the GC skew formation or the related mathematical models have been well reported, an exact impact of the replication-related elements on the genomic structure is yet actively debated, and not confirmed experimentally. However, the GC skew formation is very time consuming and challenging in the laboratory. We, therefore, used cytosine deaminase as a DNA mutator, and by monitoring the mutations during an accelerated laboratory evolution procedure with Illumina sequencing, we enabled the trial and error of the GC skew formation in high resolution. Using this technology, we succeeded in reconfirming the influence of bacterial replication machinery on the genomic structure at high resolution.

A molecular phylogeny of the cicadas (Hemiptera: Cicadidae) with a review of tribe and subfamily classification.

A molecular phylogeny and a review of family-group classification are presented for 137 species (ca. 125 genera) of the insect family Cicadidae, the true cicadas, plus two species of hairy cicadas (Tettigarctidae) and two outgroup species from Cercopidae. Five genes, two of them mitochondrial, comprise the 4992 base-pair molecular dataset. Maximum-likelihood and Bayesian phylogenetic results are shown, including analyses to address potential base composition bias. Tettigarcta is confirmed as the sister-clade of the Cicadidae and support is found for three subfamilies identified in an earlier morphological cladistic analysis. A set of paraphyletic deep-level clades formed by African genera are together named as Tettigomyiinae n. stat. Taxonomic reassignments of genera and tribes are made where morphological examination confirms incorrect placements suggested by the molecular tree, and 11 new tribes are defined (Arenopsaltriini n. tribe, Durangonini n. tribe, Katoini n. tribe, Lacetasini n. tribe, Macrotristriini n. tribe, Malagasiini n. tribe, Nelcyndanini n. tribe, Pagiphorini n. tribe, Pictilini n. tribe, Psaltodini n. tribe, and Selymbriini n. tribe). Tribe Tacuini n. syn. is synonymized with Cryptotympanini, and Tryellina n. syn. is synonymized with an expanded Tribe Lamotialnini. Tribe Hyantiini n. syn. is synonymized with Fidicinini. Tribe Sinosenini is transferred to Cicadinae from Cicadettinae, Cicadatrini is moved to Cicadettinae from Cicadinae, and Ydiellini and Tettigomyiini are transferred to Tettigomyiinae n. stat from Cicadettinae. While the subfamily Cicadinae, historically defined by the presence of timbal covers, is weakly supported in the molecular tree, high taxonomic rank is not supported for several earlier clades based on unique morphology associated with sound production.