Strong Functional Constraints Research Articles

Ultrastructure and phylogenetic significance of spermatozoa in Nemertea

AbstractSperm ultrastructure is frequently employed as a source for phylogenetic inference due to the ease of accessing spermatozoa. Despite being unicellular, sperm cells exhibit a relatively high number of diverse characters and character states. Spermatozoa are subject to strong sexual selection as they are finely tuned for maximizing male reproductive success. Given this strong functional constraint, one might anticipate the emergence of convergent characters in line with similar modes of reproduction. As a result, it might be expected that sperm cells do not possess substantial phylogenetic signal, with functional constraints overshadowing any evolutionary heritage. To test this assumption, we conducted a study on sperm ultrastructure in 11 nemertean species, representing closely related groups and major nemertean lineages. We analyzed these data for their phylogenetic signal within the context of the most recent nemertean phylogenies. Our findings demonstrate that, at all systematic levels, functional constraints imposed by sexual selection, or the mode of reproduction do not supersede the influence of historical constraints on sperm ultrastructure.

Evolution and functional role prediction of the CYP6DE and CYP6DJ subfamilies in Dendroctonus (Curculionidae: Scolytinae) bark beetles.

Dendroctonus-bark beetles are natural components and key ecological agents of coniferous forests. They spend most of their lives under the bark, where they are exposed to highly toxic terpenes present in the oleoresin. Cytochrome P450 (CYP) is a multigene family involved in the detoxification of these compounds. It has been demonstrated that CYP6DE and CYP6DJ subfamilies hydroxylate monoterpenes, whose derivatives can act as pheromone synergist compounds or be pheromones themselves in these insects. Given the diversity and functional role of CYPs, we investigated whether these cytochromes have retained their function throughout the evolution of these insects. To test this hypothesis, we performed a Bayesian phylogenetic analysis to determine phylogenetic subgroups of cytochromes in these subfamilies. Subgroups were mapped and reconciled with the Dendroctonus phylogeny. Molecular docking analyses were performed with the cytochromes of each subgroup and enantiomers of α-pinene and β-pinene, (+)-3-carene, β-myrcene and R-(+)-limonene. In addition, functional divergence analysis was performed to identify critical amino acid sites that influence changes in catalytic site conformation and/or protein folding. Three and two phylogenetic subgroups were recovered for the CYP6DE and CYP6DJ subfamilies, respectively. Mapping and reconciliation analysis showed different gain and loss patterns for cytochromes of each subgroup. Functional predictions indicated that the cytochromes analyzed are able to hydroxylate all monoterpenes; however, they showed preferential affinities to different monoterpenes. Functional divergence analyses indicated that the CYP6DE subfamily has experimented type I and II divergence, whereas the CYP6DJ subfamily has evolved under strong functional constraints. Results suggest cytochromes of the CYP6DE subfamily evolve to reinforce their detoxifying capacity hydroxylating mainly α- and β-pinene to (+) and (-)-trans-verbenol, being the negative enantiomer used as a pheromone by several Dendroctonus species; whereas cytochromes of the CYP6DJ subfamily appear to retain their original function related to the detoxification of these compounds.

A Comprehensive Analysis of the Fowleria variegata (Valenciennes, 1832) Mitochondrial Genome and Its Phylogenetic Implications within the Family Apogonidae.

Controversies surrounding the phylogenetic relationships within the family Apogonidae have persisted due to the limited molecular data, obscuring the evolution of these diverse tropical marine fishes. This study presents the first complete mitochondrial genome of Fowleria variegata, a previously unrepresented genus, using high-throughput Illumina sequencing. Through a comparative mitogenomic analysis, F. variegate was shown to exhibit a typical genome architecture and composition, including 13 protein-coding, 22 tRNA and 2 rRNA genes and a control region, consistent with studies of other Apogonidae species. Nearly all protein-coding genes started with ATG, while stop codons TAA/TAG/T were observed, along with evidence of strong functional constraints imposed via purifying selection. Phylogenetic reconstruction based on maximum likelihood and Bayesian approaches provided robust evidence that F. variegata forms a basal lineage closely related to P. trimaculatus within Apogonidae, offering novel perspectives into the molecular evolution of this family. By generating new mitogenomic resources and evolutionary insights, this study makes important headway in elucidating the phylogenetic relationships and mitogenomic characteristics of Apogonidae fishes. The findings provide critical groundwork for future investigations into the drivers of diversification, speciation patterns, and adaptive radiation underlying the extensive ecological diversity and biological success of these marine fishes using phylogenomics and population genomics approaches.

