De novo genes in Arachis hypogaea cv. Tifrunner: systematic identification, molecular evolution, and potential contributions to cultivated peanut

De novo genes are derived from non-coding sequences, and they can play essential roles in organisms. Cultivated peanut (Arachis hypogaea) is a major oil and protein crop derived from a cross between Arachis duranensis and Arachis ipaensis. However, few de novo genes have been documented in Arachis. Here, we identified 381 de novo genes in A. hypogaea cv. Tifrunner based on comparison with five closely related Arachis species. There are distinct differences in gene expression patterns and gene structures between conserved and de novo genes. The identified de novo genes originated from ancestral sequence regions associated with metabolic and biosynthetic processes, and they were subsequently integrated into existing regulatory networks. De novo paralogs and homoeologs were identified in A. hypogaea cv. Tifrunner. De novo paralogs and homoeologs with conserved expression have mismatching cis-acting elements under normal growth conditions. De novo genes potentially have pluripotent functions in responses to biotic stresses as well as in growth and development based on quantitative trait locus data. This work provides a foundation for future research examining gene birth processes and gene function in Arachis and related taxa.

Duplication of NRAMP3 gene in poplars generated two homologous transporters with distinct functions

Transition metals are essential for a wealth of metabolic reactions, but their concentrations need to be tightly controlled across cells and cell compartments, as metal excess or imbalance has deleterious effects. Metal homeostasis is achieved by a combination of metal transport across membranes and metal binding to a variety of molecules. Gene duplication is a key process in evolution, as emergence of advantageous mutations on one of the copies can confer a new function. Here, we report that the poplar genome contains two paralogues encoding NRAMP3 metal transporters localized in tandem. All Populus species analyzed had two copies of NRAMP3, whereas only one could be identified in Salix species indicating that duplication occurred when the two genera separated. Both copies are under purifying selection and encode functional transporters, as shown by expression in the yeast heterologous expression system. However, genetic complementation revealed that only one of the paralogues has retained the original function in release of metals stored in the vacuole previously characterized in A. thaliana. Confocal imaging showed that the other copy has acquired a distinct localization to the Trans Golgi Network (TGN). Expression in poplar suggested that the copy of NRAMP3 localized on the TGN has a novel function in the control of cell-to-cell transport of manganese. This work provides a clear case of neo-functionalization through change in the subcellular localization of a metal transporter as well as evidence for the involvement of the secretory pathway in cell-to-cell transport of manganese.

Size-dependent flight capacity and propensity in a range-expanding invasive insect

For capital-breeding insects, all resources available for adult metabolic needs are accumulated during larval feeding. Therefore, body size at adult eclosion represents the total energetic capacity of the individual. For female capital breeders, body size is strongly correlated with lifetime fecundity, while in males, body size, which correlates with fitness, is less understood. In capital-breeding species with wingless, flightless, or dispersal-limited females, flight potential for male Lepidoptera has important implications for mate-finding and may be correlated with body size. At low population densities, failure to mate has been identified as an important Allee effect and can drive the success or failure of invasive species at range edges and in species of conservation concern. Th capital-breeding European subspecies of Lymantria dispar (L.), was introduced to North America in 1869 and now ranges across much of eastern North America. In L. dispar, females are flightless and mate-finding is entirely performed by males. We quantified male L. dispar flight capacity and propensity relative to morphological and physiological characteristics using fixed-arm flight mills. A range of male body sizes was produced by varying the protein content of standard artificial diets while holding other dietary components constant. Wing length, a proxy for body size, relative thorax mass, and forewing aspect were all important predictors of total flight distance and maximum speed. These results have important implications for mate-finding and invasion dynamics in L. dispar and may apply broadly to other capital-breeding insects.

Phylogenetic and Spatiotemporal Analyses of the Complete Genome Sequences of Avian Coronavirus Infectious Bronchitis Virus in China During 1985-2020: Revealing Coexistence of Multiple Transmission Chains and the Origin of LX4-Type Virus

