Molecular characterization of the meq oncogene of Marek's disease virus in vaccinated Brazilian poultry farms reveals selective pressure on prevalent strains

Marek's disease virus (MDV) has become an increasingly virulent pathogen in the poultry industry despite vaccination efforts to control it. Brazil has experienced a significant rise of Marek's disease (MD) outbreaks in recent years. Our study aimed to analyze the complete meq gene sequences to understand the molecular epidemiological basis of MD outbreaks in Brazilian vaccinated layer farms. We detected a high incidence rate of visceral MD (67.74%) and multiple circulating MDV strains. The most prevalent and geographically widespread genotype presented several clinical and molecular characteristics of a highly virulent strain and evolving under positive selective pressure. Phylogenetic and phylogeographic analysis revealed a closer relationship with strains from the USA and Japan. This study sheds light on the circulation of MDV strains capable of infecting vaccinated birds. We emphasize the urgency of adopting preventive measures to manage MDV outbreaks threatening the poultry farming industry.

Selective and non-selective evolutionary signatures found in the simplest replicative biological entities

Mitoviruses, which are considered evolutionary relics of extinct alpha-proteobacteria RNA phages, represent one of the simplest self-replicating biological systems. This study aims to quantitatively describe genomes and identify potential genomic signatures that support the protein phylogenetic-based classification criterion. Genomic variables, such as mononucleotide and dinucleotide composition, codon usage bias, and minimal free energy derived from optimized predicted RNA secondary structure, were analyzed. From the values obtained, the main evolutionary pressures were discussed, indicating that natural selection plays a significant role in shaping mitovirus genomes. However, neutral evolution also makes a significant contribution. This study reveals a significant discovery of structural divergence in Kvaramitovirus. The energy minimization approach employed to study 2D folding in this study reveals a distinct spatial organization of their genomes, providing evidence for the hypothesis of a single evolutionary event of circularization in the most recent common ancestor of the lineage. This hypothesis was discussed in light of recent discoveries by other researchers that partially support the existence of mitoviruses with circular genomes. Finally, this study represents a significant advancement in the understanding of mitoviruses, as it quantitatively describes the nucleotide sequence at the family and genus taxonomic levels. Additionally, we provide hypotheses that can be experimentally validated to inspire new research and address the gaps in knowledge of this fascinating, basally divergent RNA virus lineage.

Multidrug-resistant Stenotrophomonas maltophilia in residential aged care facilities: An emerging threat

Stenotrophomonas maltophilia is a multidrug-resistant (MDR), Gram-negative bacterium intrinsically resistant to beta-lactams, including last-resort carbapenems. As an opportunistic pathogen, it can cause serious healthcare-related infections. This study assesses the prevalence, resistance profiles, and genetic diversity of S. maltophilia isolated from residential aged care facilities (RACFs). RACFs are known for their overuse and often inappropriate use of antibiotics, creating a strong selective environment that favors the development of bacterial resistance. The study was conducted on 73 S. maltophilia isolates recovered from wastewater and facility swab samples obtained from three RACFs and a retirement village. Phenotypic and genotypic assessments of the isolates revealed high carbapenem resistance, exemplifying their intrinsic beta-lactam resistance. Alarmingly, 49.3% (36/73) of the isolates were non-wild type for colistin, with minimum inhibitory concentration values of > 4 mg/L, and 11.0% (8/73) were resistant to trimethoprim-sulfamethoxazole. No resistance mechanisms were detected for either antimicrobial. Genotypic assessment of known lineages revealed isolates clustering with Sm17 and Sm18, lineages not previously reported in Australia, suggesting the potential ongoing spread of MDR S. maltophilia. Lastly, although only a few isolates were biocide tolerant (2.7%, 2/73), their ability to grow in high concentrations (64 mg/L) of triclosan is concerning, as it may be selecting for their survival and continued dissemination.

Sequencing of the Complete Mitochondrial Genome of the Big Brown Mactra Clam, Mactra grandis (Venerida: Mactridae)

Mitochondrial genomes are playing an increasingly important role in molluscan taxonomy, germplasm, and evolution studies. The first complete mitochondrial genome of the commercial big brown mactra clam, Mactra grandis, was characterized using Illumina next-generation sequencing in this study. The 17,289 bp circular genome has a typical gene organization of 13 protein-coding genes (PCGs), 2 rRNAs, and 22 tRNAs, with an obvious (A + T)-bias of 64.54%. All PCGs exhibited a homogeneous bias in nucleotide composition with a (A + T)-bias, a positive GC skew, and a negative AT skew. Results of phylogenetic analysis showed that Mactra grandis was most closely related to Mactra cygnus. The functional gene arrangement of the two species was identical but different from other Mactra species. The congeneric relationships among Mactra species were demonstrated by genetic distance analysis. Additionally, the selective pressure analysis suggested that cox1 was highly efficient for discriminating closely related species in genus Mactra, while nad2 was the most appropriate marker for population genetic analysis.

