Population transcriptomics uncover the relative roles of positive selection and differential expression in Batrachium bungei adaptation to the Qinghai-Tibetan plateau.

KEY MESSAGE: Positive selection genes are related to metabolism, while differentially expressed genes are related to photosynthesis, suggesting that genetic adaptation and expression regulation may play independent roles in different gene classes. Genome-wide investigation of the molecular mechanisms for high-altitude adaptation is an intriguing topic in evolutionary biology. The Qinghai-Tibet Plateau (QTP) with its extremely variable environments is an ideal site for studying high-altitude adaptation. Here, we used transcriptome data of 100 individuals from 20 populations collected from various altitudes on the QTP to investigate the adaptive mechanisms of the aquatic plant Batrachium bungei at both the genetic and transcriptional level. To explore genes and biological pathways that may contribute to QTP adaptation, we employed a two-step approach, in which we identified positively selected genes and differentially expressed genes using the landscape genomic and differential expression approaches. The positive selection analysis showed that genes involved in metabolic regulation played a crucial role in B. bungei adaptation to the extreme environments of the QTP, especially intense ultraviolet radiation. Altitude-based differential expression analysis suggested that B. bungei could increase the rate of energy dissipation or reduce the efficiency of light energy absorption by down regulating the expression of photosynthesis-related genes to adapt to the strong ultraviolet radiation. Weighted gene co-expression network analysis identified ribosomal genes as hubs of altitude adaptation in B. bungei. Only a small part of genes (about 10%) overlapped between positively selected genes and differentially expressed genes in B. bungei, suggesting that genetic adaptation and gene expression regulation might play relatively independent roles in different categories of functional genes. Taken together, this study enriches our understanding of the high-altitude adaptation mechanism of B. bungei on the QTP.

Comparative Analysis of the Chloroplast Genomes of the Chinese Endemic Genus Urophysa and Their Contribution to Chloroplast Phylogeny and Adaptive Evolution.

Urophysa is a Chinese endemic genus comprising two species, Urophysa rockii and Urophysa henryi. In this study, we sequenced the complete chloroplast (cp) genomes of these two species and of their relative Semiquilegia adoxoides. Illumina sequencing technology was used to compare sequences, elucidate the intra- and interspecies variations, and infer the phylogeny relationship with other Ranunculaceae family species. A typical quadripartite structure was detected, with a genome size from 158,473 to 158,512 bp, consisting of a pair of inverted repeats separated by a small single-copy region and a large single-copy region. We analyzed the nucleotide diversity and repeated sequences components and conducted a positive selection analysis by the codon-based substitution on single-copy coding sequence (CDS). Seven regions were found to possess relatively high nucleotide diversity, and numerous variable repeats and simple sequence repeats (SSR) markers were detected. Six single-copy genes (atpA, rpl20, psaA, atpB, ndhI, and rbcL) resulted to have high posterior probabilities of codon sites in the positive selection analysis, which means that the six genes may be under a great selection pressure. The visualization results of the six genes showed that the amino acid properties across each column of all species are variable in different genera. All these regions with high nucleotide diversity, abundant repeats, and under positive selection will provide potential plastid markers for further taxonomic, phylogenetic, and population genetics studies in Urophysa and its relatives. Phylogenetic analyses based on the 79 single-copy genes, the whole complete genome sequences, and all CDS sequences showed same topologies with high support, and U. rockii was closely clustered with U. henryi within the Urophysa genus, with S. adoxoides as their closest relative. Therefore, the complete cp genomes in Urophysa species provide interesting insights and valuable information that can be used to identify related species and reconstruct their phylogeny.

Kranz and single-cell forms of C4 plants in the subfamily Suaedoideae show kinetic C4 convergence for PEPC and Rubisco with divergent amino acid substitutions.

