Search

Abstract To reveal the molecular function of elongation family of very long chain fatty acids(ELO) protein in Cyrtotrachelus buqueti, we have identified 15 ELO proteins from C.buqueti genome. 15 CbuELO proteins were located on four chromosomes. Their isoelectric points ranged from 9.22 to 9.68, and they were alkaline. These CbuELO proteins were stable and hydrophobic. CbuELO proteins had transmembrane movement, and had multiple phosphorylation sites. The secondary structure of CbuELO proteins was mainly α-helix. A total of 10 conserved motifs were identified in CbuELO protein family. Phylogenetic analysis showed that molecular evolutionary relationships of ELO protein family between C. buqueti and Tribolium castaneum was the closest. Developmental transcriptome analysis indicated that CbuELO10, CbuELO13 and CbuELO02 genes were key enzyme genes that determine the synthesis of very long chain fatty acids in pupae and eggs, CbuELO6 and CbuELO7 were that in the male, and CbuELO8 and CbuELO11 were that in the larva. Transcriptome analysis under different temperature conditions indicated that CbuELO1, CbuELO5, CbuELO12 and CbuELO14 participated in regulating temperature stress responses. Transcriptome analysis at different feeding times showed CbuELO12 gene expression level in all feeding time periods was significant downregulation. The qRT-PCR experiment verified expression level changes of CbuELO gene family under different temperature and feeding time conditions. Protein-protein interaction analysis showed that 9 CbuELO proteins were related to each other, CbuELO1, CbuELO4 and CbuELO12 had more than one interaction relationship. These results lay a theoretical foundation for further studying its molecular function during growth and development of C. buqueti.

Abstract DNA binding with one finger(Dof) gene family is a class of transcription factors which play an important role on plant growth and development. Genome-wide identification results indicated that there were 45 Dof genes(ColDof) in C.oleifera genome. All 45 ColDof proteins were non-transmembrane and non-secretory proteins. Phosphorylation site analysis showed that biological function of ColDof proteins were mainly realized by phosphorylation at serine (Ser) site. The secondary structure of 44 ColDof proteins was dominated by random coil, and only one ColDof protein was dominated by α-helix. ColDof genes’ promoter region contained a variety of cis-acting elements, including light responsive regulators, gibberellin responsive regulators, abscisic acid responsive regulators, auxin responsive regulators and drought induction responsive regulators. The SSR sites analysis showed that the proportion of single nucleotide repeats and the frequency of A/T in ColDof genes were the largest. Non-coding RNA analysis showed that 45 ColDof genes contained 232 miRNAs. Transcription factor binding sites of ColDof genes showed that ColDof genes had 5793 ERF binding sites, 4381 Dof binding sites, 2206 MYB binding sites, 3702 BCR-BPC binding sites. ColDof9, ColDof39 and ColDof44 were expected to have the most TFBSs. The collinearity analysis showed that there were 40 colinear locis between ColDof proteins and AtDof proteins. Phylogenetic analysis showed that ColDof gene family was most closely related to that of Camellia sinensis var. sinensis cv.Biyun and Camellia lanceoleosa. Protein-protein interaction analysis showed that ColDof34, ColDof20, ColDof28, ColDof35, ColDof42 and ColDof26 had the most protein interactions. The transcriptome analysis of C. oleifera seeds showed that 21 ColDof genes were involved in the growth and development process of C. oleifera seeds, and were expressed in 221 C. oleifera varieties. The results of qRT-PCR experiments treated with different concentrations NaCl and PEG6000 solutions indicated that ColDof1, ColDof2, ColDof14 and ColDof36 not only had significant molecular mechanisms for salt stress tolerance, but also significant molecular functions for drought stress tolerance in C. oleifera. The results of this study provide a reference for further understanding of the function of ColDof genes in C.oleifera.

Abstract Background Insect Cytochrome P450 monooxygenase (CYPs or P450s) plays an important role in detoxifying insecticides, causing insect populations to develop resistance. However, the molecular functions of P450 gene family in Cyrtotrachelus buqueti genome are still lacking. Results In this study, 71 CbuP450 genes have been identified. The amino acids length of CbuP450 proteins was between 183 aa ~ 1041 aa. They are proteins with transmembrane domains. The main component of their secondary structure is α-helix and random coils. Phylogenetic analysis showed that C. buqueti and Rhynchophorus ferrugineus were the most closely related. This gene family has 29 high-frequency codons, which tend to use A/T bases and A/T ending codons. Gene expression analysis showed that CbuP450_23 in the female adult may play an important role on high temperature resistance, and CbuP450_17 in the larval may play an important role on low temperature tolerance. CbuP450_10, CbuP450_17, CbuP450_23, CbuP450_10, CbuP450_16, CbuP450_20, CbuP450_23 and CbuP450_ 29 may be related to the regulation of bamboo fiber degradation genes in C. buqueti. Protein interaction analysis indicates that most CbuP450 proteins are mainly divided into three aspects: encoding the biosynthesis of ecdysteroids, participating in the decomposition of synthetic insecticides, metabolizing insect hormones, and participating in the detoxification of compounds. Conclusions We systematically analyzed the gene and protein characteristics, gene expression, and protein interactions of CbuP450 gene family, revealing the key genes involved in the stress response of CbuP450 gene family in the resistance of C. buqueti to high or low temperature stress, and identified the key CbuP450 proteins involved in important life activity metabolism. These results provided a reference for further research on the function of P450 gene family in C. buqueti.

Abstract Background Myxozoa is a class of cnidarian parasites that encompasses over 2,400 species. Phylogenetic relationships among myxozoans remain highly debated, owing to both a lack of informative morphological characters and a shortage of molecular markers. Mitochondrial (mt) genomes are a common marker in phylogeny and biogeography. However, only five complete myxozoan mt genomes have been sequenced: four belonging to two closely related genera, Enteromyxum and Kudoa, and one from the genus Myxobolus. Interestingly, while cytochrome oxidase genes could be identified in Enteromyxum and Kudoa, no such genes were found in Myxobolus squamalis, and another member of the Myxobolidae (Henneguya salminicola) was found to have lost its entire mt genome. To evaluate the utility of mt genomes to reconstruct myxozoan relationships and to understand if the loss of cytochrome oxidase genes is a characteristic of myxobolids, we sequenced the mt genome of five myxozoans (Myxobolus wulii, M. honghuensis, M. shantungensis, Thelohanellus kitauei and, Sphaeromyxa zaharoni) using Illumina and Oxford Nanopore platforms. Results Unlike Enteromyxum, which possesses a partitioned mt genome, the five mt genomes were encoded on single circular chromosomes. An mt plasmid was found in M. wulii, as described previously in Kudoa iwatai. In all new myxozoan genomes, five protein-coding genes (cob, cox1, cox2, nad1, and nad5) and two rRNAs (rnl and rns) were recognized, but no tRNA. We found that Myxobolus and Thelohanellus species shared unidentified reading frames, supporting the view that these mt open reading frames are functional. Our phylogenetic reconstructions based on the five conserved mt genes agree with previously published trees based on the 18S rRNA gene. Conclusions Our results suggest that the loss of cytochrome oxidase genes is not a characteristic of all myxobolids, the ancestral myxozoan mt genome was likely encoded on a single circular chromosome, and mt plasmids exist in a few lineages. Our findings indicate that myxozoan mt sequences are poor markers for reconstructing myxozoan phylogenetic relationships because of their fast-evolutionary rates and the abundance of repeated elements, which complicates assembly.

ABSTRACTThe nuclear envelope (NE) separates translation and transcription and is the location of multiple functions, including chromatin organization and nucleocytoplasmic transport. The molecular basis for many of these functions have diverged between eukaryotic lineages. Trypanosoma brucei, a member of the early branching eukaryotic lineage Discoba, highlights many of these, including a distinct lamina and kinetochore composition. Here, we describe a cohort of proteins interacting with both the lamina and NPC, which we term lamina-associated proteins (LAPs). LAPs represent a diverse group of proteins, including two candidate NPC-anchoring pore membrane proteins (POMs) with architecture conserved with S. cerevisiae and H. sapiens, and additional peripheral components of the NPC. While many of the LAPs are Kinetoplastid specific, we also identified broadly conserved proteins, indicating an amalgam of divergence and conservation within the trypanosome NE proteome, highlighting the diversity of nuclear biology across the eukaryotes, increasing our understanding of eukaryotic and NPC evolution.