Selective Pressure and Evolution of SARS-CoV-2 Lineages BF.7 and BQ.1.1 Circulating in Italy from July to December 2022

In this work, we studied the selective pressure and evolutionary analysis on the SARS-CoV-2 BF.7 and BQ.1.1 lineages circulating in Italy from July to December 2022. Two different datasets were constructed: the first comprised 694 SARS-CoV-2 BF.7 lineage sequences and the second comprised 734 BQ.1.1 sequences, available in the Italian COVID-19 Genomic (I-Co-Gen) platform and GISAID (last access date 15 December 2022). Alignments were performed with MAFFT v.7 under the Galaxy platform. The HYPHY software was used to study the selective pressure. Four positively selected sites (two in nsp3 and two in the spike) were identified in the BF.7 dataset, and two (one in ORF8 and one in the spike gene) were identified in the BQ.1.1 dataset. Mutation analysis revealed that R408S and N440K are very common in the spike of the BF.7 genomes, as well as L452R among BQ.1.1. N1329D and Q180H in nsp3 were found, respectively, at low and rare frequencies in BF.7, while I121L and I121T were found to be rare in ORF8 for BQ.1.1. The positively selected sites may have been driven by the selection for increased viral fitness, under circumstances of defined selective pressure, as well by host genetic factors.

The rapidly evolving X-linked MIR-506 family fine-tunes spermatogenesis to enhance sperm competition

Despite rapid evolution across eutherian mammals, the X-linked MIR-506 family miRNAs are located in a region flanked by two highly conserved protein-coding genes (SLITRK2 and FMR1) on the X chromosome. Intriguingly, these miRNAs are predominantly expressed in the testis, suggesting a potential role in spermatogenesis and male fertility. Here, we report that the X-linked MIR-506 family miRNAs were derived from the MER91C DNA transposons. Selective inactivation of individual miRNAs or clusters caused no discernible defects, but simultaneous ablation of five clusters containing 19 members of the MIR-506 family led to reduced male fertility in mice. Despite normal sperm counts, motility, and morphology, the KO sperm were less competitive than wild-type sperm when subjected to a polyandrous mating scheme. Transcriptomic and bioinformatic analyses revealed that these X-linked MIR-506 family miRNAs, in addition to targeting a set of conserved genes, have more targets that are critical for spermatogenesis and embryonic development during evolution. Our data suggest that the MIR-506 family miRNAs function to enhance sperm competitiveness and reproductive fitness of the male by finetuning gene expression during spermatogenesis.

Characteristics and Comparative Analysis of the Special-Structure (Non-Single-Circle) Mitochondrial Genome of Capsicum pubescens Ruiz & Pav

Chilean peppers, cultivated from Capsicum pubescens, are globally renowned as popular vegetable and spice crops. C. pubescens belongs to the Capsicum L. (pepper) family and is one of the five pepper cultivars grown in China. In this study, we assembled and annotated the complete mt genome of C. pubescens. We investigated several aspects of its genome, including characteristics, codon usage, RNA editing sites, repeat sequences, selective pressure, gene clusters, and phylogenetic relationships. Furthermore, we compared it with other plant mt genomes. The data we obtained will provide valuable information for studying evolutionary processes in the Capsicum genus and will assist in the functional analysis of Capsicum mitogenomes.

The Rapidly Evolving X-linked miR-506 Family Finetunes Spermatogenesis to Enhance Sperm Competition

Despite rapid evolution across eutherian mammals, the X-linked miR-506 family miRNAs are located in a region flanked by two highly conserved protein-coding genes (Slitrk2 and Fmr1) on the X chromosome. Intriguingly, these miRNAs are predominantly expressed in the testis, suggesting a potential role in spermatogenesis and male fertility. Here, we report that the X-linked miR-506 family miRNAs were derived from the MER91C DNA transposons. Selective inactivation of individual miRNAs or clusters caused no discernable defects, but simultaneous ablation of five clusters containing nineteen members of the miR-506 family led to reduced male fertility in mice. Despite normal sperm counts, motility and morphology, the KO sperm were less competitive than wild-type sperm when subjected to a polyandrous mating scheme. Transcriptomic and bioinformatic analyses revealed that these X-linked miR-506 family miRNAs, in addition to targeting a set of conserved genes, have more targets that are critical for spermatogenesis and embryonic development during evolution. Our data suggest that the miR-506 family miRNAs function to enhance sperm competitiveness and reproductive fitness of the male by finetuning gene expression during spermatogenesis.

