Nonsynonymous Substitution Rates Research Articles

Host Adaptive Evolution of Avian-Origin H3N2 Canine Influenza Virus.

Since its first isolation in around 2007, the avian-origin H3N2 canine influenza virus (CIV) has become established and continues to circulate in dog populations. This virus serves as a useful model for deciphering the complex evolutionary process of interspecies transmission of influenza A virus (IAV) from one species to its subsequent circulation in another mammalian host. The present investigation is a comprehensive effort to identify and characterize genetic changes that accumulated in the avian-origin H3N2 CIV during its circulation in the dog. We revealed that H3N2 CIV experiences greater selection pressure with extremely high global non-synonymous to synonymous substitution ratios per codon (dN/dS ratio) for each gene compared to the avian reservoir viruses. A total of 54 amino acid substitutions were observed to have accumulated and become fixed in the H3N2 CIV population based on our comprehensive codon-based frequency diagram analysis. Of these substitutions, 11 sites also display high prevalence in H3N8 CIV, indicating that convergent evolution has occurred on different lineages of CIV. Notably, six substitutions, including HA-G146S, M1-V15I, NS1-E227K, PA-C241Y, PB2-K251R, and PB2-G590S, have been reported to play imperative roles in facilitating the transmission and spillover of IAVs across species barriers. Most of these substitutions were found to have become fixed in around 2015, which might have been a favorable factor that facilitating the spread of these CIV lineages from South Asia to North America and subsequent further circulation in these areas. We also detected 12 sites in six viral genes with evidence for positive selection by comparing the rates of non-synonymous and synonymous substitutions at each site. Besides, our study reports trends of enhanced ongoing adaptation of H3N2 CIV to their respective host cellular systems, based on the codon adaptation index analysis, which points toward increasing fitness for efficient viral replication. In addition, a reduction in the abundance of the CpG motif, as evident from an analysis of relative dinucleotide abundance, may contribute to the successful evasion of host immune recognition. The present study provides key insights into the adaptive changes that have accumulated in the avian-origin H3N2 viral genomes during its establishment and circulation into dog populations.

Analysis of recently duplicated TYRP1 genes and their effect on the formation of black patches in Oujiang-color common carp (Cyprinus carpio var. color).

Tyrp1 gene, as a member of the tyrosinase family, has undergone a recent duplication event during fourth-round whole genome duplication in common carp. In this research, three Tyrp1 genes were identified in Oujiang-color common carp (Cyprinus carpio var. color). The similar expression patterns and close phylogenetic relationship indicated that Tyrp1c is homologous to Tyrp1b and possibly originated from the ancient Tyrp1b. The rates of synonymous and non-synonymous substitution (Ka /Ks ) in Tyrp1 across teleost phylogeny indicated that Tyrp1a is more likely to be in the process of purifying selection. The CRISPR/Cas9 system was used to disrupt the Tyrp1 genes in zebrafish and the WB (black patches on white skin) strain of Oujiang-color common carp. The Tyrp1 loss of function variants in zebrafish and WB carp showed severe melanin deficiency in the skin. Meanwhile, inactivation of a single Tyrp1 gene did not obstruct melanin synthesis, which proved that the functional redundancy of Tyrp1 genes existed in both zebrafish and Oujiang-color common carp. Among the mosaic individuals with Tyrp1 genes in disrupted-color common carp, various mutations in Tyrp1b gene induced gray or brown phenotypes, suggesting that it may be bifunctional in Oujiang-color common carp. In addition, the phenotype of WB variants was different from that of WW (whole white skin), suggesting that Tyrp1 genes were not the key factor that caused the difference between WB and WW.

Population Genomics of the Maize Pathogen Ustilago maydis: Demographic History and Role of Virulence Clusters in Adaptation.

