Nonsynonymous Substitution Rates Research Articles

Relaxed purifying selection pressure drives accelerated and dynamic gene rearrangements in thrips (Insecta: Thysanoptera) mitochondrial genomes

Insect mitochondrial genomes (mitogenome) generally present a typical gene order, which is considered as the ancestral arrangement. All sequenced mitogenomes in the Thysanoptera display high levels of gene rearrangement. Due to limited number of thrips mitogenomes sequenced, how gene rearrangement may be shaped by evolution remain unclear. Here, we analyzed 33 thrips mitogenomes, including 14 newly sequenced. These mitogenomes were diverse in organization, nucleotides substitution and gene arrangements. We found 28 highly rearranged gene orders with the breakpoints of gene rearrangements from 25 to 33. Reconstruction of the ancestors mitochondrial gene arrangements states indicated that Tubulifera have more complex pathways than Terebrantia in the gene order evolution. Molecular calibration estimated that divergence of two suborders occurred in the middle Triassic while the radiation of thrips was associated with the arose and flourish of angiosperm. Our evolutionary hypothesis testing suggests that relaxation of selection pressure enabled the early phase of Thysanoptera evolution, followed by a stronger selective pressure fixed diversification. Our analyses found gene inversion increases the nonsynonymous substitution rates and provide an evolutionary hypothesis driving the diverse gene orders.

Divergent selection and climate adaptation fuel genomic differentiation between sister species of Sphagnum (peat moss).

New plant species can evolve through the reinforcement of reproductive isolation via local adaptation along habitat gradients. Peat mosses (Sphagnaceae) are an emerging model system for the study of evolutionary genomics and have well-documented niche differentiation among species. Recent molecular studies have demonstrated that the globally distributed species Sphagnum magellanicum is a complex of morphologically cryptic lineages that are phylogenetically and ecologically distinct. Here, we describe the architecture of genomic differentiation between two sister species in this complex known from eastern North America: the northern S. diabolicum and the largely southern S. magniae. We sampled plant populations from across a latitudinal gradient in eastern North America and performed whole genome and restriction-site associated DNA sequencing. These sequencing data were then analyzed computationally. Using sliding-window population genetic analyses we find that differentiation is concentrated within 'islands' of the genome spanning up to 400 kb that are characterized by elevated genetic divergence, suppressed recombination, reduced nucleotide diversity and increased rates of non-synonymous substitution. Sequence variants that are significantly associated with genetic structure and bioclimatic variables occur within genes that have functional enrichment for biological processes including abiotic stress response, photoperiodism and hormone-mediated signalling. Demographic modelling demonstrates that these two species diverged no more than 225 000 generations ago with secondary contact occurring where their ranges overlap. We suggest that this heterogeneity of genomic differentiation is a result of linked selection and reflects the role of local adaptation to contrasting climatic zones in driving speciation. This research provides insight into the process of speciation in a group of ecologically important plants and strengthens our predictive understanding of how plant populations will respond as Earth's climate rapidly changes.

Taxonomic Status and Phylogenetic Relationship of the Charadriidae Family Based on Complete Mitogenomes.

The Charadriiformes provide a good source for researching evolution owing to their diverse distribution, behavior, morphology, and ecology. However, in the Charadrii, family-level relationships remain understudied, and the monophyly of Charadriidae is also a subject of controversy. In the present study, we generated complete mitogenomes for two species, Charadrius leschenaultii and Charadrius mongolus, which were found to be 16,905 bp and 16,844 bp in length, respectively. Among the 13 protein codon genes, we observed variation in the rate of non-synonymous substitution rates, with the slowest rate found in COI and the fastest rate observed in ATP8. The Ka/Ks ratio for all Charadriidae species was significantly lower than one, which inferred that the protein-coding genes underwent purifying selection. Phylogenetic analysis based on the genes of Cyt b, 12S and ND2 revealed that the genus Pluvialis is the sister group of three families (Haematopodidae, Ibidorhynchidae, Recurvirostridae). However, the phylogenetic analysis based on complete mitogenomes indicated that the genus Pluvialis is within the Charadriidae family. This study highlights the importance of carefully selecting the number of genes used to obtain accurate estimates of the species tree. It also suggests that relying on partial mtDNA genes with fast-evolving rates may lead to misleading results when resolving the Pluvialis sister group. Future research should focus on sequencing more mitogenomes at different taxonomic levels to gain a better understanding of the features and phylogenetic relationships within the Charadriiformes order.