Infectious bronchitis (IB) virus (IBV) causes considerable economic losses to poultry production. The data on transmission dynamics of IBV in China are limited. The complete genome sequences of 212 IBV isolates in China during 1985–2020 were analyzed as well as the characteristics of the phylogenetic tree, recombination events, dN/dS ratios, temporal dynamics, and phylogeographic relationships. The LX4 type (GI-19) was found to have the highest dN/dS ratios and has been the most dominant genotype since 1999, and the Taiwan-I type (GI-7) and New type (GVI-1) showed an increasing trend. A total of 59 recombinants were identified, multiple recombination events between the field and vaccine strains were found in 24 isolates, and the 4/91-type (GI-13) isolates were found to be more prone to being involved in the recombination. Bayesian phylogeographic analyses indicated that the Chinese IBVs originated from Liaoning province in the early 1900s. The LX4-type viruses were traced back to Liaoning province in the late 1950s and had multiple transmission routes in China and two major transmission routes in the world. Viral phylogeography identified three spread regions for IBVs (including LX4 type) in China: Northeastern China (Heilongjiang, Liaoning, and Jilin), north and central China (Beijing, Hebei, Shanxi, Shandong, and Jiangsu), and Southern China (Guangxi and Guangdong). Shandong has been the epidemiological center of IBVs (including LX4 type) in China. Overall, our study highlighted the reasons why the LX4-type viruses had become the dominant genotype and its origin and transmission routes, providing more targeted strategies for the prevention and control of IB in China.

A Systematic Investigation Unveils High Coinfection Status of Porcine Parvovirus Types 1 through 7 in China from 2016 to 2020

Porcine parvovirus genotype 1 (PPV1) causes reproductive disorder in swine and is prevalent in China. Recently, six new genotypes of PPVs (PPV2 through PPV7) have also been detected in Chinese swine herds. However, the coinfection status of all these seven genotypes of PPVs (PPV1-7) in China was not clarified yet. In this study, we developed a panel of PPV1–7 PCR assays with satisfied specificity, sensitivity and reproducibility and then applied to the detection of PPV1–7 in 435 clinical samples collected from eight provinces of China in 2016–2020. A total of 55.40% samples (241 out of 435) were PPV positive, while PPV2 and PPV3 (both 22.53%) belonging to the genus of Tetraparvovirus were the most prevalent genotypes. Noticeably, PPV1–7 strains were more prevalent in nursery and finishing pigs than in suckling pigs. In addition, coinfection could be detected in all eight provinces and 27.36% (119/435) samples were coinfected with two to five genotypes of PPVs. Meanwhile, the coinfection of PPVs with PCV2 was 22.30% (97/435). Twenty complete genomes of representative PPV1–7 were determined, and phylogenetic analysis confirmed the genotyping results by sequence comparisons and PCR assays. Remarkably, the PPV7 HBTZ20180519-152 strain from domestic pig was recombined from parental JX15-like and JX38-like isolates from wild boars. Selective pressure analysis based on VP2 sequences of PPV1–7 showed that they were predominantly under negative selection, while few positive selection sites could be detected in VP2 of PPV7. Overall, this systematic investigation unveils high prevalence and coinfection of PPV1–7 in China from 2016 to 2020.

IMPORTANCE Porcine parvoviruses (PPVs) are prevalent in China associating with reproductive failure in swine. The coinfection of seven genotypes of PPVs (PPV1-7) might have synergistic effects on PPV1 associated SMEDI syndrome. However, the coinfection status of PPV1–7 in China is not clear yet. This study showed that PPV1–7 strains are highly prevalent (55.40%) in China and mainly in nursery and finishing pigs in recent years. In addition, the coinfections of different genotypes of PPVs (27.36%) and PPVs with PCV2 (22.30%) are common. Geographic analysis indicated that different genotypes of PPVs are widely cocirculating in China. Intriguingly, a PPV7 strain from the domestic pig was detected as a recombinant from two wild boar isolates. Selective pressure analyses showed that PPV1–7 are mainly under purifying selection. Our findings provide the first systematic investigation on the prevalence, coinfection, and evolution of PPV1 through PPV7 in Chinese swineherds from 2016 to 2020.