The tight interaction between pathogens and their hosts results in reciprocal selective forces that impact the genetic diversity of the interacting species. The footprints of this selection differ between pathosystems because of distinct life-history traits, demographic histories, or genome architectures. Here, we studied the genome-wide patterns of genetic diversity of 22 isolates of the causative agent of the corn smut disease, Ustilago maydis, originating from five locations in Mexico, the presumed center of origin of this species. In this species, many genes encoding secreted effector proteins reside in so-called virulence clusters in the genome, an arrangement that is so far not found in other filamentous plant pathogens. Using a combination of population genomic statistical analyses, we assessed the geographical, historical, and genome-wide variation of genetic diversity in this fungal pathogen.We report evidence of two partially admixed subpopulations that are only loosely associated with geographic origin. Using the multiple sequentially Markov coalescent model, we inferred the demographic history of the two pathogen subpopulations over the last 0.5 Myr. We show that both populations experienced a recent strong bottleneck starting around 10,000 years ago, coinciding with the assumed time of maize domestication. Although the genome average genetic diversity is low compared with other fungal pathogens, we estimated that the rate of nonsynonymous adaptive substitutions is three times higher in genes located within virulence clusters compared with nonclustered genes, including nonclustered effector genes. These results highlight the role that these singular genomic regions play in the evolution of this pathogen.

Evolution of MHC class I genes in Japanese and Russian raccoon dogs, Nyctereutes procyonoides (Carnivora: Canidae)

Major histocompatibility complex (MHC) genes have been widely studied to assess the immunological fitness and evolutionary adaptation of animal populations. Among the Canidae, the raccoon dog’s adventurous nature, omnivorous behavior, and high variability of intracellular pathogens make it ideal to study selection on MHC class I in a non-model canid species. Here, we examined allelic diversity and evolutionary patterns of MHC class I genes in the raccoon dog (Nyctereutes procyonoides). We identified 48 novel MHC class I alleles from 31 raccoon dogs from Japan and Russia. Some alleles were geographically restricted, whereas others were widely distributed across the species’ range. The rate of non-synonymous substitutions was greater than that of synonymous substitutions for both exon 2 and exon 3 encoding α1 and α2 domains, respectively, in the α chain of the MHC class I protein. Positively selected sites at the amino acid level were evident in both the α1 and α2 domains, and a recombination breakpoint was found in exon 3. Bayesian phylogenetic trees showed no evidence of trans-species polymorphism (TSP) with alleles from carnivoran species in other families but did detect TSP between raccoon dogs and the domestic dog, Canis familiaris, indicative of long-term balancing selection in canids. Our results indicate that the extensive allelic diversity of MHC class I in Japanese and Russian raccoon dogs has been influenced and maintained by pathogen-driven positive selection, recombination, and long-term balancing selection.Supplementary InformationThe online version contains supplementary material available at 10.1007/s13364-021-00561-y.

Complete chloroplast genome features and phylogenetic analysis of Eruca sativa (Brassicaceae).

Eruca sativa Mill. (Brassicaceae) is an important edible vegetable and a potential medicinal plant due to the antibacterial activity of its seed oil. Here, the complete chloroplast (cp) genome of E. sativa was de novo assembled with a combination of long PacBio reads and short Illumina reads. The E. sativa cp genome had a quadripartite structure that was 153,522 bp in size, consisting of one large single-copy region of 83,320 bp and one small single-copy region of 17,786 bp which were separated by two inverted repeat (IRa and IRb) regions of 26,208 bp. This complete cp genome harbored 113 unique genes: 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. Forty-nine long repetitive sequences and 69 simple sequence repeats were identified in the E. sativa cp genome. A codon usage analysis of the E. sativa cp genome showed a bias toward codons ending in A/T. The E. sativa cp genome was similar in size, gene composition, and linearity of the structural region when compared with other Brassicaceae cp genomes. Moreover, the analysis of the synonymous (Ks) and non-synonymous (Ka) substitution rates demonstrated that protein-coding genes generally underwent purifying selection pressure, expect ycf1, ycf2, and rps12. A phylogenetic analysis determined that E. sativa is evolutionarily close to important Brassica species, indicating that it may be possible to transfer favorable E. sativa alleles into other Brassica species. Our results will be helpful to advance genetic improvement and breeding of E. sativa, and will provide valuable information for utilizing E. sativa as an important resource to improve other Brassica species.

Adaptation and convergence in circadian\u2010related genes in Iberian freshwater fish

