The Nonsynonymous/Synonymous Substitution Rate Ratio versus the Radical/Conservative Replacement Rate Ratio in the Evolution of Mammalian Genes

Genome-wide identification and expression analysis of the WRKY gene family in Rhododendron henanense subsp. lingbaoense

Background

This work explored the characteristics of the WRKY transcription factor family in Rhododendron henanense subsp. lingbaoense (Rhl) and the expression patterns of these genes under abiotic stress by conducting bioinformatics and expression analyses.

Methods

RhlWRKY genes were identified from a gene library of Rhl. Various aspects of these genes were analyzed, including genetic structures, conserved sequences, physicochemical properties, cis-acting elements, and chromosomal location. RNA-seq was employed to analyze gene expression in five different tissues of Rhl: roots, stems, leaves, flowers, and hypocotyls. Additionally, qRT-PCR was used to detect changes in the expression of five RhlWRKY genes under abiotic stress.

Result

A total of 65 RhlWRKY genes were identified and categorized into three subfamilies based on their structural characteristics: Groups I, II, and III. Group II was further divided into five subtribes, with shared similar genetic structures and conserved motifs among members of the same subtribe. The physicochemical properties of these proteins varied, but the proteins are generally predicted to be hydrophilic. Most proteins are predicted to be in the cell nucleus, and distributed across 12 chromosomes. A total of 84 cis-acting elements were discovered, with many related to responses to biotic stress. Among the identified RhlWRKY genes, there were eight tandem duplicates and 97 segmental duplicates. The majority of duplicate gene pairs exhibited Ka/Ks values <1, indicating purification under environmental pressure. GO annotation analysis indicated that WRKY genes regulate biological processes and participate in a variety of molecular functions. Transcriptome data revealed varying expression levels of 66.15% of WRKY family genes in all five tissue types (roots, stems, leaves, flowers, and hypocotyls). Five RhlWRKY genes were selected for further characterization and there were changes in expression levels for these genes in response to various stresses.

Conclusion

The analysis identified 65 RhlWRKY genes, among which the expression of WRKY_42 and WRKY_17 were mainly modulated by the drought and MeJA, and WRKY_19 was regulated by the low-temperature and high-salinity conditions. This insight into the potential functions of certain genes contributes to understanding the growth regulatory capabilities of Rhl.

Genome-wide identification and expression analysis of Na+/H+antiporter (NHX) genes in tomato under salt stress

Plant Na +/H + antiporter (NHX) genes enhance salt tolerance by preventing excessive Na+ accumulation in the cytosol through partitioning of Na+ ions into vacuoles or extracellular transport across the plasma membrane. However, there is limited detailed information regarding the salt stress responsive SlNHXs in the most recent tomato genome. We investigated the role of this gene family's expression patterns in the open flower tissues under salt shock in Solanum lycopersicum using a genome-wide approach. A total of seven putative SlNHX genes located on chromosomes 1, 4, 6, and 10 were identified, but no ortholog of the NHX5 gene was identified in the tomato genome. Phylogenetic analysis revealed that these genes are divided into three different groups. SlNHX proteins with 10-12 transmembrane domains were hypothetically localized in vacuoles or cell membranes. Promoter analysis revealed that SlNHX6 and SlNHX8 are involved with the stress-related MeJA hormone in response to salt stress signaling. The structural motif analysis of SlNHX1, -2, -3, -4, and -6 proteins showed that they have highly conserved amiloride binding sites. The protein-protein network revealed that SlNHX7 and SlNHX8 interact physically with Salt Overly Sensitive (SOS) pathway proteins. Transcriptome analysis demonstrated that the SlNHX2 and SlNHX6 genes were substantially expressed in the open flower tissues. Moreover, quantitative PCR analysis indicated that all SlNHX genes, particularly SlNHX6 and SlNHX8, are significantly upregulated by salt shock in the open flower tissues. Our results provide an updated framework for future genetic research and development of breeding strategies against salt stress in the tomato.

Selection of mutant Listeria phages under food-relevant conditions can enhance application potential

Bacteriophages are viruses that infect and kill bacteria. Currently, phage products are available for the control of the pathogen Listeria monocytogenes in food products in the United States. In this study, we explore whether experimental evolution can be used to generate phages with improved abilities to function under specific food-relevant conditions. Ultra-pasteurized oat and whole milk were chosen as test matrices as they represent different food groups, yet have similar physical traits and macronutrient composition. We showed that (i) wild-type phage LP-125 infection kinetics are different in the two matrices and (ii) LP-125 has a significantly higher burst size in oat milk. From this, we attempted to evolve LP-125 to have improved infection kinetics in whole milk. Ancestral LP-125 was passaged through 10 rounds of amplification in milk conditions. Plaque-purified DNA samples from milk-selected phages were isolated and sequenced, and mutations present in the isolated phages were identified. We found two nonsynonymous substitutions in LP125_108 and LP125_112 genes, which encode putative baseplate-associated glycerophosphoryl diester phosphodiesterase and baseplate protein, respectively. Protein structural modeling showed that the substituted amino acids in the mutant phages are predicted to localize to surface-exposed helices on the corresponding structures, which might affect the surface charge of proteins and their interaction with the bacterial cell. The phage containing the LP125_112 mutation adsorbed significantly faster than the ancestral phage in both oat and whole milk. Follow-up experiments suggest that fat content may be a key factor for the expression of the phenotype of this mutation. IMPORTANCE Bacteriophages are one of the tools available to control the foodborne pathogen, Listeria monocytogenes. Phage products must work under a broad range of food conditions to be an effective control for L. monocytogenes. Here, we show that the experimental evolution of phages can be used to generate new phages with phenotypes useful under specific conditions. We used this approach to select for a mutant phage that more efficiently binds to L. monocytogenes that is grown in whole milk and oat milk. We show that the fat content of these milks is necessary for the expression of this phenotype. Our findings show that experimental evolution can be used to select for improved phages with better performance under specific conditions. This approach has the potential to support the development of condition-specific phage-based biocontrols in the food industry.

Two Complete Mitochondrial Genomes of Leptobrachium (Anura: Megophryidae: Leptobrachiinae): Characteristics, Population Divergences, and Phylogenetic Implications

The mustache toads Leptobrachium boringii and Leptobrachium liui are two attractive species in Megophryidae, in which adult males have mustache-like keratinized nuptial spines on their upper lip. However, both are under threat due to multiple factors, of which scientific studies are still very limited. In this study, two new complete mitochondrial genomes of L. boringii and L. liui were sequenced, assembled, and annotated based on next-generation sequencing. The mitogenome lengths of L. boringii and L. liui were found to be 17,100 and 17,501 bp, respectively, with both containing 13 protein coding genes, 23 tRNAs, 2 rRNAs, and 1 non-coding control region. Nucleotide diversity analyses indicate that atp8, atp6, and nad2 showed higher nucleotide diversity than cox1, cox3, and cytb. The intraspecific genetic distances among three different populations of L. boringii exceed 4%, and those between two populations of L. liui reach 7%. Phylogenetic relationships support their division into two subfamilies of Megophryidae (Leptobrachiinae and Megophryinae) as well as two species groups within Leptobrachium, corresponding to the number of keratinized nuptial spines (10-48 in the L. boringii species group vs. 2-6 in the L. liui species group). The two new mitogenomes reported in this study provide valuable data for future molecular evolutionary and conservation studies of the genus Leptobrachium and other Megophryidae toads.