The slender anole, Anolis apletophallus , is a small arboreal lizard of the rainforest understory of central and eastern Panama. This species has been the subject of numerous ecological and evolutionary studies over the past 60 years as a result of attributes that make it especially amenable to field and laboratory science. Slender anoles are highly abundant, short-lived (nearly 100% annual turnover), easy to manipulate in both the lab and field, and are ubiquitous in the forests surrounding the Smithsonian Tropical Research Institute in Panama, where researchers have access to high-quality laboratory facilities. Here, we present a high-quality genome for the slender anole, which is an important new resource for studying this model species. We assembled and annotated the slender anole genome by combining 3 technologies: Oxford Nanopore, 10× Genomics Linked-Reads, and Dovetail Omni-C. We compared this genome with the recently published brown anole (Anolis sagrei) and the canonical green anole (Anolis carolinensis) genomes. Our genome is the first assembled for an Anolis lizard from mainland Central or South America, the regions that host the majority of diversity in the genus. This new reference genome is one of the most complete genomes of any anole assembled to date and should facilitate deeper studies of slender anole evolution, as well as broader scale comparative genomic studies of both mainland and island species. In turn, such studies will further our understanding of the well-known adaptive radiation of Anolis lizards. [ABSTRACT FROM AUTHOR]

The NODULE-INCEPTION-like protein (NLP) family is a plant-specific transcription factor (TF) family involved in nitrate transport and assimilation in plants, which are essential for improving plant nitrogen use efficiency. Currently, the molecular nature and evolutionary trajectory of NLP genes in the C4 model crop foxtail millet are unknown. Therefore, we performed a comprehensive analysis of NLP and molecular evolution in foxtail millet by scanning the genomes of foxtail millet and representative species of the plant kingdom. We identified seven NLP genes in the foxtail millet genome, all of which are individually and separately distributed on different chromosomes. They were not structurally identical to each other and were mainly expressed on root tissues. We unearthed two key genes (Si5G004100.1 and Si6G248300.1) with a variety of excellent characteristics. Regarding its molecular evolution, we found that NLP genes in Gramineae mainly underwent dispersed duplication, but maize NLP genes were mainly generated via WGD events. Other factors such as base mutations and natural selection have combined to promote the evolution of NLP genes. Intriguingly, the family in plants showed a gradual expansion during evolution with more duplications than losses, contrary to most gene families. In conclusion, this study advances the use of NLP genetic resources and the understanding of molecular evolution in cereals.

Lagerstroemia indica is an important commercial tree known for the ornamental value. In this study, the complete chloroplast genome sequence of Lagerstroemia indica “Pink Velour” (Lagerstroemia “Pink Velour”) was 152,174 bp in length with a GC content of 39.50%. It contained 85 protein coding genes (PCGs), 37 tRNAs, and 8 rRNA genes. 207 simple sequence repeats (SSRs) and 31 codons with relative synonymous codon (RSCU)value > 1 were detected. Phylogenetic analysis divided 10 Lagerstroemia species into evolutionary branches of clade A and clade B. We conducted a comparative analysis of Lagerstroemia “Pink Velours” complete chloroplast genome with the genomes of six closely related Lagerstroemia species from different origins. The structural features of all seven species were similar, except for the deletion of ycf1 nucleobases at the JSA boundary. The large single-copy (LSC) and the small single-copy (SSC) had a higher sequence divergence than the IR region, and 8 genes that were highly divergent (trnK-UUU, petN, psbF, psbJ, ndhE, ndhD, ndhI, ycf1) had been identified and could be used as molecular markers in future studies. High nucleotide diversity was present in genes belonging to the photosynthesis category. Mutation of single nucleic acid was mainly influenced by codon usage. The value percentage of nonsynonymous substitutions (Ka) and synonymous substitutions (Ks) in 6 Lagerstroemia species revealed that more photosynthesis genes have Ka or Ks only in Lagerstroemia fauriei, Lagerstroemia limii, and Lagerstroemia subcostata. These advances will facilitate the breeding of closely related Lagerstroemia species and deepen understanding on climatic adaptation of Lagerstroemia plants.

Commonly used methods for inferring phylogenies were designed before the emergence of high-throughput sequencing and can generally not accommodate the challenges associated with noisy, diploid sequencing data. In many applications, diploid genomes are still treated as haploid through the use of ambiguity characters; while the uncertainty in genotype calling-arising as a consequence of the sequencing technology-is ignored. In order to address this problem, we describe two new probabilistic approaches for estimating genetic distances: distAngsd-geno and distAngsd-nuc, both implemented in a software suite named distAngsd. These methods are specifically designed for next-generation sequencing data, utilize the full information from the data, and take uncertainty in genotype calling into account. Through extensive simulations, we show that these new methods are markedly more accurate and have more stable statistical behaviors than other currently available methods for estimating genetic distances-even for very low depth data with high error rates.

• Congenital myopathies have been characterized based on structural abnormalities. • Genetic alterations are typically associated with certain disease entities. • New genes recently identified, enlarge the spectrum of congenital myopathies. • Certain congenital myopathies do not show morphological characteristics. Congenital myopathies are defined by early clinical onset, slow progression, hereditary nature and disease-specific myopathological lesions – however, with exceptions - demanding special techniques in regard to morphological diagnostic and research work-up. To identify an index disease in a family requires a muscle biopsy – and no congenital myopathy has ever been first described at autopsy. The nosographic history commenced when – in addition to special histopathological techniques in the earliest classical triad of central core disease, 1956, nemaline myopathy, 1963, and centronuclear myopathy, 1966/67, within a decade - electron microscopy and enzyme histochemistry were applied to unfixed frozen muscle tissue and, thus, revolutionized diagnostic and research myopathology. During the following years, the list of structure-defined congenital myopathies grew to some 40 conditions. Then, the introduction of immunohistochemistry allowed myopathological documentation of proteins and their abnormalities in individual congenital myopathies. Together with the diagnostic evolution of molecular genetics, many more congenital myopathies were described, without new disease-specific lesions or only already known ones. These were nosographically defined by individual mutations in hitherto congenital myopathies-unrelated genes. This latter development may also affect the nomenclature of congenital myopathies in that the mutant gene needs to be attached to the individually identified congenital myopathies with or without the disease-specific lesion, such as CCD- RYR1 or CM -RYR1. This principle is similar to that of the nomenclature of Congenital Disorders of Glycosylation. Retroactive molecular characterization of originally and first described congenital myopathies has only rarely been achieved. [ABSTRACT FROM AUTHOR]

The transfer of DNA among distantly related organisms is relatively common in bacteria but less prevalent in eukaryotes. Among fungi and plants, few of these events have been reported. Two segments of fungal mitochondrial DNA have been recently discovered in the mitogenome of orchids. Here, we build on that work to understand the timing of those transfer events, which orchids retain the fungal DNA and the fate of the foreign DNA during orchid evolution. We update the content of the large DNA fragment and establish that it was transferred to the most recent common ancestor of a highly diverse clade of epidendroid orchids that lived ~28–43 Mya. Also, we present hypotheses of the origin of the small transferred fragment. Our findings deepen the knowledge of these interesting DNA transfers among organelles and we formulate a probable mechanism for these horizontal gene transfer events. [ABSTRACT FROM AUTHOR]