BackgroundThe circadian clock is a biological timing system that improves the ability of organisms to deal with environmental fluctuations. At the molecular level it consists of a network of transcription-translation feedback loops, involving genes that activate (bmal and clock – positive loop) and repress expression (cryptochrome (cry) and period (per) – negative loop). This is regulated by daily alternations of light but can also be affected by temperature. Fish, as ectothermic, depend on the environmental temperature and thus are good models to study its integration within the circadian system. Here, we studied the molecular evolution of circadian genes in four Squalius freshwater fish species, distributed across Western Iberian rivers affected by two climatic types with different environmental conditions (e.g., light and temperature). S. carolitertii and S. pyrenaicus inhabit the colder northern region under Atlantic climate type, while S. torgalensis, S. aradensis and some populations of S. pyrenaicus inhabit the warmer southern region affected by summer droughts, under Mediterranean climate type.ResultsWe identified 16 circadian-core genes in the Squalius species using a comparative transcriptomics approach. We detected evidence of positive selection in 12 of these genes using methods based on dN/dS. Positive selection was mainly found in cry and per genes of the negative loop, with 55 putatively adaptive substitutions, 16 located on protein domains. Evidence for positive selection is predominant in southern populations affected by the Mediterranean climate type. By predicting protein features we found that changes at sites under positive selection can impact protein thermostability by changing their aliphatic index and isoelectric point. Additionally, in nine genes, the phylogenetic clustering of species that belong to different clades but inhabit southern basins with similar environmental conditions indicated evolutionary convergence. We found evidence for increased nonsynonymous substitution rate in convergent lineages, likely due to positive selection at 27 sites, mostly in cry genes.ConclusionsOur results support that temperature may be a selective pressure driving the evolution of genes involved in the circadian system. By integrating sequence-based functional protein prediction with dN/dS-based methods to detect selection we uncovered adaptive convergence in the southern populations, probably related to their similar thermal conditions.

Molecular Evolution of rbcL in Orthotrichales (Bryophyta): Site Variation, Adaptive Evolution, and Coevolutionary Patterns of Amino Acid Replacements.

Molecular evolution of the large subunit of the RuBisCO enzyme is understudied in early diverging land plants. These groups show morphological and eco-physiological adaptations to the uneven and intermittent distribution of water in the terrestrial environment. This might have prompted a continuous fine-tuning of RuBisCO under a selective pressure modifying the species-specific optima for photosynthesis in contrasting microdistributions and environmental niches. To gain a better insight into the molecular evolution of RuBisCO large subunits, the aim of this study was to assess the pattern of evolutionary change in the amino acid residues in a monophyletic group of Bryophyta (Orthotrichaceae). Tests for positive, neutral, or purifying selection at the amino acid level were assessed by comparing rates (ω) of non-synonymous (dN) and synonymous (dS) nucleotide substitutions along a Maximum Likelihood phylogenetic tree. Molecular adaptation tests using likelihood ratio tests, reconstruction of ancestral amino acid sites, and intra-protein coevolution analyses were performed. Variable amino acid sites (39) were unevenly distributed across the LSU. The residues are located on rbcL sites that are highly variable in higher plants and close to key regions implying dimer-dimer (L2L2), RuBisCO-activase interactions, and conformational functions during catalysis. Ten rbcL sites (32, 33, 91, 230, 247, 251, 255, 424, 449 and 475) have been identified by the Bayesian Empirical Bayes inference to be under positive selection and under adaptive evolution under the M8 model. The pattern of amino acid variation suggests that it is not lineage specific, but rather representative of a case of convergent evolution, suggesting recurrent changes that potentially favor the same amino acid substitutions that are likely optimized the RuBisCO activity.

Two Korean Endemic Clematis Chloroplast Genomes: Inversion, Reposition, Expansion of the Inverted Repeat Region, Phylogenetic Analysis, and Nucleotide Substitution Rates.