Quasispecies of SARS-CoV-2 revealed by single nucleotide polymorphisms (SNPs) analysis

New SARS-CoV-2 mutants have been continuously indentified with enhanced transmission ever since its outbreak in early 2020. As an RNA virus, SARS-CoV-2 has a high mutation rate due to the low fidelity of RNA polymerase. To study the single nucleotide polymorphisms (SNPs) dynamics of SARS-CoV-2, 158 SNPs with high confidence were identified by deep meta-transcriptomic sequencing, and the most common SNP type was C > T. Analyses of intra-host population diversity revealed that intra-host quasispecies' composition varies with time during the early onset of symptoms, which implicates viral evolution during infection. Network analysis of co-occurring SNPs revealed the most abundant non-synonymous SNP 22,638 in the S glycoprotein RBD region and 28,144 in the ORF8 region. Furthermore, SARS-CoV-2 variations differ in an individual's respiratory tissue (nose, throat, BALF, or sputum), suggesting independent compartmentalization of SARS-CoV-2 populations in patients. The positive selection analysis of the SARS-CoV-2 genome uncovered the positive selected amino acid G251V on ORF3a. Alternative allele frequency spectrum (AAFS) of all variants revealed that ORF8 could bear alternate alleles with high frequency. Overall, the results show the quasispecies' profile of SARS-CoV-2 in the respiratory tract in the first two months after the outbreak.

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.

Solids Residence Time Impacts Carbon Dynamics and Bioenergy Feedstock Potential in Phototrophic Wastewater Treatment Systems

The use of wastewater-grown microalgae has the potential to reduce the cost of algae-derived biofuels while simultaneously advancing nutrient recovery at water resource recovery facilities (WRRFs). However, a significant barrier has been the low yield and high protein content of phototrophic biomass. Here, we examine the use of solids residence time (SRT) as a selective pressure in driving biochemical composition, yield, biofuel production, and WRRF nutrient management cost. We cultivated mixed phototrophic communities in controlled, laboratory-scale photobioreactors on the local WRRF secondary effluent to link SRT with biochemical composition and techno-economic analysis to yield insights into biomass composition and downstream processing effects on minimum fuel selling price. SRT significantly impacted biochemical composition, with total and dynamic carbohydrates the highest at low SRT (total carbohydrates being 0.60 and 0.32 mg-carbohydrate·mg-protein^-1 at SRT 5 and 15 days, respectively). However, there were distinct differences between extant, steady-state performance and intrinsic potential, and longer SRT communities were able to accumulate significant fractions (51% on an ash-free dry weight basis, AFDW %) of carbohydrate reserves under nutrient starvation. Overall, hydrothermal liquefaction (HTL) was found to be more suitable than lipid extraction for hydrotreating (LEH) and combined algal processing (CAP) for conversion of biomass to fuels, but LEH and CAP became more competitive when intrinsic carbon storage potential was realized. The results suggest that the use of algae for nutrient recovery could reduce the nutrient management cost at WRRFs through revenue from algal biofuels, with HTL resulting in a net revenue.

Comparative phylogenetic analysis of SARS-CoV-2 spike protein-possibility effect on virus spillover

Coronavirus disease 2019 has developed into a dramatic pandemic with tremendous global impact. The receptor-binding motif (RBM) region of the causative virus, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), binds to host angiotensin-converting enzyme 2 (ACE2) receptors for infection. As ACE2 receptors are highly conserved within vertebrate species, SARS-CoV-2 can infect significant animal species as well as human populations. An analysis of SARS-CoV-2 genotypes isolated from human and significant animal species was conducted to compare and identify mutation and adaptation patterns across different animal species. The phylogenetic data revealed seven distinct phylogenetic clades with no significant relationship between the clades and geographical locations. A high rate of variation within SARS-CoV-2 mink isolates implies that mink populations were infected before human populations. Positions of most single-nucleotide polymorphisms (SNPs) within the spike (S) protein of SARS-CoV-2 genotypes from the different hosts are mostly accumulated in the RBM region and highlight the pronounced accumulation of variants with mutations in the RBM region in comparison with other variants. These SNPs play a crucial role in viral transmission and pathogenicity and are keys in identifying other animal species as potential intermediate hosts of SARS-CoV-2. The possible roles in the emergence of new viral strains and the possible implications of these changes, in compromising vaccine effectiveness, deserve urgent considerations.