An improved assembly of the "Cascade" hop (Humulus lupulus) genome uncovers signatures of molecular evolution and refines time of divergence estimates for the Cannabaceae family

We present a chromosome-level assembly of the Cascade hop (Humulus lupulus L. var. lupulus) genome. The hop genome is large (2.8 Gb) and complex, and early attempts at assembly were fragmented. Recent advances have made assembly of the hop genome more tractable, transforming the extent of investigation that can occur. The chromosome-level assembly of Cascade was developed by scaffolding the previously reported Cascade assembly generated with PacBio long-read sequencing and polishing with Illumina short-read DNA sequencing. We developed gene models and repeat annotations and used a controlled bi-parental mapping population to identify significant sex-associated markers. We assessed molecular evolution in gene sequences, gene family expansion and contraction, and time of divergence from Cannabis sativa and other closely related plant species using Bayesian inference. We identified the putative sex chromosome in the female genome based on significant sex-associated markers from the bi-parental mapping population. While the estimate of repeat content (~64%) is similar to the estimate for the hemp genome, syntenic blocks in hop contain a greater percentage of LTRs. Hop is enriched for disease resistance-associated genes in syntenic gene blocks and expanded gene families. The Cascade chromosome-level assembly will inform cultivation strategies and serve to deepen our understanding of the hop genomic landscape, benefiting hop researchers and the Cannabaceae genomics community.

Genomic profile of eGFP-expressing canine distemper virus that undergoes serial plaque-to-plaque transfers

Canine distemper virus (CDV) is classified into the genus Morbillivirus in the family Paramyxoviridae. This virus has a single-stranded genomic RNA with negative polarity. The wild-type CDV genome is generally composed of 15 690 nucleotides. We previously rescued an enhanced green fluorescence protein (eGFP)-tagged recombinant CDV (rCDV-eGFP) using reverse genetics. In this study, the rCDV-eGFP at passage-7 was subjected to 38 serial plaque-to-plaque transfers (or bottleneck passages) and two extra common passages in cells. In theory, the effect of Muller’s ratchet may fix deleterious mutations in a single viral population after consecutive plaque-to-plaque transfers. In order to uncover a mutated landscape of the rCDV-eGFP under the circumstances of bottleneck passages, the passage-47 progeny was collected for the in-depth analysis via next-generation sequencing. The result revealed a total of nine single-nucleotide mutations (SNMs) in the viral antigenome. Out of them, SNMs at nt 1832, 5022, 5536, 5580, 5746, 6913 and 8803 were identified as total single-nucleotide substitution, i.e., 100% of mutation frequency. The result suggested no notable formation of viral quasispecies in the rCDV-eGFP population after consecutive plaque-to-plaque transfers.

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 (d_N) to synonymous distance (d_S). Under the assumption that synonymous substitutions in the coding region are selectively neutral, the d_N/d_S ratio can be used to statistically detect the adaptive evolution (or purifying selection) if d_N/d_S > 1 (or d_N/d_S < 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 d_N/d_S < 1. Consequently, the statistical power for testing d_N/d_S = 1 may be insufficient to distinguish between different selection modes. In this paper, we propose a more powerful test, called d_N/d_S-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 d_N/d_S-H model predicts d_N/d_S = 1-H for neutral evolution, d_N/d_S < 1-H for nearly neutral selection, and d_N/d_S > 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 C_N/C_S-H test).

Evolution of γ chain cytokines: mechanisms, methods and applications

The common γ chain family of cytokines and their receptors play fundamental roles in the immune system. Evolutionary studies of γ chain cytokines have elegantly illustrated how the immune system adapts to ever-changing environmental conditions. Indeed, these studies have revealed the uniqueness of cytokine evolution, which exhibits strong positive selection pressure needed to adapt to rapidly evolving threats whilst still conserving their receptor binding capabilities. In this review, we summarise the evolutionary mechanisms that gave rise to the characteristically diverse family of γ chain cytokines. We also speculate on the benefits of studying cytokine evolution, which may provide alternative ways to design novel cytokine therapeutic strategies. Additionally, we discuss current evolutionary models that elucidate the emergence of distinct cytokines (IL-4 and IL-13) and cytokine receptors (IL-2Rα and IL-15Rα). Finally, we address and reflect on the difficulties associated with evolutionary studies of rapidly evolving genes and describe a variety of computational methods that have revealed numerous aspects of cytokine evolution.

Genome-Wide Insights Into the Organelle Translocation of Photosynthetic NDH-1 Genes During Evolution

Translocation of chloroplast-located genes to mitochondria or nucleus is considered to be a safety strategy that impedes mutation of photosynthetic genes and maintains their household function during evolution. The organelle translocation strategy is also developed in photosynthetic NDH-1 (pNDH-1) genes but its understanding is still far from complete. Here, we found that the mutation rate of the conserved pNDH-1 genes was gradually reduced but their selection pressure was maintained at a high level during evolution from cyanobacteria to angiosperm. By contrast, oxygenic photosynthesis-specific (OPS) pNDH-1 genes had an opposite trend, explaining the reason why they were transferred from the reactive oxygen species (ROS)-enriched chloroplast to the ROS-barren nucleus. Further, genome-wide sequence analysis supported the possibility that all conserved pNDH-1 genes lost in chloroplast genomes of Chlorophyceae and Pinaceae were transferred to the ROS-less mitochondrial genome as deduced from their truncated pNDH-1 gene fragments. Collectively, we propose that the organelle translocation strategy of pNDH-1 genes during evolution is necessary to maintain the function of the pNDH-1 complex as an important antioxidant mechanism for efficient photosynthesis.

Novel Salmonella Phage, vB_Sen_STGO-35-1, Characterization and Evaluation in Chicken Meat

Salmonellosis is one of the most frequently reported zoonotic foodborne diseases worldwide, and poultry is the most important reservoir of Salmonella enterica serovar Enteritidis. The use of lytic bacteriophages (phages) to reduce foodborne pathogens has emerged as a promising biocontrol intervention for Salmonella spp. Here, we describe and evaluate the newly isolated Salmonella phage STGO-35-1, including: (i) genomic and phenotypic characterization, (ii) an analysis of the reduction of Salmonella in chicken meat, and (iii) genome plasticity testing. Phage STGO-35-1 represents an unclassified siphovirus, with a length of 47,483 bp, a G + C content of 46.5%, a headful strategy of packaging, and a virulent lifestyle. Phage STGO-35-1 reduced S. Enteritidis counts in chicken meat by 2.5 orders of magnitude at 4 °C. We identified two receptor-binding proteins with affinity to LPS, and their encoding genes showed plasticity during an exposure assay. Phenotypic, proteomic, and genomic characteristics of STGO-35-1, as well as the Salmonella reduction in chicken meat, support the potential use of STGO-35-1 as a targeted biocontrol agent against S. Enteritidis in chicken meat. Additionally, computational analysis and a short exposure time assay allowed us to predict the plasticity of genes encoding putative receptor-binding proteins.

Genome-wide identification and analysis of WRKY gene family in maize provide insights into regulatory network in response to abiotic stresses

BACKGROUND

The WRKY transcription factor family plays significant roles in biotic and abiotic stress responses, which has been associated with various biological processes in higher plants. However, very little is known regarding the structure and function of WRKY genes in maize.

RESULTS