Tissue-Predisposition to Cancer Driver Mutations

Driver mutations are considered the cornerstone of cancer initiation. They are defined as mutations that convey a competitive fitness advantage, and hence, their mutation frequency in premalignant tissue is expected to exceed the basal mutation rate. In old terms, that translates to "the survival of the fittest" and implies that a selective process underlies the frequency of cancer driver mutations. In that sense, each tissue is its own niche that creates a molecular selective pressure that may favor the propagation of a mutation or not. At the heart of this stands one of the biggest riddles in cancer biology: the tissue-predisposition to cancer driver mutations. The frequency of cancer driver mutations among tissues is non-uniform: for instance, mutations in APC are particularly frequent in colorectal cancer, and 99% of chronic myeloid leukemia patients harbor the driver BCR-ABL1 fusion mutation, which is rarely found in solid tumors. Here, we provide a mechanistic framework that aims to explain how tissue-specific features, ranging from epigenetic underpinnings to the expression of viral transposable elements, establish a molecular basis for selecting cancer driver mutations in a tissue-specific manner.

Evolutionary Analysis of a Parrot Bornavirus 2 Detected in a Sulphur-Crested Cockatoo (Cacatua galerita) Suggests a South American Ancestor

Parrot bornavirus (PaBV) is an RNA virus that causes Proventricular Dilatation Disease (PDD), neurological disorders, and death in Psittaciformes. Its diversity in South America is poorly known. We examined a Cacatua galerita presenting neuropathies, PDD, and oculopathies as the main signs. We detected PaBV through reverse transcription polymerase chain reaction (RT-PCR) and partial sequencing of the nucleoprotein (N) and matrix (M) genes. Maximum likelihood and Bayesian phylogenetic inferences classified it as PaBV-2. The nucleotide identity of the sequenced strain ranged from 88.3% to 90.3% against genotype PaBV-2 and from 80.2% to 84.4% against other genotypes. Selective pressure analysis detected signs of episodic diversifying selection in both the N and M genes. No recombination events were detected. Phylodynamic analysis estimated the time to the most recent common ancestor (TMRCA) as the year 1758 for genotype PaBV-2 and the year 1049 for the Orthobornavirus alphapsittaciforme species. Substitution rates were estimated at 2.73 × 10-4 and 4.08 × 10-4 substitutions per year per site for N and M, respectively. The analysis of population dynamics showed a progressive decline in the effective population size during the last century. Timescale phylogeographic analysis revealed a potential South American ancestor as the origin of genotypes 1, 2, and 8. These results contribute to our knowledge of the evolutionary origin, diversity, and dynamics of PaBVs in South America and the world. Additionally, it highlights the importance of further studies in captive Psittaciformes and the potential impact on endangered wild birds.

Global Prevalence and Hemagglutinin Evolution of H7N9 Avian Influenza Viruses from 2013 to 2022

H7N9 avian influenza viruses have caused severe harm to the global aquaculture industry and human health. For further understanding of the characteristics of prevalence and hemagglutinin evolution of H7N9 avian influenza viruses, we generated the global epidemic map of H7N9 viruses from 2013 to 2022, constructed a phylogenetic tree, predicted the glycosylation sites and compared the selection pressure of the hemagglutinin. The results showed that although H7N9 avian influenza appeared sporadically in other regions worldwide, China had concentrated outbreaks from 2013 to 2017. The hemagglutinin genes were classified into six distinct lineages: A, B, C, D, E and F. After 2019, H7N9 viruses from the lineages B, E and F persisted, with the lineage B being the dominant. The hemagglutinin of highly pathogenic viruses in the B lineage has an additional predicted glycosylation site, which may account for their persistent pandemic, and is under more positive selection pressure. The most recent ancestor of the H7N9 avian influenza viruses originated in September 1991. The continuous evolution of hemagglutinin has led to an increase in virus pathogenicity in both poultry and humans, and sustained human-to-human transmission. This study provides a theoretical basis for better prediction and control of H7N9 avian influenza.

Tracking the Selective Pressure Profile and Gene Flow of SARS-CoV-2 Delta Variant in Italy from April to October 2021 and Frequencies of Key Mutations from Three Representative Italian Regions