Complete Mitochondrial Genome Sequence and Phylogenetic Analysis of Procambarus clarkii and Cambaroides dauricus from China.

To enhance the management and protection of crayfish genetic diversity and germplasm resources in Cambaroides dauricus (C. dauricus), a common species of Procambarus clarkii (P. clarkii) was used as a control group to compare the whole mitochondrial genome sequence using Illumina sequencing technology. This study found that the mitochondrial genome of C. dauricus is 15580 bp in length, with a base composition of A (31.84%), G (17.66%), C (9.42%), and T (41.08%) and a C + G content of 27.08%. The C + G in the D-loop is rich in 17.06%, indicating a significant preference. The mitochondrial genome of C. dauricus contains 13 protein-coding genes, 22 tRNA genes, and 2 rRNA genes, with most of the genes labeled in the negative direction, except for a few genes that are labeled in the positive direction. The start codons of the ten coding sequences are ATG, and the quintessential TAA and TAG are the stop codons. This study also found that the Ka/Ks ratios of most protein-coding genes in the mitochondria of both shrimps are lower than 1, indicating weak natural selection, except for nad 2, nad 5, and cox 1. The Ka/Ks ratio of cox 3 is the lowest (less than 0.1), indicating that this protein-coding gene bears strong natural selection pressure and functional constraint in the process of mitochondrial genetic evolution of both shrimps. Furthermore, we constructed phylogenetic analyses based on the entire sequence, which effectively distinguishes the high body from other shrimp species of the genus based on the mitochondrial genome. This study provides molecular genetic data for the diversity investigation and protection of fishery resources with Chinese characteristics and a scientific reference for the evolutionary study of Procambarus.

DN/dS-H, a New Test to Distinguish Different Selection Modes in Protein Evolution and Cancer Evolution.

One of the most popular measures in the analysis of protein sequence evolution is the ratio of nonsynonymous distance (dN) to synonymous distance (dS). Under the assumption that synonymous substitutions in the coding region are selectively neutral, the dN/dS ratio can be used to statistically detect the adaptive evolution (or purifying selection) if dN/dS > 1 (or dN/dS < 1) significantly. However, due to strong structural constraints and/or variable functional constraints imposed on amino acid sites, most encoding genes in most species have demonstrated dN/dS < 1. Consequently, the statistical power for testing dN/dS = 1 may be insufficient to distinguish between different selection modes. In this paper, we propose a more powerful test, called dN/dS-H, in which a new parameter H, a relative measure of rate variation among sites, was introduced. Given the condition of strong purifying selections at some sites, the dN/dS-H model predicts dN/dS = 1-H for neutral evolution, dN/dS < 1-H for nearly neutral selection, and dN/dS > 1-H for adaptive evolution. The potential of this new method for resolving the neutral-adaptive debates is illustrated by the protein sequence evolution in vertebrates, Drosophila and yeasts, as well as somatic cancer evolution (specialized as the CN/CS-H test).

Fossorial adaptations in African mole-rats (Bathyergidae) and the unique appendicular phenotype of naked mole-rats

Life underground has constrained the evolution of subterranean mammals to maximize digging performance. However, the mechanisms modulating morphological change and development of fossorial adaptations in such taxa are still poorly known. We assessed the morpho-functional diversity and early postnatal development of fossorial adaptations (bone superstructures) in the appendicular system of the African mole-rats (Bathyergidae), a highly specialized subterranean rodent family. Although bathyergids can use claws or incisors for digging, all genera presented highly specialized bone superstructures associated with scratch-digging behavior. Surprisingly, Heterocephalus glaber differed substantially from other bathyergids, and from fossorial mammals by possessing a less specialized humerus, tibia and fibula. Our data suggest strong functional and developmental constraints driving the selection of limb specializations in most bathyergids, but more relaxed pressures acting on the limbs of H. glaber. A combination of historical, developmental and ecological factors in Heterocephalus are hypothesized to have played important roles in shaping its appendicular phenotype.