In this study, a total of 140 ZmWRKY proteins encoded by 125 ZmWRKY genes were eventually identified in maize. On the basis of features of molecular structure and a comparison of phylogenetic relationships of WRKY transcription factor families from Arabidopsis, rice and maize, all 140 ZmWRKY proteins in maize were divided into three main groups (Groups I, II and III) and the Group II was further classified into five subgroups. The characteristics of exon-intron structure of these putative ZmWRKY genes and conserved protein motifs of their encoded ZmWRKY proteins were also presented respectively, which was in accordance with the group classification results. Promoter analysis suggested that ZmWRKY genes shared many abiotic stress-related elements and hormone-related elements. Gene duplication analysis revealed that the segmental duplication and purifying selection might play a significant role during the evolution of the WRKY gene family in maize. Using RNA-seq data, transcriptome analysis indicated that most of ZmWRKY genes displayed differential expression patterns at different developmental stages of maize. Further, by quantitative real-time PCR analysis, twenty-one ZmWRKY genes were confirmed to respond to two different abiotic stress treatments, suggesting their potential roles in various abiotic stress responses. In addition, RNA-seq dataset was used to conduct weighted gene co-expression network analysis (WGCNA) in order to recognize gene subsets possessing similar expression patterns and highly correlated with each other within different metabolic networks. Further, subcellular localization prediction, functional annotation and interaction analysis of ZmWRKY proteins were also performed to predict their interactions and associations involved in potential regulatory network.

CONCLUSIONS

Taken together, the present study will serve to present an important theoretical basis for further exploring function and regulatory mechanism of ZmWRKY genes in the growth, development, and adaptation to abiotic stresses in maize.

The mitochondrial genome and phylogenetic characteristics of the Thick-billed Green-Pigeon, Treron curvirostra: the first sequence for the genus

Members of the genus Treron (Columbidae) are widely distributed in southern Asia and the Indo-Malayan Region but their relationships are poorly understood. Better knowledge of the systematic status of this genus may help studies of historical biogeography and taxonomy. The complete mitochondrial genome of T.curvirostra was characterized, a first for the genus. It is 17,414 base pairs in length, containing two rRNAs, 22 tRNAs, 13 protein coding genes (PCGs), and one D-loop with a primary structure that is similar to that found in most members of Columbidae. Most PCGs start with the common ATG codon but are terminated by different codons. The highest value of the Ka/Ks ratio within 13 PCGs was found in ATP8 with 0.1937, suggesting that PCGs of the mitochondrial genome tend to be conservative in Columbidae. Moreover, the phylogenetic relationships within Columbidae, which was based on sequences of 13 PCGs, showed that (T.curvirostra + Hemiphaganovaeseelandiae) were clustered in one clade, suggesting a potentially close relationship between Treron and Hemiphaga. However, the monophyly of the subfamilies of Columbidae recognized by the Interagency Taxonomic Information System could not be corroborated. Hence, the position of the genus Treron in the classification of Columbidae may have to be revised.

Combining liquid biopsy and radiomics for personalized treatment of lung cancer patients. State of the art and new perspectives

Lung cancer has become a paradigm for precision medicine in oncology, and liquid biopsy (LB) together with radiomics may have a great potential in this scenario. They are both minimally invasive, easy to perform, and can be repeated during patient's follow-up. Also, increasing evidence suggest that LB and radiomics may provide an efficient way to screen and diagnose tumors at an early stage, including the monitoring of any change in the tumor molecular profile. This could allow treatment optimization, improvement of patients' quality of life, and healthcare-related costs reduction. Latest reports on lung cancer patients suggest a combination of these two strategies, along with cutting-edge data analysis, to decode valuable information regarding tumor type, aggressiveness, progression, and response to treatment. The approach seems more compatible with clinical practice than the current standard, and provides new diagnostic companions being able to suggest the best treatment strategy compared to conventional methods. To implement radiomics and liquid biopsy directly into clinical practice, an artificial intelligence (AI)-based system could help to link patients' clinical data together with tumor molecular profiles and imaging characteristics. AI could also solve problems and limitations related to LB and radiomics methodologies. Further work is needed, including new health policies and the access to large amounts of high-quality and well-organized data, allowing a complementary and synergistic combination of LB and imaging, to provide an attractive choice e in the personalized treatment of lung cancer.

Long-Lived Species of Bivalves Exhibit Low MT-DNA Substitution Rates

Bivalves represent valuable taxonomic group for aging studies given their wide variation in longevity (from 1–2 to >500 years). It is well known that aging is associated to the maintenance of Reactive Oxygen Species homeostasis and that mitochondria phenotype and genotype dysfunctions accumulation is a hallmark of these processes. Previous studies have shown that mitochondrial DNA mutation rates are linked to lifespan in vertebrate species, but no study has explored this in invertebrates. To this end, we performed a Bayesian Phylogenetic Covariance model of evolution analysis using 12 mitochondrial protein-coding genes of 76 bivalve species. Three life history traits (maximum longevity, generation time and mean temperature tolerance) were tested against 1) synonymous substitution rates (dS), 2) conservative amino acid replacement rates (Kc) and 3) ratios of radical over conservative amino acid replacement rates (Kr/Kc). Our results confirm the already known correlation between longevity and generation time and show, for the first time in an invertebrate class, a significant negative correlation between dS and longevity. This correlation was not as strong when generation time and mean temperature tolerance variations were also considered in our model (marginal correlation), suggesting a confounding effect of these traits on the relationship between longevity and mtDNA substitution rate. By confirming the negative correlation between dS and longevity previously documented in birds and mammals, our results provide support for a general pattern in substitution rates.

Analysis of selection in protein-coding sequences accounting for common biases

The evolution of protein-coding genes is usually driven by selective processes, which favor some evolutionary trajectories over others, optimizing the subsequent protein stability and activity. The analysis of selection in this type of genetic data is broadly performed with the metric nonsynonymous/synonymous substitution rate ratio (dN/dS). However, most of the well-established methodologies to estimate this metric make crucial assumptions, such as lack of recombination or invariable codon frequencies along genes, which can bias the estimation. Here, we review the most relevant biases in the dN/dS estimation and provide a detailed guide to estimate this metric using state-of-the-art procedures that account for such biases, along with illustrative practical examples and recommendations. We also discuss the traditional interpretation of the estimated dN/dS emphasizing the importance of considering complementary biological information such as the role of the observed substitutions on the stability and function of proteins. This review is oriented to help evolutionary biologists that aim to accurately estimate selection in protein-coding sequences.

Complete Genome Sequence of a Serotype 7 Listeria monocytogenes Strain, FSL R9-0915

Listeria monocytogenes serotype 7 lacks glycosidic constituents in wall teichoic acids. Here, we present the complete genome sequence of L. monocytogenes serotype 7 strain FSL R9-0915 and an analysis of genes known to affect L. monocytogenes antigenicity. This strain is used as a control strain in Listeria phage host range analyses.

Oxytocin amino acid variation within Neotropical primates: new genetic variants in hormone and receptor sequences and evidence for evolutionary forces driving this unexpected diversity

Oxytocin is a mammalian neuropeptide hormone that mediates behaviours important to reproduction. Despite almost universal amino acid sequence conservation across most groups of mammals, several unique forms have been reported across Neotropical primates. To explore sequence diversity, we investigated the genes encoding oxytocin and its receptor across the Atelidae, which was known to contain at least three unique oxytocin sequences. Additionally, we included the genus Cebus, within the Cebidae, to further explore the ubiquity of the Pro8 variant in this family. We found a novel amino acid variant (Val3) within the Atelidae radiation, bringing the total number of oxytocin sequences within Neotropical primates to seven. Analyses of physicochemical properties revealed conservative substitutions that are likely tolerated within the selective constraints imposed by receptor binding. Furthermore, we report radical substitutions at the eighth codon and evidence for co-evolution between Pro8 and a ligand-binding region of the oxytocin receptor in the Atelidae, supporting the notion that this variant may affect binding specificity. Overall, we suggest that selective constraint on binding specificity may maintain proper oxytocin function and that the diversification of amino acid sequence is likely due to a variety of processes such as relaxed constraint, neutral mutation, positive selection and coevolution.

