Characteristics of the mitochondrial genome of Rana omeimontis and related species in Ranidae: Gene rearrangements and phylogenetic relationships

The Omei wood frog (Rana omeimontis), endemic to central China, belongs to the family Ranidae. In this study, we achieved detail knowledge about the mitogenome of the species. The length of the genome is 20,120 bp, including 13 protein‐coding genes (PCGs), 22 tRNA genes, two rRNA genes, and a noncoding control region. Similar to other amphibians, we found that only nine genes (ND6 and eight tRNA genes) are encoded on the light strand (L) and other genes on the heavy strand (H). Totally, The base composition of the mitochondrial genome included 27.29% A, 28.85% T, 28.87% C, and 15.00% G, respectively. The control regions among the Rana species were found to exhibit rich genetic variability and A + T content. R. omeimontis was clustered together with R. chaochiaoensis in phylogenetic tree. Compared to R. amurensis and R. kunyuensi, it was more closely related to R. chaochiaoensis, and a new way of gene rearrangement (ND6‐trnE‐Cytb‐D‐loop‐trnL2 (CUN)‐ND5‐D‐loop) was also found in the mitogenome of R. amurensis and R. kunyuensi. Our results about the mitochondrial genome of R. omeimontis will contribute to the future studies on phylogenetic relationship and the taxonomic status of Rana and related Ranidae species.     

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

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

Mitochondrial genomes of five Hyphessobrycon tetras and their phylogenetic implications

To date, the taxonomic status and phylogenetic affinities within Hyphessobrycon, even among other genera in Characidae, remain unclear. Here, we determined five new mitochondrial genomes (mitogenomes) of Hyphessobrycon species (H. elachys, H. flammeus, H. pulchripinnis, H. roseus, and H. sweglesi). The mitogenomes were all classical circular structures, with lengths ranging from 16,008 to 17,224 bp. The type of constitutive genes and direction of the coding strand that appeared in the mitogenomes were identical to those of other species in Characidae. The highest value of the Ka/Ks ratio within 13 protein‐coding genes (PCGs) was found in ND2 with 0.83, suggesting that they were subject to purifying selection in the Hyphessobrycon genus. Comparison of the control region sequences among seven Hyphessobrycon fish revealed that repeat units differ in length and copy number across different species, which led to sharp differences in mitogenome sizes. Phylogenetic trees based on the 13 PCGs did not support taxonomic relationships, as the Hyphessobrycon fish mixed with those from other genera. These data were combined to explore higher level relationships within Characidae and could aid in the understanding of the evolution of this group.     

Host plant species and carbohydrate supplements affecting rate of multiplication of redlegged earth mite

The effects of young plants of different legume species on the rate of multiplication of Halotydeus destructor (Tucker)(Acarina: Penthaleidae) were examined in one, two or four week experiments in the laboratory. Rates of multiplication of H. destructor were substantially greater on Vicia sativa (common vetch), Phaseolus vulgaris (dwarf bean) and Pisum sativum (snow pea), than on subterranean clover, Trifolium subterraneum; cultivars Dalkeith, Geraldton and Junee. Levels of soluble carbohydrate and soluble nitrogen were greater in vetch, bean and pea, than in the subclovers. Concentrations of 20 amino acids were not associated with differences in mite performance; pea contained higher levels of glutamine and homoserine than the other plants, which may have adversely affected mite multiplication.In further experiments a carbohydrate supplement caused a fourfold increase in rate of multiplication on subclover cv Junee, but had no effect on rate of multiplication on vetch cv Blanchefleur. The carbohydrate supplement did not cause the rate of multiplication on subclover to equal that on vetch. Subclover is not an optimum host plant for H. destructor nutritionally. The nutritional value of host plants is one factor affecting mite multiplication.     

Integrative analysis provides multi‐omics evidence for the pathogenesis of placenta percreta