Adaptation to graviportality in Rhinocerotoidea? An investigation through the long bone shape variation in their hindlimb

AbstractWeight support is a strong functional constraint modelling limb bones in heavy quadrupeds. However, the complex relations between bone shape, mass, size and body proportions have been poorly explored. Rhinocerotoidea is one of the groups showing the highest body mass reached by terrestrial mammals through time. Here, we explore the evolutionary variation of shape in hindlimb stylopod and zeugopod bones and its relationship with mass, size and gracility in this superfamily. Our results show that bones undergo a general increase in robustness towards high masses, associated with reinforcements of the main muscle insertions. The shape of the femur, carrying a marked phylogenetic signal, varies conjointly with mass, size and gracility, whereas that of the tibia appears related to gracility and mass only. The shape of the fibula does not vary according to that of the tibia. Moreover, congruent variation of shape between the distal part of the femur and the complete tibia underlines the potentially strong covariation of the elements constituting the knee joint. These results, coupled with those previously obtained from forelimb study, allow a better comprehension of the relationship between bone shape and mass among Rhinocerotoidea, and a refining of the concept of ‘graviportality’ in this superfamily.

Shared evolutionary trajectories of three independent neo-sex chromosomes in Drosophila.

Dosage compensation (DC) on the X Chromosome counteracts the deleterious effects of gene loss on the Y Chromosome. However, DC is not efficient if the X Chromosome also degenerates. This indeed occurs in Drosophila miranda, in which both the neo-Y and the neo-X are under accelerated pseudogenization. To examine the generality of this pattern, we investigated the evolution of two additional neo-sex chromosomes that emerged independently in D. albomicans and D. americana and reanalyzed neo-sex chromosome evolution in D. miranda. Comparative genomic and transcriptomic analyses revealed that the pseudogenization rate on the neo-X is also accelerated in D. albomicans and D. americana although to a lesser extent than in D. miranda. In males, neo-X-linked genes whose neo-Y-linked homologs are pseudogenized tended to be up-regulated more than those whose neo-Y-linked homologs remain functional. Moreover, genes under strong functional constraint and genes highly expressed in the testis tended to remain functional on the neo-X and neo-Y, respectively. Focusing on the D. miranda and D. albomicans neo-sex chromosomes that emerged independently from the same autosome, we further found that the same genes tend to become pseudogenized in parallel on the neo-Y. These genes include Idgf6 and JhI-26, which may be unnecessary or even harmful in males. Our results indicate that neo-sex chromosomes in Drosophila share a common evolutionary trajectory after their emergence, which may prevent sex chromosomes from being an evolutionary dead end.

Phenotypic regionalization of the vertebral column in the thorny skate Amblyraja radiata: Stability and variation.

Regionalization of the vertebral column occurred early during vertebrate evolution and has been extensively investigated in mammals. However, less data are available on vertebral regions of crown gnathostomes. This is particularly true for batoids (skates, sawfishes, guitarfishes, and rays) whose vertebral column has long been considered to be composed of the same two regions as in teleost fishes despite the presence of a synarcual. However, the numerous vertebral units in chondrichthyans may display a more complex regionalization pattern than previously assumed and the intraspecific variation of such pattern deserves a thorough investigation. In this study, we use micro-computed tomography (µCT) scans of vertebral columns of a growth series of thorny skates Amblyraja radiata to provide the first fine-scale morphological description of vertebral units in a batoids species. We further investigate axial regionalization using a replicable clustering analysis on presence/absence of vertebral elements to decipher the regionalization of the vertebral column of A. radiata. We identify four vertebral regions in this species. The two anteriormost regions, named synarcual and thoracic, may undergo strong developmental or functional constraints because they display stable patterns of shapes and numbers of vertebral units across all growth stages. The third region, named hemal transitional, is characterized by high inter-individual morphological variation and displays a transition between the monospondylous (one centrum per somite) to diplospondylous (two centra per somite) conditions. The posteriormost region, named caudal, is subdivided into three sub-regions with shapes changing gradually along the anteroposterior axis. These regionalized patterns are discussed in light of ecological habits of skates.

Documentation‐based functional constraint generation for library methods