Argonaute (AGO) proteins are central players in RNA interference in eukaryotes. They associate with small RNAs (sRNA) and lead to transcriptional or posttranscriptional silencing of targets, thereby regulating diverse biological processes. The molecular and biological functions of AGO proteins have been extensively characterized, particularly in a few angiosperm species, leading to the recognition that the AGO family has expanded to accommodate diverse sRNAs thereby performing diverse biological functions. However, understanding of the expansion of AGO proteins in plants is still limited, due to a dearth of knowledge of AGO proteins in green algal groups. Here, we identified more than 2900 AGO proteins from 244 plant species, including green algae, and performed a large-scale phylogenetic analysis. The phylogeny shows that the plant AGO family gave rise to four clades after the emergence of hydrobiontic algae and prior to the emergence of land plants. Subsequent parallel expansion in ferns and angiosperms resulted in eight main clades in angiosperms: AGO2, AGO7, AGO6, AGO4, AGO1, AGO10a, AGO10b and AGO5. On the basis of this phylogeny, we identified two novel AGO4 orthologs that Arabidopsis does not have, and redefined AGO10, which is composed of AGO10a and AGO10b. Finally, we propose a hypothetical evolutionary model of AGO proteins in plants. Our studies provide a deeper understanding of the phylogenetic relationships of AGO family members in the green lineage, which would help to further reveal their roles as RNAi effectors.

The XBB.1.16 SARS-CoV-2 variant, also known as Arcturus, is a recent descendant lineage of the recombinant XBB (nicknamed Gryphon). Compared to its direct progenitor, XBB.1, XBB.1.16 carries additional spike mutations in key antigenic sites, potentially conferring an ability to evade the immune response compared to other circulating lineages. In this context, we conducted a comprehensive genome-based survey to gain a detailed understanding of the evolution and potential dangers of the XBB.1.16 variant, which became dominant in late June. Genetic data indicates that the XBB.1.16 variant exhibits an evolutionary background with limited diversification, unlike dangerous lineages known for rapid changes. The evolutionary rate of XBB.1.16, which amounts to 3.95 × 10−4 subs/site/year, is slightly slower than that of its direct progenitors, XBB and XBB.1.5, which have been circulating for several months. A Bayesian Skyline Plot reconstruction suggests that the peak of genetic variability was reached in early May 2023, and currently, it is in a plateau phase with a viral population size similar to the levels observed in early March. Structural analyses indicate that, overall, the XBB.1.16 variant does not possess structural characteristics markedly different from those of the parent lineages, and the theoretical affinity for ACE2 does not seem to change among the compared variants. In conclusion, the genetic and structural analyses of SARS-CoV-2 XBB.1.16 do not provide evidence of its exceptional danger or high expansion capability. Detected differences with previous lineages are probably due to genetic drift, which allows the virus constant adaptability to the host, but they are not necessarily connected to a greater danger. Nevertheless, continuous genome-based monitoring is essential for a better understanding of its descendants and other lineages. [ABSTRACT FROM AUTHOR]

Recent advances in malaria genetics and genomics have transformed many aspects of malaria research in areas of molecular evolution, epidemiology, transmission, host–parasite interaction, drug resistance, pathogenicity, and vaccine development. Here, in addition to introducing some background information on malaria parasite biology, parasite genetics/genomics, and genotyping methods, we discuss some applications of genetic and genomic approaches in vaccine development and in studying interactions with microbiota. Genetic and genomic data can be used to search for novel vaccine targets, design an effective vaccine strategy, identify protective antigens in a whole-organism vaccine, and evaluate the efficacy of a vaccine. Microbiota has been shown to influence disease outcomes and vaccine efficacy; studying the effects of microbiota in pathogenicity and immunity may provide information for disease control. Malaria genetics and genomics will continue to contribute greatly to many fields of malaria research. [ABSTRACT FROM AUTHOR]

Bitter taste perception plays vital roles in animal behavior and fitness. By signaling the presence of toxins in foods, particularly noxious defense compounds found in plants, it enables animals to avoid exposure. In vertebrates, bitter perception is initiated by TAS2Rs, a family of G protein-coupled receptors expressed on the surface of taste buds. There, oriented toward the interior of the mouth, they monitor the contents of foods, drinks and other substances as they are ingested. When bitter compounds are encountered, TAS2Rs respond by triggering neural pathways leading to sensation. The importance of this role placed TAS2Rs under selective pressures in the course of their evolution, leaving signatures in patterns of gene gain and loss, sequence polymorphism, and population structure consistent with vertebrates' diverse feeding ecologies. The protective value of bitter taste is reduced in modern humans because contemporary food supplies are safe and abundant. However, this is not always the case. Some crops, particularly in the developing world, retain surprisingly high toxicity and bitterness remains an important measure of safety. Bitter perception also shapes health through its influence on preference driven behaviors such as diet choice, alcohol intake and tobacco use. Further, allelic variation in TAS2Rs is extensive, leading to individual differences in taste sensitivity that drive these behaviors, shaping susceptibility to disease. Thus, bitter taste perception occupies a critical intersection between ancient evolutionary processes and modern human health.

Human severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is most closely related, by average genetic distance, to two coronaviruses isolated from bats, RaTG13 and RmYN02. However, there is a segment of high amino acid similarity between human SARS-CoV-2 and a pangolin-isolated strain, GD410721, in the receptor-binding domain (RBD) of the spike protein, a pattern that can be caused by either recombination or by convergent amino acid evolution driven by natural selection. We perform a detailed analysis of the synonymous divergence, which is less likely to be affected by selection than amino acid divergence, between human SARS-CoV-2 and related strains. We show that the synonymous divergence between the bat-derived viruses and SARS-CoV-2 is larger than between GD410721 and SARS-CoV-2 in the RBD, providing strong additional support for the recombination hypothesis. However, the synonymous divergence between pangolin strain and SARS-CoV-2 is also relatively high, which is not consistent with a recent recombination between them, instead, it suggests a recombination into RaTG13. We also find a 14-fold increase in the dN /dS ratio from the lineage leading to SARS-CoV-2 to the strains of the current pandemic, suggesting that the vast majority of nonsynonymous mutations currently segregating within the human strains have a negative impact on viral fitness. Finally, we estimate that the time to the most recent common ancestor of SARS-CoV-2 and RaTG13 or RmYN02 based on synonymous divergence is 51.71 years (95% CI, 28.11-75.31) and 37.02 years (95% CI, 18.19-55.85), respectively.

BackgroundA mechanistic understanding of the spread of SARS-CoV-2 and diligent tracking of ongoing mutagenesis are of key importance to plan robust strategies for confining its transmission. Large numbers of available sequences and their dates of transmission provide an unprecedented opportunity to analyze evolutionary adaptation in novel ways. Addition of high-resolution structural information can reveal the functional basis of these processes at the molecular level. Integrated systems biology-directed analyses of these data layers afford valuable insights to build a global understanding of the COVID-19 pandemic.ResultsHere we identify globally distributed haplotypes from 15,789 SARS-CoV-2 genomes and model their success based on their duration, dispersal, and frequency in the host population. Our models identify mutations that are likely compensatory adaptive changes that allowed for rapid expansion of the virus. Functional predictions from structural analyses indicate that, contrary to previous reports, the Asp614Gly mutation in the spike glycoprotein (S) likely reduced transmission and the subsequent Pro323Leu mutation in the RNA-dependent RNA polymerase led to the precipitous spread of the virus. Our model also suggests that two mutations in the nsp13 helicase allowed for the adaptation of the virus to the Pacific Northwest of the USA. Finally, our explainable artificial intelligence algorithm identified a mutational hotspot in the sequence of S that also displays a signature of positive selection and may have implications for tissue or cell-specific expression of the virus.ConclusionsThese results provide valuable insights for the development of drugs and surveillance strategies to combat the current and future pandemics.

We introduce a systematic method of approximating finite-time transition probabilities for continuous-time insertion-deletion models on sequences. The method uses automata theory to describe the action of an infinitesimal evolutionary generator on a probability distribution over alignments, where both the generator and the alignment distribution can be represented by pair hidden Markov models (HMMs). In general, combining HMMs in this way induces a multiplication of their state spaces; to control this, we introduce a coarse-graining operation to keep the state space at a constant size. This leads naturally to ordinary differential equations for the evolution of the transition probabilities of the approximating pair HMM. The TKF91 model emerges as an exact solution to these equations for the special case of single-residue indels. For the more general case of multiple-residue indels, the equations can be solved by numerical integration. Using simulated data, we show that the resulting distribution over alignments, when compared to previous approximations, is a better fit over a broader range of parameters. We also propose a related approach to develop differential equations for sufficient statistics to estimate the underlying instantaneous indel rates by expectation maximization. Our code and data are available at https://github.com/ihh/trajectory-likelihood.