Sex chromosome heteromorphism and the Fast-X effect in poeciliids.

Fast-X evolution has been observed in a range of heteromorphic sex chromosomes. However, it remains unclear how early in the process of sex chromosome differentiation the Fast-X effect becomes detectible. Recently, we uncovered an extreme variation in sex chromosome heteromorphism across poeciliid fish species. The common guppy, Poecilia reticulata, Endler's guppy, P. wingei, swamp guppy, P. picta and para guppy, P. parae, appear to share the same XY system and exhibit a remarkable range of heteromorphism. Species outside this group lack this sex chromosome system. We combined analyses of sequence divergence and polymorphism data across poeciliids to investigate X chromosome evolution as a function of hemizygosity and reveal the causes for Fast-X effects. Consistent with the extent of Y degeneration in each species, we detect higher rates of divergence on the X relative to autosomes, a signal of Fast-X evolution, in P. picta and P. parae, species with high levels of X hemizygosity in males. In P. reticulata, which exhibits largely homomorphic sex chromosomes and little evidence of hemizygosity, we observe no change in the rate of evolution of X-linked relative to autosomal genes. In P. wingei, the species with intermediate sex chromosome differentiation, we see an increase in the rate of nonsynonymous substitutions on the older stratum of divergence only. We also use our comparative approach to test for the time of origin of the sex chromosomes in this clade. Taken together, our study reveals an important role of hemizygosity in Fast-X evolution.

Genome-wide investigation of NLP gene family members in alfalfa (Medicago sativa L.): evolution and expression profiles during development and stress

BackgroundNIN-like protein (NLP) transcription factors (TFs) compose a plant-specific gene family whose members play vital roles in plant physiological processes, especially in the regulation of plant growth and the response to nitrate-nitrogen. However, no systematic identification or analysis of the NLP gene family has been reported in alfalfa. The recently completed whole-genome sequence of alfalfa has allowed us to investigate genome-wide characteristics and expression profiles.Results53 MsNLP genes were identified from alfalfa and renamed according to their respective chromosome distributions. Phylogenetic analysis demonstrated that these MsNLPs can be classified into three groups on the basis of their conserved domains. Gene structure and protein motif analyses showed that closely clustered MsNLP genes were relatively conserved within each subgroup. Synteny analysis revealed four fragment duplication events of MsNLPs in alfalfa. The ratios of nonsynonymous (Ka) and synonymous (Ks) substitution rates of gene pairs indicated that the MsNLP genes underwent purifying selection during evolution. Examination of the expression patterns of different tissues revealed specific expression patterns of the MsNLP genes in the leaves, indicating that these genes are involved in plant functional development. Prediction of cis-acting regulatory elements and expression profiles further demonstrated that the MsNLP genes might play important roles in the response to abiotic stress and in phytohormone signal transduction processes.ConclusionThis study represents the first genome-wide characterization of MsNLP in alfalfa. Most MsNLPs are expressed mainly in leaves and respond positively to abiotic stresses and hormonal treatments. These results provide a valuable resource for an improved understanding of the characteristics and biological roles of the MsNLP genes in alfalfa.

Structural Characteristics and Phylogenetic Analysis of the Mitochondrial Genomes of Four Krisna Species (Hemiptera: Cicadellidae: Iassinae).

Krisna species are insects that have piercing-sucking mouthparts and belong to the Krisnini tribe in the Iassinae subfamily of leafhoppers in the Cicadellidae family. In this study, we sequenced and compared the mitochondrial genomes (mitogenomes) of four Krisna species. The results showed that all four mitogenomes were composed of cyclic double-stranded molecules and contained 13 protein-coding genes (PCGs) and 22 and 2 genes coding for tRNAs and rRNAs, respectively. Those mitogenomes exhibited similar base composition, gene size, and codon usage patterns for the protein-coding genes. The analysis of the nonsynonymous substitution rate (Ka)/synonymous substitution rate (Ks) showed that evolution occurred the fastest in ND4 and the slowest in COI. 13 PCGs that underwent purification selection were suitable for studying phylogenetic relationships within Krisna. ND2, ND6, and ATP6 had highly variable nucleotide diversity, whereas COI and ND1 exhibited the lowest diversity. Genes or gene regions with high nucleotide diversity can provide potential marker candidates for population genetics and species delimitation in Krisna. Analyses of parity and neutral plots showed that both natural selection and mutation pressure affected the codon usage bias. In the phylogenetic analysis, all subfamilies were restored to a monophyletic group; the Krisnini tribe is monophyletic, and the Krisna genus is paraphyletic. Our study provides novel insights into the significance of the background nucleotide composition and codon usage patterns in the CDSs of the 13 mitochondrial PCGs of the Krisna genome, which could enable the identification of a different gene organization and may be used for accurate phylogenetic analysis of Krisna species.