Molecular Evolution in Small Steps under Prevailing Negative Selection: A Nearly Universal Rule of Codon Substitution

The widely accepted view that evolution proceeds in small steps is based on two premises: 1) negative selection acts strongly against large differences and 2) positive selection favors small-step changes. The two premises are not biologically connected and should be evaluated separately. We now extend a previous approach to studying codon evolution in the entire genome. Codon substitution rate is a function of the physicochemical distance between amino acids (AAs), equated with the step size of evolution. Between nine pairs of closely related species of plants, invertebrates, and vertebrates, the evolutionary rate is strongly and negatively correlated with a set of AA distances (Δ_U, scaled to [0, 1]). Δ_U, a composite measure of evolutionary rates across diverse taxa, is influenced by almost all of the 48 physicochemical properties used here. The new analyses reveal a crucial trend hidden from previous studies: Δ_U is strongly correlated with the evolutionary rate (R² > 0.8) only when the genes are predominantly under negative selection. Because most genes in most taxa are strongly constrained by negative selection, Δ_U has indeed appeared to be a nearly universal measure of codon evolution. In conclusion, molecular evolution at the codon level generally takes small steps due to the prevailing negative selection. Whether positive selection may, or may not, follow the small-step rule is addressed in a companion study.

Molecular Evolution in Large Steps-Codon Substitutions under Positive Selection

Molecular evolution is believed to proceed in small steps. The step size can be defined by a distance reflecting physico-chemical disparities between amino acid (AA) pairs that can be exchanged by single 1-bp mutations. We show that AA substitution rates are strongly and negatively correlated with this distance but only when positive selection is relatively weak. We use the McDonald and Kreitman test to separate the influences of positive and negative selection. While negative selection is indeed stronger on AA substitutions generating larger changes in chemical properties of AAs, positive selection operates by different rules. For 65 of the 75 possible pairs, positive selection is comparable in strength regardless of AA distance. However, the ten pairs under the strongest positive selection all exhibit large leaps in chemical properties. Five of the ten pairs are shared between Drosophila and Hominoids, thus hinting at a common but modest biochemical basis of adaptation across taxa. The hypothesis that adaptive changes often take large functional steps will need to be extensively tested. If validated, molecular models will need to better integrate positive and negative selection in the search for adaptive signal.

An evolutionary model motivated by physicochemical properties of amino acids reveals variation among proteins

Motivation

The relative rates of amino acid interchanges over evolutionary time are likely to vary among proteins. Variation in those rates has the potential to reveal information about constraints on proteins. However, the most straightforward model that could be used to estimate relative rates of amino acid substitution is parameter-rich and it is therefore impractical to use for this purpose.

Results

A six-parameter model of amino acid substitution that incorporates information about the physicochemical properties of amino acids was developed. It showed that amino acid side chain volume, polarity and aromaticity have major impacts on protein evolution. It also revealed variation among proteins in the relative importance of those properties. The same general approach can be used to improve the fit of empirical models such as the commonly used PAM and LG models.

Availability and implementation

Perl code and test data are available from https://github.com/ebraun68/sixparam.

Supplementary information

Supplementary data are available at Bioinformatics online.

Two Phages of the Genera Felixounavirus Subjected to 12 Hour Challenge on Salmonella Infantis Showed Distinct Genotypic and Phenotypic Changes

Salmonella Infantis is considered in recent years an emerging Salmonella serovar, as it has been associated with several outbreaks and multidrug resistance phenotypes. Phages appear as a possible alternative strategy to control Salmonella Infantis (SI). The aims of this work were to characterize two phages of the Felixounavirus genus, isolated using the same strain of SI, and to expose them to interact in challenge assays to identify genetic and phenotypic changes generated from these interactions. These two phages have a shared nucleotide identity of 97% and are differentiated by their host range: one phage has a wide host range (lysing 14 serovars), and the other has a narrow host range (lysing 6 serovars). During the 12 h challenge we compared: (1) optical density of SI, (2) proportion of SI survivors from phage-infected cultures, and (3) phage titer. Isolates obtained through the assays were evaluated by efficiency of plating (EOP) and by host-range characterization. Genomic modifications were characterized by evaluation of single nucleotide polymorphisms (SNPs). The optical density (600 nm) of phage-infected SI decreased, as compared to the uninfected control, by an average of 0.7 for SI infected with the wide-host-range (WHR) phage and by 0.3 for SI infected with the narrow-host-range (NHR) phage. WHR phage reached higher phage titer (7 × 1011 PFU/mL), and a lower proportion of SI survivor was obtained from the challenge assay. In SI that interacted with phages, we identified SNPs in two genes (rfaK and rfaB), which are both involved in lipopolysaccharide (LPS) polymerization. Therefore, mutations that could impact potential phage receptors on the host surface were selected by lytic phage exposure. This work demonstrates that the interaction of Salmonella phages (WHR and NHR) with SI for 12 h in vitro leads to emergence of new phenotypic and genotypic traits in both phage and host. This information is crucial for the rational design of phage-based control strategies.

Mapping Selection within Drosophila melanogaster Embryo's Anatomy

We present a survey of selection across Drosophila melanogaster embryonic anatomy. Our approach integrates genomic variation, spatial gene expression patterns, and development with the aim of mapping adaptation over the entire embryo's anatomy. Our adaptation map is based on analyzing spatial gene expression information for 5,969 genes (from text-based annotations of in situ hybridization data directly from the BDGP database, Tomancak et al. 2007) and the polymorphism and divergence in these genes (from the project DGRP, Mackay et al. 2012).The proportion of nonsynonymous substitutions that are adaptive, neutral, or slightly deleterious are estimated for the set of genes expressed in each embryonic anatomical structure using the distribution of fitness effects-alpha method (Eyre-Walker and Keightley 2009). This method is a robust derivative of the McDonald and Kreitman test (McDonald and Kreitman 1991). We also explore whether different anatomical structures differ in the phylogenetic age, codon usage, or expression bias of the genes they express and whether genes expressed in many anatomical structures show more adaptive substitutions than other genes.We found that: 1) most of the digestive system and ectoderm-derived structures are under selective constraint, 2) the germ line and some specific mesoderm-derived structures show high rates of adaptive substitution, and 3) the genes that are expressed in a small number of anatomical structures show higher expression bias, lower phylogenetic ages, and less constraint.

Extreme Deviations from Expected Evolutionary Rates in Archaeal Protein Families