Previous studies on the chloroplast genome in Clematis focused on the chloroplast structure within Anemoneae. The chloroplast genomes of Cleamtis were sequenced to provide information for studies on phylogeny and evolution. Two Korean endemic Clematis chloroplast genomes (Clematis brachyura and C. trichotoma) range from 159,170 to 159,532 bp, containing 134 identical genes. Comparing the coding and non-coding regions among 12 Clematis species revealed divergent sites, with carination occurring in the petD-rpoA region. Comparing other Clematis chloroplast genomes suggested that Clematis has two inversions (trnH-rps16 and rps4), reposition (trnL-ndhC), and inverted repeat (IR) region expansion. For phylogenetic analysis, 71 protein-coding genes were aligned from 36 Ranunculaceae chloroplast genomes. Anemoneae (Anemoclema, Pulsatilla, Anemone, and Clematis) clades were monophyletic and well-supported by the bootstrap value (100%). Based on 70 chloroplast protein-coding genes, we compared nonsynonymous (dN) and synonymous (dS) substitution rates among Clematis, Anemoneae (excluding Clematis), and other Ranunculaceae species. The average synonymoussubstitution rates (dS)of large single copy (LSC), small single copy (SSC), and IR genes in Anemoneae and Clematis were significantly higher than those of other Ranunculaceae species, but not the nonsynonymous substitution rates (dN). This study provides fundamental information on plastid genome evolution in the Ranunculaceae.

Kinetochore proteins and microtubule-destabilizing factors are fast evolving in eutherian mammals.

Centromeres have central functions in chromosome segregation, but centromeric DNA and centromere-binding proteins evolve rapidly in most eukaryotes. The selective pressure(s) underlying the fast evolution of centromere-binding proteins are presently unknown. An attractive possibility is that selfish centromeres promote their preferential inclusion in the oocyte and centromeric proteins evolve to suppress meiotic drive (centromere drive hypothesis). We analysed the selective patterns of mammalian genes that encode kinetochore proteins and microtubule (MT)-destabilizing factors. We show that several of these proteins evolve at the same rate or faster than proteins with a role in centromere specification. Elements of the kinetochore that bind MTs or that bridge the interaction between MTs and the centromere represented the major targets of positive selection. These data are in line with the possibility that the genetic conflict fuelled by meiotic drive extends beyond genes involved in centromere specification. However, we cannot exclude that different selective pressures underlie the rapid evolution of MT-destabilizing factors and kinetochore components. Whatever the nature of such pressures, they must have been constant during the evolution of eutherian mammals, as we found a surprisingly good correlation in dN/dS (ratio of the rate of nonsynonymous and synonymous substitutions) across orders/clades. Finally, when phylogenetic relationships were accounted for, we found little evidence that the evolutionary rates of these genes change with testes size, a proxy for sperm competition. Our data indicate that, in analogy to centromeric proteins, kinetochore components are fast evolving in mammals. This observation may imply that centromere drive plays out at multiple levels or that these proteins adapt to lineage-specific centromeric features.

Comparative genomic analysis reveals evolutionary and structural attributes of MCM gene family in Arabidopsis thaliana and Oryza sativa

The mini-chromosome maintenance (MCM) family, a large and functionally diverse protein family belonging to the AAA+ superfamily, is essential for DNA replication in all eukaryotic organisms. The MCM 2-7 form a hetero-hexameric complex which serves as licensing factor necessary to ensure the proper genomic DNA replication during the S phase of cell cycle. MCM 8-10 are also associated with the DNA replication process though their roles are particularly unclear. In this study, we report an extensive in silico analysis of MCM gene family (MCM 2-10) in Arabidopsis and rice. Comparative analysis of genomic distribution across eukaryotes revealed conservation of core MCMs 2-7 while MCMs 8-10 are absent in some taxa. Domain architecture analysis underlined MCM 2-10 subfamily specific features. Phylogenetic analyses clustered MCMs into 9 clades as per their subfamily. Duplication events are prominent in plant MCM family, however no duplications are observed in Arabidopsis and rice MCMs. Synteny analysis among Arabidopsis thaliana, Oryza sativa, Glycine max and Zea mays MCMs demonstrated orthologous relationships and duplication events. Further, estimation of synonymous and non-synonymous substitution rates illustrated evolution of MCM family under strong constraints. Expression profiling using available microarray data and qRT-PCR revealed differential expression under various stress conditions, hinting at their potential use to develop stress resilient crops. Homology modeling of Arabidopsis and rice MCM 2-7 and detailed comparison with yeast MCMs identified conservation of eukaryotic specific insertions and extensions as compared to archeal MCMs. Protein-protein interaction analysis revealed an extensive network of putative interacting partners mainly involved in DNA replication and repair. The present study provides novel insights into the MCM family in Arabidopsis and rice and identifies unique features, thus opening new perspectives for further targeted analyses.

