The complete mitochondrial genome (mitogenome) of Neoperla bimaculata was sequenced and annotated in this study. We found that the mitogenome of N. bimaculata is 15,774 bp in length with an A + T content of 64.3%. It exhibits the classic structure of a mitogenome. Most protein-coding genes (PCGs) of the mitogenome initiate with the standard start codon ATN. Ten PCGs use the standard stop codon TAA/TAG, while the COI, COII, and ND5 genes terminate with a single T nucleotide. Phylogenetic analyses suggested that N. bimaculata, along with two unpublished Neoperla species, formed a cluster within the phylogenetic tree. Our results indicated that the genus Neoperla and Neoperlops were sister groups. Meanwhile, the monophyly of Perlinae and Acroneuriinae was supported in the mitochondrial phylogeny.

Bulbophyllum is one of the largest genera and presents some of the most intricate taxonomic problems in the family Orchidaceae, including species of ornamental and medical importance. The lack of knowledge regarding the characterization of Bulbophyllum chloroplast (cp) genomes has imposed current limitations on our study. Here, we report the complete cp genomes of seven Bulbophyllum species, including B. ambrosia, B. crassipes, B. farreri, B. hamatum, B. shanicum, B. triste, and B. violaceolabellum, and compared with related taxa to provide a better understanding of their genomic information on taxonomy and phylogeny. A total of 28 Bulbophyllum cp genomes exhibit typical quadripartite structures with lengths ranging from 145,092 bp to 165,812 bp and a GC content of 36.60% to 38.04%. Each genome contained 125-132 genes, encompassing 74-86 protein-coding genes, 38 tRNA genes, and eight rRNA genes. The genome arrangements, gene contents, and length were similar, with differences observed in ndh gene composition. It is worth noting that there were exogenous fragment insertions in the IR regions of B. crassipes. A total of 18-49 long repeats and 38-80 simple sequence repeats (SSRs) were detected and the single nucleotide (A/T) was dominant in Bulbophyllum cp genomes, with an obvious A/T preference. An analysis of relative synonymous codon usage (RSCU) revealed that leucine (Leu) was the most frequently used codon, while cysteine (Cys) was the least used. Six highly variable regions (rpl32-trnLUAG > trnTUGU-trnLUAA > trnFGAA-ndhJ > rps15-ycf1 > rbcL-accD > psbI-trnSGCU) and five coding sequences (ycf1 > rps12 > matK > psbK > rps15) were identified as potential DNA markers based on nucleotide diversity. Additionally, 31,641 molecular diagnostic characters (MDCs) were identified in complete cp genomes. A phylogenetic analysis based on the complete cp genome sequences and 68 protein-coding genes strongly supported that 28 Bulbophyllum species can be divided into four branches, sects. Brachyantha, Cirrhopetalum, and Leopardinae, defined by morphology, were non-monophyly. Our results enriched the genetic resources of Bulbophyllum, providing valuable information to illustrate the complicated taxonomy, phylogeny, and evolution process of the genus.

Rubus irritans Focke is a type of tonifying kidney-essence herb used in China. We present the complete chloroplast (cp) genome of R. irritans, a member of the genus Rubus. The complete cp genome of R. irritans was 155,286 bp long and consisted of an 84,613 bp long large single-copy (LSC) region, an 18,697 bp long SSC region, and a pair of 25,988 bp long inverted repeats (IR). Furthermore, the plastid genome contained 130 genes, including 85 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. The overall GC content of the genome was 37.29%. Based on the complete cp genome, phylogenetic analysis revealed that R. irritans is closely related to R. amabilis.

The calmodulin binding transcription activator (CAMTA) is a transcription factor that is widely present in eukaryotes with conserved structure. It contributes to the response to biotic and abiotic stresses and promotes the growth and development of plants. Although previous studies have investigated the number and function of CAMTAs in some species, there is still a lack of comprehensive understanding of the evolutionary process, phylogenetic relationship, expression patterns, and functions of CAMTAs in plants. Here we identified 465 CMATA genes from 112 plants and systematically studied the origin of CAMTA family, gene expansion, functional differentiation, gene structure, and conservative motif distribution. Based on these analyses, we presented the evidence that CAMTA family was originated from chlorophyta, and we speculated that CAMTA might experience obvious structure variation during its early evolution, and that the number of CAMTA genes might gradually increase in higher plants. To reveal potential functions of CAMTA genes, we analyzed the expression patterns of 12 representative species and found significant species specificity, tissue specificity, and developmental stage specificity of CAMTAs. The results also indicated that the CAMTA genes might promote the maturation and senescence. The expression levels and regulatory networks of CAMTAs revealed that CAMTAs could enhance cold tolerance of rice by regulating carbohydrate metabolism-related genes to accumulate carbohydrates or by modulating target genes together with other transcription factors. Our study provides an insight into the molecular evolution of CAMTA family and lays a foundation for further study of related biological functions.