Origin of new biological functions is a complex phenomenon ranging from single-nucleotide substitutions to the gain of new genes via horizontal gene transfer or duplication. Neofunctionalization and subfunctionalization of proteins is often attributed to the emergence of paralogs that are subject to relaxed purifying selection or positive selection and thus evolve at accelerated rates. Such phenomena potentially could be detected as anomalies in the phylogenies of the respective gene families. We developed a computational pipeline to search for such anomalies in 1,834 orthologous clusters of archaeal genes, focusing on lineage-specific subfamilies that significantly deviate from the expected rate of evolution. Multiple potential cases of neofunctionalization and subfunctionalization were identified, including some ancient, house-keeping gene families, such as ribosomal protein S10, general transcription factor TFIIB and chaperone Hsp20. As expected, many cases of apparent acceleration of evolution are associated with lineage-specific gene duplication. On other occasions, long branches in phylogenetic trees correspond to horizontal gene transfer across long evolutionary distances. Significant deceleration of evolution is less common than acceleration, and the underlying causes are not well understood; functional shifts accompanied by increased constraints could be involved. Many gene families appear to be “highly evolvable,” that is, include both long and short branches. Even in the absence of precise functional predictions, this approach allows one to select targets for experimentation in search of new biology.

Developmental Constraints on Genome Evolution in Four Bilaterian Model Species

Developmental constraints on genome evolution have been suggested to follow either an early conservation model or an "hourglass" model. Both models agree that late development strongly diverges between species, but debate on which developmental period is the most conserved. Here, based on a modified "Transcriptome Age Index" approach, that is, weighting trait measures by expression level, we analyzed the constraints acting on three evolutionary traits of protein coding genes (strength of purifying selection on protein sequences, phyletic age, and duplicability) in four species: Nematode worm Caenorhabditis elegans, fly Drosophila melanogaster, zebrafish Danio rerio, and mouse Mus musculus. In general, we found that both models can be supported by different genomic properties. Sequence evolution follows an hourglass model, but the evolution of phyletic age and of duplicability follow an early conservation model. Further analyses indicate that stronger purifying selection on sequences in the middle development are driven by temporal pleiotropy of these genes. In addition, we report evidence that expression in late development is enriched with retrogenes, which usually lack efficient regulatory elements. This implies that expression in late development could facilitate transcription of new genes, and provide opportunities for acquisition of function. Finally, in C. elegans, we suggest that dosage imbalance could be one of the main factors that cause depleted expression of high duplicability genes in early development.

Radical amino acid mutations persist longer in the absence of sex

Harmful mutations are ubiquitous and inevitable, and the rate at which these mutations are removed from populations is a critical determinant of evolutionary fate. Closely related sexual and asexual taxa provide a particularly powerful setting to study deleterious mutation elimination because sexual reproduction should facilitate mutational clearance by reducing selective interference between sites and by allowing the production of offspring with different mutational complements than their parents. Here, we compared the rate of removal of conservative (i.e., similar biochemical properties) and radical (i.e., distinct biochemical properties) nonsynonymous mutations from mitochondrial genomes of sexual versus asexual Potamopyrgus antipodarum, a New Zealand freshwater snail characterized by coexisting and ecologically similar sexual and asexual lineages. Our analyses revealed that radical nonsynonymous mutations are cleared at higher rates than conservative changes and that sexual lineages eliminate radical changes more rapidly than asexual counterparts. These results are consistent with reduced efficacy of purifying selection in asexual lineages allowing harmful mutations to remain polymorphic longer than in sexual lineages. Together, these data illuminate some of the population-level processes contributing to mitochondrial mutation accumulation and suggest that mutation accumulation could influence the outcome of competition between sexual and asexual lineages.

Population size may shape the accumulation of functional mutations following domestication

BACKGROUND

Population genetics theory predicts an important role of differences in the effective population size (N e ) among species on shaping the accumulation of functional mutations by regulating the selection efficiency. However, this correlation has never been tested in domesticated animals.

RESULTS

Here, we synthesized 62 whole genome data in eight domesticated species (cat, dog, pig, goat, sheep, chicken, cattle and horse) and compared domesticates with their wild (or ancient) relatives. Genes with significantly different selection pressures (revealed by nonsynonymous/synonymous substitution rate ratios, Ka/Ks or ω) between domesticated (Dω) and wild animals (Wω) were determined by likelihood-ratio tests. Species-level effective population sizes (N e ) were evaluated by the pairwise sequentially Markovian coalescent (PSMC) model, and Dω/Wω were calculated for each species to evaluate the changes in accumulation of functional mutations after domestication relative to pre-domestication period. Correlation analysis revealed that the most recent (~ 10.000 years ago) N e (s) are positively correlated with Dω/Wω. This result is consistent with the corollary of the nearly neutral theory, that higher N e could boost the efficiency of positive selection, which might facilitate the overall accumulation of functional mutations. In addition, we also evaluated the accumulation of radical and conservative mutations during the domestication transition as: Dradical/Wradical and Dconservative/Wconservative, respectively. Surprisingly, only Dradical/Wradical ratio exhibited a positive correlation with N e (p < 0.05), suggesting that domestication process might magnify the accumulation of radical mutations in species with larger N e .

CONCLUSIONS

Our results confirm the classical population genetics theory prediction and highlight the important role of species' N e in shaping the patterns of accumulation of functional mutations, especially radical mutations, in domesticated animals. The results aid our understanding of the mechanisms underlying the accumulation of functional mutations after domestication, which is critical for understanding the phenotypic diversification associated with this process.

Integrated transcriptomic analysis of Trichosporon Asahii uncovers the core genes and pathways of fluconazole resistance

Trichosporon asahii (T. asahii) has emerged as a dangerous pathogen that causes rare but life-threatening infections. Its resistance to certain antifungal agents makes it difficult to treat, especially for patients undergoing long-term antibiotic therapy. In this study, we performed a series of fluconazole (FLC) perturbation experiments for two T. asahii strains, a clinical isolate stain CBS 2479 (T2) and an environmental isolate strain CBS 8904 (T8), to uncover potential genes and pathways involved in FLC resistance. We achieved 10 transcriptomes of T2 and T8 that were based on dose and time series of FLC perturbations. Systematic comparisons of the transcriptomes revealed 32 T2 genes and 25 T8 genes that are highly sensitive to different FLC perturbations. In both T2 and T8 strains with the phenotype of FLC resistance, the processes of oxidation-reduction and transmembrane transport were detected to be significantly changed. The antifungal susceptibility testing of FLC and penicillin revealed their resistance pathways are merged. Accumulated mutations were found in 564 T2 and 225 T8 genes, including four highly mutated genes that are functionally related to the target of rapamycin complex (TOR). Our study provides abundant data towards genome-wide understanding of the molecular basis of FLC resistance in T. asahii.

Mitochondrial dual-coding genes in Trypanosoma brucei

Trypanosoma brucei is transmitted between mammalian hosts by the tsetse fly. In the mammal, they are exclusively extracellular, continuously replicating within the bloodstream. During this stage, the mitochondrion lacks a functional electron transport chain (ETC). Successful transition to the fly, requires activation of the ETC and ATP synthesis via oxidative phosphorylation. This life cycle leads to a major problem: in the bloodstream, the mitochondrial genes are not under selection and are subject to genetic drift that endangers their integrity. Exacerbating this, T. brucei undergoes repeated population bottlenecks as they evade the host immune system that would create additional forces of genetic drift. These parasites possess several unique genetic features, including RNA editing of mitochondrial transcripts. RNA editing creates open reading frames by the guided insertion and deletion of U-residues within the mRNA. A major question in the field has been why this metabolically expensive system of RNA editing would evolve and persist. Here, we show that many of the edited mRNAs can alter the choice of start codon and the open reading frame by alternative editing of the 5’ end. Analyses of mutational bias indicate that six of the mitochondrial genes may be dual-coding and that RNA editing allows access to both reading frames. We hypothesize that dual-coding genes can protect genetic information by essentially hiding a non-selected gene within one that remains under selection. Thus, the complex RNA editing system found in the mitochondria of trypanosomes provides a unique molecular strategy to combat genetic drift in non-selective conditions.