Integrating massive RNA-seq data to elucidate transcriptome dynamics in Drosophila melanogaster.

The volume of ribonucleic acid (RNA)-seq data has increased exponentially, providing numerous new insights into various biological processes. However, due to significant practical challenges, such as data heterogeneity, it is still difficult to ensure the quality of these data when integrated. Although some quality control methods have been developed, sample consistency is rarely considered and these methods are susceptible to artificial factors. Here, we developed MassiveQC, an unsupervised machine learning-based approach, to automatically download and filter large-scale high-throughput data. In addition to the read quality used in other tools, MassiveQC also uses the alignment and expression quality as model features. Meanwhile, it is user-friendly since the cutoff is generated from self-reporting and is applicable to multimodal data. To explore its value, we applied MassiveQC to Drosophila RNA-seq data and generated a comprehensive transcriptome atlas across 28 tissues from embryogenesis to adulthood. We systematically characterized fly gene expression dynamics and found that genes with high expression dynamics were likely to be evolutionarily young and expressed at late developmental stages, exhibiting high nonsynonymous substitution rates and low phenotypic severity, and they were involved in simple regulatory programs. We also discovered that human and Drosophila had strong positive correlations in gene expression in orthologous organs, revealing the great potential of the Drosophila system for studying human development and disease.

The rate of adaptive molecular evolution in wild and domesticated Saccharomyces cerevisiae populations.

Through its fermentative capacities, Saccharomyces cerevisiae was central in the development of civilisation during the Neolithic period, and the yeast remains of importance in industry and biotechnology, giving rise to bona fide domesticated populations. Here, we conduct a population genomic study of domesticated and wild populations of S. cerevisiae. Using coalescent analyses, we report that the effective population size of yeast populations decreased since the divergence with S. paradoxus. We fitted models of distributions of fitness effects to infer the rate of adaptive ( ) and non-adaptive ( ) non-synonymous substitutions in protein-coding genes. We report an overall limited contribution of positive selection to S. cerevisiae protein evolution, albeit with higher rates of adaptive evolution in wild compared to domesticated populations. Our analyses revealed the signature of background selection and possibly Hill-Robertson interference, as recombination was found to be negatively correlated with and positively correlated with . However, the effect of recombination on was found to be labile, as it is only apparent after removing the impact of codon usage bias on the synonymous site frequency spectrum and disappears if we control for the correlation with , suggesting that it could be an artefact of the decreasing population size. Furthermore, the rate of adaptive non-synonymous substitutions is significantly correlated with the residue solvent exposure, a relation that cannot be explained by the population's demography. Together, our results provide a detailed characterisation of adaptive mutations in protein-coding genes across S. cerevisiae populations.

Rate accelerations in plastid and mitochondrial genomes of Cyperaceae occur in the same clades

Cyperaceae, the second largest family in the monocot order Poales, comprises >5500 species and includes the genus Eleocharis with ∼ 250 species. A previous study of complete plastomes of two Eleocharis species documented extensive structural heteroplasmy, gene order changes, high frequency of dispersed repeats along with gene losses and duplications. To better understand the phylogenetic distribution of gene and intron content as well as rates and patterns of sequence evolution within and between mitochondrial and plastid genomes of Eleocharis and Cyperaceae, an additional 29 Eleocharis organelle genomes were sequenced and analyzed. Eleocharis experienced extensive gene loss in both genomes while loss of introns was mitochondria-specific. Eleocharis has higher rates of synonymous (dS) and nonsynonymous (dN) substitutions in the plastid and mitochondrion than most sampled angiosperms, and the pattern was distinct from other eudicot lineages with accelerated rates. Several clades showed higher dS and dN in mitochondrial genes than in plastid genes. Furthermore, nucleotide substitution rates of mitochondrial genes were significantly accelerated on the branch leading to Cyperaceae compared to most angiosperms. Mitochondrial genes of Cyperaceae exhibited dramatic loss of RNA editing sites and a negative correlation between RNA editing and dS values was detected among angiosperms. Mutagenic retroprocessing and dysfunction of DNA replication, repair and recombination genes are the most likely cause of striking rate accelerations and loss of edit sites and introns in Eleocharis and Cyperaceae organelle genomes.