SummaryAlthough software libraries promote code reuse and facilitate software development, they increase the complexity of programme analysis tasks. To effectively analyse programmes built on top of software libraries, it is essential to have specifications for the library methods that can be easily processed by analysis tools. However, the availability of such specifications is seriously limited at the moment. Manually writing the specifications can be prohibitively expensive and error‐prone, while existing automated approaches to inferring the specifications seldom produce results that are strong enough to be used in programme analysis. In this work, we propose the Doc2smt approach to generating strong functional constraints in SMT for library methods based on their documentations. Doc2smt first applies natural language processing (NLP) techniques and a set of rules to translate a method's natural language documentation into a large number of candidate constraint clauses in OCL. Then, it utilizes a manually enhanced domain model to identify OCL candidate constraint clauses that comply with the problem domain in static validation, translates well‐formed OCL constraints into the SMT‐LIB format, and checks whether each SMB‐LIB constraint rightly abstracts the functionalities of the method under consideration via testing in dynamic validation. In the end, it reports the first functional constraint that survives both validations to the user as the result. We have implemented the approach into a supporting tool with the same name. In experiments conducted on 451 methods from the Java Collections Framework and the Java IO library, Doc2smt generated correct constraints for 309 methods, with the average generation time for each correct constraint being merely 2.7 min. We have also applied the generated constraints to facilitate symbolic‐execution‐based test generation with the Symbolic Java PathFinder (SPF) tool. For 24 utility methods manipulating Java container and IO objects, SPF with access to the generated constraints produced 51.2 times more test cases than SPF without the access.

Ultra-conserved sequences in the genomes of highly diverse Anopheles mosquitoes, with implications for malaria vector control.

DNA sequences that are exactly conserved over long evolutionary time scales have been observed in a variety of taxa. Such sequences are likely under strong functional constraint and they have been useful in the field of comparative genomics for identifying genome regions with regulatory function. A potential new application for these ultra-conserved elements (UCEs) has emerged in the development of gene drives to control mosquito populations. Many gene drives work by recognizing and inserting at a specific target sequence in the genome, often imposing a reproductive load as a consequence. They can therefore select for target sequence variants that provide resistance to the drive. Focusing on highly conserved, highly constrained sequences lowers the probability that variant, gene drive-resistant alleles can be tolerated. Here, we search for conserved sequences of 18 bp and over in an alignment of 21 Anopheles genomes, spanning an evolutionary timescale of 100 million years, and characterize the resulting sequences according to their location and function. Over 8000 UCEs were found across the alignment, with a maximum length of 164 bp. Length-corrected gene ontology analysis revealed that genes containing Anopheles UCEs were over-represented in categories with structural or nucleotide-binding functions. Known insect transcription factor binding sites were found in 48% of intergenic Anopheles UCEs. When we looked at the genome sequences of 1142 wild-caught mosquitoes, we found that 15% of the Anopheles UCEs contained no polymorphisms. Our list of Anopheles UCEs should provide a valuable starting point for the selection and testing of new targets for gene-drive modification in the mosquitoes that transmit malaria.

Major antigenic site B of human influenza H3N2 viruses has an evolving local fitness landscape

Antigenic drift of influenza virus hemagglutinin (HA) is enabled by facile evolvability. However, HA antigenic site B, which has become immunodominant in recent human H3N2 influenza viruses, is also evolutionarily constrained by its involvement in receptor binding. Here, we employ deep mutational scanning to probe the local fitness landscape of HA antigenic site B in six different human H3N2 strains spanning from 1968 to 2016. We observe that the fitness landscape of HA antigenic site B can be very different between strains. Sequence variants that exhibit high fitness in one strain can be deleterious in another, indicating that the evolutionary constraints of antigenic site B have changed over time. Structural analysis suggests that the local fitness landscape of antigenic site B can be reshaped by natural mutations via modulation of the receptor-binding mode. Overall, these findings elucidate how influenza virus continues to explore new antigenic space despite strong functional constraints.

Towards an integrative understanding of stamen–corolla tube modifications and floral architecture in Boraginaceae s.s. (Boraginales)