The SARS-CoV-2 Delta variant of concern (VOC) was often associated with serious clinical course of the COVID-19 disease. Herein, we investigated the selective pressure, gene flow and evaluation on the frequencies of mutations causing amino acid substitutions in the Delta variant in three Italian regions. A total of 1500 SARS-CoV-2 Delta genomes, collected in Italy from April to October 2021 were investigated, including a subset of 596 from three Italian regions. The selective pressure and the frequency of amino acid substitutions and the prediction of their possible impact on the stability of the proteins were investigated. Delta variant dataset, in this study, identified 68 sites under positive selection: 16 in the spike (23.5%), 11 in nsp2 (16.2%) and 10 in nsp12 (14.7%) genes. Three of the positive sites in the spike were located in the receptor-binding domain (RBD). In Delta genomes from the three regions, 6 changes were identified as very common (>83.7%), 4 as common (>64.0%), 21 at low frequency (2.1%-25.0%) and 29 rare (≤2.0%). The detection of positive selection on key mutations may represent a model to identify recurrent signature mutations of the virus.

New Generations of Tyrosine Kinase Inhibitors in Treating NSCLC with Oncogene Addiction: Strengths and Limitations

Tyrosine kinase inhibitors (TKIs) revolutionized the treatment of patients with advanced or metastatic non-small cell lung cancer (NSCLC) harboring most driver gene alterations. Starting from the first generation, research rapidly moved to the development of newer, more selective generations of TKIs, obtaining improved results in terms of disease control and survival. However, the use of novel generations of TKIs is not without limitations. We reviewed the main results obtained, as well as the ongoing clinical trials with TKIs in oncogene-addicted NSCLC, together with the biology underlying their potential strengths and limitations. Across driver gene alterations, novel generations of TKIs allowed delayed resistance, prolonged survival, and improved brain penetration compared to previous generations, although with different toxicity profiles, that generally moved their use from further lines to the front-line treatment. However, the anticipated positioning of novel generation TKIs leads to abolishing the possibility of TKI treatment sequencing and any role of previous generations. In addition, under the selective pressure of such more potent drugs, resistant clones emerge harboring more complex and hard-to-target resistance mechanisms. Deeper knowledge of tumor biology and drug properties will help identify new strategies, including combinatorial treatments, to continue improving results in patients with oncogene-addicted NSCLC.

The evolution of bitter taste receptor gene in primates: Gene duplication and selection

Bitter taste perception plays an important role in preventing animals from digesting poisonous and harmful substances. In primates, especially the Cercopithecidae species, most species feed on plants; thus, it is reasonable to speculate that most of the bitter taste receptor genes (T2Rs) of primates are under purifying selection to maintain the functional stability of bitter taste perception. Gene duplication has happened in T2Rs frequently, and what will be the fate of T2Rs copies is another question we are concerned about. To answer these questions, we selected the T2Rs of primates reported in another study and conducted corresponding selective pressure analyses to determine what kind of selective pressure was acting on them. Further, we carried out selective pressure analyses on gene copies and their corresponding ancestors by considering several possible situations. The results showed that among the 25 gene groups examined here, 15 groups are subject to purifying selection and others are under relaxed selection, with many positively selected sites detected. Gene copies existed in several groups, but only some groups (clade1_a1-b2, clade1_c-c2, clade1_d1-d3, clade1_f1-f2, T2R10, T2R13, and T2R42) have positively selected sites, inferring that they may have some relation to functional divergence. Taken together, T2Rs in primates are under diverse selective pressures, and most gene copies are subject to the same selective pressures. In such cases, the copies may be just to keep the function conservative, and more copies can increase the quantity of the bitter taste receptor, raise the efficiency of bitter substance recognition, and finally enhance the fitness of feeding during the evolutionary course of primates. This study can improve our understanding of T2Rs evolution in primates.

Antigenic and molecular characterization of isolates of the Brazilian genotype BR-I (GI-11) of infectious bronchitis virus supports its recognition as BR-I serotype