Comparative Analyses of the Chloroplast Genomes of Patchouli Plants and Their Relatives in Pogostemon (Lamiaceae).

Pogostemon Desf., the largest genus of the tribe Pogostemoneae (Lamiaceae), consists of ca. 80 species distributed mainly from South and Southeast Asia to China. The genus contains many patchouli plants, which are of great economic importance but taxonomically difficult. Therefore, it is necessary to characterize more chloroplast (cp) genomes for infrageneric phylogeny analyses and species identification of Pogostemon, especially for patchouli plants. In this study, we newly generated four cp genomes for three patchouli plants (i.e., Pogostemon plectranthoides Desf., P. septentrionalis C. Y. Wu et Y. C. Huang, and two cultivars of P. cablin (Blanoco) Benth.). Comparison of all samples (including online available cp genomes of P. yatabeanus (Makino) Press and P. stellatus (Lour.) Kuntze) suggested that Pogostemon cp genomes are highly conserved in terms of genome size and gene content, with a typical quadripartite circle structure. Interspecific divergence of cp genomes has been maintained at a relatively low level, though seven divergence hotspot regions were identified by stepwise window analysis. The nucleotide diversity (Pi) value was correlated significantly with gap proportion (indels), but significantly negative with GC content. Our phylogenetic analyses based on 80 protein-coding genes yielded high-resolution backbone topologies for the Lamiaceae and Pogostemon. For the overall mean substitution rates, the synonymous (dS) and nonsynonymous (dN) substitution rate values of protein-coding genes varied approximately threefold, while the dN values among different functional gene groups showed a wider variation range. Overall, the cp genomes of Pogostemon will be useful for phylogenetic reconstruction, species delimitation and identification in the future.

Nucleotide substitution rates of diatom plastid encoded protein genes are positively correlated with genome architecture

Diatoms are the largest group of heterokont algae with more than 100,000 species. As one of the single-celled photosynthetic organisms that inhabit marine, aquatic and terrestrial ecosystems, diatoms contribute ~ 45% of global primary production. Despite their ubiquity and environmental significance, very few diatom plastid genomes (plastomes) have been sequenced and studied. This study explored patterns of nucleotide substitution rates of diatom plastids across the entire suite of plastome protein-coding genes for 40 taxa representing the major clades. The highest substitution rate was lineage-specific within the araphid 2 taxon Astrosyne radiata and radial 2 taxon Proboscia sp. Rate heterogeneity was also evident in different functional classes and individual genes. Similar to land plants, proteins genes involved in photosynthetic metabolism have lower synonymous and nonsynonymous substitutions rates than those involved in transcription and translation. Significant positive correlations were identified between substitution rates and measures of genomic rearrangements, including indels and inversions, which is a similar result to what was found in legume plants. This work advances the understanding of the molecular evolution of diatom plastomes and provides a foundation for future studies.

Transcriptome analysis and comparison reveal divergence between the Mediterranean and the greenhouse whiteflies

Both the Mediterranean (MED) species of the Bemisia tabaci whitefly complex and the greenhouse whitefly (Trialeurodes vaporariorum, TV) are important agricultural pests. The two species of whiteflies differ in many aspects such as morphology, geographical distribution, host plant range, plant virus transmission, and resistance to insecticides. However, the molecular basis underlying their differences remains largely unknown. In this study, we analyzed the genetic divergences between the transcriptomes of MED and TV. In total, 2,944 pairs of orthologous genes were identified. The average identity of amino acid sequences between the two species is 93.6%. The average nonsynonymous (Ka) and synonymous (Ks) substitution rates and the ratio of Ka/Ks of the orthologous genes are 0.0389, 2.23 and 0.0204, respectively. The low average Ka/Ks ratio indicates that orthologous genes tend to be under strong purified selection. The most divergent gene classes are related to the metabolisms of xenobiotics, cofactors, vitamins and amino acids, and this divergence may underlie the different biological characteristics between the two species of whiteflies. Genes of differential expression between the two species are enriched in carbohydrate metabolism and regulation of autophagy. These findings provide molecular clues to uncover the biological and molecular differences between the two species of whiteflies.