An Update on Trichoderma Mitogenomes: Complete De Novo Mitochondrial Genome of the Fungal Biocontrol Agent Trichoderma harzianum (Hypocreales, Sordariomycetes), an Ex-Neotype Strain CBS 226.95, and Tracing the Evolutionary Divergences of Mitogenomes in Trichoderma

Members of the genus Trichoderma (Hypocreales), widely used as biofungicides, biofertilizers, and as model fungi for the industrial production of CAZymes, have actively been studied for the applications of their biological functions. Recently, the study of the nuclear genomes of Trichoderma has expanded in the directions of adaptation and evolution to gain a better understanding of their ecological traits. However, Trichoderma's mitochondria have received much less attention despite mitochondria being the most necessary element for sustaining cell life. In this study, a mitogenome of the fungus Trichoderma harzianum CBS 226.95 was assembled de novo. A 27,632 bp circular DNA molecule was revealed with specific features, such as the intronless of all core PCGs, one homing endonuclease, and a putative overlapping tRNA, on a closer phylogenetic relationship with T. reesei among hypocrealean fungi. Interestingly, the mitogenome of T. harzianum CBS 226.95 was predicted to have evolved earlier than those of other Trichoderma species and also assumed with a selection pressure in the cox3. Considering the bioavailability, both for the ex-neotype strain of the T. harzianum species complex and the most globally representative commercial fungal biocontrol agent, our results on the T. harzianum CBS 226.95 mitogenome provide crucial information which will be helpful criteria in future studies on Trichoderma.

A Codon-Pair Bias Associated With Network Interactions in Influenza A, B, and C Genomes

A new codon-pair bias present in the genomes of different types of influenza virus is described. Codons with fewer network interactions are more frequency paired together than other codon-pairs in influenza A, B, and C genomes. A shared feature among three different influenza types suggests an evolutionary bias. Codon-pair preference can affect both speed of protein translation and RNA structure. This newly identified bias may provide insight into drivers of virus evolution.

Ongoing evolution of Chlamydia trachomatis lymphogranuloma venereum: exploring the genomic diversity of circulating strains

Lymphogranuloma venereum (LGV), the invasive infection of the sexually transmissible infection (STI) Chlamydia trachomatis, is caused by strains from the LGV biovar, most commonly represented by ompA-genotypes L2b and L2. We investigated the diversity in LGV samples across an international collection over seven years using typing and genome sequencing. LGV-positive samples (n=321) from eight countries collected between 2011 and 2017 (Spain n=97, Netherlands n=67, Switzerland n=64, Australia n=53, Sweden n=37, Hungary n=31, Czechia n=30, Slovenia n=10) were genotyped for pmpH and ompA variants. All were found to contain the 9 bp insertion in the pmpH gene, previously associated with ompA-genotype L2b. However, analysis of the ompA gene shows ompA-genotype L2b (n=83), ompA-genotype L2 (n=180) and several variants of these (n=52; 12 variant types), as well as other/mixed ompA-genotypes (n=6). To elucidate the genomic diversity, whole genome sequencing (WGS) was performed from selected samples using SureSelect target enrichment, resulting in 42 genomes, covering a diversity of ompA-genotypes and representing most of the countries sampled. A phylogeny of these data clearly shows that these ompA-genotypes derive from an ompA-genotype L2b ancestor, carrying up to eight SNPs per isolate. SNPs within ompA are overrepresented among genomic changes in these samples, each of which results in an amino acid change in the variable domains of OmpA (major outer membrane protein, MOMP). A reversion to ompA-genotype L2 with the L2b genomic backbone is commonly seen. The wide diversity of ompA-genotypes found in these recent LGV samples indicates that this gene is under immunological selection. Our results suggest that the ompA-genotype L2b genomic backbone is the dominant strain circulating and evolving particularly in men who have sex with men (MSM) populations.