AbstractMorphological studies rarely address floral organ modifications or their integration into floral architecture and floral function. Boraginaceae show two prominent types of stamen–corolla tube modifications: faucal and basal scales. Both types, especially faucal scales, are widely used in classification. Here, the ontogeny and morphology of faucal and basal scales are studied in 29 species from eight tribes and all three subfamilies of Boraginaceae s.s. (=Boraginoideae sensu APG IV) using scanning electron microscopy. Integration into floral architecture is visualized with micro-computed tomography (µCT). Faucal and basal scales are present in 18 and 27 species, respectively. Both types of scales develop late in flower ontogeny, but with variable timing. Faucal scales are morphologically far more variable than basal scales. Faucal scales are located close to the anthers and sometimes are involved in anther cone formation. Basal scales cover either the gynoecial disc nectary or the entire ovary. The different scale morphologies identified here enclose complex internal spaces, but they show no obvious phylogenetic patterns. This probably indicates strong functional constraints and adaptive pressures. This is a first in-depth study of the morphology and development of stamen–corolla tube modifications in Boraginaceae s.s., demonstrating that in situ three-dimensional visualization of floral architecture with µCT provides unprecedented insights into the evolution and functional integration of stamen–corolla tube modifications in Boraginales.

An Update on the Evolution of Glucosyltransferase (Gtf) Genes in Streptococcus.

In many caries-promoting Streptococcus species, glucosyltransferases (Gtfs) are recognized as key enzymes contributing to the modification of biofilm structures, disruption of homeostasis of healthy microbiota community and induction of caries development. It is therefore of great interest to investigate how Gtf genes have evolved in Streptococcus. In this study, we conducted a comprehensive survey of Gtf genes among 872 streptococci genomes of 37 species and identified Gtf genes from 364 genomes of 18 species. To clarify the relationships of these Gtf genes, 45 representative sequences were used for phylogenic analysis, which revealed two clear clades. Clade I included 12 Gtf genes from nine caries-promoting species of the Mutans and Downei groups, which produce enzymes known to synthesize sticky, water-insoluble glucans (WIG) that are critical for modifying biofilm structures. Clade II primarily contained Gtf genes responsible for synthesizing water-soluble glucans (WSG) from all 18 species, and this clade further diverged into three subclades (IIA, IIB, and IIC). An analysis of 16 pairs of duplicated Gtf genes revealed high divergence levels at the C-terminal repeat regions, with ratios of the non-synonymous substitution rate (dN) to synonymous substitution rate (dS) ranging from 0.60 to 1.03, indicating an overall relaxed constraint in this region. However, among the clade I Gtf genes, some individual repeat units possessed strong functional constraints by the same criterion. Structural variations in the repeat regions were also observed, with detection of deletions or recent duplications of individual repeat units. Overall, by establishing an updated phylogeny and further elucidating their evolutionary patterns, this work enabled us to gain a greater understanding of the origination and divergence of Gtf genes in Streptococcus.

Genetic conflicts with Plasmodium parasites and functional constraints shape the evolution of erythrocyte cytoskeletal proteins

Plasmodium parasites exerted a strong selective pressure on primate genomes and mutations in genes encoding erythrocyte cytoskeleton proteins (ECP) determine protective effects against Plasmodium infection/pathogenesis. We thus hypothesized that ECP-encoding genes have evolved in response to Plasmodium-driven selection. We analyzed the evolutionary history of 15 ECP-encoding genes in primates, as well as of their Plasmodium-encoded ligands (KAHRP, MESA and EMP3). Results indicated that EPB42, SLC4A1, and SPTA1 evolved under pervasive positive selection and that episodes of positive selection tended to occur more frequently in primate species that host a larger number of Plasmodium parasites. Conversely, several genes, including ANK1 and SPTB, displayed extensive signatures of purifying selection in primate phylogenies, Homininae lineages, and human populations, suggesting strong functional constraints. Analysis of Plasmodium genes indicated adaptive evolution in MESA and KAHRP; in the latter, different positively selected sites were located in the spectrin-binding domains. Because most of the positively selected sites in alpha-spectrin localized to the domains involved in the interaction with KAHRP, we suggest that the two proteins are engaged in an arms-race scenario. This observation is relevant because KAHRP is essential for the formation of “knobs”, which represent a major virulence determinant for P. falciparum.

In-silico approach to determine the possibility of a recombinant broad-spectrum vaccine for bluetongue disease