Thyroid cancer surgery has evolved dramatically with advances in our understanding of the biological behaviour of WDTC. Molecular profiling is shedding light on the subset that may behave aggressively. In an era when thyroid cancer management is becoming increasingly conservative, decision making regarding the extent of surgery must be objectively guided by molecular markers. The aim of the present article is to summarise the current published literature and provide possible practice recommendations. An online search for relevant published articles was performed using several databases. Title, abstract, and full-text screening, along with data extraction, was performed by two independent reviewers after the inclusion and exclusion criteria were defined. A total of 1241 articles were identified, and 82 relevant articles were extracted and scrutinised. BRAF V600E and TERT promoter mutations were found to be associated with an increased risk of disease recurrence and distant metastases. Several other mutations have been identified that enhance disease aggressiveness (such as RET/PTC, PTEN, and TP53). One of the most important determinants of the outcome in WDTC is the extent of surgical resection. The evolution of molecular testing has reached a stage of personalised incorporation into surgical practice. Guidelines for molecular testing and surgery in WDTC will need to be clearly defined, arguably representing the next chapter in the management of the disease. [ABSTRACT FROM AUTHOR]

While the relationship of protective human leukocyte antigen (HLA) class I alleles and HIV progression is well defined, the interaction of HLA-mediated protection and CD8 T-cell exhaustion is less well characterized. To gain insight into the influence of HLA-B*57:01 on the deterioration of CD8 T-cell responses during HIV infection in the absence of antiretroviral treatment, we compared HLA-B*57:01-restricted HIV-specific CD8 T-cell responses to responses restricted by other HLA class I alleles longitudinally after control of peak viremia. Detailed characterization of polyfunctionality, differentiation phenotypes, transcription factor, and inhibitory receptor expression revealed progression of CD8 T-cell exhaustion over the course of the infection in both patient groups. However, early effects on the phenotype of the total CD8 T-cell population were apparent only in HLA-B*57-negative patients. The HLA-B*57:01-restricted, HIV epitope-specific CD8 T-cell responses showed beneficial functional patterns and significantly lower frequencies of inhibitory receptor expression, i.e., PD-1 and coexpression of PD-1 and TIGIT, within the first year of infection. Coexpression of PD-1 and TIGIT was correlated with clinical markers of disease progression and declining percentages of the T-bethi Eomesdim CD8 T-cell population. In accordance with clinical and immunological deterioration in the HLA-B*57:01 group, the difference in PD-1 and TIGIT receptor expression did not persist to later stages of the disease.IMPORTANCE Given the synergistic nature of TIGIT and PD-1, the coexpression of those inhibitory receptors should be considered when evaluating T-cell pathogenesis, developing immunomodulatory therapies or vaccines for HIV, and when using immunotherapy or vaccination for other causes in HIV-infected patients. HIV-mediated T-cell exhaustion influences the patient´s disease progression, immune system and subsequently non-AIDS complications, and efficacy of vaccinations against other pathogens. Consequently, the possibilities of interfering with exhaustion are numerous. Expanding the use of immunomodulatory therapies to include HIV treatment depends on information about possible targets and their role in the deterioration of the immune system. Furthermore, the rise of immunotherapies against cancer and elevated cancer incidence in HIV-infected patients together increase the need for detailed knowledge of T-cell exhaustion and possible interactions. A broader approach to counteract immune exhaustion to alleviate complications and improve efficacy of other vaccines also promises to increase patients' health and quality of life.

Large population-genomic sequencing studies can enable highly-powered analyses of sequence signatures of natural selection. Genome repositories now available for Saccharomyces yeast make it a premier model for studies of the molecular mechanisms of adaptation. We mined the genomes of hundreds of isolates of the sister species S. cerevisiae and S. paradoxus to identify sequence hallmarks of adaptive divergence between the two. From the top hits we focused on a set of genes encoding membrane proteins of the peroxisome, an organelle devoted to lipid breakdown and other specialized metabolic pathways. In-depth population- and comparative-genomic sequence analyses of these genes revealed striking divergence between S. cerevisiae and S. paradoxus And from transcriptional profiles we detected non-neutral, directional cis-regulatory variation at the peroxisome membrane genes, with overall high expression in S. cerevisiae relative to S. paradoxus Taken together, these data support a model in which yeast species have differentially tuned the expression of peroxisome components to boost their fitness in distinct niches.

There are many compelling examples of molecular convergence at individual genes. However, the prevalence and the relative importance of adaptive genome-wide convergence remain largely unknown. Many recent works have reported striking examples of excess genome-wide convergence, but some of these studies have been called into question because of the use of inappropriate null models. Here, we sequenced and compared the genomes of 12 species of anole lizards that have independently converged on suites of adaptive behavioral and morphological traits. Despite extensive searches for a genome-wide signature of molecular convergence, we found no evidence supporting molecular convergence at specific amino acids either at individual genes or at genome-wide comparisons; we also uncovered no evidence supporting an excess of adaptive convergence in the rates of amino acid substitutions within genes. Our findings indicate that comprehensive phenotypic convergence is not mirrored at genome-wide protein-coding levels in anoles, and therefore, that adaptive phenotypic convergence is likely not constrained by the evolution of many specific protein sequences or structures.

Comparative genomics has begun to elucidate the genomic basis of social life in insects, but insight into the genomic basis of spider sociality has lagged behind. To begin, to characterize genomic signatures associated with the evolution of social life in spiders, we performed one of the first spider comparative genomics studies including five solitary species and two social species, representing two independent origins of sociality in the genus Stegodyphus. We found that the two social spider species had a large expansion of gene families associated with transport and metabolic processes and an elevated genome-wide rate of molecular evolution compared with the five solitary spider species. Genes that were rapidly evolving in the two social species relative to the five solitary species were enriched for transport, behavior, and immune functions, whereas genes that were rapidly evolving in the solitary species were enriched for energy metabolism processes. Most rapidly evolving genes in the social species Stegodyphus dumicola were broadly expressed across four tissues and enriched for transport functions, but 12 rapidly evolving genes showed brain-specific expression and were enriched for social behavioral processes. Altogether, our study identifies putative genomic signatures and potential candidate genes associated with spider sociality. These results indicate that future spider comparative genomic studies, including broader sampling and additional independent origins of sociality, can further clarify the genomic causes and consequences of social life.

Many biomolecular machines need to be both fast and efficient. How has evolution optimized these machines along the tradeoff between speed and efficiency? We explore this question using optimizable dynamical models along coordinates that are plausible evolutionary degrees of freedom. Data on 11 motors and ion pumps are consistent with the hypothesis that evolution seeks an optimal balance of speed and efficiency, where any further small increase in one of these quantities would come at great expense to the other. For FoF1-ATPases in different species, we also find apparent optimization of the number of subunits in the c-ring, which determines the number of protons pumped per ATP synthesized. Interestingly, these ATPases appear to more optimized for efficiency than for speed, which can be rationalized through their key role as energy transducers in biology. The present modeling shows how the dynamical performance properties of biomolecular motors and pumps may have evolved to suit their corresponding biological actions.