Genes and sites under adaptation at the phylogenetic scale also exhibit adaptation at the population-genetic scale

Adaptation in protein-coding sequences can be detected from multiple sequence alignments across species or alternatively by leveraging polymorphism data within a population. Across species, quantification of the adaptive rate relies on phylogenetic codon models, classically formulated in terms of the ratio of nonsynonymous over synonymous substitution rates. Evidence of an accelerated nonsynonymous substitution rate is considered a signature of pervasive adaptation. However, because of the background of purifying selection, these models are potentially limited in their sensitivity. Recent developments have led to more sophisticated mutation-selection codon models aimed at making a more detailed quantitative assessment of the interplay between mutation, purifying, and positive selection. In this study, we conducted a large-scale exome-wide analysis of placental mammals with mutation-selection models, assessing their performance at detecting proteins and sites under adaptation. Importantly, mutation-selection codon models are based on a population-genetic formalism and thus are directly comparable to the McDonald and Kreitman test at the population level to quantify adaptation. Taking advantage of this relationship between phylogenetic and population genetics analyses, we integrated divergence and polymorphism data across the entire exome for 29 populations across 7 genera and showed that proteins and sites detected to be under adaptation at the phylogenetic scale are also under adaptation at the population-genetic scale. Altogether, our exome-wide analysis shows that phylogenetic mutation-selection codon models and the population-genetic test of adaptation can be reconciled and are congruent, paving the way for integrative models and analyses across individuals and populations.

HexSE: Simulating evolution in overlapping reading frames.

Gene overlap occurs when two or more genes are encoded by the same nucleotides. This phenomenon is found in all taxonomic domains, but is particularly common in viruses, where it may provide a mechanism to increase the information content of compact genomes. The presence of overlapping reading frames (OvRFs) can skew estimates of selection based on the rates of non-synonymous and synonymous substitutions, since a substitution that is synonymous in one reading frame may be non-synonymous in another and vice versa. To understand the impact of OvRFs on molecular evolution, we implemented a versatile simulation model of nucleotide sequence evolution along a phylogeny with any distribution of open reading frames in linear or circular genomes. We use a custom data structure to track the substitution rates at every nucleotide site, which is determined by the stationary nucleotide frequencies, transition bias and the distribution of selection biases (dN/dS) in the respective reading frames. Our simulation model is implemented in the Python scripting language. All source code is released under the GNU General Public License version 3 and are available at https://github.com/PoonLab/HexSE.

Evolutionary Patterns of the Chloroplast Genome in Vanilloid Orchids (Vanilloideae, Orchidaceae).

The Vanilloideae (vanilloids) is one of five subfamilies of Orchidaceae and is composed of fourteen genera and approximately 245 species. In this study, the six new chloroplast genomes (plastomes) of vanilloids (two Lecanorchis, two Pogonia, and two Vanilla species) were decoded, and then the evolutionary patterns of plastomes were compared to all available vanilloid plastomes. Pogonia japonica has the longest plastome, with 158,200 bp in genome size. In contrast, Lecanorchis japonica has the shortest plastome with 70,498 bp in genome size. The vanilloid plastomes have regular quadripartite structures, but the small single copy (SSC) region was drastically reduced. Two different tribes of Vanilloideae (Pogonieae and Vanilleae) showed different levels of SSC reductions. In addition, various gene losses were observed among the vanilloid plastomes. The photosynthetic vanilloids (Pogonia and Vanilla) showed signs of stage 1 degradation and had lost most of their ndh genes. The other three species (one Cyrotsia and two Lecanorchis), however, had stage 3 or stage 4 degradation and had lost almost all the genes in their plastomes, except for some housekeeping genes. The Vanilloideae were located between the Apostasioideae and Cypripedioideae in the maximum likelihood tree. A total of ten rearrangements were found among ten Vanilloideae plastomes when compared to the basal Apostasioideae plastomes. The four sub-regions of the single copy (SC) region shifted into an inverted repeat (IR) region, and the other four sub-regions of the IR region shifted into the SC regions. Both the synonymous (dS) and nonsynonymous (dN) substitution rates of IR in-cooperated SC sub-regions were decelerated, while the substitution rates of SC in-cooperated IR sub-regions were accelerated. A total of 20 protein-coding genes remained in mycoheterotrophic vanilloids. Almost all these protein genes show accelerated base substitution rates compared to the photosynthetic vanilloids. Two of the twenty genes in the mycoheterotrophic species faced strong "relaxed selection" pressure (p-value < 0.05).

