Protein-coding Gene Sequences Research Articles

Phylogenetic and divergence analysis of Pentatomidae, with a comparison of the mitochondrial genomes of two related species (Hemiptera, Pentatomidae).

Pentatomidae, the most diverse family of Pentatomoidea, is found worldwide. Currently, the phylogenetic relationships among Pentatomidae tribes remain unstable, and subfamily divergence has not been estimated. Here, we sequenced and analyzed the complete mitochondrial genomes of two species of Lelia, and studied the phylogenetic relationships among Pentatominae tribes. We also selected three available fossil as the calibration points in the family, and preliminarily discussed the divergence time of Pentatomidae. Trees of Pentatomidae were reconstructed using the Bayesian inference method. Divergence times of Pentatominae were estimated based on the nucleotide sequences of protein-coding genes with a relaxed clock log-normal model in BEASTv.1.8.2. The results showed that the gene arrangements, nucleotide composition, and codon preferences were highly conserved in Lelia. Further, a phylogenetic analysis recovered Eysarcorini, Strachiini, Phyllocephalini, and Menidini as monophyletic with strong support, however, the monophyly of Antestiini, Nezarini, Carpocorini, Pentatomini and Cappaeini were rejected. Moreover, Pentatominae diverged from Pentatomidae soon after the origin of the Cretaceous Period, at approximately 110.38 Ma. This study enriches the mitochondrial genome database of Pentatomidae and provides a reference for further phylogenetic studies, and provides a more accurate estimate of divergence time.

Genome-wide characterization of extrachromosomal circular DNA in SLE and functional analysis reveal their association with apoptosis

Extrachromosomal circular DNA (eccDNA) derived from linear chromosomes, are showed typical nucleosomal ladder pattern in agarose gel which as a known feature of apoptosis and demonstrated to be immunogenicity. In systemic lupus erythematosus (SLE) patients, elevated levels of cell-free DNA (cfDNA) can be found in either linear forms or circular forms, while circular ones are much less common and harder to detect. The molecular characteristics and function of circular forms in plasma SLE patients remains elusive. Herein, we characterized the hallmarks of plasma eccDNA in SLE patients, including the lower normalized number and GC content of eccDNA in SLE plasma than in the healthy, and SLE eccDNA number positively correlated with C3 and negatively with anti-dsDNA antibodies. The differential eccGenes (eccDNAs carrying the protein coding gene sequence) of SLE was significantly enriched in apoptosis-related pathways. The artificially synthesized eccDNA with sequences of the PRF1 exon region could promote transcriptional expression of PRF1, IFNA and IFIT3 and inhibit early-stage apoptosis. Plasma eccDNA can serve as a novel autoantigen in the pathogenesis of SLE.

The complete mitochondrial genome and phylogenetic analysis of Cyanopterus ninghais (Hymenoptera: Braconidae)

Cyanopterus ninghais is an important gregarious ectoparasitoid during the larval stage of Monochamus alternatus, a key vector for pine wilt disease in Asia. In this study, the complete mitochondrial genome of C. ninghais was sequenced and analyzed. The mitochondrial genome of C. ninghais is 15,386 bp in length, comprising 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), and 2 ribosomal RNA genes (rRNAs). The nucleotide composition is 41.32% A, 8.29% G, 6.06% C, and 44.33% T. Phylogenetic trees of Braconidae were constructed using 13 PCG sequences via Bayesian inference (BI) and maximum likelihood (ML) analyses to determine their phylogenetic position. Both ML and BI analyses revealed that C. ninghais is closely related to Euurobracon yokahamae, Virgulibracon endoxylaphagus, and Habrobracon hebetor.

Comparative Analysis of Plastome Sequences of Seven Tulipa L. (Liliaceae Juss.) Species from Section Kolpakowskianae Raamsd. Ex Zonn and Veldk.