The two carboxylation reactions performed by phosphoenolpyruvate carboxylase (PEPC) and ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) are vital in the fixation of inorganic carbon for C4 plants. The abundance of PEPC is substantially elevated in C4 leaves, while the location of Rubisco is restricted to one of two chloroplast types. These differences compared with C3 leaves have been shown to result in convergent enzyme optimization in some C4 species. Investigation into the kinetic properties of PEPC and Rubisco from Kranz C4, single cell C4, and C3 species in Chenopodiaceae s. s. subfamily Suaedoideae showed that these major carboxylases in C4 Suaedoideae species lack the same mutations found in other C4 systems which have been examined; but still have similar convergent kinetic properties. Positive selection analysis on the N-terminus of PEPC identified residues 364 and 368 to be under positive selection with a posterior probability >0.99 using Bayes empirical Bayes. Compared with previous analyses on other C4 species, PEPC from C4 Suaedoideae species have different convergent amino acids that result in a higher K m for PEP and malate tolerance compared with C3 species. Kinetic analysis of Rubisco showed that C4 species have a higher catalytic efficiency of Rubisco (k catc in mol CO2 mol(-1) Rubisco active sites s(-1)), despite lacking convergent substitutions in the rbcL gene. The importance of kinetic changes to the two-carboxylation reactions in C4 leaves related to amino acid selection is discussed.

Positive selection of Kranz and non-Kranz C4 phosphoenolpyruvate carboxylase amino acids in Suaedoideae (Chenopodiaceae).

In subfamily Suaedoideae, four independent gains of C4 photosynthesis are proposed, which includes two parallel origins of Kranz anatomy (sections Salsina and Schoberia) and two independent origins of single-cell C4 anatomy (Bienertia and Suaeda aralocaspica). Additional phylogenetic support for this hypothesis was generated from sequence data of the C-terminal portion of the phosphoenolpyruvate carboxylase (PEPC) gene used in C4 photosynthesis (ppc-1) in combination with previous sequence data. ppc-1 sequence was generated for 20 species in Suaedoideae and two outgroup Salsola species that included all types of C4 anatomies as well as two types of C3 anatomies. A branch-site test for positively selected codons was performed using the software package PAML. From labelling of the four branches where C4 is hypothesized to have developed (foreground branches), residue 733 (maize numbering) was identified to be under positive selection with a posterior probability >0.99 and residue 868 at the >0.95 interval using Bayes empirical Bayes (BEB). When labelling all the branches within C4 clades, the branch-site test identified 13 codons to be under selection with a posterior probability >0.95 by BEB; this is discussed considering current information on functional residues. The signature C4 substitution of an alanine for a serine at position 780 in the C-terminal end (which is considered a major determinant of affinity for PEP) was only found in four of the C4 species sampled, while eight of the C4 species and all the C3 species have an alanine residue; indicating that this substitution is not a requirement for C4 function.

Identification and adaptive evolution analysis of glutaredoxin genes in Populus spp.

Glutaredoxin (GRX) is a class of small redox proteins widely involved in cellular redox homeostasis and the regulation of various cellular processes. The role of GRX gene in the differentiation of Populus spp. is rarely reported. We compared the similarities and differences of GRX genes among four sections of poplar using bioinformatics, corrected the annotations of some GRX genes, and focused on analysing their transcript profiling and adaptive evolution in Populus spp. A total of 219 GRX genes were identified in four sections of poplar, among which annotations for 13 genes were corrected. Differences in GRX genes were found between sect. Turanga, represented by P. euphratica, and other poplar sections. Most notably, P. euphratica had the smallest number of duplication events for GRX genes (n = 9) and no tandem duplications, whereas there were >25 duplication events for all other poplars. Furthermore, we detected 18 pairs of GRX genes under positive selection pressure in various sections of poplar, and identified two groups of GRX genes in the Salicaceae that potentially underwent positive selection. Expression profiling results showed that the PtrGRX34 and its orthologous genes were upregulated under stress treatments. In summary, the GRX gene family underwent expansion during poplar differentiation, and some genes underwent rapid evolution during this process, which may be beneficial for Populus spp. to adapt to environmental changes. This study may provide more insights into the molecular mechanisms of Populus spp. adaptation to environmental changes and the adaptive evolution of GRX genes.

Whole-genome sequencing and evolutionary analysis of the wild edible mushroom, Morchella eohespera.