Pernicious placenta previa with placenta percreta (PP) is a catastrophic condition during pregnancy. However, the underlying pathogenesis remains unclear. In the present study, the placental tissues of normal cases and PP tissues of pernicious placenta previa cases were collected to determine the expression profile of protein‐coding genes, miRNAs, and lncRNAs through sequencing. Weighted gene co‐expression network analysis (WGCNA), accompanied by miRNA target prediction and correlation analysis, were employed to select potential hub protein‐coding genes and lncRNAs. The expression levels of selected protein‐coding genes, Wnt5A and MAPK13, were determined by quantitative PCR and immunohistochemical staining, and lncRNA PTCHD1‐AS and PAPPA‐AS1 expression levels were determined by quantitative PCR and fluorescence in situ hybridization. The results indicated that 790 protein‐coding genes, 382 miRNAs, and 541 lncRNAs were dysregulated in PP tissues, compared with normal tissues. WGCNA identified coding genes in the module (ME) black and ME turquoise modules that may be involved in the pathogenesis of PP. The selected potential hub protein‐coding genes, Wnt5A and MAPK13, were down‐regulated in PP tissues, and their expression levels were positively correlated with the expression levels of PTCHD1‐AS and PAPPA‐AS1. Further analysis demonstrated that PTCHD1‐AS and PAPPA‐AS1 regulated Wnt5A and MAPK13 expression by interacting with specific miRNAs. Collectively, our results provided multi‐omics data to better understand the pathogenesis of PP and help identify predictive biomarkers and therapeutic targets for PP.     

Reproduction of ectoparasitic mites in a coevolved system: Varroa spp.—Eastern honey bees,Apis cerana

Parasite host shifts can impose a high selective pressure on novel hosts. Even though the coevolved systems can reveal fundamental aspects of host–parasite interactions, research often focuses on the new host–parasite relationships. This holds true for two ectoparasitic mite species, Varroa destructor and Varroa jacobsonii, which have shifted hosts from Eastern honey bees, Apis cerana, to Western honey bees, Apis mellifera, generating colony losses of these pollinators globally. Here, we study infestation rates and reproduction of V. destructor and V. jacobsonii haplotypes in 185 A. cerana colonies of six populations in China and Thailand to investigate how coevolution shaped these features. Reproductive success was mostly similar and low, indicating constraints imposed by hosts and/or mite physiology. Infestation rates varied between mite haplotypes, suggesting distinct local co‐evolutionary scenarios. The differences in infestation rates and reproductive output between haplotypes did not correlate with the virulence of the respective host‐shifted lineages suggesting distinct selection scenarios in novel and original host. The occasional worker brood infestation was significantly lower than that of drone brood, except for the V. destructor haplotype (Korea) from which the invasive lineage derived. Whether mites infesting and reproducing in atypical intraspecific hosts (i.e., workers and queens) actually predisposes for and may govern the impact of host shifts on novel hosts should be determined by identifying the underlying mechanisms. In general, the apparent gaps in our knowledge of this coevolved system need to be further addressed to foster the adequate protection of wild and managed honey bees from these mites globally.     

Evolutionary patterns of nucleotide substitution rates in plastid genomes of Quercus

Molecular evolution, including nucleotide substitutions, plays an important role in understanding the dynamics and mechanisms of species evolution. Here, we sequenced whole plastid genomes (plastomes) of Quercus fabri, Quercus semecarpifolia, Quercus engleriana, and Quercus phellos and compared them with 14 other Quercus plastomes to explore their evolutionary relationships using 67 shared protein‐coding sequences. While many previously identified evolutionary relationships were found, our findings do not support previous research which retrieve Quercus subg. Cerris sect. Ilex as a monophyletic group, with sect. Ilex found to be polyphyletic and composed of three strongly supported lineages inserted between sections Cerris and Cyclobalanposis. Compared with gymnosperms, Quercus plastomes showed higher evolutionary rates (Dn/Ds = 0.3793). Most protein‐coding genes experienced relaxed purifying selection, and the high Dn value (0.1927) indicated that gene functions adjusted to environmental changes effectively. Our findings suggest that gene interval regions play an important role in Quercus evolution. We detected greater variation in the intergenic regions (trnH‐psbA, trnK_UUU‐rps16, trnfM_CAU‐rps14, trnS_GCU‐trnG_GCC, and atpF‐atpH), intron losses (petB and petD), and pseudogene loss and degradation (ycf15). Additionally, the loss of some genes suggested the existence of gene exchanges between plastid and nuclear genomes, which affects the evolutionary rate of the former. However, the connective mechanism between these two genomes is still unclear.     