Tulipa L. is a genus of significant economic, environmental, and cultural importance in several parts of the world. The exact number of species in the genus remains uncertain due to inherent taxonomic challenges. We utilized next-generation sequencing technology to sequence and assemble the plastid genomes of seven Tulipa species collected in Kazakhstan and conducted a comparative analysis. The total number of annotated genes was 136 in all seven studied Tulipa species, 114 of which were unique, including 80 protein-coding, 30 tRNA, and 4 rRNA genes. Nine regions (petD, ndhH, ycf2-ycf3, ndhA, rpl16, clpP, ndhD-ndhF, rpoC2, and ycf1) demonstrated significant nucleotide variability, suggesting their potential as molecular markers. A total of 1388 SSRs were identified in the seven Tulipa plastomes, with mononucleotide repeats being the most abundant (60.09%), followed by dinucleotide (34.44%), tetranucleotide (3.90%), trinucleotide (1.08%), pentanucleotide (0.22%), and hexanucleotide (0.29%). The Ka/Ks values of the protein-coding genes ranged from 0 to 3.9286, with the majority showing values <1. Phylogenetic analysis based on a complete plastid genome and protein-coding gene sequences divided the species into three major clades corresponding to their subgenera. The results obtained in this study may contribute to understanding the phylogenetic relationships and molecular taxonomy of Tulipa species.

Comparison and phylogenetic analysis of the mitochondrial genomes of Synodontis eupterus and Synodontis polli

We aimed to distinguish Synodontis eupterus and Synodontis polli. We performed sequencing and bioinformatic analysis of their mitochondrial genomes and constructed a phylogenetic tree of Mochokidae fish using maximum likelihood and Bayesian methods based on protein-coding gene (PCG) sequences of 14 Mochokidae species. The total length of the S. eupterus mitochondrial genome was 16,579 bp, including 13 (PCGs), 22 tRNA genes, two rRNA genes, and one D-loop, with an AT-biased nucleotide composition (56.0%). The total length of the S. polli mitochondrial genome was 16,544 bp, including 13 PCGs, 22 tRNA genes, two rRNA genes, and one D-loop, with an AT-biased nucleotide composition (55.0%). In both species, except for COI, PCGs use ATG as the starting codon, the vast majority use TAG or TAA as the ending codon, and a few use incomplete codons (T - or TA -) as the ending codon. Phylogenetic analysis showed that S. eupterus and Synodontis clarias converged into one branch, S. polli and Synodontis petricola converged into one branch, Mochokiella paynei, Mochokus brevis, and nine species of the genus Synodontis converged into one branch, and M. paynei clustered with the genus Synodontis. This study lays a foundation for rebuilding a clearer Mochokidae fish classification system.

Mitochondrial genome of the fluke pond snail, Austropeplea cf. brazieri (Gastropoda: Lymnaeidae)

BackgroundLymnaeid snails of the genus Austropeplea are an important vector of the liver fluke (Fasciola hepatica), contributing to livestock production losses in Australia and New Zealand. However, the species status within Austropeplea is ambiguous due to heavy reliance on morphological analysis and a relative lack of genetic data. This study aimed to characterise the mitochondrial genome of A. cf. brazieri, an intermediate host of liver fluke in eastern Victoria.MethodsThe mitochondrial genome was assembled and annotated from a combination of second- and third-generation sequencing data. For comparative purposes, we performed phylogenetic analyses of the concatenated nucleotide sequences of the mitochondrial protein-coding genes, cytochrome c oxidase subunit 1 and 16S genes.ResultsThe assembled mt genome was 13,757 base pairs and comprised 37 genes, including 13 protein-coding genes, 22 transfer RNA genes and 2 ribosomal RNA genes. The mt genome length, gene order and nucleotide compositions were similar to related species of lymnaeids. Phylogenetic analyses of the mt nucleotide sequences placed A. cf. brazieri within the same clade as Orientogalba ollula with strong statistical supports. Phylogenies of the cox1 and 16S mt sequences were constructed due to the wide availability of these sequences representing the lymnaeid taxa. As expected in both these phylogenies, A. cf. brazieri clustered with other Austropeplea sequences, but the nodal supports were low.ConclusionsThe representative mt genome of A. cf. brazieri should provide a useful resource for future molecular, epidemiology and parasitological studies of this socio-economically important lymnaeid species.Graphical abstract