Similar to other RNA viruses, the emergence of Betacoronavirus relies on cross-species viral transmission, which requires careful health surveillance monitoring of protein-coding information as well as genome-wide analysis. Although the evolutionary jump from natural reservoirs to humans may be mainly traced-back by studying the effect that hotspot mutations have on viral proteins, it is largely unexplored if other impacts might emerge on the structured RNA genome of Betacoronavirus. In this survey, the protein-coding and viral genome architecture were simultaneously studied to uncover novel insights into cross-species horizontal transmission events. We analysed 1,252,952 viral genomes of SARS-CoV, MERS-CoV, and SARS-CoV-2 distributed across the world in bats, intermediate animals, and humans to build a new landscape of changes in the RNA viral genome. Phylogenetic analyses suggest that bat viruses are the most closely related to the time of most recent common ancestor of Betacoronavirus, and missense mutations in viral proteins, mainly in the S protein S1 subunit: SARS-CoV (G > T; A577S); MERS-CoV (C > T; S746R and C > T; N762A); and SARS-CoV-2 (A > G; D614G) appear to have driven viral diversification. We also found that codon sites under positive selection on S protein overlap with non-compensatory mutations that disrupt secondary RNA structures in the RNA genome complement. These findings provide pivotal factors that might be underlying the eventual jumping the species barrier from bats to intermediate hosts. Lastly, we discovered that nearly half of the Betacoronavirus genomes carry highly conserved RNA structures, and more than 90% of these RNA structures show negative selection signals, suggesting essential functions in the biology of Betacoronavirus that have not been investigated to date. Further research is needed on negatively selected RNA structures to scan for emerging functions like the potential of coding virus-derived small RNAs and to develop new candidate antiviral therapeutic strategies.

Abstract Background Speciation genes contribute disproportionately to species divergence, but few examples exist, especially in vertebrates. Here we test whether Zan, which encodes the sperm acrosomal protein zonadhesin that mediates species-specific adhesion to the egg’s zona pellucida, is a speciation gene in placental mammals. Results Genomic ontogeny reveals that Zan arose by repurposing of a stem vertebrate gene that was lost in multiple lineages but retained in Eutheria on acquiring a function in egg recognition. A 112-species Zan sequence phylogeny, representing 17 of 19 placental Orders, resolves all species into monophyletic groups corresponding to recognized Orders and Suborders, with

Increasing evidence indicates that wide distributed bryophyte taxa with homogeneous morphology may represent separate evolutionary lineages. The evolutionary histories of these cryptic lineages may be related to historical factors, such as the climatic oscillations in the Quaternary. Thus, the post-glacial demographic signatures paired with cryptic speciation may result in complex phylogeographic patterns. This research has two aims: to determine whether the widespread moss Racomitrium lanuginosum represents cryptic molecular taxa across the Northern Hemisphere and to infer the effects of Quaternary glaciations on spatial genetic diversity. We used the internal transcribed spacer (ITS) marker to resolve the phylogeographic history of the species and single nucleotide polymorphisms (genotyping-by-sequencing) to infer the genetic structure and demographic history. Finally, we assessed the historical changes in the distribution range using species distribution models. Racomitrium lanuginosum comprises distinct molecular lineages sympatrically distributed in the Northern Hemisphere. We also uncovered long-distance dispersal from eastern North America to Scandinavia and potential in situ survival in northern Scandinavia. Due to the genetic signatures, the Alaska Peninsula could be considered a glacial refugium. The species experienced post-glacial expansion northwards in the Northern Hemisphere, mainly from the Alaska Peninsula. Our results exemplify the complex phylogeographic history in cold environments and contribute to recognizing evolutionary patterns in the Northern Hemisphere. [ABSTRACT FROM AUTHOR]

Codon usage bias (CUB) refers to an unequal occurrence of specific synonymous codons in a genome with variations within and among species. Findings suggest that CUB is significant in interpreting evolutionary trends of species and/or genes. However, information of CUB remained unclear for organisms in the family Potamididae. This study examined compositional features and codon usage patterns of coI gene among potamidids to identify factors shaping CUB. Quantification of CUB and the identification of causative factors were performed by assessing nucleotide composition, effective number of codons (ENC), relative synonymous codon usage (RSCU), parity rule (PR2) analysis and neutrality plot analysis as well as correspondence analysis (CoA). Results from neutrality plot suggested that natural selection is the dominant evolutionary factor for CUB in the coI gene among the different genera of potamidids. In contrast, CoA, PR2 and ENC analyses indicated the existence of mutational pressure with minor contribution. Nucleotide composition investigations indicated that most codons are AT rich; with T nucleotides count being highly presented. Thus, compositional restraint was also among the reasons behind the pattern of codon usage in coI gene of species in the potamidids. Although, the present study involved one mitochondrial gene, coI, the findings have demonstrated that various genera of potamidids are differentially adapted in the environment, and that coI gene may be differentially expressed and remain to be elucidated.

Pustular microbial mats in Shark Bay, Western Australia, are modern analogs of microbial systems that colonized peritidal environments before the evolution of complex life. To understand how these microbial communities evolved to grow and metabolize in the presence of various environmental stresses, the horizontal gene transfer (HGT) detection tool, MetaCHIP, was used to identify the horizontal transfer of genes related to stress response in 83 metagenome-assembled genomes from a Shark Bay pustular mat. Subsequently, maximum-likelihood phylogenies were constructed using these genes and their most closely related homologs from other environments in order to determine the likelihood of these HGT events occurring within the pustular mat. Phylogenies of several stress-related genes—including those involved in response to osmotic stress, oxidative stress and arsenic toxicity—indicate a potentially long history of HGT events and are consistent with these transfers occurring outside of modern pustular mats. The phylogeny of a particular osmoprotectant transport gene reveals relatively recent adaptations and suggests interactions between Planctomycetota and Myxococcota within these pustular mats. Overall, HGT phylogenies support a potentially broad distribution in the relative timing of the HGT events of stress-related genes and demonstrate ongoing microbial adaptations and evolution in these pustular mat communities.

Glutamate decarboxylase (GAD) pathway (GDP) is a major acid resistance mechanism enabling microorganisms’ survival in low pH environments. We aimed to study the molecular evolution and population genetics of GDP in Lactic Acid Bacteria (LAB) to understand evolutionary processes shaping adaptation to acidic environments comparing species where the GDP genes are organized in an operon structure (Levilactobacillus brevis) versus lack of an operon structure (Lactiplantibacillus plantarum). Within species molecular population genetic analyses of GDP genes in L. brevis and L. plantarum sampled from diverse fermented food and other environments showed abundant synonymous and non-synonymous nucleotide diversity, mostly driven by low frequency changes, distributed throughout the coding regions for all genes in both species. GAD genes showed higher level of replacement polymorphism compared to transporter genes (gadC and YjeM) for both species, and GAD genes that are outside of an operon structure showed even higher level of replacement polymorphism. Population genetic tests suggest negative selection against replacement changes in all genes. Molecular structure and amino acid characteristics analyses showed that in none of the GDP genes replacement changes alter 3D structure or charge distribution supporting negative selection against non-conservative amino acid changes. Phylogenetic and between species divergence analyses suggested adaptive protein evolution on GDP genes comparing phylogenetically distant species, but conservative evolution comparing closely related species. GDP genes within an operon structure showed slower molecular evolution and higher conservation. All GAD and transporter genes showed high codon usage bias in examined LAB species suggesting high expression and utilization of acid resistance genes. Substantial discordances between species, GAD, and transporter gene tree topologies were observed suggesting molecular evolution of GDP genes do not follow speciation events. Distribution of operon structure on the species tree suggested multiple independent gain or loss of operon structure in LABs. In conclusion, GDP genes in LABs exhibit a dynamic molecular evolutionary history shaped by gene loss, gene transfer, negative and positive selection to maintain its active role in acid resistance mechanism, and enable organisms to thrive in acidic environments.

Understanding the molecular underpinnings of the evolution of complex (multi-part) systems is a fundamental topic in biology. One unanswered question is to what the extent do similar or different genes and regulatory interactions underlie similar complex systems across species? Animal eyes and phototransduction (light detection) are outstanding systems to investigate this question because some of the genetics underlying these traits are well characterized in model organisms. However, comparative studies using non-model organisms are also necessary to understand the diversity and evolution of these traits. Here, we compare the characteristics of photoreceptor cells, opsins, and phototransduction cascades in diverse taxa, with a particular focus on cnidarians. In contrast to the common theme of deep homology, whereby similar traits develop mainly using homologous genes, comparisons of visual systems, especially in non-model organisms, are beginning to highlight a "deep diversity" of underlying components, illustrating how variation can underlie similar complex systems across taxa. Although using candidate genes from model organisms across diversity was a good starting point to understand the evolution of complex systems, unbiased genome-wide comparisons and subsequent functional validation will be necessary to uncover unique genes that comprise the complex systems of non-model groups to better understand biodiversity and its evolution. [ABSTRACT FROM AUTHOR]