Herbicides Tolerance in a Pseudomonas Strain Is Associated With Metabolic Plasticity of Antioxidative Enzymes Regardless of Selection

Agriculture uses many food production chains, and herbicides participate in this process by eliminating weeds through different biochemical strategies. However, herbicides can affect non-target organisms such as bacteria, which can suffer damage if there is no efficient control of reactive oxygen species. It is not clear, according to the literature, whether the efficiency of this control needs to be selected by the presence of xenobiotics. Thus, the Pseudomonas sp. CMA 6.9 strain, collected from biofilms in an herbicide packaging washing tank, was selected for its tolerance to pesticides and analyzed for activities of different antioxidative enzymes against the herbicides Boral^®, absent at the isolation site, and Heat^®, present at the site; both herbicides have the same mode of action, the inhibition of the enzyme protoporphyrinogen oxidase. The strain showed tolerance to both herbicides in doses up to 45 times than those applied in agriculture. The toxicity of these herbicides, which is greater for Boral^®, was assessed by means of oxidative stress indicators, growth kinetics, viability, and amounts of peroxide and malondialdehyde. However, the studied strain showed two characteristic antioxidant response systems for each herbicide: glutathione-s-transferase acting to control malondialdehyde in treatments with Boral^®; and catalase, ascorbate peroxidase, and guaiacol peroxidase in the control of peroxide induced by Heat^®. It is possible that this modulation of the activity of different enzymes independent of previous selection characterizes a system of metabolic plasticity that may be more general in the adaptation of microorganisms in soil and water environments subjected to chemical contaminants. This is relevant to the impact of pesticides on the diversity and abundance of microbial species as well as a promising line of metabolic studies in microbial consortia for use in bioremediation.

Population genetic structure of sharpbelly Hemiculter leucisculus (Basilesky, 1855) and morphological diversification along climate gradients in China

Sharpbelly Hemiculter leucisculus (Basilewski, 1855) is a small, widespread, and native cyprinid fish with prominent habitat suitability and high invasive potential and is becoming the dominant species in freshwater ecosystems under intensified environmental disturbances. But how H. leucisculus acclimates to extremely heterogeneous environments remains unclear. In current study, the genetic structure of H. leucisculus was analyzed using Bayesian phylogenetic inference, haplotype network, and STRUCTURE base on cytb gene across 18 populations spanning 20 degrees of latitude and 18 degrees of longitude in China. The morphological diversification of body size and shape for H. leucisculus along the climate gradient was studied. The results showed that the 18 H. leucisculus populations were divided into 3 clusters: one cluster mainly from Huanghe River Basin, another cluster mainly from Yangzi River Basin, and H cluster containing Hainan and Beihai populations. The fish from southern populations were deeper bodied while individuals from northern populations were more slender. Inland individuals were more streamlined while coastal individuals were of deeper body. The partial Mantel test predicts that the potential mechanism underlining the intraspecies morphological diversification along climate gradients is primarily the divergent selection pressures among different environments, while genetic variation had less contribution to morphological differentiation. The formation of the Nanling Mountain Range could drive genetic differentiation between Beihai population and those from Yangzi River Basin. The present results highlight strong selective pressures of climate on widespread species and enrich morphological differentiation basis of acclimation for species with high habitat suitability and invasive potential.

Epidemiology, Genetic Characterization, and Evolution of Hunnivirus Carried by Rattus norvegicus and Rattus tanezumi: The First Epidemiological Evidence from Southern China