Highly accurate classification and discovery of microbial protein-coding gene functions using FunGeneTyper: an extensible deep learning framework.

High-throughput DNA sequencing technologies decode tremendous amounts of microbial protein-coding gene sequences. However, accurately assigning protein functions to novel gene sequences remain a challenge. To this end, we developed FunGeneTyper, an extensible framework with two new deep learning models (i.e., FunTrans and FunRep), structured databases, and supporting resources for achieving highly accurate (Accuracy > 0.99, F1-score > 0.97) and fine-grained classification of antibiotic resistance genes (ARGs) and virulence factor genes. Using an experimentally confirmed dataset of ARGs comprising remote homologous sequences as the test set, our framework achieves by-far-the-best performance in the discovery of new ARGs from human gut (F1-score: 0.6948), wastewater (0.6072), and soil (0.5445) microbiomes, beating the state-of-the-art bioinformatics tools and sequence alignment-based (F1-score: 0.0556-0.5065) and domain-based (F1-score: 0.2630-0.5224) annotation approaches. Furthermore, our framework is implemented as a lightweight, privacy-preserving, and plug-and-play neural network module, facilitating its versatility and accessibility to developers and users worldwide. We anticipate widespread utilization of FunGeneTyper (https://github.com/emblab-westlake/FunGeneTyper) for precise classification of protein-coding gene functions and the discovery of numerous valuable enzymes. This advancement will have a significant impact on various fields, including microbiome research, biotechnology, metagenomics, and bioinformatics.

Unraveling the phylogenetics of genetically closely related species, Haemaphysalis japonica and Haemaphysalis megaspinosa, using entire tick mitogenomes and microbiomes

Ticks have a profound impact on public health. Haemaphysalis is one of the most widespread genera in Asia, including Japan. The taxonomy and genetic differentiation of Haemaphysalis spp. is challenging. For instance, previous studies struggled to distinguish Haemaphysalis japonica and Haemaphysalis megaspinosa due to the dearth of nucleotide sequence polymorphisms in widely used barcoding genes. The classification of H. japonica japonica and its related sub-species Haemaphysalis japonica douglasi or Haemaphysalis jezoensis is also confused due to their high morphological similarity and a lack of molecular data that support the current classification. We used mitogenomes and microbiomes of H. japonica and H. megaspinosa to gain deeper insights into the phylogenetic relationships and genetic divergence between two species. Phylogenetic analyses of concatenated nucleotide sequences of protein-coding genes and ribosomal DNA genes distinguished H. japonica and H. megaspinosa as monophyletic clades, with further subdivision within the H. japonica clade. The 16S rRNA and NAD5 genes were valuable markers for distinguishing H. japonica and H. megaspinosa. Population genetic structure analyses indicated that genetic variation within populations accounted for a large proportion of the total variation compared to variation between populations. Microbiome analyses revealed differences in alpha and beta diversity between H. japonica and H. megaspinosa: H. japonica had the higher diversity. Coxiella sp., a likely endosymbiont, was found in both Haemaphysalis species. The abundance profiles of likely endosymbionts, pathogens, and commensals differed between H. japonica and H. megaspinosa: H. megaspinosa was more diverse.

Sequencing and Analysis of the Complete Mitochondrial Genome of Lentipes ikeae.

We sequenced and analyzed the complete mitochondrial genome of Lentipes ikeae and explored the phylogenetic relationships among Sicydiinae based on mitochondrial genome sequences. The complete mitochondrial genome sequence of L. ikeae was determined using the Illumina HiSeq X Ten sequencing platform, and the gene structural characteristics and base composition were analyzed. Based on the mitochondrial genome sequences of 28 Sicydiinae species published in GenBank and mitochondrial protein-coding genes (PCGs), Acanthogobius flavimanus (Gobionellinae) was selected as an outgroup to construct phylogenetic trees of Sicydiinae using the maximum likelihood and Bayesian inference methods. The mitochondrial genome of L. ikeae (GenBank number: OP764680) has a total length of 16,498 bp and encodes 13 PCGs, 22 transfer RNA genes, two ribosomal RNA genes, and a D-loop (control) region. Gene rearrangement is not observed. The mitochondrial genome of L. ikeae exhibits an AT preference, with AT skew > 0 and GC skew < 0 across the entire genome. The phylogenetic relationships of Sicydiinae based on 13 mitochondrial PCG sequences are Sicydium + (Stiphodon + (Sicyopus + Lentipes)) + Sicyopterus, indicating that Sicydium, Sicyopterus, Lentipes, and Stiphodon are all monophyletic groups.

Compositional Features and Codon Usage Pattern of Genes Associated with Parkinson's Disease.

Codon usage bias (CUB) is the phenomenon of non-uniform usage of synonymous codons in which some codons are more used than others and it helps in understanding the molecular organization of genome. Bioinformatic approach was used to analyze the protein-coding sequences of genes associated with Parkinson's disease (PD) to explore compositional features and codon usage pattern as no details work was reported yet. The average improved effective number of codons (Nc) and Nc prime were 42.74 and 44.26 respectively, indicated that CUB was low in these genes. In most of the genes, the overall GC content was almost 50% and GC content at the 1st codon position was the highest while GC content at the 2nd codon position was lowest. Relative synonymous codon usage (RSCU) analysis elucidated over-represented (p > 1.6) and under-represented codons (p < 0.6). The GTG (Val) is the only codon over-represented in all genes. Over-represented codons except (GTG) were A or T ending while under-represented codons (except ACT) were G or C ending. The codons namely TTA (Leu), CTA (Leu), ATC (Ile), ATA (Ile), AGT (Ser), AAC (Asn), TGT (Cys), TGC (Cys), CGC (Arg), AGA (Arg), and AGG (Arg) were absent in SNCA1 to SNCA8 genes. The codon TCG (Ser) was absent in all genes except UCHL1 and PINK1. Correspondence analysis (COA) revealed that the pattern of codon usage differs among genes associated with PD. Neutrality plot analysis indicated some of the points are diagonal distribution suggested that mutation pressure influenced the CUB in genes associated with PD.

Comparative Mitogenome Analyses of Fifteen Ramshorn Snails and Insights into the Phylogeny of Planorbidae (Gastropoda: Hygrophila).

Ramshorn snails from the family Planorbidae are important freshwater snails due to their low trophic level, and some of them act as intermediate hosts for zoonotic trematodes. There are about 250 species from 40 genera of Planorbidae, but only 14 species from 5 genera (Anisus, Biomphalaria, Bulinus, Gyraulus, and Planorbella) have sequenced complete mitochondrial genomes (mitogenomes). In this study, we sequenced and assembled a high-quality mitogenome of a ramshorn snail, Polypylis sp. TS-2018, which represented the first mitogenome of the genus. The mitogenome of Polypylis sp. TS-2018 is 13,749 bp in length, which is shorter than that of most gastropods. It contains 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, and 2 ribosomal RNA (rRNA). We compared mitogenome characteristics, selection pressure, and gene rearrangement among all of the available mitogenomes of ramshorn snails. We found that the nonsynonymous and synonymous substitution rates (Ka/Ks) of most PCGs indicated purifying and negative selection, except for atp8 of Anisus, Biomphalaria, and Gyraulus, which indicated positive selection. We observed that transpositions and reverse transpositions occurred on 10 tRNAs and rrnS, which resulted in six gene arrangement types. We reconstructed the phylogenetic trees using the sequences of PCGs and rRNAs and strongly supported the monophyly of each genus, as well as three tribes in Planorbidae. Both the gene rearrangement and phylogenetic results suggested that Polypylis had a close relationship with Anisus and Gyraulus, while Bulinus was the sister group to all of the other genera. Our results provide useful data for further investigation of species identification, population genetics, and phylogenetics among ramshorn snails.

Complete mitochondrial genome of larged-eye pygmy goby Trimma macrophthalmus (Teleostei, Gobiidae) and its phylogenetic implications

Trimma macrophthalmus is a colorful pygmy goby that is widely associated with coral and distributed in the Indo-Pacific Ocean. Here, the complete mitochondrial genome of T. macrophthalmus was first sequenced and annotated. The entire mitogenome is a typical circular molecule (16,943 bp in total length) containing 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, and 1 control region (CR). The overall base composition is 29.6% for A, 15.6% for G, 27.9% for C and 26.9% for T. Phylogenetic tree was reconstructed using concatenated 13 protein coding genes sequences of the 38 related taxa. The phylogenetic analysis indicated that T. macrophthalmus was most closely related to T. okinawae. These findings will contribute for further genetic studies such as evolution, DNA barcoding and molecular phylogeny of genus Trimma and related gobiid fishes.

The complete mitochondrial genome of Poraniopsis inflata (Asteroidea: Valvatida: Poraniidae) from Dokdo Island, Korea

This study presents the complete mitochondrial genome sequence of Poraniopsis inflata, providing valuable information on its genetic and taxonomic studies. Through next-generation sequencing, we successfully obtained the complete mitogenome of P. inflata, spanning a length of 16,322 bp. This genome structure encompasses 13 protein-coding genes (PCGs), 22 transfer RNA genes, and two ribosomal RNA genes. The phylogenetic analysis, based on a dataset of 13 PCG sequences, illuminated the phylogenetic relationships of P. inflata with other species of class Asteroidea and a species of echinoderm classes. The maximum likelihood phylogenetic tree showed that P. inflata closely clustered with Linckia laevigata. By revealing its mitochondrial genome and positioning it within the Asteroidea lineage, this study provides insights into the phylogenetic context of P. inflata.

Complete mitochondrial genome of Guigarracailaoensis Wang, Chen & Zheng, 2022 (Cypriniformes, Cyprinidae) and its phylogenetic implications.

Guigarracailaoensis is a member of family Cyprinidae, subfamily Labeoninae (Cypriniformes) which was recently discovered in southwestern China. Following its initial description, additional information on this species has remained notably scarce. In the current study, we assemble the complete mitochondrial genome (mitogenome) of G.cailaoensis using the Illumina sequencing platform. The mitogenome is identified as a circular, double-stranded DNA sequence of 16,593 base pairs, encompassing 13 protein-coding genes (PCGs), 22 transfer RNA genes, two ribosomal RNA genes, and a putative control region. Maximum-likelihood and Bayesian-inference approaches were used to construct phylogenetic trees for three datasets: (i) PCG sequences of the complete mitogenome (dataset 1); (ii) PCG sequences of the complete mitogenome combined with nuclear DNA (ncDNA) (Rag1) sequence (dataset 2); and (iii) ncDNA (Rag1) sequences (dataset 3). Phylogenetic analyses position G.cailaoensis as a sister taxon to the lineage consisting of Paraqianlabeolineatus Zhao, Sullivan, Zhang & Peng, 2014 and Pseudogyrinocheilusprochilus Fang, 1933 in dataset 1, and to Pseudogyrinocheilusprochilus in dataset 2, species lacking an oral disc on the lower lip. However, G.cailaoensis showed a close relationship to the lineage consisting of Discogobio and Discocheilus in dataset 3, species possessing an oral disc on the lower lip. Nonetheless, a variety of species with an oral disc on the lower lip are clustered into different lineages across the three datasets that may indicate that the development of the oral disc is homoplastic within the subfamily Labeoninae. The outcomes of this study have the potential to support conservation efforts for this species and to enrich our understanding of genetic resources in the area.

WGCCRR: a web-based tool for genome-wide screening of convergent indels and substitutions of amino acids.

Genome-wide analyses of proteincoding gene sequences are being employed to examine the genetic basis of adaptive evolution in many organismal groups. Previous studies have revealed that convergent/parallel adaptive evolution may be caused by convergent/parallel amino acid changes. Similarly, detailed analysis of lineage-specific amino acid changes has shown correlations with certain lineage-specific traits. However, experimental validation remains the ultimate measure of causality. With the increasing availability of genomic data, a streamlined tool for such analyses would facilitate and expedite the screening of genetic loci that hold potential for adaptive evolution, while alleviating the bioinformatic burden for experimental biologists. In this study, we present a user-friendly web-based tool called WGCCRR (Whole Genome Comparative Coding Region Read) designed to screen both convergent/parallel and lineage-specific amino acid changes on a genome-wide scale. Our tool allows users to replicate previous analyses with just a few clicks, and the exported results are straightforward to interpret. In addition, we have also included amino acid indels that are usually neglected in previous work. Our website provides an efficient platform for screening candidate loci for downstream experimental tests. The tool is available at: https://fishevo.xmu.edu.cn/.

Complete mitochondrial genomes of three vulnerable cave bat species and their phylogenetic relationships within the order Chiroptera.

The IUCN Red List of Threatened Species contains 175 Brazilian bat species that are threatened by extinction in some degree. From this perspective, it is essential to expand the knowledge about the genetic diversity of vulnerable bats. Genomic sequencing can be useful to generate robust and informative genetic references, increasing resolution when analyzing relationships among populations, species, or higher taxonomic levels. In this study, we sequenced and characterized in detail the first complete mitochondrial genomes of Furipterus horrens, Lonchorhina aurita, and Natalus macrourus, and investigated their phylogenetic position based on amino acid sequences of protein-coding genes (PCGs). The mitogenomes of these species are 16,516, 16,697, and 16,668 bp in length, respectively, and each comprises 13 PCGs, 22 tRNA genes, two rRNA genes, and a putative control region (CR). In the three species, genes were arranged similarly to all other previously described bat mitogenomes, and nucleotide composition was also consistent with the reported range. The length and arrangement of rrnS and rrnL were also consistent with those of other bat species, showing a positive AT-skew and a negative GC-skew. Except for trnS1, for which we did not observe the DHU arm, all other tRNAs showed the cloverleaf secondary structure in the three species. In addition, the mitogenomes showed minor differences in start and stop codons, and in all PCGs, codons ending in adenine were more common compared to those ending in guanine. We found that PCGs of the three species use multiple codons to encode each amino acid, following the previously documented pattern. Furthermore, all PCGs are under purifying selection, with atp8 experiencing the most relaxed purifying selection. Considering the phylogenetic reconstruction, F. horrens was recovered as sister to Noctilio leporinus, L. aurita and Tonatia bidens shared a node within Phyllostomidae, and N. macrourus appeared as sister to Molossidae and Vespertilionidae.

The complete chloroplast genome of an Arctic lycophyte, Huperzia arctica (Tolm.) Sipliv. 1973 (Lycopodiaceae)

Huperzia arctica (Tolm.) Sipliv. 1973 is a lycophyte species belonging to the Lycopodiaceae family, which is widely distributed in the Arctic region of Svalbard, Norway. To determine its taxonomic position, we sequenced the chloroplast genome of H. arctica. The complete chloroplast genome of H. arctica is 153,956 bp in length with 122 annotated genes, including 87 protein-coding genes, 31 tRNA genes, and 4 rRNA genes. To evaluate its evolutionary position, we performed phylogenetic analysis using 36 conserved orthologous protein-coding gene sequences from the chloroplast genomes of H. arctica and publicly available data from other Lycopodiaceae. H. arctica formed a monophyletic group with four other Huperzia species: H. lucidula, H. serrata, H. crispata, and H. javanica. However, it appeared as a separate species with a highly supported bootstrap value.

Overlaps Between CDS Regions of Protein-Coding Genes in the Human Genome: A Case Study on the NR1D1-THRA Gene Pair.

For several decades, it has been known that a substantial number of genes within human DNA exhibit overlap; however, the biological and evolutionary significance of these overlaps remain poorly understood. This study focused on investigating specific instances of overlap where the overlapping DNA region encompasses the coding DNA sequences (CDSs) of protein-coding genes. The results revealed that proteins encoded by overlapping CDSs exhibit greater disorder than those from nonoverlapping CDSs. Additionally, these DNA regions were identified as GC-rich. This could be partially attributed to the absence of stop codons from two distinct reading frames rather than one. Furthermore, these regions were found to harbour fewer single-nucleotide polymorphism (SNP) sites, possibly due to constraints arising from the overlapping state where mutations could affect two genes simultaneously.While elucidating these properties, the NR1D1-THRA gene pair emerged as an exceptional case with highly structured proteins and a distinctly conserved sequence across eutherian mammals. Both NR1D1 and THRA are nuclear receptors lacking a ligand-binding domain at their C-terminus, which is the region where these gene pairs overlap. The NR1D1 gene is involved in the regulation of circadian rhythm, while the THRA gene encodes a thyroid hormone receptor, and both play crucial roles in various physiological processes. This study suggests that, in addition to their well-established functions, the specifically overlapping CDS regions of these genes may encode protein segments with additional, yet undiscovered, biological roles.

Microbial symbionts are shared between ants and their associated beetles.

The transmission of microbial symbionts across animal species could strongly affect their biology and evolution, but our understanding of transmission patterns and dynamics is limited. Army ants (Formicidae: Dorylinae) and their hundreds of closely associated insect guest species (myrmecophiles) can provide unique insights into interspecific microbial symbiont sharing. Here, we compared the microbiota of workers and larvae of the army ant Eciton burchellii with those of 13 myrmecophile beetle species using 16S rRNA amplicon sequencing. We found that the previously characterized specialized bacterial symbionts of army ant workers were largely absent from ant larvae and myrmecophiles, whose microbial communities were usually dominated by Rickettsia, Wolbachia, Rickettsiella and/or Weissella. Strikingly, different species of myrmecophiles and ant larvae often shared identical 16S rRNA genotypes of these common bacteria. Protein-coding gene sequences confirmed the close relationship of Weissella strains colonizing army ant larvae, some workers and several myrmecophile species. Unexpectedly, these strains were also similar to strains infecting dissimilar animals inhabiting very different habitats: trout and whales. Together, our data show that closely interacting species can share much of their microbiota, and some versatile microbial species can inhabit and possibly transmit across a diverse range of hosts and environments.

Characterization and Molecular Phylogenetic Analysis of Subfamily Erebinae (Lepidoptera: Noctuoidea: Erebidae) Using Five Complete Mitochondrial Genomes.

In this study, the complete mitogenomes of Sympis rufibasis, Lacera noctilio, Oxyodes scrobiculata, Mocis undata, and Artena dotata were newly sequenced to bring up-to-date the database using the next-generation sequencing methods. The gene order of all sequenced mitogenomes was identical consisting of 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs), and a non-coding A+T-rich region, which were common to other Lepidopteran insects. All protein-coding genes (PCGs) initiated with a canonical ATN codon and ended with TAN or an incomplete stop codon, single T. The A+T-rich region of S. rufibasis, L. noctilio, O. scrobiculata, M. undata, and A. dotata are 406bp, 462bp, 372bp, 410bp, and 406bp long, respectively, containing number of characteristics that are distinctive to Noctuoidea moths. We analyzed concatenated amino acid sequences of protein-coding genes not including rRNAs, using Maximum Likelihood and Bayesian Inference methods. The phylogenetic analyses indicated that the tribe relationships within Erebinae were reconstructed as (Sypnini+((Erebini 1+Poaphilini 1)+((Euclidiini+Catocalini+(Hypopyrini+Erebini 2))+((Hulodini+(Poaphilini 2+Ophiusini))))). Phylogenetic analyses supported and confirmed the monophyly of the subfamilies' relationships as follows: (Hypeninae+Lymantriinae)+((Scoliopterginae+((Calpinae+Erebinae)+((Herminiinae+Aganainae)+Arctiinae)))) within Erebidae.