The A. cajennense tick complex has been thought to be the main vector of Rickettsia rickettsii in Central and South America. Studies in Colombia have determined the presence of species A. patinoi and A. mixtum of the A. cajennense complex, but it is unknown which species of this complex exist in northwestern Colombia. Our aim was to identify the species of the A. cajennense complex that are present in northwestern Colombia. We sampled ticks of A. cajennense sensu lato infesting equids. Females identified according to the morphology of their genital pore were selected for genetic confirmation. Specimens from each locality were selected to perform molecular and genetic analysis. Specimens were analyzed from five departments (Antioquia, Bolívar, Córdoba, Magdalena, and Sucre). Morphologically 65 specimens were identified as A. patinoi and 5 as A. mixtum. Molecular analysis allowed to confirm the morphological identification of 27 specimens. In this study A. patinoi was widely distributed in the departments of Antioquia, Bolívar, and Córdoba with allopatric and sympatric distribution in some places. These two species in the region could have unexpected effects on the epidemiology of rickettsiosis.

To better understand the diversity and phylogeny of Perlodidae, we sequenced and annotated the complete mitochondrial genome (mitogenome) of Arcynopteryx dichroa. This mitogenome was 16,215 bp long and encoded 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), 2 ribosomal RNA unit genes (rRNAs), and a control region like other plecopteran. The nucleotide composition of A. dichroa mitochondrial genome obviously biases toward A and T. The A + T content of the whole mitochondrial genome, the PCGs, tRNAs, rRNAs, and the control region were: 69.3%, 67.6%, 69.8%, 71.3%, and 78.7%. Phylogenetic relationship based on the concatenated sequences of 13 PCGs and two ribosomal RNAs showed that the family Perlodidae and Chloroperlidae are sister relationship, family Perlidae being the sister group to the clade (Perlodidae + Chloroperlidae). The monophyly of the family Perlodidae is well supported but the subfamily Perlodinae and Isoperlinae are not monophyletic groups.

BACKGROUND: Hepatitis E virus (HEV) is a common cause of acute viral hepatitis with significant morbidity and mortality, particularly in pregnant women. There are four major genotypes which can cause disease in humans. Genotypes 1 and 2 are usually associated with outbreaks and spread via facal/oral route or contaminated water. Genotypes 3 and 4 are zoonotic and usually associated with handling of pigs or consumption of contaminated pork. The strains circulating in Australia have never been characterized.
UNASSIGNED: The aims for this project are to identify the HEV genotypes found in Australia and link them to possible sources of transmission by phylogenetic analysis.
METHODS: Between 2015 and 2020, 91 HEV isolates were sequenced and genotyped using an in-house PCR. Sixty-six of these were also sequenced by using the international HEVnet primers. Genotypes were determined using the BLASTn program. Relatedness to other strains in Australia was determined by phylogenetic analyses of the HEVnet sequences. Isolates were also stratified by state of origin, gender, age, predisposing factors and travel history (if known).
RESULTS: Of the 91 HEV isolates sequenced, 55 (60.4%) were genotype 1. There were 34 (37.4%) genotype 3 strains and two genotype 4 (2.2%). At least 20 of the genotype 1 strains have been linked to travel in India, and another three with Pakistan. Five of the \"Indian\" strains were closely related and are suspected to have originated in Gujarat. Phylogenetic analysis also showed that 12 genotype 3 strains were genetically related and potentially acquired in/from New South Wales, Australia. The two genotype 4 strains may have originated in China.
CONCLUSIONS: This is the first study to describe the HEV isolates identified in Australia. The results infer that HEV may be acquired during overseas travel as well as locally, presumably from consumption of pork or pork-related products. The phylogenetic analyses also reveal clusters of infection originating from India and Pakistan. This study provides some insight into the source and epidemiology of HEV infection in Australia which may be used to guide public health procedure and enable the implementation of measures to deal with potential outbreaks of infection.