Hunnivirus is a novel member of the family Picornaviridae. A single species, Hunnivirus A, is currently described. However, there is limited information on the identification of Hunnivirus to date, and thereby the circulation of Hunnivirus is not fully understood. Thus, the objective of this study was to investigate the prevalence, genomic characteristics, and evolution of rat hunnivirus in southern China. A total of 404 fecal samples were subjected to detection of Hunnivirus from urban rats (Rattus norvegicus and Rattus tanezumi) using PCR assay based on specific primers targeted to partial 3D regions, with the prevalence of 17.8% in Rattus norvegicus and 15.6% in Rattus tanezumi. An almost full-length rat hunnivirus sequence (RatHuV/YY12/CHN) and the genome structure were acquired in the present study. Phylogenetic analysis of the P1 coding regions suggested the RatHuV/YY12/CHN sequence was found to be within the genotype of Hunnivirus A4. The negative selection was further identified based on analysis of non-synonymous to synonymous substitution rates. The present findings suggest that hunniviruses are common in urban rats. Further research is needed for increased surveillance and awareness of potential risks to human health.

Positive selection drives the evolution of a primate bitter taste receptor gene

Bitter taste perception is important in food choice of animals, and it is mediated by bitter taste receptor (T2R) containing three regions: extracellular regions (ECs), transmembrane regions (TMs), and intracellular regions (ICs). It is hypothesized that ECs, TMs, and ICs are under different selective pressures, with ECs being unstable while TMs and ICs being constrained. To test this hypothesis, we examined the selective pressures on one of the bitter taste receptor genes-T2R1 and its different areas from 35 primates and found that T2R1 was under neutral evolution but with some positively selected sites in it. Further analyses suggested that TMs, ICs, and the concatenated transmembrane region TM1237 were under purifying selection; in contrast, extracellular regions, the first and second extracellular loop (EL1, EL2), were subject to positive selection with several positively selected sites in them. Therefore, this study supported the above-mentioned hypothesis. The reason why EL1 and EL2 of T2R1 have positively selected sites is probably due to their participation in forming the cap-like structure involved in ligand binding. Positive selection acts as a driving force of the T2R1 functional differentiation and confers the ability to discern various bitter substances for primates, which could help them to adapt to the changing environment during the evolutionary course.

Evolutionary Analysis of Pre-S/S Mutations in HBeAg-Negative Chronic Hepatitis B With HBsAg < 100 IU/ml

Background: Hepatitis B surface antigen (HBsAg) and viral load are important clinical indicators for antiviral therapy. Few studies have evaluated viral sequence biomarkers predicting the risk of hepatocellular carcinoma (HCC) in the stage, which show a low serological response (HBsAg < 100 IU/ml) and high viral levels (HBV DNA > 2,000 IU/ml). This study aims to determine the trend of the biological prevalence within the pre-S/S regions of special model of inactive CHB infection. Methods: We used Sanger sequencing, quantitative HBV serology (HBeAg and HBsAg), and liver function index to identify whether HBV genome sequences are associated with long-term risk of further HCC progression in special inactive CHB infection. Results: HBV sequencing analysis of 28 CHB patients with special infectious pattern showed higher genetic diversity among four opening reading frames (ORFs) (p < 0.001). However, dN/dS ratios of HBsAg and pre-C/C regions in the experimental group showed no significantly different from those in the HCC group (p = 0.06), while significantly lower in polymerase and HBxAg regions of the experimental group (p < 0.001). In addition, seven positively selected sites were identified in pre-S1, five in pre-S2, and four in S, in which five sites (128H/135Q/135R/139L/141P) were among "α" determinant. Conclusions: These mutations in the pre-S/S region might be associated with the HCC phenotype of low HBsAg expression, with the P region possibly impacting high viral loads. Increased viral diversity across the HBV genome is also associated with low levels of HBsAg. The cumulative evolutionary changes in the HBV pre-S/S regions shows that facilitate immune evasion should be monitored individually. Due to the similarity of evolutionary characteristics in HCC, low serological responses and high viremia may be associated with the risk of further disease progression.

Breeding selection imposed a differential selective pressure on the wheat root-associated microbiome