In African trypanosomes, many of the mitochondrial mRNAs require extensive RNA editing before they can be translated. During this process, each edited transcript can undergo hundreds of cleavage/ligation events as U-residues are inserted or deleted to generate a translatable open reading frame. A major paradox has been why this incredibly metabolically expensive process would evolve and persist. In this work, we show that many of the mitochondrial genes in trypanosomes are dual-coding, utilizing different reading frames to potentially produce two very different proteins. Access to both reading frames is made possible by alternative editing of the 5’ end of the transcript. We hypothesize that dual-coding genes may work to protect the mitochondrial genes from mutations during growth in the mammalian host, when many of the mitochondrial genes are not being used. Thus, the complex RNA editing system may be maintained because it provides a unique molecular strategy to combat genetic drift.

Limited locomotive ability relaxed selective constraints on molluscs mitochondrial genomes

Mollusca are the second largest phylum in the animal kingdom with different types of locomotion. Some molluscs are poor-migrating, while others are free-moving or fast-swimming. Most of the energy required for locomotion is provided by mitochondria via oxidative phosphorylation. Here, we conduct a comparative genomic analysis of 256 molluscs complete mitochondrial genomes and evaluate the role of energetic functional constraints on the protein-coding genes, providing a new insight into mitochondrial DNA (mtDNA) evolution. The weakly locomotive molluscs, compared to strongly locomotive molluscs, show significantly higher Ka/Ks ratio, which suggest they accumulated more nonsynonymous mutations in mtDNA and have experienced more relaxed evolutionary constraints. Eleven protein-coding genes (CoxI, CoxII, ATP6, Cytb, ND1-6, ND4L) show significant difference for Ka/Ks ratios between the strongly and weakly locomotive groups. The relaxation of selective constraints on Atp8 arise in the common ancestor of bivalves, and the further relaxation occurred in marine bivalves lineage. Our study thus demonstrates that selective constraints relevant to locomotive ability play an essential role in evolution of molluscs mtDNA.

Genome-wide analysis of the CCCH zinc finger family identifies tissue specific and stress responsive candidates in chickpea (Cicer arietinum L.)

The CCCH zinc finger is a group of proteins characterised by a typical motif consisting of three cysteine residues and one histidine residue. These proteins have been reported to play important roles in regulation of plant growth, developmental processes and environmental responses. In the present study, genome wide analysis of the CCCH zinc finger gene family was carried out in the available chickpea genome. Various bioinformatics tools were employed to predict 58 CCCH zinc finger genes in chickpea (designated CarC3H1-58), which were analysed for their physio-chemical properties. Phylogenetic analysis classified the proteins into 12 groups in which members of a particular group had similar structural organization. Further, the numbers as well as the types of CCCH motifs present in the CarC3H proteins were compared with those from Arabidopsis and Medicago truncatula. Synteny analysis revealed valuable information regarding the evolution of this gene family. Tandem and segmental duplication events were identified and their Ka/Ks values revealed that the CarC3H gene family in chickpea had undergone purifying selection. Digital, as well as real time qRT-PCR expression analysis was performed which helped in identification of several CarC3H members that expressed preferentially in specific chickpea tissues as well as during abiotic stresses (desiccation, cold, salinity). Moreover, molecular characterization of an important member CarC3H45 was carried out. This study provides comprehensive genomic information about the important CCCH zinc finger gene family in chickpea. The identified tissue specific and abiotic stress specific CCCH genes could be potential candidates for further characterization to delineate their functional roles in development and stress.

The impacts of drift and selection on genomic evolution in insects

Genomes evolve through a combination of mutation, drift, and selection, all of which act heterogeneously across genes and lineages. This leads to differences in branch-length patterns among gene trees. Genes that yield trees with the same branch-length patterns can be grouped together into clusters. Here, we propose a novel phylogenetic approach to explain the factors that influence the number and distribution of these gene-tree clusters. We apply our method to a genomic dataset from insects, an ancient and diverse group of organisms. We find some evidence that when drift is the dominant evolutionary process, each cluster tends to contain a large number of fast-evolving genes. In contrast, strong negative selection leads to many distinct clusters, each of which contains only a few slow-evolving genes. Our work, although preliminary in nature, illustrates the use of phylogenetic methods to shed light on the factors driving rate variation in genomic evolution.

Incorporating tree-thinking and evolutionary time scale into developmental biology

Phylogenetic approaches are indispensable in any comparative molecular study involving multiple species. These approaches are in increasing demand as the amount and availability of DNA sequence information continues to increase exponentially, even for organisms that were previously not extensively studied. Without the sound application of phylogenetic concepts and knowledge, one can be misled when attempting to infer ancestral character states as well as the timing and order of evolutionary events, both of which are frequently exerted in evolutionary developmental biology. The ignorance of phylogenetic approaches can also impact non-evolutionary studies and cause misidentification of the target gene or protein to be examined in functional characterization. This review aims to promote tree-thinking in evolutionary conjecture and stress the importance of a sense of time scale in cross-species comparisons, in order to enhance the understanding of phylogenetics in all biological fields including developmental biology. To this end, molecular phylogenies of several developmental regulatory genes, including those denoted as "cryptic pan-vertebrate genes", are introduced as examples.

Studies of nosocomial outbreaks of hepatitis B in nursing homes in Germany suggest a major role of hepatitis B e antigen expression in disease severity and progression

Hepatitis B virus (HBV) causes acute or chronic hepatitis B. Local outbreaks of HBV infections in skilled nursing facilities is a matter of growing concern in developed countries. Here, we investigated two outbreaks of hepatitis B that recently occurred in nursing homes in Germany. The outbreak at location A was associated with acute fulminant hepatitis with fatal outcome in several cases, while individuals infected at location B developed asymptomatic or mild hepatitis B. Sequence analysis of viruses involved in these outbreaks revealed different, but unique HBV strains for each location. Each of the strains produced high viremia of more than 10(9) virions/mL serum. We found that the mild course of hepatitis B at location B was caused by a circulating wild-type HBV genotype A2 strain, which is commonly found in Central Europe. Complete genome sequences of isolates obtained from infected patients revealed nearly 100% sequence identity at the nucleotide level as well as expression of HBV e protein (HBeAg), a known T cell tolerogen in the incubation or chronic phases of HBV infection. By contrast, the outbreak at location A was associated with an HBV genotype D2 variant that lacked HBeAg expression, suggesting that immunopathology and selection of specific HBV variants played a major role in the severe (or even fulminant) acute hepatitis observed at location A. Importantly, all patients were diagnosed with type 2 diabetes mellitus, a known risk factor for healthcare-associated transmission of HBV. The study leads us to suggest that, besides strict adherence to hygiene standards, additional efforts are required to reduce the risk of HBV transmission and fulminant disease progression in healthcare settings and nursing homes. In this context, a general screening for HBsAg and active hepatitis B vaccination should be considered for people living in nursing homes, especially for those with diagnosed diabetes or other predisposing factors for HBV transmission.

Genome Biology and the Evolution of Cell-Size Diversity

Cell size is highly variable among different species across the Tree of Life. For decades, biologists have generated hypotheses to explain this variation, in many cases, drawing on the correlations that exist among cell size, genome size, nucleus size, and various physiological and developmental parameters. In recent years, our understanding of the molecular processes that generate variation in genome size over evolutionary time, as well as the processes that maintain homeostasis in cell size over ontogenetic time, has increased dramatically. The goal of this article is to highlight how information from these fields can be integrated to generate new hypotheses to explain cell-size diversity.