WRKY proteins belong to one of the largest families of transcription factors. They have important functions in plant growth and development, signal transduction and stress responses. However, little information is available regarding the WRKY family in drumstick (Moringa oleifera Lam.). In the present study, we identified 54 MoWRKY genes in this species using genomic data. On the basis of structural features of the proteins they encode, the MoWRKY genes were classified into three main groups, with the second group being further divided into five subgroups. Phylogenetic trees constructed from the sequences of WRKY domains and overall amino acid compositions derived from drumstick and Arabidopsis were similar; the results indicated that the WRKY domain was the main evolutionary unit of WRKY genes. Gene structure and conserved motif analysis showed that genes with similar structures and proteins with similar motif compositions were usually clustered in the same class. Selective pressure analysis indicated that although neutral evolution and positive selection have happened in several MoWRKY genes, most have evolved under strong purifying selection. Moreover, different subgroups had evolved at different rates. The levels of expression of MoWRKY genes in response to five different abiotic stresses (salt, heat, drought, H2O2, cold) were evaluated by reverse transcription polymerase chain reaction (RT-PCR) and quantitative RT-PCR (qRT-PCR), with the results indicating that these genes had different expression levels and that some may be involved in abiotic stress responses. Our results will provide a foundation for cloning genes with specific functions for use in further research and applications.

Bacillus velezensis strain WRN014 was isolated from banana fields in Hainan, China. Bacillus velezensis is an important member of the plant growth-promoting rhizobacteria (PGPR) which can enhance plant growth and control soil-borne disease. The complete genome of Bacillus velezensis WRN014 was sequenced by combining Illumina Hiseq 2500 system and Pacific Biosciences SMRT high-throughput sequencing technologies. Then, the genome of Bacillus velezensis WRN014, together with 45 other completed genome sequences of the Bacillus velezensis strains, were comparatively studied. The genome of Bacillus velezensis WRN014 was 4,063,541bp in length and contained 4,062 coding sequences, 9 genomic islands and 13 gene clusters. The results of comparative genomic analysis provide evidence that (i) The 46 Bacillus velezensis strains formed 2 obviously closely related clades in phylogenetic trees. (ii) The pangenome in this study is open and is increasing with the addition of new sequenced genomes. (iii) Analysis of single nucleotide polymorphisms (SNPs) revealed local diversification of the 46 Bacillus velezensis genomes. Surprisingly, SNPs were not evenly distributed throughout the whole genome. (iv) Analysis of gene clusters revealed that rich gene clusters spread over Bacillus velezensis strains and some gene clusters are conserved in different strains. This study reveals that the strain WRN014 and other Bacillus velezensis strains have potential to be used as PGPR and biopesticide.

BACKGROUND: G-protein coupled receptors (GPCRs) are the largest family of transmembrane receptors in fungi, where they play important roles in signal transduction. Among them, the Pth11-related GPCRs form a large and divergent protein family, and are only found in fungi in Pezizomycotina. However, the evolutionary process and potential functions of Pth11-related GPCRs remain largely unknown.
RESULTS: Twenty genomes of fungi in Pezizomycotina covering different nutritional strategies were mined for putative Pth11-related GPCRs. Phytopathogens encode much more putative Pth11-related GPCRs than symbionts, saprophytes, or entomopathogens. Based on the phylogenetic tree, these GPCRs can be divided into nine clades, with each clade containing fungi in different taxonomic orders. Instead of fungi from the same order, those fungi with similar nutritional strategies were inclined to share orthologs of putative Pth11-related GPCRs. Most of the CFEM domain-containing Pth11-related GPCRs, which were only included in two clades, were detected in phytopathogens. Furthermore, many putative Pth11-related GPCR genes of phytopathogens were upregulated during invasive plant infection, but downregulated under biotic stress. The expressions of putative Pth11-related GPCR genes of saprophytes and entomopathogens could be affected by nutrient conditions, especially the carbon source. The gene expressions revealed that Pth11-related GPCRs could respond to biotic/abiotic stress and invasive plant infection with different expression patterns.
CONCLUSIONS: Our results indicated that the Pth11-related GPCRs existed before the diversification of Pezizomycotina and have been gained and/or lost several times during the evolutionary process. Tandem duplications and trophic variations have been important factors in this evolution.