Genome-wide identification and characterization of HSP70 gene family in four species of cotton

Heat shock proteins (HSPs) are important elements of the cellular group of molecular chaperones. Specifically, HSP70 proteins protect cells from being damaged when plants are exposed to environmental stresses. These proteins are catalysts that manage the correct folding of other proteins, and they play a key role in the development of tolerance against biotic and abiotic stresses. In the present study, 113 HSP70 genes were retrieved from the available genome assemblies of four cotton species, including Gossypium hirsutum, G. barbadense, G. arboreum, and G. raimondii. The HSP70 genes were clustered into 11 subfamilies based on phylogeny. One hundred and nine (109) gene duplications were found across these four species. Localization of genes revealed that several HSP70 genes reside in the cytoplasm. Synonymous and non-synonymous substitution rates revealed that functional segregation of HSP70 genes in cotton is due to purifying selection. Furthermore, HSP70 genes in cotton are expressed constitutively during developmental stages. These findings are valuable to understand the complex mechanism of HSP70 gene regulation that occurs in signaling pathways in response to plant stress.

Molecular Evolution and Developmental Expression of Melanin Pathway Genes in Lepidoptera

Pigmentation is involved in a wide array of biological functions across insect orders, including body patterning, thermoregulation, and immunity. The melanin pathway, in particular, has been characterized in several species. However, molecular evolution of the genes involved in this pathway is poorly explored. We traced the molecular evolution of six melanin pathway genes in 53 species of Lepidoptera covering butterflies and moths, and representing over 100 million years of diversification. We compared the rates of synonymous and nonsynonymous substitutions within and between these genes to study signatures of selection at the level of individual sites, genes, and branches of the gene tree. We found that molecular evolution of all six genes was governed by strong purifying selection. Yet, a number of sites showed signs of being under positive selection, including in the highly conserved domain regions of three genes. Further, we traced the expression of these genes across developmental stages, tissues, and sexes in the Papilio polytes butterfly using a developmental transcriptome dataset. We observed that the expression patterns of the genes in P. polytes largely reflected their known tissue-specific function in other species. The expression of sequentially acting genes in the melanin pathway was correlated. Interestingly, melanin pathway genes also showed a sexually dimorphic pattern of developmental heterochrony; i.e., females showed prominent upregulation of melanin pathway genes in pre-pupal stage compared to males, while males showed prominent upregulation in 9-day pupal wings compared to females. Our evolutionary and developmental analyses suggest that the vast diversity of wing patterning and pigmentation in Lepidoptera may have evolved despite largely constrained sequence evolution, with potential contribution from differential developmental expression of genes in a highly conserved pathway.

Intra-Protein Coevolution Is Increasingly Functional with Greater Proximity to Fertilization

Intramolecular coevolution of amino acid sites has repeatedly been studied to improve predictions on protein structure and function. Thereby, the focus was on bacterial proteins with available crystallographic data. However, intramolecular coevolution has not yet been compared between protein sets along a gradient of functional proximity to fertilization. This is especially true for the potential effect of external selective forces on intraprotein coevolution. In this study, we investigated both aspects in equally sized sets of mammalian proteins representing spermatozoa, testis, entire body, and liver. For coevolutionary analyses, we derived the proportion of covarying sites per protein from amino acid alignments of 10 mammalian orthologues each. In confirmation of the validity of our coevolution proxy, we found positive associations with the nonsynonymous or amino acid substitution rate in all protein sets. However, our coevolution proxy negatively correlated with the number of protein interactants (node degree) in male reproductive protein sets alone. In addition, a negative association of our coevolution proxy with protein hydrophobicity was significant in sperm proteins only. Accordingly, the restrictive effect of protein interactants was most pronounced in male reproductive proteins, and the tendency of sperm proteins to form internal structures decreased the more coevolutionary sites they had. Both aspects illustrate that the share of outward and thus functional coevolution increases with greater proximity to fertilization. We found this conclusion confirmed by additional comparisons within sperm proteins. Thus, sperm proteins with high hydrophobicity had the lowest proportions of covarying sites and, according to gene annotations, localized more frequently to internal cellular structures. They should therefore be less exposed to postcopulatory forms of sexual selection. Their counterparts with low hydrophobicity had larger proportions of covarying sites and more often resided at the cell membrane or were secreted. At the cellular level, they are thus closer to externally induced forces of postcopulatory selection which are known for their potential to increase substitution rates. In addition, we show that the intermediary status of the testicular protein set in correlation analyses is probably due to a special combination of reproductive and somatic involvements.