Differential gene expression of the honey bees Apis mellifera and A. cerana induced by Varroa destructor infection

Varroa destructor mite is currently the most serious threat to the world bee industry. Differences in mite tolerance are reported between two honey bee species Apis mellifera and Apis cerana. Differential gene expression of two honey bee species induced by V. destructor infection was investigated by constructing two suppression subtractive hybridization (SSH) libraries, as first steps toward elucidating molecular mechanisms of Varroa tolerance. From the SSH libraries, we obtained 289 high quality sequences which clustered into 132 unique sequences grouped in 26 contigs and 106 singlets where 49 consisted in A. cerana subtracted library and 83 in A. mellifera. Using BLAST, we found that 85% sequences had counterpart known genes whereas 15% were undescribed. A Gene Ontology analysis classified 51 unique sequences into different functional categories. Eight of these differentially expressed genes, representative of different regulation patterns, were confirmed by qRT-PCR. Upon the mite induction, the differentially expressed genes from both bee species were different, except hex 110 gene, which was up-regulated in A. cerana but down-regulated in A. mellifera, and Npy-r gene, which was down-regulated in both species. In general, most of the differential expression genes were involved in metabolic processes and nerve signaling. The results provide information on the molecular response of these two bee species to Varroa infection.     

A method to test compounds for feeding deterrence towards redlegged earth mite (Acarina: Penthaleidae)

A bioassay, based on a membrane sachet technique, has been developed to identify antifeeding compounds affecting the redlegged earth mite, Halotydeus destructor (Acari: Penthaleidae). The method consists of counting H. destructor numbers on membrane sachets in choice experiments, which was quicker and more efficient than weighing the mites. Five per cent aqueous glucose solution was used as a feeding stimulant, with Tween 80® at 5% concentration as a solubilising agent for water-insoluble compounds. (+)-Catechin, rutin, biochanin A, formononetin, chlorogenic acid, and gramine acted as feeding deterrents at 1% concentration. Quercetin (1 %) and azadirachtin (100 ppm) had no significant effect. At lower concentrations (0.01%), compounds showed antifeeding (gramine), phagostimulating (quercetin and chlorogenic acid), or no effects on mite numbers. Dose-dependent deterrent effects of plant extracts were demonstrated with the bioassay, which could be used for other mites.     

Structure and evolution of the mitochondrial genome of Exorista sorbillans: the Tachinidae (Diptera: Calyptratae) perspective

The first complete mitochondrial genome (mitogenome) of Tachinidae Exorista sorbillans (Diptera) is sequenced by PCR-based approach. The circular mitogenome is 14,960?bp long and has the representative mitochondrial gene (mt gene) organization and order of Diptera. All protein-coding sequences are initiated with ATN codon; however, the only exception is Cox I gene, which has a 4-bp ATCG putative start codon. Ten of the thirteen protein-coding genes have a complete termination codon (TAA), but the rest are seated on the H strand with incomplete codons. The mitogenome of E. sorbillans is biased toward A+T content at 78.4?%, and the strand-specific bias is in reflection of the third codon positions of mt genes, and their T/C ratios as strand indictor are higher on the H strand more than those on the L strand pointing at any strain of seven Diptera flies. The length of the A+T-rich region of E. sorbillans is 106?bp, including a tandem triple copies of a13-bp fragment. Compared to Haematobia irritans, E. sorbillans holds distant relationship with Drosophila. Phylogenetic topologies based on the amino acid sequences, supporting that E. sorbillans (Tachinidae) is clustered with strains of Calliphoridae and Oestridae, and superfamily Oestroidea are polyphyletic groups with Muscidae in a clade.     

The complete mitochondrial genome of Eriocheir japonica sinensis (Decapoda: Varunidae) and its phylogenetic analysis