Evaluating the Potential Fitness Effects of Chinook Salmon (Oncorhynchus tshawytscha) Aquaculture Using Non-Invasive Population Genomic Analyses of MHC Nucleotide Substitution Spectra

Genetic diversity plays a vital role in the adaptability of salmon to changing environmental conditions that can introduce new selective pressures on populations. Variability among local subpopulations may increase the chance that certain advantageous genes are passed down to future generations to mitigate susceptibility to novel diseases, warming oceans, loss of genetic stocks, and ocean acidification. Class I and II genes of the major histocompatibility complex (MHC) are crucial for the fitness of Chinook salmon due to the role they play in disease and pathogen resistance. The objective of this study was to assess the DNA sequence variability among wild and hatchery populations of Alaskan Chinook salmon at the class I α1 and class II β1 exons of the MHC. We hypothesized that the 96 wild samples taken from the Deshka River would display greater levels of observed heterozygosity (Ho) relative to expected heterozygosity (He) in suggesting that individuals with similar phenotypes mate with one another more frequently than would be expected under random mating patterns. Conversely, since no mate selection occurs in the William Jack Hernandez Sport Fish hatchery, we would not expect to see this discrepancy (He = Ho) in the 96 hatchery fish tested in this study. Alternatively, we hypothesized that post-mating selection is driving higher levels of observed heterozygosity as opposed to mate selection. If this is the case, we will observe higher than expected levels of heterozygosity among hatchery salmon. Both populations displayed higher levels of observed heterozygosity than expected heterozygosity at the Class I and II loci but genetic differentiation between the spatially distinct communities was minimal. Class I sequences showed evidence of balancing selection, despite high rates of non-synonymous substitutions observed, specifically at the peptide binding regions of both MHC genes.

PANAS: Pipeline and a Case Study to Obtain Synonymous and Nonsynonymous Substitution Rates in Genes of Platyhelminthes

Substitution rates in genes can provide information about their evolutionary dynamics and divergence among either species or protein groups. This work presents PANAS, a computational pipeline developed to facilitate the obtention of substitution rates in large amounts of genomic data from different taxa. We tested PANAS by analyzing the evolution of the mitochondrial protein-coding genes (PCGs) from 18 species of platyhelminths. Our results suggest that the cox1, cox2, cytb, and nad1 genes are under strong selection pressure at the nucleotide and protein level in platyhelminths.

PlantTribes2: Tools for comparative gene family analysis in plant genomics.

Plant genome-scale resources are being generated at an increasing rate as sequencing technologies continue to improve and raw data costs continue to fall; however, the cost of downstream analyses remains large. This has resulted in a considerable range of genome assembly and annotation qualities across plant genomes due to their varying sizes, complexity, and the technology used for the assembly and annotation. To effectively work across genomes, researchers increasingly rely on comparative genomic approaches that integrate across plant community resources and data types. Such efforts have aided the genome annotation process and yielded novel insights into the evolutionary history of genomes and gene families, including complex non-model organisms. The essential tools to achieve these insights rely on gene family analysis at a genome-scale, but they are not well integrated for rapid analysis of new data, and the learning curve can be steep. Here we present PlantTribes2, a scalable, easily accessible, highly customizable, and broadly applicable gene family analysis framework with multiple entry points including user provided data. It uses objective classifications of annotated protein sequences from existing, high-quality plant genomes for comparative and evolutionary studies. PlantTribes2 can improve transcript models and then sort them, either genome-scale annotations or individual gene coding sequences, into pre-computed orthologous gene family clusters with rich functional annotation information. Then, for gene families of interest, PlantTribes2 performs downstream analyses and customizable visualizations including, (1) multiple sequence alignment, (2) gene family phylogeny, (3) estimation of synonymous and non-synonymous substitution rates among homologous sequences, and (4) inference of large-scale duplication events. We give examples of PlantTribes2 applications in functional genomic studies of economically important plant families, namely transcriptomics in the weedy Orobanchaceae and a core orthogroup analysis (CROG) in Rosaceae. PlantTribes2 is freely available for use within the main public Galaxy instance and can be downloaded from GitHub or Bioconda. Importantly, PlantTribes2 can be readily adapted for use with genomic and transcriptomic data from any kind of organism.