Bluetongue (BT) disease caused by BT virus (BTV) is an acute haemorrhagic disease of domesticated and wild ruminants. 29 BTV serotypes have been identified so far throughout the world with 24 circulating in India conferring no cross protective immunity. Hence, the objective of the current study was to investigate the possibility of developing a recombinant subunit protein based vaccine for BT that is cross protective to all the known serotypes of BTV. VP2 is a major serotype determining protein and vaccination with VP2 was known to induce neutralizing antibody response. Since BTV circulates as 29 immunologically distinct serotypes, the current study aimed to identify conserved region of VP2. The amino acid sequence alignment showed that although VP2 was very variable among serotypes, several regions were relatively more conserved between serotypes; most conserved region (cVP2) of which was evident from 338-383aa. Moreover, selection pressure analysis revealed that all the codons in cVP2 region were under strong negative selection (dN/dS=0.16) suggesting strong structural and functional constraints on cVP2. In addition, two MHCI, MHCII each, and three B cell epitopes were predicted based on the physicochemical properties of the amino acids in cVP2 region. This suggests strong immunogenicity of cVP2 making it a suitable vaccine candidate. Taken together, cVP2 appears to be a promising target for a broad-spectrum subunit vaccine moiety for BT that can confer protection against the available serotypes of BTV. Since all the structural proteins of BTV will not be used in such subunit vaccines, observations from current study can also be extended to develop a marker vaccine for DIVA against BT.

Low levels of genetic diversity associated with evidence of negative selection on the Babesia bovis apical membrane antigen 1 from parasite populations in Thailand

Babesia bovis, a parasite infecting cattle and buffalo, continues to spread throughout the developing world. The babesial vaccine was developed to be a sustainable alternative treatment to control the parasite. However, genetic diversity is a major obstacle for designing and developing a safe and effective vaccine. The apical membrane antigen 1 (AMA-1) is considered to be a potential vaccine candidate antigen among immunogenic genes of B. bovis. To gain a more comprehensive understanding of B. bovis AMA-1 (BbAMA-1), three B. bovis DNA samples were randomly selected to characterize in order to explore genetic diversity and natural selection and to predict the antigen epitopes. The sequence analysis revealed that BbAMA-1 has a low level of polymorphism and is highly conserved (95.46–99.94%) among Thai and global isolates. The majority of the polymorphic sites were observed in domains I and III. Conversely, domain II contained no polymorphic sites. We report the first evidence of strong negative or purifying selection across the full length of the gene, especially in domain I, by demonstrating a significant excess of the average number of synonymous (dS) over the non-synonymous (dN) substitutions. Finally, we also predict the linear and conformational B-cell epitope. The predicted B-cell epitopes appeared to be involved with the amino acid changes. Collectively, the results suggest that the conserved BbAMA-1 may be used to detect regional differences in the B. bovis parasite. Importantly, the limitation of BbAMA-1 diversity under strong negative selection indicates strong functional constraints on this gene. Thus, the gene could be a valuable target vaccine candidate antigen.

HIV Progression Depends on Codon and Amino Acid Usage Profile of Envelope Protein and Associated Host-Genetic Influence.

Acquired immune deficiency syndrome (AIDS) is a spectrum of conditions caused by infection with the human immunodeficiency virus (HIV). Two types of HIV have been characterized: HIV-1 and HIV-2. The present study investigated whether evolutionary selection pressure differs between rapid progressor (RP), slow progressor (SP), and long-term non-progressor (LTNP) of HIV-I infected individuals. An unexpected association between the evolutionary rate of substitution in envelope (env) gene and disease progression is observed. Our present study suggests that env genes of LTNP are subject to unusually strong functional constraint with respect to RP. We also observed that the three categories of env genes i.e., RP, SP, and LTNP, had their own characteristic pattern of amino acid usage and SP and LTNP sequences shared similar patterns of amino acid usage different from RP sequences and evolutionary rate significantly influenced the amino acid usage pattern of the three different types of env gene sequences. It was also noted that the evolutionary rate for the glycosylation sites of LTNP and SP sequences were even significantly less than the RP sequences. Comparative analysis on the influence of human host on the three categories of env genes are well correlated with the rates of disease progression suggesting the adaptive strategies of the viruses for successful residence and infection. Host associated selective constraints appeared most relaxed on the RP sequences and strongest in LTNP sequences. The present study clearly portrays how evolutionary selection pressure differs between three categories of env genes i.e., RP, SP, and LTNP. The env genes, coding for the env glycoproteins, experience severe selection constraints from the host due to their constant exposure to the host immune system. In this perspective it might be suggested that env gene evolution occurs mainly by negative selection with the occurrence of mutation that might not reach fixation in the viral population. This work also confers a deeper insight into the crucial effects of host factors that govern the overall progression of HIV infection.