Plants-microbiome associations are the result of millions of years of co-evolution. Due to breeding-accelerated plant evolution in non-native and highly managed soil, plant-microbe links could have been lost. We hypothesized that post-domestication breeding of wheat changed the root-associated microbiome. To test this, we analyzed root-associated fungal and bacterial communities shortly after emergence of seedlings representing a transect of wheat evolution including modern wheat, landraces and ancestors. Numbers of observed microbial taxa were highest in landraces bred in low-input agricultural systems, and lowest in ancestors that had evolved in native soils. The microbial communities of modern cultivars were different from those of landraces and ancestors. Old wheat accessions enriched Acidobacteria and Actinobacteria, while modern cultivars enriched OTUs from Candidatus Saccharibacteria, Verrucomicrobia and Firmicutes. The fungal pathogens Fusarium, Neoascochyta and Microdochium enriched in modern cultivars. Both bacterial and fungal communities followed a neutral assembly model when bulk soil was considered as the source community, but accessions of the ancient Triticum turgidum and T. monococcum created a more isolated environment in their roots. In conclusion, wheat root-associated microbiomes have dramatically changed through a transect of breeding history.

A Transcriptomic Approach to the Recruitment of Venom Proteins in a Marine Annelid

The growing number of known venomous marine invertebrates indicates that chemical warfare plays an important role in adapting to diversified ecological niches, even though it remains unclear how toxins fit into the evolutionary history of these animals. Our case study, the Polychaeta Eulalia sp., is an intertidal predator that secretes toxins. Whole-transcriptome sequencing revealed proteinaceous toxins secreted by cells in the proboscis and delivered by mucus. Toxins and accompanying enzymes promote permeabilization, coagulation impairment and the blocking of the neuromuscular activity of prey upon which the worm feeds by sucking pieces of live flesh. The main neurotoxins ("phyllotoxins") were found to be cysteine-rich proteins, a class of substances ubiquitous among venomous animals. Some toxins were phylogenetically related to Polychaeta, Mollusca or more ancient groups, such as Cnidaria. Some toxins may have evolved from non-toxin homologs that were recruited without the reduction in molecular mass and increased specificity of other invertebrate toxins. By analyzing the phylogeny of toxin mixtures, we show that Polychaeta is uniquely positioned in the evolution of animal venoms. Indeed, the phylogenetic models of mixed or individual toxins do not follow the expected eumetazoan tree-of-life and highlight that the recruitment of gene products for a role in venom systems is complex.

Inference of Nipah virus evolution, 1999-2015

Despite near-annual human outbreaks of Nipah virus (NiV) disease in Bangladesh, typically due to individual spillover events from the local bat population, only twenty whole-genome NiV sequences exist from humans and ten from bats. NiV whole-genome sequences from annual outbreaks have been challenging to generate, primarily due to the low viral load in human throat swab and serum specimens. Here, we used targeted enrichment with custom NiV-specific probes and generated thirty-five additional unique full-length genomic sequences directly from human specimens and viral isolates. We inferred the temporal and geographic evolutionary history of NiV in Bangladesh and expanded a tool to visualize NiV spatio-temporal spread from a Bayesian continuous diffusion analysis. We observed that strains from Bangladesh segregated into two distinct clades that have intermingled geographically in Bangladesh over time and space. As these clades expanded geographically and temporally, we did not observe evidence for significant branch and site-specific selection, except for a single site in the Henipavirus L polymerase. However, the Bangladesh 1 and 2 clades are differentiated by mutations initially occurring in the polymerase, with additional mutations accumulating in the N, G, F, P, and L genes on external branches. Modeling the historic geographical and temporal spread demonstrates that while widespread, NiV does not exhibit significant genetic variation in Bangladesh. Thus, future public health measures should address whether NiV within in the bat population also exhibits comparable genetic variation, if zoonotic transmission results in a genetic bottleneck and if surveillance techniques are detecting only a subset of NiV.

Insights Into the Molecular Evolution of AT-Hook Motif Nuclear Localization Genes in Brassica napus