Sequence diversity and molecular evolutionary rates between buffalo and cattle

Identification of genes of importance regarding production traits in buffalo is impaired by a paucity of genomic resources. Choice to fill this gap is to exploit data available for cow. The cross-species application of comparative genomics tools is potential gear to investigate the buffalo genome. However, this is dependent on nucleotide sequences similarity. In this study, gene diversity between buffalo and cattle was determined using 86 gene orthologues. There was approximately 3% difference in all genes in terms of nucleotide diversity and 0.267 ± 0.134 in amino acids, indicating the possibility for successfully using cross-species strategies for genomic studies. There were significantly higher non-synonymous substitutions both in cattle and buffalo; however, there was similar difference in terms of dN- dS (4.414 versus 4.745) in buffalo and cattle, respectively. Higher rate of non-synonymous substitutions at similar level in buffalo and cattle indicated a similar positive selection pressure. Results for relative rate test were assessed with the chi-squared test. There was no significance difference on unique mutations between cattle and buffalo lineages at synonymous sites. However, there was a significance difference on unique mutations for non-synonymous sites, indicating ongoing mutagenic process that generates substitutional mutation at approximately the same rate at silent sites. Moreover, despite of common ancestry, our results indicate a different divergent time among genes of cattle and buffalo. This is the first demonstration that variable rates of molecular evolution may be present within the family Bovidae.

Subgenotypes and mutations in the s and polymerase genes of hepatitis B virus carriers in the West Bank, palestine

The mutation rate and genetic variability of hepatitis B virus (HBV) are crucial factors for efficient treatment and successful vaccination against HBV. Until today, genetic properties of this virus among the Palestinian population remain unknown. Therefore, we performed genetic analysis of the overlapping S and polymerase genes of HBV, isolated from 40 Palestinian patients' sera. All patients were HBsAg positive and presented with a viral load above 105 HBV genome copies/ml. The genotyping results of the S gene demonstrated that HBV D1 was detected in 90% of the samples representing the most prominent subgenotype among Palestinians carrying HBV. Various mutations existed within the S gene; in five patients four known escape mutations including the common G145R and D144E were found. Furthermore, a ratio of 4.25 of non-synonymous to synonymous mutations in the S gene indicated a strong selection pressure on the HBs antigen loops of HBV strains circulating in those Palestinian patients. Although all patients were treatment-naïve, with the exception of one, several mutations were found in the HBV polymerase gene, but none pointed to drug resistance. The study presented here is the first report to address subgenotypes and mutation analyses of HBV S and polymerase genes in Palestine.

High prevalence of non-synonymous substitutions in mtDNA of cichlid fishes from Lake Victoria

When a population size is reduced, genetic drift may fix slightly deleterious mutations, and an increase in nonsynonymous substitution is expected. It has been suggested that past aridity has seriously affected and decreased the populations of cichlid fishes in Lake Victoria, while geographical studies have shown that the water levels in Lake Tanganyika and Lake Malawi have remained fairly constant. The comparably stable environments in the latter two lakes might have kept the populations of cichlid fishes large enough to remove slightly deleterious mutations. The difference in the stability of cichlid fish population sizes between Lake Victoria and the Lakes Tanganyika and Malawi is expected to have caused differences in the nonsynonymous/synonymous ratio, ω (=dN/dS), of the evolutionary rate. Here, we estimated ω and compared it between the cichlids of the three lakes for 13 mitochondrial protein-coding genes using maximum likelihood methods. We found that the lineages of the cichlids in Lake Victoria had a significantly higher ω for several mitochondrial loci. Moreover, positive selection was indicated for several codons in the mtDNA of the Lake Victoria cichlid lineage. Our results indicate that both adaptive and slightly deleterious molecular evolution has taken place in the Lake Victoria cichlids' mtDNA genes, whose nonsynonymous sites are generally conserved.

Genomic characterization of coxsackievirus type A24 strains associated with acute flaccid paralysis and rarely identified Hopkins syndrome

The full-length genome sequence analysis of four coxsackievirus A24 (CV-A24) strains, detected in three paralytic and one post-asthmatic paralytic (Hopkins syndrome) cases, is reported here for the first time. A phylogenetic tree constructed on the basis of entire genomes displayed topology similar to that of the full-VP1 tree, classifying the study strains in genogroup CV-A24vGIV along with their temporal counterparts in strains from non-paralytic cases. The strains of the study formed a single genetic cluster C4 within CV-A24vGIV and showed 3.5-19.4 % nucleotide sequence divergence, with 2-4 novel nucleotide mutations in the 5'NCR and 3-8 unique amino acid substitutions in the polyprotein, with respect to the CV-A24 strains associated with non-paralytic cases. Among the nucleotide mutations, A299U was identified in the 5'NCRs of all of the study strains. CV-A24v strains of the same genogroup with few genomic variations but different disease manifestations need to be explored to investigate the molecular basis of evolution of neurovirulence.

Efficient purging of deleterious mutations in plants with haploid selfing

In diploid organisms, selfing reduces the efficiency of selection in removing deleterious mutations from a population. This need not be the case for all organisms. Some plants, for example, undergo an extreme form of selfing known as intragametophytic selfing, which immediately exposes all recessive deleterious mutations in a parental genome to selective purging. Here, we ask how effectively deleterious mutations are removed from such plants. Specifically, we study the extent to which deleterious mutations accumulate in a predominantly selfing and a predominantly outcrossing pair of moss species, using genome-wide transcriptome data. We find that the selfing species purge significantly more nonsynonymous mutations, as well as a greater proportion of radical amino acid changes which alter physicochemical properties of amino acids. Moreover, their purging of deleterious mutation is especially strong in conserved regions of protein-coding genes. Our observations show that selfing need not impede but can even accelerate the removal of deleterious mutations, and do so on a genome-wide scale.

Urotensin-II system in genetic control of blood pressure and renal function

Urotensin-II controls ion/water homeostasis in fish and vascular tone in rodents. We hypothesised that common genetic variants in urotensin-II pathway genes are associated with human blood pressure or renal function. We performed family-based analysis of association between blood pressure, glomerular filtration and genes of the urotensin-II pathway (urotensin-II, urotensin-II related peptide, urotensin-II receptor) saturated with 28 tagging single nucleotide polymorphisms in 2024 individuals from 520 families; followed by an independent replication in 420 families and 7545 unrelated subjects. The expression studies of the urotensin-II pathway were carried out in 97 human kidneys. Phylogenetic evolutionary analysis was conducted in 17 vertebrate species. One single nucleotide polymorphism (rs531485 in urotensin-II gene) was associated with adjusted estimated glomerular filtration rate in the discovery cohort (p = 0.0005). It showed no association with estimated glomerular filtration rate in the combined replication resource of 8724 subjects from 6 populations. Expression of urotensin-II and its receptor showed strong linear correlation (r = 0.86, p<0.0001). There was no difference in renal expression of urotensin-II system between hypertensive and normotensive subjects. Evolutionary analysis revealed accumulation of mutations in urotensin-II since the divergence of primates and weaker conservation of urotensin-II receptor in primates than in lower vertebrates. Our data suggest that urotensin-II system genes are unlikely to play a major role in genetic control of human blood pressure or renal function. The signatures of evolutionary forces acting on urotensin-II system indicate that it may have evolved towards loss of function since the divergence of primates.