The antigenic and molecular characteristics of BR-I infectious bronchitis viruses (IBVs) isolated from Brazil are reported. IBVs isolated from commercial flocks with different clinical manifestations between 2003 and 2019 were submitted to antigenic and molecular characterization. The complete S1 glycoprotein gene of 11 field isolates was amplified and sequenced. The virus neutralization (VN) test showed 94.75% neutralization with a BR-I isolate and 30% or less against other worldwide reference strains. The nucleotide and amino acid sequence analyses revealed 84.3-100% and 83.5-100% identity among them, respectively. The identity values ranged from 57.1 to 82.6% for nucleotides and from 46.6-84.4% for amino acids compared with those of other genotypes. By phylogenetic tree analysis, the Brazilian isolates were branched into the BR-I genotype (lineage GI-11), which was differentiated from foreign reference strains. Selective pressure analyses of BR-I IBVs revealed evolution under purifying selection (negative pressure) for the complete S1 gene but four specific sites (87, 121, 279, and 542) under diversifying selection (positive pressure). Profiles of cleavage sites and potential N-glycosylation sites differed from those of other genotypes. The low molecular relationship among the Brazilian viruses and foreign serotypes was concordant with the VN test results. The low antigenic relatedness (ranging from 5.3-30% between Brazilian genotype BR-I and reference IBV serotypes of North America, Europe, and Asia) indicates that the BR-I genotype is a different serotype, referred to for the first time and hereafter as serotype BR-I. RESEARCH HIGHLIGHTSStrains of the BR-I genotype presented robust antigenic and molecular similarity.BR-I strains evolved under purifying selection mode (negative pressure).The BR-I genotype originated in Brazil and dispersed to other countries.BR-I genotype viruses can be referred to as the BR-I serotype.

Influence of Different Evolutive Forces on GDF5 Gene Variability

The GDF5 gene is involved in the development of skeletal elements, synovial joint formation, tendons, ligaments, and cartilage. Several polymorphisms are present within the gene, and two of them, rs143384 and 143383, were reported to be correlated with osteoarticular disease or muscle flexibility. The aim of this research is to verify if the worldwide distribution of the rs143384 polymorphism among human populations was shaped by selective pressure, or if it was the result of random genetic drift events. Ninety-four individuals of both the male and female sexes, 18-28 years old, from Sardinia were analyzed. We observed the following genotype frequencies: 28.72% of AA homozygotes, 13.83% of GG homozygotes, and 57.45% of AG heterozygotes. The allele frequencies were 0.574 for allele A and 0.426 for allele G. The relationships between the populations were verified via Multidimensional Scaling (MDS). Our data show (i) a clear heterogeneity within the African populations; (ii) a strong differentiation between the African populations and the other populations; and that (iii) the Sardinian population is placed within the European cluster. To reveal possible traces of selective pressure, the Population Branch Statistic (PBS) was calculated; both the rs143384 and 143383 SNPs have low PBS values, suggesting that there are no signals of selective pressure in those areas of the gene.

Atovaquone exposure and Pneumocystis jirovecii cytochrome b mutations: French data and review of the literature

Pneumocystis jirovecii is a transmissible fungus responsible for severe pneumonia (Pneumocystis pneumonia [PCP]) in immunocompromised patients. Missense mutations due to atovaquone selective pressure have been identified on cytochrome b (CYB) gene of P. jirovecii. It was recently shown that atovaquone prophylaxis can lead to the selection of specific P. jirovecii CYB mutants potentially resistant to atovaquone among organ transplant recipients. In this context, our objectives were to provide data on P. jirovecii CYB mutants and the putative selective pressure exerted by atovaquone on P. jirovecii organisms in France. A total of 123 patients (124 P. jirovecii specimens) from four metropolitan hospitals and two overseas hospitals were retrospectively enrolled. Fourteen patients had prior exposure to atovaquone, whereas 109 patients did not at the time of P. jirovecii detection. A 638 base-pair fragment of the CYB gene of P. jirovecii was amplified and sequenced. A total of 10 single nucleotide polymorphisms (SNPs) were identified. Both missense mutations C431T (Ala144Val) and C823T (Leu275Phe), located at the Qo active site of the enzyme, were significantly associated with prior atovaquone exposure, these mutations being conversely incidental in the absence of prior atovaquone exposure (P < 0.001). Considering that the aforementioned hospitals may be representative of the national territory, these findings suggest that the overall presence of P. jirovecii CYB mutants remains low in France.

ACE2 receptor polymorphism in humans and animals increases the risk of the emergence of SARS-CoV-2 variants during repeated intra- and inter-species host-switching of the virus