AT-hook motif nuclear localization (AHL) proteins belong to a family of transcription factors, and play important roles in plant growth and development and response to various stresses through protein-DNA and protein-protein interactions. To better understand the Brassica napus AHL gene family, AHL genes in B. napus and related species were analyzed. Using Arabidopsis as a reference, 122 AHL gene family members were first identified in B. napus. According to the phylogenetic tree and gene organization, the BnaAHLs were classified into two clades (Clade-A and Clade-B) and three types (Type-I, Type-II, and Type-III). Gene organization and motif distribution analysis suggested that the AHL gene family is relatively conserved during evolution. These BnaAHLs are unevenly distributed on 38 chromosomes and expanded by whole-genome duplication (WGD) or segmental duplication. And large-scale loss events have also occurred in evolution. All types of BnaAHLs are subject to purification or neutral selection, while some positive selection sites are also identified in Type-II and Type-III groups. At the same time, the purification effect of Type-I members are stronger than that of the others. In addition, RNA-seq data and cis-acting element analysis also suggested that the BnaAHLs play important roles in B. napus growth and development, as well as in response to some abiotic and biotic stresses. Protein-protein interaction analysis identified some important BnaAHL-binding proteins, which also play key roles in plant growth and development. This study is helpful to fully understand the origin and evolution of the AHL gene in B. napus, and lays the foundation for their functional studies.

Antimicrobial Resistance as a Hidden Menace Lurking Behind the COVID-19 Outbreak: The Global Impacts of Too Much Hygiene on AMR

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a new coronavirus that was recently discovered in 2019. While the world is working hard to overcome and control the coronavirus disease 2019 (COVID-19) pandemic, it is also crucial to be prepared for the great impacts of this outbreak on the development of antimicrobial resistance (AMR). It is predicted that inappropriate and too much use of antibiotics, biocides, and disinfectants during this pandemic may raise disastrous effects on antibiotic stewardship programs and AMR control all around the world. Furthermore, the use of certain antibiotics alone or in combination with antiviral agents or other medications for the treatment of secondary bacterial infections among COVID-19 patients may be regarded as a major factor that negatively affects host immune response by disrupting mitochondrial function and activity. Herein, we suggest that the current management strategies to control AMR and prioritize antibiotic stewardship schemes should be extremely highlighted in relation to the COVID-19 outbreak. The rising concerns about excessive use of antimicrobials and biocides and taking too much hygiene also need to be addressed during this pandemic due to their impacts on AMR, public health, and the environment.

Efflux Pump-Driven Antibiotic and Biocide Cross-Resistance in Pseudomonas aeruginosa Isolated from Different Ecological Niches: A Case Study in the Development of Multidrug Resistance in Environmental Hotspots

Pseudomonas aeruginosa is an opportunistic pathogen displaying high intrinsic antimicrobial resistance and the ability to thrive in different ecological environments. In this study, the ability of P. aeruginosa to develop simultaneous resistance to multiple antibiotics and disinfectants in different natural niches were investigated using strains collected from clinical samples, veterinary samples, and wastewater. The correlation between biocide and antimicrobial resistance was determined by employing principal component analysis. Molecular mechanisms linking biocide and antimicrobial resistance were interrogated by determining gene expression using RT-qPCR and identifying a potential genetic determinant for co- and cross-resistance using whole-genome sequencing. A subpopulation of P. aeruginosa isolates belonging to three sequence types was resistant against the common preservative benzalkonium chloride and showed cross-resistance to fluoroquinolones, cephalosporins, aminoglycosides, and multidrug resistance. Of these, the epidemiological high-risk ST235 clone was the most abundant. The overexpression of the MexAB-OprM drug efflux pump resulting from amino acid mutations in regulators MexR, NalC, or NalD was the major contributing factor for cross-resistance that could be reversed by an efflux pump inhibitor. This is the first comparison of antibiotic-biocide cross-resistance in samples isolated from different ecological niches and serves as a confirmation of laboratory-based studies on biocide adapted isolates. The isolates from wastewater had a higher incidence of multidrug resistance and biocide-antibiotic cross-resistance than those from clinical and veterinary settings.