Antennal transcriptome profiles of anopheline mosquitoes reveal human host olfactory specialization in Anopheles gambiae

BACKGROUND

Two sibling members of the Anopheles gambiae species complex display notable differences in female blood meal preferences. An. gambiae s.s. has a well-documented preference for feeding upon human hosts, whereas An. quadriannulatus feeds on vertebrate/mammalian hosts, with only opportunistic feeding upon humans. Because mosquito host-seeking behaviors are largely driven by the sensory modality of olfaction, we hypothesized that hallmarks of these divergent host seeking phenotypes will be in evidence within the transcriptome profiles of the antennae, the mosquito's principal chemosensory appendage.

RESULTS

To test this hypothesis, we have sequenced antennal mRNA of non-bloodfed females from each species and observed a number of distinct quantitative and qualitative differences in their chemosensory gene repertoires. In both species, these gene families show higher rates of sequence polymorphisms than the overall rates in their respective transcriptomes, with potentially important divergences between the two species. Moreover, quantitative differences in odorant receptor transcript abundances have been used to model potential distinctions in volatile odor receptivity between the two sibling species of anophelines.

CONCLUSION

This analysis suggests that the anthropophagic behavior of An. gambiae s.s. reflects the differential distribution of olfactory receptors in the antenna, likely resulting from a co-option and refinement of molecular components common to both species. This study improves our understanding of the molecular evolution of chemoreceptors in closely related anophelines and suggests possible mechanisms that underlie the behavioral distinctions in host seeking that, in part, account for the differential vectorial capacity of these mosquitoes.

Pseudoplusia includens single nucleopolyhedrovirus: genetic diversity, phylogeny and hypervariability of the pif-2 gene

The soybean looper (Pseudoplusia includens Walker, 1857) has become a major pest of soybean crops in Brazil. In order to determine the genetic diversity and phylogeny of variants of Pseudoplusia includens single nucleopolyhedrovirus (PsinSNPV-IA to -IG), partial sequences of the genes lef-8, lef-9, pif-2, phr and polh were obtained following degenerate PCR and phylogenetic trees constructed using maximum parsimony and Bayesian methods. The aligned sequences showed polymorphisms among the isolates, where the pif-2 gene was by far the most variable and is predicted to be under positive selection. Furthermore, some of the pif-2 DNA sequence mutations are predicted to result in significant amino acid substitutions, possibly leading to changes in oral infectivity of this baculovirus. Cladistic analysis revealed two closely related monophyletic groups, one containing PsinNPV isolates IB, IC and ID and another containing isolates IA, IE, IF and IG. The phylogeny of PsinSNPV in relation to 56 other baculoviruses was also determined from the concatenated partial LEF-8, LEF-9, PIF-2 and POLH/GRAN deduced amino acid sequences, using maximum-parsimony and Bayesian methods. This analysis clearly places PsinSNPV with the Group II Alphabaculovirus, where PsinSNPV is most closely related to Chrysodeixis chalcites NPV and Trichoplusia ni SNPV.

Small open reading frames associated with morphogenesis are hidden in plant genomes

It is likely that many small ORFs (sORFs; 30-100 amino acids) are missed when genomes are annotated. To overcome this limitation, we identified ∼8,000 sORFs with high coding potential in intergenic regions of the Arabidopsis thaliana genome. However, the question remains as to whether these coding sORFs play functional roles. Using a designed array, we generated an expression atlas for 16 organs and 17 environmental conditions among 7,901 identified coding sORFs. A total of 2,099 coding sORFs were highly expressed under at least one experimental condition, and 571 were significantly conserved in other land plants. A total of 473 coding sORFs were overexpressed; ∼10% (49/473) induced visible phenotypic effects, a proportion that is approximately seven times higher than that of randomly chosen known genes. These results indicate that many coding sORFs hidden in plant genomes are associated with morphogenesis. We believe that the expression atlas will contribute to further study of the roles of sORFs in plants.

Significance of population size on the fixation of nonsynonymous mutations in genes under varying levels of selection pressure

Previous studies observed a higher ratio of divergences at nonsynonymous and synonymous sites (ω = dN/dS) in species with a small population size compared to that estimated for those with a large population size. Here we examined the theoretical relationship between ω, effective population size (Ne), and selection coefficient (s). Our analysis revealed that when purifying selection is high, ω of species with small Ne is much higher than that of species with large Ne. However the difference between the two ω reduces with the decline in selection pressure (s → 0). We examined this relationship using primate and rodent genes and found that the ω estimated for highly constrained genes of primates was up to 2.9 times higher than that obtained for their orthologous rodent genes. Conversely, for genes under weak purifying selection the ω of primates was only 17% higher than that of rodents. When tissue specificity was used as a proxy for selection pressure we found that the ω of broadly expressed genes of primates was up to 2.1-fold higher than that of their rodent counterparts and this difference was only 27% for tissue specific genes. Since most of the nonsynonymous mutations in constrained or broadly expressed genes are deleterious, fixation of these mutations is influenced by Ne. This results in a higher ω of these genes in primates compared to those from rodents. Conversely, the majority of nonsynonymous mutations in less-constrained or tissue-specific genes are neutral or nearly neutral and therefore fixation of them is largely independent of Ne, which leads to the similarity of ω in primates and rodents.

Molecular evolutionary analysis of the duck MYOD gene family and its differential expression pattern in breast muscle development

1. The objective of the research was to investigate the molecular evolutionary relationships between the duck myogenic determination factors (MYOD) gene family members and their roles in muscle development. 2. The four members of the duck MYOD gene family were cloned using RT-PCR, and their relative mRNA expression during duck muscle development was measured using qRT-PCR. 3. The results showed that MyoD and Myf5 clustered together, as did MyoG and MRF4 based on their complete amino acid sequence and the basic helix-loop-helix domain. Results of the evolutionary level analysis were consistent with that of the differential expression patterns during duck breast muscle development. As determined by qRT-PCR, MyoD and Myf5 were highly expressed in 22-day embryos, while MyoG and MRF4 expression was high in 14-day embryos. 4. We conclude that the entire MYOD gene family in the duck originated from a common ancestral gene and evolved after two duplication events. The roles of the MYOD gene family members in duck muscle development are similar to those in mammals.

Independent domestication of Asian rice followed by gene flow from japonica to indica

Results from studies on the domestication process of Asian rice Oryza sativa have been controversial because of its complicated evolutionary history. Previous studies have yielded two alternative hypotheses about the origin(s) of the two major groups of O. sativa: japonica and indica. One study proposes a single common wild ancestor, whereas the other suggests that there were multiple domestication events of different types of wild rice. Here, we provide clear evidence of the independent domestication of japonica and indica obtained via high-throughput sequencing and a large-scale comparative analysis of two wild rice accessions (W1943 and W0106) and two cultivars (a japonica cultivar called "Nipponbare" and an indica cultivar called "Guangluai-4"). The different domestication processes of the two cultivar groups appear to have led to distinct patterns of molecular evolution in protein-coding regions. The intensity of purifying selection was relaxed only in the japonica group, possibly because of a bottleneck effect. Moreover, a genome-wide comparison between Nipponbare, Guangluai-4, and another indica cultivar (93-11) suggests multiple hybridization events between japonica and indica, both before and after the divergence of the indica cultivars. We found that a large amount of genomic DNA, including domestication-related genes, was transferred from japonica to indica, which might have been important in the development of modern rice. Our study provides an overview of the dynamic process of Asian rice domestication, including independent domestication events and subsequent gene flow.