Like other coronaviruses, SARS-CoV-2 has ability to spread through human-to-human transmission and to circulate from humans to animals and from animals to humans. A high frequency of SARS-CoV-2 mutations has been observed in the viruses isolated from both humans and animals, suggesting a genetic fitness under positive selection in both ecological niches. The most documented positive selection force driving SARS-CoV-2 mutations is the host-specific immune response. However, after electrostatic interactions with lipid rafts, the first contact between the virus and host proteins is the viral spike-cellular receptor binding. Therefore, it is likely that the first level of selection pressure impacting viral fitness relates to the virus's affinity for its receptor, the angiotensin I converting enzyme 2 (ACE2). Although sufficiently conserved in a huge number of species to support binding of the viral spike with enough affinity to initiate fusion, ACE2 is highly polymorphic both among species and within a species. Here, we provide evidence suggesting that when the viral spike-ACE2 receptor interaction is not optimal, due to host-switching, mutations can be selected to improve the affinity of the spike for the ACE2 expressed by the new host. Notably, SARS-CoV-2 is mutation-prone in the spike receptor binding domain (RBD), allowing a better fit for ACE2 orthologs in animals. It is possibly that this may also be true for rare human alleles of ACE2 when the virus is spreading to billions of people. In this study, we present evidence that human subjects expressing the rare E₃₂₉G allele of ACE2 with higher allele frequencies in European populations exhibit a improved affinity for the SARS-CoV-2 spike N₅₀₁Y variant of the virus. This may suggest that this viral N₅₀₁Y variant emerged in the human population after SARS-CoV-2 had infected a human carrying the rare E₃₂₉G allele of ACE2. In addition, this viral evolution could impact viral replication as well as the ability of the adaptive humoral response to control infection with RBD-specific neutralizing antibodies. In a shifting landscape, this ACE2-driven genetic drift of SARS-CoV-2 which we have named the 'boomerang effect', could complicate the challenge of preventing COVID with a SARS-CoV-2 spike-derived vaccine.

Antibiotic, Heavy Metal, and Biocide Concentrations in a Wastewater Treatment Plant and Its Receiving Water Body Exceed PNEC Limits: Potential for Antimicrobial Resistance Selective Pressure

Although the rise in antimicrobial resistance has been attributed mainly to the extensive and indiscriminate use of antimicrobials such as antibiotics and biocides in humans, animals and on plants, studies investigating the impact of this use on water environments in Africa are minimal. This study quantified selected antibiotics, heavy metals, and biocides in an urban wastewater treatment plant (WWTP) and its receiving water body in Kwazulu-Natal, South Africa, in the context of the predicted no-effect concentrations (PNEC) for the selection of antimicrobial resistance (AMR). Water samples were collected from the WWTP effluent discharge point and upstream and downstream from this point. Heavy metals were identified and quantified using the United States Environmental Protection Agency (US EPA) method 200.7. Biocides and antibiotic residues were determined using validated ultra-high-performance liquid chromatography with tandem mass spectrometry-based methods. The overall highest mean antibiotic, metal and biocide concentrations were observed for sulfamethoxazole (286.180 µg/L), neodymium (Nd; 27.734 mg/L), and benzalkonium chloride (BAC 12) (7.805 µg/L), respectively. In decreasing order per sampling site, the pollutant concentrations were effluent > downstream > upstream. This implies that the WWTP significantly contributed to the observed pollution in the receiving water. Furthermore, most of the pollutants measured recorded values exceeding the recommended predicted no-effect concentration (PNEC) values, suggesting that the microbes in such water environments were at risk of developing resistance due to the selection pressure exerted by these antimicrobials. Further studies are required to establish such a relationship.

Molecular Phylogeny and Evolution of the Tuerkayana (Decapoda: Brachyura: Gecarcinidae) Genus Based on Whole Mitochondrial Genome Sequences

Tuerkayana is of particular interest because it has been separated, in recent years, from Cardisoma and Discoplax but studies of its taxonomic status, especially from a whole mitochondrial genome perspective, have been lacking. In this study, the mitogenomes of four species (Tuerkayana magnum, Tuerkayana rotundum, Tuerkayana hirtipes, and Tuerkayana celeste) of Tuerkayana are sequenced and contrasted with other species in Brachyura for the first time. The phylogenetic tree of Brachyura, which includes 206 crab species (189 species of Brachyuran and 17 Anomura species) with a complete mitogenome, was constructed to evaluate the phylogenetic position of Tuerkayana and Gecarcinidae within Brachyuran, and explore the monophyly of Gecarcinidae. Furthermore, two single gene trees based on cox1 and 16SrRNA separately within interspecies of Gecarcinidae were reconstructed, providing molecular evidence for Tuerkayana and further clarifying the division of genera in Gecarcinidae. Based on the mitogenome dataset of 206 crabs, the branch-site model was utilized to explore selective pressure in individual codons with CodeML. The strong selective pressure shown in nad6 indicates that it may have played a significant role in the evolution of Gecarcinidae.