Link Your Sites (LYS) Scripts: Automated Search of Protein Structures and Mapping of Sites Under Positive Selection Detected by PAML

The visualization of the molecular context of an amino acid mutation in a protein structure is crucial for the assessment of its functional impact and the understanding of its evolutionary implications. Currently, searches for fast evolving amino acid positions using codon substitution models like those implemented in PAML (Yang and Nielsen in Estimating synonymous and nonsynonymous substitution rates under realistic evolutionary models. Mol Biol Evol 17(1):32–43, 2000; Zhang et al. in Evaluation of an improved branch-site likelihood method for detecting positive selection at the molecular level. Mol Biol Evol 22(12):2472–2479, 2005) are done in almost complete proteomes, generating large numbers of candidate proteins making the analysis of individual protein structures and models very time-consuming. Here we present the package Link Your Sites (LYS) that can be used to reduce the number of analysed targets to those for which structural information can be retrieved. LYS consists of two python wrapper scripts, where the first one (i) mines the RCSB database (Berman et al. in The protein data bank. Nucleic Acids Res 28(1):235–242, 2000) using the BLAST alignment tool to find the best matching homologous sequences, (ii) fetches their domain positions by using Prosites (Hamelryck and Manderick in Pdb file parser and structure class implemented in python. Bioinformatics 19(17):2308–2310, 2003; Sigrist et al. in Prosite: a documented database using patterns and profiles as motif descriptors. Brief Bioinf 3(3):265–274, 2002; Sigrist et al. in New and continuing developments at prosite. Nucleic Acids Res 41(D1):D344–D347, 2012), (iii) parses the output of PAML extracting the positional information of fast-evolving sites and transforms them into the coordinate system of the protein structure, (iv) outputs one file per gene with the equivalence among the positions in the input sequence and homologous structure. The second script produces figures to be used in publications highlighting the positively selected sites mapped on regions that are known to have functional relevance.

Molecular evolution and genetic variations of \u2009V and W proteins derived by RNA editing in Avian Paramyxoviruses

The newly assigned subfamily Avulavirinae in the family Paramyxoviridae includes avian paramyxoviruses (APMVs) isolated from a wide variety of avian species across the globe. Till date, 21 species of APMVs are reported and their complete genome sequences are available in GenBank. The APMV genome comprises of a single stranded, negative sense, non-segmented RNA comprising six transcriptional units (except APMV-6 with seven units) each coding for a structural protein. Additionally, by co-transcriptional RNA editing of phosphoprotein (P) gene, two mRNAs coding for accessory viral proteins, V and W, are generated along with unedited P mRNA. However, in APMV-11, the unedited mRNA codes for V protein while +2 edited mRNA translates to P protein, similar to members of subfamily Rubulavirinae in the same family. Such RNA editing in paramyxoviruses enables maximizing the coding capacity of their smaller genome. The three proteins of P gene: P, V and W, share identical N terminal but varied C terminal sequences that contribute to their unique functions. Here, we analyzed the P gene editing site, V and W sequences of all 21 APMV species known so far (55 viruses) by using bioinformatics and report their genetic variations and molecular evolution. The variations observed in the sequence and hexamer phase positions of the P gene editing sites is likely to influence the levels and relative proportions of P, V and W proteins’ expressions which could explain the differences in the pathogenicity of APMVs. The V protein sequences of APMVs had conserved motifs similar to V proteins of other paramyxoviruses including the seven cysteine residues involved in MDA5 interference, STAT1 degradation and interferon antagonism. Conversely, W protein sequences of APMVs were distinct. High sequence homology was observed in both V and W proteins between strains of the same species than between species except in APMV-3 which was the most divergent APMV species. The estimates of synonymous and non-synonymous substitution rates suggested negative selection pressure on the V and W proteins within species indicating their low evolution rate. The molecular clock analysis revealed higher conservation of V protein sequence compared to W protein indicating the important role played by V protein in viral replication, pathogenesis and immune evasion. However, we speculate the genetic diversity of W proteins could impact the degree of pathogenesis, variable interferon antagonistic activity and the wide host range exhibited by APMV species. Phylogenetically, V proteins of APMVs clustered into three groups similar to the recent classification of APMVs into three new genera while no such pattern could be deciphered in the analysis of W proteins except that strains of same species grouped together. This is the first comprehensive study describing in detail the genetic variations and the molecular evolution of P gene edited, accessory viral proteins of Avian paramyxoviruses.