Abstract Background The genus Tulotis has been classified into the genus Platanthera in the present taxonomic studies since the morphological characteristics of this genus is very similar to that of Platanthera. Platanthera ussuriensis, formerly named as Tulotis ussuriensis, is a small terrestrial orchid species and has been listed as wild plant under State protection (category II) in China. An improved understanding of the genomic information will enable future applications of conservation strategy as well as phylogenetic studies for this rare orchid species. The objective of this research was to characterize and compare the chloroplast genome of P. ussuriensis with other closely related species of Orchidaceae. Results The chloroplast genome sequence of P. ussuriensis is 155,016 bp in length, which included a pair of inverted repeats (IRs) of 26,548 bp that separated a large single copy (LSC) region of 83,984 bp and a small single copy (SSC) region of 17,936 bp. The annotation contained a total of 132 genes, including 86 protein-coding genes, 38 tRNA genes and 8 rRNA genes. The simple sequence repeat (SSR) analysis showed that there were 104 SSRs in the chloroplast genome of P. ussuriensis. RNA editing sites recognition indicated 72 RNA editing events occurred, and all codon changes were C to T conversions. Comparative genomics showed that the chloroplast sequence of Platanthera related species were relatively conserved, while there were still some high variation regions that could be used as molecular markers. Moreover, Platanthera related species showed similar IR/SSC and IR/LSC borders. The phylogenetic analysis suggested that P. ussuriensis had a closer evolutionary relationship with P. japonica followed by the remaining Platanthera species. Conclusion Orchidaceae is a key group of biodiversity protection and also a hot spot group in the plant taxonomy and evolution studies due to their characteristics of high specialization and rapid evolution. This research determined the complete chloroplast genome of P. ussuriensis for the first time, and compared the sequence with other closely related orchid species. These results provide a foundation for future genomic and molecular evolution of the Orchidaceae species, and provide insights into the development of conservation strategy for Platanthera species.

As a member of the forkhead box L gene family, foxl2 plays a significant role in gonadal development and the regulation of reproduction. During the evolution of deuterostome, whole genome duplication (WGD)-enriched lineage diversifications and regulation mechanisms occurs. However, only limited research exists on foxl2 duplication in teleost or other vertebrate species. In this study, two foxl2 paralogs, foxl2 and foxl2l, were identified in the transcriptome of spotted knifejaw (Oplegnathus punctatus), which had varying expressions in the gonads. The foxl2 was expressed higher in the ovary, while foxl2l was expressed higher in the testis. Phylogenetic reconstruction, synteny analysis, and the molecular evolution test confirmed that foxl2 and foxl2l likely originated from the first two WGD. The expression patterns test using qRT-PCR and ISH as well as motif scan analysis revealed evidence of potentially functional divergence between the foxl2 and foxl2l paralogs in spotted knifejaw. Our results indicate that foxl2 and foxl2l may originate from the first two WGD, be active in transcription, and have undergone functional divergence. These results shed new light on the evolutionary trajectories of foxl2 and foxl2l and highlights the need for further detailed functional analysis of these two duplicated paralogs.

Dkks have inhibitory effects on the Wnt signaling pathway, which is involved in the development of skin and its appendages and the regulation of hair growth. The nucleotide sequences were compared and analyzed to further investigate the relationship between the structure and function of the Dkk gene family and vertebrate epidermal hair. The analysis of the molecular evolution of the Dkk family revealed that the evolution rate of the genes changed significantly after speciation, with the Aves and Reptilia branches showing accelerated evolution. Additionally, positive selection was observed at specific sites. The tertiary structure of the protein was also predicted. The analysis of the functional divergence of the Dkk family revealed that the functional divergence coefficient of each gene was greater than 0, with most of the functional divergence sites were located in the Cys-2 domain and a few in the Cys-1 domain. This suggests that the amino acid and functional divergence sites may play a role in regulating the binding of the Dkk family to LRP5/6, and thus affect the inhibition of Wnt signaling, leading to different functions of Dkk1, Dkk2, and Dkk4 in the development of skin hair follicles. In addition, the Dkk families of Aves and Reptilia may have undergone adaptive evolution and functional divergence.

The XBB.1.16 SARS-CoV-2 variant, also known as Arcturus, is a recent descendant lineage of the recombinant XBB (nicknamed Gryphon). Compared to its direct progenitor, XBB.1, XBB.1.16 carries additional spike mutations in key antigenic sites, potentially conferring an ability to evade the immune response compared to other circulating lineages. In this context, we conducted a comprehensive genome-based survey to gain a detailed understanding of the evolution and potential dangers of the XBB.1.16 variant, which became dominant in late June. Genetic data indicates that the XBB.1.16 variant exhibits an evolutionary background with limited diversification, unlike dangerous lineages known for rapid changes. The evolutionary rate of XBB.1.16, which amounts to 3.95 × 10−4 subs/site/year, is slightly slower than that of its direct progenitors, XBB and XBB.1.5, which have been circulating for several months. A Bayesian Skyline Plot reconstruction suggests that the peak of genetic variability was reached in early May 2023, and currently, it is in a plateau phase with a viral population size similar to the levels observed in early March. Structural analyses indicate that, overall, the XBB.1.16 variant does not possess structural characteristics markedly different from those of the parent lineages, and the theoretical affinity for ACE2 does not seem to change among the compared variants. In conclusion, the genetic and structural analyses of SARS-CoV-2 XBB.1.16 do not provide evidence of its exceptional danger or high expansion capability. Detected differences with previous lineages are probably due to genetic drift, which allows the virus constant adaptability to the host, but they are not necessarily connected to a greater danger. Nevertheless, continuous genome-based monitoring is essential for a better understanding of its descendants and other lineages.

Thermal tolerance is a key determinant of species distribution. Despite much study, the genetic basis of adaptive evolution of thermal tolerance, including the relative contributions of transcriptional regulation versus protein evolution, remains unclear. Populations of the intertidal copepod Tigriopus californicus are adapted to local thermal regimes across their broad geographic range. Upon thermal stress, adults from a heat tolerant southern population, San Diego (SD), upregulate several heat shock proteins (HSPs) to higher levels than those from a less tolerant northern population, Santa Cruz (SC). Suppression of a specific HSP, HSPB1, significantly reduces T. californicus survival following acute heat stress. Sequencing of HSPB1 revealed population specific nucleotide substitutions in both promoter and coding regions of the gene. HSPB1 promoters from heat tolerant populations contain two canonical heat shock elements (HSEs), the binding sites for heat shock transcription factor (HSF), whereas less tolerant populations have mutations in these conserved motifs. Allele specific expression of HSPB1 in F1 hybrids between tolerant and less tolerant populations showed significantly biased expression favoring alleles from tolerant populations and supporting the adaptive divergence in these cis-regulatory variants. The functional impact of population-specific nonsynonymous substitutions in HSPB1 coding sequences was tested by assessing the thermal stabilization properties of SD versus SC HSPB1 protein variants. Recombinant HSPB1 from the southern SD population showed greater capacity for protecting protein structure under elevated temperature. Our results indicate that both regulatory and protein coding sequence evolution within a single gene appear to contribute to thermal tolerance phenotypes and local adaptation among conspecific populations.

Understanding genetic variation in human immunodeficiency virus (HIV) is clinically and immunologically important for patient treatment and vaccine development. We investigated the longitudinal intra-host genetic variation of HIV in over 3,000 individuals in the US National HIV Surveillance System with at least four reported HIV-1 polymerase (pol) sequences. In this population, we identified 149 putative instances of superinfection (i.e. an individual sequentially infected with genetically divergent, polyphyletic viruses). Unexpectedly, we discovered a group of 240 individuals with consecutively sampled viral strains that were >0.015 substitutions/site divergent, despite remaining monophyletic in the phylogeny. Viruses in some of these individuals had a maximum genetic divergence approaching that found between two random, unrelated HIV-1 subtype-B pol sequences within the US population. Individuals with these highly divergent viruses tended to be diagnosed nearly a decade earlier in the epidemic than people with superinfection or virus with less intra-host genetic variation, and they had distinct transmission risk factor profiles. To better understand this genetic variation in cases with extremely divergent, monophyletic viruses, we performed molecular clock phylogenetic analysis. Our findings suggest that, like Hepatitis C virus, extremely divergent HIV lineages can be maintained within an individual and reemerge over a period of years.