Mitochondrial phylogeography and molecular evolution of the rhodopsin visual pigment in troglobitic populations of Astyanax mexicanus (De Filippi, 1853).

Cave-adapted animals provide a unique opportunity to study the evolutionary mechanisms underlying phenotypic, metabolic, behavioral, and genetic evolution in response to cave environments. The Mexican tetra ( Astyanax mexicanus) is considered a unique model system as it shows both surface and cave-dwelling morphs. To date, at least 33 different cave populations have been identified, with phylogenetic studies suggesting an origin from at least two independent surface lineages, thereby providing a unique opportunity to study parallel evolution. In the present study, we carried out the most exhaustive phylogeographic study of A. mexicanus to date, including cave and surface localities, using two mitochondrial markers (cytochrome b (cyt b) and cytochrome c oxidase subunit I ( COI)) and nuclear rhodopsin visual pigment ( rho). Additionally, we inferred the molecular evolution of rho within the two contrasting environments (cave and surface) and across three geographic regions (Sierra de El Abra, Sierra de Guatemala, and Micos). In total, 267 individuals were sequenced for the two mitochondrial fragments and 268 individuals were sequenced for the rho visual pigment from 22 cave and 46 surface populations. Phylogeographic results based on the mitochondrial data supported the two-lineage hypothesis, except for the Pachón and Chica caves, whose introgression has been largely documented. The Sierra de El Abra region depicted the largest genetic diversity, followed by the Sierra de Guatemala region. Regarding the phylogeographic patterns of rho, we recovered exclusive haplogroups for the Sierra de El Abra (Haplogroup I) and Sierra de Guatemala regions (Haplogroup IV). Moreover, a 544 bp deletion in the rho gene was observed in the Escondido cave population from Sierra de Guatemala, reducing the protein from seven to three intramembrane domains. This change may produce a loss-of-function (LOF) but requires further investigation. Regarding nonsynonymous ( dN) and synonymous ( dS) substitution rates (omega values ω), our results revealed the prevailing influence of purifying selection upon the rho pigment for both cave and surface populations (ω<1), but relaxation at the El Abra region. Notably, in contrast to the other two regions, we observed an increase in the number of dN mutations for Sierra de El Abra. However, given that a LOF was exclusively identified in the Sierra de Guatemala region, we cannot dismiss the possibility of a pleiotropic effect on the Rho protein.

The associated evolution among the extensive RNA editing, GC‐biased mutation, and PPR family expansion in the organelle genomes of Selaginellaceae

AbstractExtensive C‐to‐U editing has been reported from plastid genomes (plastomes) and mitochondrial genomes (mitogenomes) of spikemoss. While “reverse” U‐to‐C editing was recorded in other seed‐free plants such as hornworts, quillworts, and ferns, it was not observed in spikemosses. However, no comprehensive study on the association between RNA editing and other genomic features was conducted for the organelle genomes of spikemosses. Here, we report thousands of C‐to‐U editing sites from plastomes and mitogenomes of two species: 1767 and 2394 edits in Selaginella remotifolia, and 4091 and 2786 edits in Selaginella nipponica, respectively. Comparative analyses revealed two different editing frequencies among plastomes, but one similar frequency in mitogenomes. The different editing frequency in the Selaginella organelle genomes is related to the nonsynonymous substitution rate and the genome structural complexity. The high guanine and cytosine (GC) content caused by GC‐biased mutations in organelle genomes might be related to the absence of U‐to‐C editing in Selaginellaceae. Using RNA‐seq and whole‐genome data, we screened the pentatricopeptide repeat (PPR) family and discovered that the number of aspartic acid–tyrosine–tryptophan (DYW) domain‐containing PPR proteins corresponded roughly to the editing abundance in the Selaginella organelle genomes. Consequently, we hypothesize that associated evolution among RNA editing, GC‐biased mutation in organelle genomes, and the PPR protein family encoded in the nuclear genome, is probably triggered by the aberrant DNA repair system in Selaginellaceae. Our study provides new insights into the association between organelle and nuclear genomes in Selaginellaceae, which would contribute to understanding the evolution of post‐transcriptional modifications of organelle genomes in land plants.