Amino Acid Substitutions in Positions 385 and 393 of the Hydrophobic Region of VP4 May Be Associated with Rotavirus Attenuation and Cell Culture Adaptation.

Rotaviruses (RVs) are the leading cause of the acute viral gastroenteritis in young children and livestock animals worldwide. Although live attenuated vaccines have been applied to control RV infection for many years, the underlying mechanisms of RV attenuation following cell culture adaption are unknown. To study these mechanisms at the genomic level, we have sequenced and conducted a comparative analysis of two virulent human (Wa, G1P[8] and M, G3P[8]) and two virulent porcine (Gottfried, G4P[6] and OSU, G5P[7]) RV strains maintained in gnotobiotic piglets for 22, 11, 12 and 9 serial passages, respectively, with their attenuated counterparts serially passaged in MA-104 cell cultures for 25, 43, 54 and 43 passages, respectively. We showed that most of the mutations were clustered in the VP4 gene, with a relatively high nonsynonymous substitution rate (81.2%). Moreover, two amino acid substitutions observed in the VP4 gene were conserved between two or more strain pairs. D385N substitution was found in M, Wa and Gottfried strains, and another one, S471H/L was present in Wa and Gottfried strains. Importantly, D385 was reported previously in another study and may be involved in regulation of virus entry. Of interest, although no 385 substitution was found in OSU strains, the attenuated OSU strain contained a unique D393H substitution within the same VP4 hydrophobic domain. Collectively, our data suggest that the VP4 hydrophobic region may play an important role in RV attenuation and aa385 and aa393 may represent potential targets for RV vaccine development using reverse genetics and site-specific mutagenesis.

Dissecting Fertility Functions of Drosophila Y Chromosome Genes with CRISPR.

Gene-poor, repeat-rich regions of the genome are poorly understood and have been understudied due to technical challenges and the misconception that they are degenerating "junk." Yet multiple lines of evidence indicate these regions may be an important source of variation that could drive adaptation and species divergence, particularly through regulation of fertility. The ∼40 Mb Y chromosome of Drosophila melanogaster contains only 16 known protein-coding genes, and is highly repetitive and entirely heterochromatic. Most of the genes originated from duplication of autosomal genes and have reduced nonsynonymous substitution rates, suggesting functional constraint. We devised a genetic strategy for recovering and retaining stocks with sterile Y-linked mutations and combined it with CRISPR to create mutants with deletions that disrupt three Y-linked genes. Two genes, PRY and FDY, had no previously identified functions. We found that PRY mutant males are subfertile, but FDY mutant males had no detectable fertility defects. FDY, the newest known gene on the Y chromosome, may have fertility effects that are conditional or too subtle to detect. The third gene, CCY, had been predicted but never formally shown to be required for male fertility. CRISPR targeting and RNA interference of CCY caused male sterility. Surprisingly, however, our CCY mutants were sterile even in the presence of an extra wild-type Y chromosome, suggesting that perturbation of the Y chromosome can lead to dominant sterility. Our approach provides an important step toward understanding the complex functions of the Y chromosome and parsing which functions are accomplished by genes vs. repeat elements.