Morels (Morchella, Ascomycota) are an extremely desired group of edible mushrooms with worldwide distribution. Morchella eohespera is a typical black morel species, belonging to the Elata clade of Morchella species. The biological and genetic studies of this mushroom are rare, largely hindering the studies of molecular breeding and evolutionary aspects. In this study, we performed de novo sequencing and assembly of the M. eohespera strain m200 genome using the third-generation nanopore sequencing platform. The whole-genome size of M. eohespera was 53.81 Mb with a contig N50 of 1.93 Mb, and the GC content was 47.70%. A total of 9,189 protein-coding genes were annotated. Molecular dating showed that M. eohespera differentiated from its relative M. conica at ~19.03 Mya (million years ago) in Burdigalian. Evolutionary analysis showed that 657 gene families were contracted and 244 gene families expanded in M. eohespera versus the related morel species. The non-coding RNA prediction results showed that there were 336 tRNAs, 76 rRNAs, and 45 snRNAs in the M. eohespera genome. Interestingly, there was a high degree of repetition (20.93%) in the M. eohespera genome, and the sizes of long interspersed nuclear elements, short interspersed nuclear elements, and long terminal repeats were 0.83 Mb, 0.009 Mb, and 4.56 Mb, respectively. Additionally, selection pressure analysis identified that a total of 492 genes in the M. eohespera genome have undergone signatures of positive selection. The results of this study provide new insights into the genome evolution of M. eohespera and lay the foundation for in-depth research into the molecular biology of the genus Morchella in the future.

Chromosome-Level Genome Assembly and Multi-Omics Dataset Provide Insights into Isoflavone and Puerarin Biosynthesis in Pueraria lobata (Wild.) Ohwi.

Pueraria lobata (wild.) Ohwi is a leguminous plant and one of the traditional Chinese herbal medicines. Its puerarin extract is widely used in the pharmaceutical industry. This study reported a chromosome-level genome assembly for P. lobata and its characteristics. The genome size was ~939.2 Mb, with a contig N50 of 29.51 Mbp. Approximately 97.82% of the assembled sequences were represented by 11 pseudochromosomes. We identified that the repetitive sequences accounted for 63.50% of the P. lobata genome. A total of 33,171 coding genes were predicted, of which 97.34% could predict the function. Compared with other species, P. lobata had 757 species-specific gene families, including 1874 genes. The genome evolution analysis revealed that P. lobata was most closely related to Glycine max and underwent two whole-genome duplication (WGD) events. One was in a gamma event shared by the core dicotyledons at around 65 million years ago, and another was in the common ancestor shared by legume species at around 25 million years ago. The collinearity analysis showed that 61.45% of the genes (54,579 gene pairs) in G. max and P. lobata had collinearity. In this study, six unique PlUGT43 homologous genes were retrieved from the genome of P. lobata, and no 2-hydroxyisoflavanone 8-C-glucoside was found in the metabolites. This also revealed that the puerarin synthesis was mainly from the glycation of daidzein. The combined transcriptome and metabolome analysis suggested that two bHLHs, six MYBs and four WRKYs were involved in the expression regulation of puerarin synthesis structural genes. The genetic information obtained in this study provided novel insights into the biological evolution of P. lobata and leguminous species, and it laid the foundation for further exploring the regulatory mechanism of puerarin synthesis.

The early SARS-CoV-2 epidemic in Senegal was driven by the local emergence of B.1.416 and the introduction of B.1.1.420 from Europe.

Molecular surveillance of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is growing in west Africa, especially in the Republic of Senegal. Here, we present a molecular epidemiology study of the early waves of SARS-CoV-2 infections in this country based on Bayesian phylogeographic approaches. Whereas the first wave in mid-2020 was characterized by a significant diversification of lineages and predominance of B.1.416, the second wave in late 2020 was composed primarily of B.1.1.420. Our results indicate that B.1.416 originated in Senegal and was exported mainly to Europe. In contrast, B.1.1.420 was introduced from Italy, gained fitness in Senegal, and then spread worldwide. Since both B.1.416 and B.1.1.420 lineages carry several positive selected mutations in the spike and nucleocapsid genes, each of which may explain their local dominance, their mutation profiles should be carefully monitored. As the pandemic continues to evolve, molecular surveillance in all regions of Africa will play a key role in stemming its spread.