The complete mitochondrial genome (mitogenome) can provide novel insights into understanding the mechanisms underlying mitogenome evolution. In this study, the complete mitogenome of Eriocheir japonica sinensis (Decapoda: Varunidae) was determined to be 16,378 bp, including 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and a D-loop region. The AT skew of the E. j. sinensis mitogenome was slightly negative (−0.016), indicating a higher number of T compared with A nucleotides. The nucleotide composition of the mitogenome was also biased toward A + T nucleotides (71.6%). All PCGs were initiated by ATN codons. Eight of the 13 PCGs harbored the incomplete termination codon by T, or TA. All other tRNA genes displayed a typical clover-leaf structure of mitochondrial tRNA. The D-loop region of the E. j. sinensis mitogenome was 918 bp in length. Based on 13 PCGs, phylogenetic analysis confirmed the placement of E. j. sinensis within the Varunidae.     

The complete mitochondrial genome sequence of the world's largest fish,the whale shark (Rhincodon typus), and its comparison with those of related shark species

The whale shark (Rhincodon typus) is the largest extant species of fish, belonging to the order Orectolobiformes. It is listed as a “vulnerable” species on the International Union for Conservation of Nature (IUCN)'s Red List of Threatened Species, which makes it an important species for conservation efforts. We report here the first complete sequence of the mitochondrial genome (mitogenome) of the whale shark obtained by next-generation sequencing methods. The assembled mitogenome is a 16,875 bp circle, comprising of 13 protein-coding genes, two rRNA genes, 22 tRNA genes and a control region. We also performed comparative analysis of the whale shark mitogenome to the available mitogenome sequences of 17 other shark species, four from the order Orectolobiformes, five from Lamniformes and eight from Carcharhiniformes. The nucleotide composition, number and arrangement of the genes in whale shark mitogenome are the same as found in the mitogenomes of the other members of the order Orectolobiformes and its closest orders Lamniformes and Carcharhiniformes, although the whale shark mitogenome had a slightly longer control region. The availability of mitogenome sequence of whale shark will aid studies of molecular systematics, biogeography, genetic differentiation, and conservation genetics in this species.     

Field associations of first generation densities of the pest mites <Emphasis Type="Italic">Halotydeus destructor</Emphasis> and <Emphasis Type="Italic">Penthaleus major</Emphasis> in pasture

Halotydeus destructor and Penthaleus major are species of earth mite commonly found at high densities in agricultural fields in Australia and other parts of the world. These mites pose a risk to a range of winter crops and pastures when seedlings emerge in autumn. In order to predict likely mite pressure, we investigated whether autumn densities in pastures can be determined from agronomic and environmental field variables. For H. destructor, field densities showed little association with a range of vegetation variables but could largely be explained using the variable field type, with high densities present when fields had mixtures of grass, clover and weeds. For P. major, we found a regional effect. In the region where most data were available, P. major field densities were associated with grass abundance, whereas an association with field type was significant but different to that found for H. destructor. For both species, densities were not associated with rainfall, but there was a weak association with soil moisture capacity. We discuss how these results can help in managing these important pest mites.     

Loss of the IR region in conifer plastomes: Changes in the selection pressure and substitution rate of protein‐coding genes

Plastid genomes (plastomes) have a quadripartite structure, but some species have drastically reduced or lost inverted repeat (IR) regions. IR regions are important for genome stability and the evolution rate. In the evolutionary process of gymnosperms, the typical IRs of conifers were lost, possibly affecting the evolutionary rate and selection pressure of genomic protein‐coding genes. In this study, we selected 78 gymnosperm species (51 genera, 13 families) for evolutionary analysis. The selection pressure analysis results showed that negative selection effects were detected in all 50 common genes. Among them, six genes in conifers had higher ω values than non‐conifers, and 12 genes had lower ω values. The evolutionary rate analysis results showed that 9 of 50 common genes differed between conifers and non‐conifers. It is more obvious that in non‐conifers, the rates of psbA (trst, trsv, ratio, dN, dS, and ω) were 2.6‐ to 3.1‐fold of conifers. In conifers, trsv, ratio, dN, dS, and ω of ycf2 were 1.2‐ to 3.6‐fold of non‐conifers. In addition, the evolution rate of ycf2 in the IR was significantly reduced. psbA is undergoing dynamic change, with an abnormally high evolution rate as a small portion of it enters the IR region. Although conifers have lost the typical IR regions, we detected no change in the substitution rate or selection pressure of most protein‐coding genes due to gene function, plant habitat, or newly acquired IRs.     

食线虫真菌洛斯里被毛孢线粒体基因组的再分析