The Microbiome of Things: Appliances, Machines, and Devices Hosting Artificial Niche-Adapted Microbial Communities

As it is the case with natural substrates, artificial surfaces of man-made devices are home to a myriad of microbial species. Artificial products are not necessarily characterized by human-associated microbiomes; instead, they can present original microbial populations shaped by specific environmental-often extreme-selection pressures. This review provides a detailed insight into the microbial ecology of a range of artificial devices, machines, and appliances, which we argue are specific microbial niches that do not necessarily fit in the "build environment" microbiome definition. Instead, we propose here the Microbiome of Things (MoT) concept analogous to the Internet of Things (IoT) because we believe it may be useful to shed light on human-made, but not necessarily human-related, unexplored microbial niches.

Mimicking the last step of gene elongation: hints from the bacterial hisF gene

Gene elongation consists in an in-tandem duplication of a gene and divergence and fusion of the two copies, resulting in a gene constituted by two divergent paralogous modules. Many present-day proteins show internal repeats of amino acid sequences, generated by gene elongation events; however, gene elongation is still a poorly studied evolutionary molecular mechanism. The most documented case is that of the histidine biosynthetic genes hisA and hisF, which derive from the gene elongation of an ancestral gene half the size of the extant ones. The aim of this work was to experimentally simulate the possible last step of the gene elongation event occurred during hisF gene evolution under selective pressure conditions. Azospirillum brasilense hisF gene, carrying a single nucleotide mutation that generates a stop codon between the two halves of the gene, was used to transform the histidine-auxotrophic Escherichia coli strain FB182 (hisF892). The transformed strain was subjected to selective pressure (i.e., low concentration/absence of histidine in the growth medium) and the obtained mutants were characterized. The restoration of prototrophy was strongly dependent on the time of incubation and on the strength of the selective pressure. The mutations involved the introduced stop codon with a single base substitution and none of the mutants restored the wild-type codon. Possible correlations between the different mutations and i) E. coli codon usage, ii) three-dimensional structures of the mutated HisF proteins, and iii) growth ability of the mutants were investigated. On the contrary, when the experiment was repeated by mutating a more conserved codon, only a synonymous substitution was obtained. Thus, experiments performed in this study allowed to mimic a possible gene elongation event occurred during the evolution of hisF gene, evidencing the ability of bacterial cells to modify their genome in short times under selective conditions.

Attenuation of Avian Flavivirus by Rewiring the Leucine and Serine Codons of Its E-NS1 Protein toward Stop Mutation To Redirect Virus Evolution

Recently, a new strategy for attenuating RNA viruses by redirecting their evolution in sequence space was confirmed for Enterovirus and Influenza viruses. Using avian flavivirus as a model, the 69 serine and 53 leucine codons on the E-NS1 genes were modified to change evolutionary direction of the viral sequence space. This means that all codons encoding serine or leucine residues were substituted with codons that are only one base different from the three stop codons, resulting in the initial position of the virus genome in sequence space being closer to the detrimental areas to achieve attenuation by reducing viral adaptability. The growth curve and plaque size of CQW1-one-to-stop (CQW1-OTS) were similar to those of CQW1-wild type (CQW1-WT) in vitro, but attenuated proliferation was detected when treated with a mutagenic reagent (ribavirin). However, comparably high CQW1-OTS and CQW1-WT lethality rates were detected in 9-day-old duck embryos and 5-day-old ducklings, suggesting that this strategy works but with limitations. With that in mind, homologous hosts in nonsensitive age (25-day-old ducks) and heterologous hosts (3-week-old Kunming mice) were employed to investigate if CQW1-OTS was attenuated under host selection pressure. Minimal attenuation of CQW1-OTS in elder ducks and apparent attenuation in mice were reported, providing reduced viral titers, mild clinical signs, and lower specific infectivity. Collectively, we experimentally demonstrate that the attenuation strategy of redirecting virus evolution in sequence space works for flavivirus. Redirection of the virus is attenuated only under some outside pressure, such as heterologous hosts or antiviral drugs treatment, limiting its usage in flaviviruses. IMPORTANCE Flaviviruses are medically important arboviruses that threaten public health, but no approved treatments are currently available. Vaccines prevent flavivirus infection. We employed duck Tembusu virus (TMUV), a mosquito-borne flavivirus, to evaluate virus redirection. TMUV is native to birds and could infect mice by intracerebral injection, making it an experimental animal model to study flavivirus characteristics in vivo. The 69 serine and 53 leucine codons on the E-NS1 proteins of CQW1 were synonymously substituted to change evolutionary direction of the virus in sequence space. In vitro mutagen reagent treatment suppressed CQW1-OTS viral multiplication, but in vivo attenuation depended on host selective pressure. CQW1-OTS viral attenuation was observed in older ducks but not sensitive ducklings; considerable attenuation was also observed in heterogenous host (mice), which provides more selective pressure on viruses. Collectively, these data indicated that there are very important preconditions for application of evaluating whether this strategy shows application prospects in novel flavivirus vaccine development.