Nucleic and Amino Acid Sequences Support Structure-Based Viral Classification.

ABSTRACTViral capsids ensure viral genome integrity by protecting the enclosed nucleic acids. Interactions between the genome and capsid and between individual capsid proteins (i.e., capsid architecture) are intimate and are expected to be characterized by strong evolutionary conservation. For this reason, a capsid structure-based viral classification has been proposed as a way to bring order to the viral universe. The seeming lack of sufficient sequence similarity to reproduce this classification has made it difficult to reject structural convergence as the basis for the classification. We reinvestigate whether the structure-based classification for viral coat proteins making icosahedral virus capsids is in fact supported by previously undetected sequence similarity. Since codon choices can influence nascent protein folding cotranslationally, we searched for both amino acid and nucleotide sequence similarity. To demonstrate the sensitivity of the approach, we identify a candidate gene for the pandoravirus capsid protein. We show that the structure-based classification is strongly supported by amino acid and also nucleotide sequence similarities, suggesting that the similarities are due to common descent. The correspondence between structure-based and sequence-based analyses of the same proteins shown here allow them to be used in future analyses of the relationship between linear sequence information and macromolecular function, as well as between linear sequence and protein folds.IMPORTANCE Viral capsids protect nucleic acid genomes, which in turn encode capsid proteins. This tight coupling of protein shell and nucleic acids, together with strong functional constraints on capsid protein folding and architecture, leads to the hypothesis that capsid protein-coding nucleotide sequences may retain signatures of ancient viral evolution. We have been able to show that this is indeed the case, using the major capsid proteins of viruses forming icosahedral capsids. Importantly, we detected similarity at the nucleotide level between capsid protein-coding regions from viruses infecting cells belonging to all three domains of life, reproducing a previously established structure-based classification of icosahedral viral capsids.

Evolution and expression of the fructokinase gene family in Saccharum

BackgroundSugarcane is an important sugar crop contributing up to about 80% of the world sugar production. Efforts to characterize the genes involved in sugar metabolism at the molecular level are growing since increasing sugar content is a major goal in the breeding of new sugarcane varieties. Fructokinases (FRK) are the main fructose phosphorylating enzymes with high substrate specificity and affinity.ResultsIn this study, by combining comparative genomics approaches with BAC resources, seven fructokinase genes were identified in S. spontaneum. Phylogenetic analysis based on representative monocotyledon and dicotyledon plant species suggested that the FRK gene family is ancient and its evolutionary history can be traced in duplicated descending order: SsFRK4, SsFRK6/SsFRK7,SsFRK5, SsFRK3 and SsFRK1/SsFRK2. Among the close orthologs, the number and position of exons in FRKs were conserved; in contrast, the size of introns varied among the paralogous FRKs in Saccharum. Genomic constraints were analyzed within the gene alleles and between S. spontaneum and Sorghum bicolor, and gene expression analysis was performed under drought stress and with exogenous applications of plant hormones. FRK1, which was under strong functional constraint selection, was conserved among the gene allelic haplotypes, and displayed dominant expression levels among the gene families in the control conditions, suggesting that FRK1 plays a major role in the phosphorylation of fructose. FRK3 and FRK5 were dramatically induced under drought stress, and FRK5 was also found to increase its expression levels in the mature stage of Saccharum. Similarly, FRK3 and FRK5 were induced in response to drought stress in Saccharum. FRK2 and FRK7 displayed lower expression levels than the other FRK family members; FRK2 was under strong genomic selection constraints whereas FRK7 was under neutral selection. FRK7 may have become functionally redundant in Saccharum through pseudogenization. FRK4 and FRK6 shared the most similar expression pattern: FRK4 was revealed to have higher expression levels in mature tissues than in premature tissues of Saccharum, and FRK6 presented a slight increase under drought stress.ConclusionsOur study presents a comprehensive genomic study of the entire FRK gene family in Saccharum, providing the foundations for approaches to characterize the molecular mechanism regulated by the SsFRK family in sugarcane.