【目的】鉴定洛斯里被毛孢OWVT-1菌株的线粒体基因组,验证公布的USA-87-5菌株线粒体基因组中的错误,对洛斯里被毛孢正确的线粒体基因组序列进行注释并开展不同被毛孢物种间的比较线粒体基因组学分析。【方法】借助DNA高通量测序数据并通过必要的Sanger测序组装OWVT-1的线粒体基因组。通过PCR验证OWVT-1与公布的USA-87-5线粒体基因组序列差异的真实性。利用多种生物信息方法分析和注释洛斯里被毛孢的线粒体基因组。【结果】公布的洛斯里被毛孢USA-87-5菌株的线粒体基因组存在几处序列错误,包括3处长片段的插入缺失和多处短片段的插入缺失。实际上,洛斯里被毛孢USA-87-5与OWVT-1菌株的线粒体基因组序列完全相同。该菌的线粒体基因组全长62949 bp,在7个基因中共插入13个内含子,部分内含子和基因间区显现出序列退化的特征。洛斯里被毛孢、明尼苏达被毛孢、线虫被毛孢的线粒体基因组具有较强的共线性关系。除一些独立的ORF外,核心蛋白编码基因、rRNA基因和tRNA基因的排列顺序非常保守。基因间区的长短是影响3种被毛孢线粒体基因组大小最主要的因素。【结论】公布的洛斯里被毛孢USA-87-5菌株线粒体基因组中存在序列错误。本文新报道了OWVT-1菌株的线粒体基因组,并进行注释和比较线粒体基因组学分析。     

Chromosome‐level genome assembly for the horned‐gall aphid provides insights into interactions between gall‐making insect and its host plant

The aphid Schlechtendalia chinensis is an economically important insect that can induce horned galls, which are valuable for the medicinal and chemical industries. Up to now, more than twenty aphid genomes have been reported. Most of the sequenced genomes are derived from free‐living aphids. Here, we generated a high‐quality genome assembly from a galling aphid. The final genome assembly is 271.52 Mb, representing one of the smallest sequenced genomes of aphids. The genome assembly is based on contig and scaffold N50 values of the genome sequence are 3.77 Mb and 20.41 Mb, respectively. Nine‐seven percent of the assembled sequences was anchored onto 13 chromosomes. Based on BUSCO analysis, the assembly involved 96.9% of conserved arthropod and 98.5% of the conserved Hemiptera single‐copy orthologous genes. A total of 14,089 protein‐coding genes were predicted. Phylogenetic analysis revealed that S. chinensis diverged from the common ancestor of Eriosoma lanigerum approximately 57 million years ago (MYA). In addition, 35 genes encoding salivary gland proteins showed differentially when S. chinensis forms a gall, suggesting they have potential roles in gall formation and plant defense suppression. Taken together, this high‐quality S. chinensis genome assembly and annotation provide a solid genetic foundation for future research to reveal the mechanism of gall formation and to explore the interaction between aphids and their host plants.     

Genetic markers in Andean Puya species (Bromeliaceae) with implications on plastome evolution and phylogeny

The Andean plant endemic Puya is a striking example of recent and rapid diversification from central Chile to the northern Andes, tracking mountain uplift. This study generated 12 complete plastomes representing nine Puya species and compared them to five published plastomes for their features, genomic evolution, and phylogeny. The total size of the Puya plastomes ranged from 159,542 to 159,839 bp with 37.3%–37.4% GC content. The Puya plastomes were highly conserved in organization and structure with a typical quadripartite genome structure. Each of the 17 consensus plastomes harbored 133 genes, including 87 protein‐coding genes, 38 tRNA (transfer RNA) genes, and eight rRNA (ribosomal RNA) genes; we found 69–78 tandem repeats, 45–60 SSRs (simple sequence repeats), and 8–22 repeat structures among 13 species. Four protein‐coding genes were identified under positive site‐specific selection in Puya. The complete plastomes and hypervariable regions collectively provided pronounced species discrimination in Puya and a practical tool for future phylogenetic studies. The reconstructed phylogeny and estimated divergence time for the lineage suggest that the diversification of Puya is related to Andean orogeny and Pleistocene climatic oscillations. This study provides plastome resources for species delimitation and novel phylogenetic and biogeographic studies.