It has recently been described that the Japanese plum "Santa Rosa" bud sport series contains variations in ripening pattern: climacteric, suppressed-climacteric and non-climacteric types. This provides an interesting model to study the role of ethylene and other key mechanisms governing fruit ripening, softening and senescence. The aim of the current study was to investigate such differences at the genomic level, using this series of plum bud sports, with special reference to genes involved in ethylene biosynthesis, signal transduction, and sugar metabolism. Genomic DNA, isolated from leaf samples of six Japanese plum cultivars ("Santa Rosa", "July Santa Rosa", "Late Santa Rosa", "Sweet Miriam", "Roysum", and "Casselman"), was used to construct paired-end standard Illumina libraries. Sequences were aligned to the Prunus persica genome, and genomic variations (SNPs, INDELS, and CNV's) were investigated. Results determined 12 potential candidate genes with significant copy number variation (CNV), being associated with ethylene perception and signal transduction components. Additionally, the Maximum Likelihood (ML) phylogenetic tree showed two sorbitol dehydrogenase genes grouping into a distinct clade, indicating that this natural group is well-defined and presents high sequence identity among its members. In contrast, the ethylene group, which includes ACO1, ACS1, ACS4, ACS5, CTR1, ERF1, ERF3, and ethylene-receptor genes, was widely distributed and clustered into 10 different groups. Thus, ACS, ERF, and sorbitol dehydrogenase proteins potentially share a common ancestor for different plant genomes, while the expansion rate may be related to ancestral expansion rather than species-specific events. Based on the distribution of the clades, we suggest that gene function diversification for the ripening pathway occurred prior to family extension. We herein report all the frameshift mutations in genes involved in sugar transport and ethylene biosynthesis detected as well as the gene CNV implicated in ripening differences.

Interdisciplinary syntheses are needed to scale up discovery of the environmental drivers and molecular basis of adaptation in nature. Here we integrated novel approaches using whole genome sequences, satellite remote sensing, and transgenic experiments to study natural loss-of-function alleles associated with drought histories in wild Arabidopsis thaliana. The genes we identified exhibit population genetic signatures of parallel molecular evolution, selection for loss-of-function, and shared associations with flowering time phenotypes in directions consistent with longstanding adaptive hypotheses seven times more often than expected by chance. We then confirmed predicted phenotypes experimentally in transgenic knockout lines. These findings reveal the importance of drought timing to explain the evolution of alternative drought tolerance strategies and further challenge popular assumptions about the adaptive value of genetic loss-of-function in nature. These results also motivate improved species-wide sequencing efforts to better identify loss-of-function variants and inspire new opportunities for engineering climate resilience in crops.

Abstract Background Methicillin-resistant Staphylococcus aureus (MRSA) is a major nosocomial pathogen subdivided into lineages termed sequence types (STs). Since the 1950s, successive waves of STs have appeared and replaced previously dominant lineages. One such event has been occurring in China since 2013, with community-associated (CA-MRSA) strains including ST59 largely replacing the previously dominant healthcare-associated (HA-MRSA) ST239. We previously showed that ST59 isolates tend to have a competitive advantage in growth experiments against ST239. However, the underlying genomic and phenotypic drivers of this replacement event are unclear. Methods Here, we investigated the replacement of ST239 using whole-genome sequencing data from 204 ST239 and ST59 isolates collected in Chinese hospitals between 1994 and 2016. We reconstructed the evolutionary history of each ST and considered two non-mutually exclusive hypotheses for ST59 replacing ST239: antimicrobial resistance (AMR) profile and/or ability to colonise and persist in the environment through biofilm formation. We also investigated the differences in cytolytic activity, linked to higher virulence, between STs. We performed an association study using the presence and absence of accessory virulence genes. Results ST59 isolates carried fewer AMR genes than ST239 and showed no evidence of evolving towards higher AMR. Biofilm production was marginally higher in ST59 overall, though this effect was not consistent across sub-lineages so is unlikely to be a sole driver of replacement. Consistent with previous observations of higher virulence in CA-MRSA STs, we observed that ST59 isolates exhibit significantly higher cytolytic activity than ST239 isolates, despite carrying on average fewer putative virulence genes. Our association study identified the chemotaxis inhibitory protein (chp) as a strong candidate for involvement in the increased virulence potential of ST59. We experimentally validated the role of chp in increasing the virulence potential of ST59 by creating Δchp knockout mutants, confirming that ST59 can carry chp without a measurable impact on fitness. Conclusions Our results suggest that the ongoing replacement of ST239 by ST59 in China is not associated to higher AMR carriage or biofilm production. However, the increase in ST59 prevalence is concerning since it is linked to a higher potential for virulence, aided by the carriage of the chp gene.

Abstract The association between plants and their pollinators is essential for increasing the diversity in angiosperms. Morphological and physiological traits, mainly floral scent, can influence the pollination dynamics and select pollinators for each plant species. In this work, we studied two proteins involved in producing volatile organic compounds in plants, conyferyl alcohol acyltransferase (CFAT) and benzoyl-CoA:benzyl alcohol/phenyl ethanol benzoyl transferase (BPBT) genes. We aimed to understand these proteins with respect to evolutionary and structural aspects and functions in Solanaceae using phylogenetic methods and comparative molecular modeling. We used Bayesian inference to describe the proteins’ evolutionary history using Petunia x hybrida as a query to search for homologs in the Solanaceae family. Theoretical 3D models were obtained for both proteins using Panicum virgatum as a template. The phylogenetic tree included several different enzymes with diverse biological roles in Solanaceae, displaying the transferase domain. We identified only one sequence of CFAT in the databases, which belongs to Petunia x hybrida, and found several BPBT sequences from the genera Nicotiana, Solanum, and Capsicum. The 3D structures of CFAT and BPBT have two different domains, and we have identified the amino acid residues essential for the enzymatic activity and interaction with substrates.

Abstract Here we sequenced and characterized the complete plastome of Bactris riparia, a species closely related to B. gasipaes and widely distributed in Western Amazonia. We performed a comparative genomic analysis with B. riparia and the other four Bactridinae species retrieved from GenBank. The plastome of B. riparia was 156,715 bp with a quadripartite structure. Gene content included 86 protein-coding genes (CDS), 38 tRNAs, and 8 rRNAs. Bactris riparia has 69 more base pairs than B. gasipaes, with identical numbers in IR, and more in LSC and SSC. The comparative analysis indicated that structure, collinearity, and IR/SSC borders of plastomes within subtribe Bactridinae are, in general, conserved. We predicted 69 SSRs in B. riparia plastome. Among them, ~80% consisted of A/T homopolymers. Among the 52 variable CDS, rbcL showed the highest non-synonymous rate, while the rps15 gene had the highest synonymous rate. Three genes (ccsA, cemA, and rpoC1) presented evidence of positive selection and 22 genes showed evidence of purifying selection. The phylogenetic tree based on plastome sequences set Bactris as more closely related to Astrocaryum than to Acrocomia. These new plastome data of B. riparia will contribute to studies about the diversity, evolutionary history, and conservation of palms.

Plant phospholipase D (PLD) exerts important roles in various biological processes, such as intracellular signaling and morphological development. Our knowledge about early land plant PLDs is still underdeveloped. In this study, we identified 84 PLD genes in six mosses, i.e., Physcomitrella patens, Ceratodon purpureus, Fontinalis antipyretica, Pleurozium schreberi, Sphagnum magellanicum, and Sphagnum fallax. These PLDs were classified into four clades (I–IV). We showed that PLD underwent rapid expansion in mosses. A total of six conserved domains and two core HKD motifs were detected. Structure analysis uncovered that the moss PLDs from within a clade generally exhibited similar exon-intron organization. Cis-elements prediction and expression analyses indicated that P. patens PLDs had key roles in stress responsiveness and plant development. Particularly, about half of the P. patens PLDs (e.g., PpPLD1, PpPLD2, and PpPLD5) were differentially expressed under biotic and abiotic stresses. We also determined the expression pattern of P. patens PLD genes in various tissues and at different stages of development. Although the moss, clubmoss, liverwort, and fern PLDs evolved largely under functional constraints, we found episodic positive selection in the moss PLDs, e.g., C. purpureus PLD2 and P. patens PLD11. We infer that the evolutionary force acting on the PLDs may have facilitated moss colonization of land. Our work provides valuable insights into the diversification of moss PLD genes, and can be used for future studies of their functions.

Background: After domestication, the evolution of phenotypically-varied sheep breeds has generated rich biodiversity. This wide phenotypic variation arises as a result of hidden genomic changes that range from a single nucleotide to several thousands of nucleotides. Thus, it is of interest and significance to reveal and understand the genomic changes underlying the phenotypic variation of sheep breeds in order to drive selection towards economically important traits. Review: Various traits contribute to the emergence of variation in sheep phenotypic characteristics, including coat color, horns, tail, wool, ears, udder, vertebrae, among others. The genes that determine most of these phenotypic traits have been investigated, which has generated knowledge regarding the genetic determinism of several agriculturally-relevant traits in sheep. In this review, we discuss the genomic knowledge that has emerged in the past few decades regarding the phenotypic traits in sheep, and our ultimate aim is to encourage its practical application in sheep breeding. In addition, in order to expand the current understanding of the sheep genome, we shed light on research gaps that require further investigation. Conclusions: Although significant research efforts have been conducted in the past few decades, several aspects of the sheep genome remain unexplored. For the full utilization of the current knowledge of the sheep genome, a wide practical application is still required in order to boost sheep productive performance and contribute to the generation of improved sheep breeds. The accumulated knowledge on the sheep genome will help advance and strengthen sheep breeding programs to face future challenges in the sector, such as climate change, global human population growth, and the increasing demand for products of animal origin. [ABSTRACT FROM AUTHOR]

Simple Summary: Infectious diseases place an economic burden on aquaculture and a limitation to its growth. This state-of-the-art review describes the application of genetics and genomics as novel tools to control infectious disease in aquaculture. Infectious diseases place an economic burden on aquaculture and a limitation to its growth. An innovative approach to mitigate their impact on production is breeding for disease resistance: selection for domestication, family-based selection, marker-assisted selection, and more recently, genomic selection. Advances in genetics and genomics approaches to the control of infectious diseases are key to increasing aquaculture efficiency, profitability, and sustainability and to reducing its environmental footprint. Interaction and co-evolution between a host and pathogen can, however, turn breeding to boost infectious disease resistance into a potential driver of pathogenic change. Parallel molecular characterization of the pathogen and its virulence and antimicrobial resistance genes is therefore essential to understand pathogen evolution over time in response to host immunity, and to apply appropriate mitigation strategies. [ABSTRACT FROM AUTHOR]

Lettuce mosaic virus (LMV) is one of the most important and destructive members of Potyviruses in the family Potyviridae that is transmitted by aphid and seed which cause economic damage and serious yield losses to different crops around the world. To explore phylogenetic relationships, population evolution, and the effect of selection forces on the complete coat protein (CP) gene of LMV isolates, 36 complete coat protein gene sequences retrieved from the GenBank database in several different countries in four continents: Central and East Asia (China, Taiwan, South Korea, and Turkey), Europe (France), Africa (Tunisia), and America (Chile and Brazil). The phylogenetic tree of LMV isolates was grouped into independent clades with the significance of FST values (>0.27). The ratio of dN/dS is calculated less than one and showed that the LMV-CP has been under negative selection. Statistical tests (Tajima’s D, Fu and Li’s D* and F*) were used to estimate non-significantly negative values for all clades and geographic populations except for the Brazilian population and Clade III in a phylogenetic group. The negatives values revealed that there is less polymorphism estimation. All geographic populations in four phylogroups of LMV seem to be at equilibrium because all neutrality test statistics were non-significant. The findings suggested that the dynamics of LMV molecular evolution may be dependent on mutation, recombination, and negative selection. No recombination events were observed in this part of the LMV genome. Therefore, this study provides the first time evolution and differentiation between populations of LMV isolates from around the world and suggested that they may be occurred by the transmission of the virus among lettuce plants by types of aphid species, migration in different geographical areas by infected seeds and plant material, and broad host range which seems that these events have played an important role in shaping the LMV population structure.

Abstract Background In plants, Calcium (Ca2+) acts as a universal messenger in various signal transduction pathways, including responses to biotic and abiotic stresses and regulation of cellular and developmental processes. The Ca2+/cation antiporter (CaCA) superfamily proteins play vital roles in the transport of Ca2+ and/or other cations. However, the characteristics of these superfamily members in Saccharum and their evolutionary and functional implications have remained unclear. Results A total of 34 CaCA genes in Saccharum spontaneum, 5 CaCA genes in Saccharum spp. R570, and 14 CaCA genes in Sorghum bicolor were identified and characterized. These genes consisted of the H+/cation exchanger (CAX), cation/Ca2+ exchanger (CCX), EF-hand / CAX (EFCAX), and Mg2+/H+ exchanger (MHX) families, among which the CCX and EFCAX could be classified into three groups while the CAX could be divided into two groups. The exon/intron structures and motif compositions suggested that the members in the same group were highly conserved. Synteny analysis of CaCAs established their orthologous and paralogous relationships among the superfamily in S. spontaneum, R570, and S. bicolor. The results of protein-protein interactions indicated that these CaCA proteins had direct or indirect interactions. Quantitative reverse transcription polymerase chain reaction (qRT-PCR) analysis demonstrated that most members of Saccharum CaCA genes exhibited a similar expression pattern in response to hormonal (abscisic acid, ABA) treatment but played various roles in response to biotic (Sporisorium scitamineum) and abiotic (cold) stresses. Furthermore, ScCAX4, a gene encoding a cytoplasm, plasma membrane and nucleus positioning protein, was isolated from sugarcane. This gene was constitutively expressed in different sugarcane tissues and its expression was only induced at 3 and 6 h time points after ABA treatment, however was inhibited and indued in the whole process under cold and S. scitamineum stresses, respectively. Conclusions This study systematically conducted comparative analyses of CaCA superfamily genes among S. spontaneum, R570, and S. bicolor, delineating their sequence and structure characteristics, classification, evolutionary history, and putative functions. These results not only provided rich gene resources for exploring the molecular mechanism of the CaCA superfamily genes but also offered guidance and reference for research on other gene families in Saccharum.

Oil-tea (Camellia oleifera) is a woody oil crop whose nectar includes galactose derivatives that are toxic to honey bees. Interestingly, some mining bees of the genus Andrena can entirely live on the nectar (and pollen) of oil-tea and are able to metabolize these galactose derivatives. We present the first next-generation genomes for five and one Andrena species that are, respectively, specialized and non-specialized oil-tea pollinators and, combining these with the published genomes of six other Andrena species which did not visit oil-tea, we performed molecular evolution analyses on the genes involved in the metabolizing of galactose derivatives. The six genes (NAGA, NAGA-like, galM, galK, galT, and galE) involved in galactose derivatives metabolism were identified in the five oil-tea specialized species, but only five (with the exception of NAGA-like) were discovered in the other Andrena species. Molecular evolution analyses revealed that NAGA-like, galK, and galT in oil-tea specialized species appeared under positive selection. RNASeq analyses showed that NAGA-like, galK, and galT were significantly up-regulated in the specialized pollinator Andrena camellia compared to the non-specialized pollinator Andrena chekiangensis. Our study demonstrated that the genes NAGA-like, galK, and galT have played an important role in the evolutionary adaptation of the oil-tea specialized Andrena species.