Genomic characterization of colistin resistance in Klebsiella spp. under the pressure of colistin

Introduction. Carbapenem-resistant Klebsiella pneumoniae (CRKP) has become a serious threat to global public health. Colistin is regarded as the last-resort antibiotic for CRKP infections, but colistin resistance among CRKP is increasingly being reported, making clinical treatment for CRKP infections more difficult.Hypothesis/Gap Statement. The molecular mechanisms of colistin resistance in Klebsiella spp. under the pressure of colistin have not been fully investigated.Aim. We aimed to investigate the phenotypic and genetic variation in two colistin-susceptible Klebsiella spp. strains under selective pressure of colistin.Methodology. One hundred microlitres of overnight cultures of the CRKP clinical strain CRKP12-130 and of ATCC 700603 was spread on five Mueller-Hinton Agar (MHA) plates with colistin concentrations of 2, 4, 8, 16 and 32 µg ml^-1, and growth of colonies was observed for five consecutive days. Colonies collected from plates were passaged daily for 10 days on MHA plates without colistin and susceptibility testing of colistin was performed by broth microdilution. Thirty-four colistin-resistant strains randomly selected were submitted to whole genome sequencing (WGS). Transcriptional levels of genes involved in colistin resistance (mgrB, phoP, phoQ, pmrA, pmrB, pmrD, pmrE and pmrK) were measured by quantitative real-time PCR.Results. A total of 114 and 119 colistin-resistant colonies were obtained from CRKP12-130 and ATCC 700603 in this study, among which 16 and 18 colonies were submitted to WGS, respectively. Among these 34 sequenced isolates, mutation in phoQ (13/16, 81.25 %) was the main genetic factor mediating colistin resistance in strains from CRKP12-130, while for strains from ATCC 700603, mutation associated with mgrB (8/18, 44.44 %) was found to be the commonest. Mutation of mgrB led to a significant increase in the MIC for colistin (from 64 to >128 µg ml^-1), and a novel mutation C28R in mgrB was first reported in this study.Conclusion. Colistin-resistant Klebsiella spp. could be easily selected under pressure of different concentrations of colistin. Mutations of mgrB, phoP, phoQ and pmrB genes were the main mechanisms leading to chromosomally mediated colistin resistance in Klebsiella spp.

Complete Mitogenome Analysis of Five Leafhopper Species of Idiocerini (Hemiptera: Cicadellidae)

Insect mitochondrial genomes (mitogenomes) are of great interest in exploring molecular evolution, phylogenetics, and biogeography. So far, only 12 mitogenomes of the leafhopper tribe Idiocerini have been released in GenBank, although the tribe comprises 488 known species including some agricultural, forestry, and horticultural pests. In order to compare and analyze the mitochondrial genome structure of Idiocerini and even the selective pressure of 13 protein-coding genes (PCGs) of the family Cicadellidae, the complete mitogenomes of five species including Nabicerus dentimus, Sahlbergotettix salicicola, Podulmorinus opacus, Podulmorinus consimilis, and a new species of a new genus were determined by next-generation sequencing. The size of the newly determined mitogenomes ranged from 14,733 bp to 15,044 bp, comprising the standard set of 13 PCGs, 22 transfer RNA genes, two ribosomal RNA genes, and a long non-coding control region (CR). The extent of purifying selection presented different pictures in the tribe and the family. The less pronounced genes (0.5 < dN/dS < 1) were nad5 and nad4l in Idiocerin, whereas in the family Cicadellidae including the sequences of Idiocerin, nad1-nad6 and cox1 genes were less pronounced. The codon encoding leucine was the most common in all species, and the codon encoding serine 1 was the most common in all species except for P. opacus. Interestingly, in P. opacus, another of the most common codons is that encoding serine 2. Among the 17 examined species of the Idiocerini, 14 species contained the tandem repeats, and 11 species of them contained the motif "TTATA". These findings will promote research on the structure and evolution of the mitochondrial genome and highlight the need for more mitogenomes in Cicadellidae.