Genomic Characterization of Jumbo Salmonella Phages That Effectively Target United Kingdom Pig-Associated Salmonella Serotypes.

A common cause of human food poisoning is through ingestion of pork products contaminated with Salmonella spp. Worryingly multi-drug resistant (MDR) Salmonella strains have been isolated from pigs, which motivates the need for alternative antimicrobials. In this study isolation and characterization of 21 lytic Salmonella phages is described. All 21 phages, labeled as SPFM phages were shown to efficiently infect MDR Salmonella strains isolated from United Kingdom pigs and phages SPFM1, SPFM3, SPFM10, SPFM14, SPFM15, SPFM17, and SPFM19 could lyse 100% of strains tested. The phage genome sizes range from 233 to 242 Kb, which qualifies them as jumbo phages. All SPFM phage genomes are approximately 95% similar to each other by average nucleotide identity, they encode between 258-307 coding sequences and share 188 core genes. Phylogenetic analysis shows these phages are most similar to phages of the genus Seoulvirus and to further characterize phages within the genus, genes under positive selection were identified. Several of the genes under evolutionary selection pressure were predicted to encode for proteins that interact with bacteria. We describe the phenotypic and genetic characterization of this novel Salmonella phage set. As the phages efficiently kill MDR Salmonella strains, they may offer a promising alternative to antibiotics.

Positive selection analysis highlights key positions in plant PP2A regulatory subunits.

A versatile hub for cellular control, the eukaryotic protein phosphatase 2A (PP2A) enzyme family is thought to achieve specificity through combinatorial complexity. Phylogenetic analysis has revealed that expansion of PP2A gene families resulted from whole genome duplications followed by non-random gene loss, and selection analysis suggests that retention of B56/PPP2R5 gene family members after genome duplication events was driven by functional diversification. Here we identify the sites at which positive selection is detected in the plant B56 gene family, and discuss the significance of selection at these positions in the context of PP2A holoenzyme structure. The pattern of positive selection observed in the B11 subclade is distinctive, and suggests selective pressure on interactions with substrates and the enzymatic core.

Genetic characterization of seasonal influenza A (H3N2) viruses in Ontario during 2010-2011 influenza season: high prevalence of mutations at antigenic sites.

BACKGROUND: The direct effect of antigenic site mutations in influenza viruses on antigenic drift and vaccine effectiveness is poorly understood. OBJECTIVE: To investigate the genetic and antigenic characteristics of human influenza A (H3N2) viruses circulating in Ontario during the early 2010-2011 winter season. STUDY DESIGN: We sequenced the hemagglutinin (HA) and neuraminidase (NA) genes from 41 A(H3N2) viruses detected in nasopharyngeal specimens. Strain typing was performed by hemagglutination inhibition (HI) assay. Molecular and phylogenetic tree analyses were conducted. RESULTS: HA and NA genes showed high similarity to the 2010-2011 vaccine strain, A/Perth/16/2009 (H3N2)-like virus (97·7-98·5% and 98·7-99·5% amino acid (AA) identity, respectively). Compared to A/Perth/16/2009 strain, HA gene mutations were documented at 28 different AA positions across all five H3 antigenic sites, with a range of 5-11 mutations in individual viruses. Thirty-six (88%) viruses had 8 AA substitutions in common; none of these had reduced HI titer. Among Ontario isolates, 11 antigenic site AAs were positively selected with an increase in glycosylation sites. CONCLUSION: The presence of antigenic site mutations with high frequency among 2010-2011 influenza H3N2 isolates confirms ongoing adaptive H3N2 evolution. These may represent early phylogenetic changes that could cause antigenic drift with further mutations. Clinical relevance of antigenic site mutations not causing drift in HI assays is unknown and requires further investigation. In addition, viral sequencing information will assist with vaccine strain planning and may facilitate early detection of vaccine escape.