Complete mitochondrial genome of the branching octocoral Paramuricea grayi (Johnson, 1861), phylogenetic relationships and divergence analysis

The Gray’s sea fan, Paramuricea grayi (Johnson, 1861), typically inhabits deep littoral and circalittoral habitats of the eastern temperate and tropical Atlantic Ocean. Along the Iberian Peninsula, where P. grayi is a dominant constituent of circalittoral coral gardens, two segregating lineages (yellow and purple morphotypes) were recently identified using single-copy nuclear orthologues. The mitochondrial genomes of 9 P. grayi individuals covering both color morphotypes were assembled from RNA-seq data, using samples collected at three sites in southern (Sagres and Tavira) and western (Cape Espichel) Portugal. The complete circular mitogenome is 18,668 bp in length, has an A + T-rich base composition (62.5%) and contains the 17 genes typically found in Octocorallia: 14 protein-coding genes (atp6, atp8, cob, cox1-3, mt-mutS, nad1-6, and nad4L), the small and large subunit rRNAs (rns and rnl), and one transfer RNA (trnM). The mitogenomes were nearly identical for all specimens, though we identified a noteworthy polymorphism (two SNPs 9 bp apart) in the mt-mutS of one purple individual that is shared with the sister species P. clavata. The mitogenomes of the two species have a pairwise sequence identity of 99.0%, with nad6 and mt-mutS having the highest rates of non-synonymous substitutions.

Complete Mitochondrial Genomes and Phylogenetic Positions of Two Longicorn Beetles, Anoplophora glabripennis and Demonax pseudonotabilis (Coleoptera: Cerambycidae).

Anoplophora glabripennis (Motschulsky, 1854) and Demonax pseudonotabilis Gressitt & Rondon, 1970 are two commonly found longicorn beetles from China. However, the lack of sufficient molecular data hinders the understanding of their evolution and phylogenetic relationships with other species of Cerambycidae. This study sequenced and assembled the complete mitochondrial genomes of the two species using the next-generation sequencing method. The mitogenomes of A. glabripennis and D. pseudonotabilis are 15,622 bp and 15,527 bp in length, respectively. The mitochondrial gene content and gene order of A. glabripennis and D. pseudonotabilis are highly conserved with other sequenced longicorn beetles. The calculation of nonsynonymous (Ka) and synonymous (Ks) substitution rates in PCGs indicated the existence of purifying selection in the two longicorn beetles. The phylogenetic analysis was conducted using the protein-coding gene sequences from available mitogenomes of Cerambycidae. The two species sequenced in this study are, respectively, grouped with their relatives from the same subfamily. The monophyly of Cerambycinae, Dorcasominae, Lamiinae, and Necydalinae was well-supported, whereas Lepturinae, Prioninae, and Spondylidinae were recovered as paraphyletic.

The extremely reduced, diverged and reconfigured plastomes of the largest mycoheterotrophic orchid lineage

BackgroundPlastomes of heterotrophic plants have been greatly altered in structure and gene content, owing to the relaxation of selection on photosynthesis-related genes. The orchid tribe Gastrodieae is the largest and probably the oldest mycoheterotrophic clade of the extant family Orchidaceae. To characterize plastome evolution across members of this key important mycoheterotrophic lineage, we sequenced and analyzed the plastomes of eleven Gastrodieae members, including representative species of two genera, as well as members of the sister group Nervilieae.ResultsThe plastomes of Gastrodieae members contain 20 protein-coding, four rRNA and five tRNA genes. Evolutionary analysis indicated that all rrn genes were transferred laterally and together, forming an rrn block in the plastomes of Gastrodieae. The plastome GC content of Gastrodia species ranged from 23.10% (G. flexistyla) to 25.79% (G. javanica). The plastome of Didymoplexis pallens contains two copies each of ycf1 and ycf2. The synonymous and nonsynonymous substitution rates were very high in the plastomes of Gastrodieae among mycoheterotrophic species in Orchidaceae and varied between genes.ConclusionsThe plastomes of Gastrodieae are greatly reduced and characterized by low GC content, rrn block formation, lineage-specific reconfiguration and gene content, which might be positively selected. Overall, the plastomes of Gastrodieae not only serve as an excellent model for illustrating the evolution of plastomes but also provide new insights into plastome evolution in parasitic plants.