Higher Rates Of Nonsynonymous Substitution Research Articles

Pervasive relaxed selection in termite genomes.

Genetic changes that enabled the evolution of eusociality have long captivated biologists. More recently, attention has focussed on the consequences of eusociality on genome evolution. Studies have reported higher molecular evolutionary rates in eusocial hymenopteran insects compared with their solitary relatives. To investigate the genomic consequences of eusociality in termites, we analysed nine genomes, including newly sequenced genomes from three non-eusocial cockroaches. Using a phylogenomic approach, we found that termite genomes have experienced lower rates of synonymous substitutions than those of cockroaches, possibly as a result of longer generation times. We identified higher rates of non-synonymous substitutions in termite genomes than in cockroach genomes, and identified pervasive relaxed selection in the former (24-31% of the genes analysed) compared with the latter (2-4%). We infer that this is due to reductions in effective population size, rather than gene-specific effects (e.g. indirect selection of caste-biased genes). We found no obvious signature of increased genetic load in termites, and postulate efficient purging of deleterious alleles at the colony level. Additionally, we identified genomic adaptations that may underpin caste differentiation, such as genes involved in post-translational modifications. Our results provide insights into the evolution of termites and the genomic consequences of eusociality more broadly.

Integrating massive RNA-seq data to elucidate transcriptome dynamics in Drosophila melanogaster.

The volume of ribonucleic acid (RNA)-seq data has increased exponentially, providing numerous new insights into various biological processes. However, due to significant practical challenges, such as data heterogeneity, it is still difficult to ensure the quality of these data when integrated. Although some quality control methods have been developed, sample consistency is rarely considered and these methods are susceptible to artificial factors. Here, we developed MassiveQC, an unsupervised machine learning-based approach, to automatically download and filter large-scale high-throughput data. In addition to the read quality used in other tools, MassiveQC also uses the alignment and expression quality as model features. Meanwhile, it is user-friendly since the cutoff is generated from self-reporting and is applicable to multimodal data. To explore its value, we applied MassiveQC to Drosophila RNA-seq data and generated a comprehensive transcriptome atlas across 28 tissues from embryogenesis to adulthood. We systematically characterized fly gene expression dynamics and found that genes with high expression dynamics were likely to be evolutionarily young and expressed at late developmental stages, exhibiting high nonsynonymous substitution rates and low phenotypic severity, and they were involved in simple regulatory programs. We also discovered that human and Drosophila had strong positive correlations in gene expression in orthologous organs, revealing the great potential of the Drosophila system for studying human development and disease.

Temporal dynamics of mildly deleterious nonsynonymous substitutions in mitochondrial gene sequences in rodents and moles.

We have previously estimated the evolutionary rate (number of substitutions/site/million years) of mitochondrial cytochrome b gene (Cytb) sequences in rodents and moles to be about 0.11 at more recent divergence times of tens of thousands of years, and to decrease rapidly to about 0.03 at more distant divergence times. Because this time dependency is thought to be caused by the removal of mildly deleterious substitutions in later generations, we focused in this study on the abundance of nonsynonymous substitutions. We collected 23 haplogroups of Cytb with signals of late Quaternary population expansion events from rodents and moles and categorized them into three groups for comparison based on predicted expansion start time: 5,000-15,000 years ago (Group I), ca. 53,000 years ago (Group II) and 130,000-230,000 years ago (Group III). We counted the numbers of nonsynonymous and synonymous substitutions in all haplogroups. The rates of nonsynonymous substitutions were lowest in Groups II and III (0.08-0.22), whereas those in Group I varied markedly. We further classified Group I into two subgroups based on high (0.29-0.43) and low (0.09-0.20) nonsynonymous substitution rates, which were likely to be associated with the start of the expansion within 10,000 years and at around 15,000 years ago, respectively. The Group II and III networks had two- or three-step star-shaped structures and tended to exhibit frequent and less frequent nonsynonymous substitutions on exterior and interior branches, respectively. Based on temporal dynamics, nonsynonymous mitochondrial DNA (mtDNA) substitutions in small mammals accounted for at most 40% of all substitutions during the early evolutionary stage and then rapidly declined, dropping to approximately 15%. The results of this study provide a good explanation of the time-dependent trend in the mtDNA evolution rate predicted in previous work.

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.

Exceptionally high rates of positive selection on the rbcL gene in the genus Ilex (Aquifoliaceae)

BackgroundThe genus Ilex (Aquifoliaceae) has a near-cosmopolitan distribution in mesic habitats from tropical to temperate lowlands and in alpine forests. It has a high rate of hybridization and plastid capture, and comprises four geographically structured plastid groups. A previous study showed that the plastid rbcL gene, coding for the large subunit of Rubisco, has a particularly high rate of non-synonymous substitutions in Ilex, when compared with other plant lineages. This suggests a strong positive selection on rbcL, involved in yet unknown adaptations. We therefore investigated positive selection on rbcL in 240 Ilex sequences from across the global range.ResultsThe rbcL gene shows a much higher rate of positive selection in Ilex than in any other plant lineage studied so far (> 3000 species) by tests in both PAML and SLR. Most positively selected residues are on the surface of the folded large subunit, suggesting interaction with other subunits and associated chaperones, and coevolution between positively selected residues is prevalent, indicating compensatory mutations to recover molecular stability. Coevolution between positively selected sites to restore global stability is common.ConclusionsThis study has confirmed the predicted high incidence of positively selected residues in rbcL in Ilex, and shown that this is higher than in any other plant lineage studied so far. The causes and consequences of this high incidence are unclear, but it is probably associated with the similarly high incidence of hybridization and introgression in Ilex, even between distantly related lineages, resulting in large cytonuclear discordance in the phylogenies.

Sex-biased gene expression and evolution of candidate reproductive transcripts in adult stages of salmon lice (Lepeophtheirus salmonis)

The salmon louse Lepeophtheirus salmonis (Krøyer 1837) displays numerous sexually dimorphic characteristics. Insights into their underlying molecular components have only recently been explored, which serve to better understand both the basic biology of the louse, associated impacts on drug sensitivity, and evolution of resistance. Expression of 16 L. salmonis genes putatively involved in sexual dimorphism and reproduction were used to determine differences between sexes and better understand responses to mating using RT-qPCR of pre-adult and adult lice. Analysis of these genes revealed the dynamic nature of sex-biased expression across stages. However, female reception of a spermatophore did not appear to impact the expression of these particular genes. Furthermore six of these transcripts and 84 others identified previously in a large-scale louse transcriptomics experiment were used to estimate differences in evolutionary rates and codon-usage bias of sex-related genes using phylogenetic analysis by maximum likelihood (PAML) and codonW. Results suggest male-biased genes are evolving at significantly greater rates than female-biased and unbiased genes as evidenced by higher rates of non-synonymous substitutions and lower codon-usage bias in these genes. These analyses expand our understanding of interactions of sex-biased expression across the pre-adult and adult life stages and provide foundations for better understanding evolutionary differences in sex-biased genes of L. salmonis.

Evolution of Gustatory Receptor Gene Family Provides Insights into Adaptation to Diverse Host Plants in Nymphalid Butterflies.

The host plant range of herbivorous insects is a major aspect of insect–plant interaction, but the genetic basis of host range expansion in insects is poorly understood. In butterflies, gustatory receptor genes (GRs) play important roles in host plant selection by ovipositing females. Since several studies have shown associations between the repertoire sizes of chemosensory gene families and the diversity of resource use, we hypothesized that the increase in the number of genes in the GR family is associated with host range expansion in butterflies. Here, we analyzed the evolutionary dynamics of GRs among related species, including the host generalist Vanessa cardui and three specialists. Although the increase of the GR repertoire itself was not observed, we found that the gene birth rate of GRs was the highest in the lineage leading to V. cardui compared with other specialist lineages. We also identified two taxon-specific subfamilies of GRs, characterized by frequent lineage-specific duplications and higher non-synonymous substitution rates. Together, our results suggest that frequent gene duplications in GRs, which might be involved in the detection of plant secondary metabolites, were associated with host range expansion in the V. cardui lineage. These evolutionary patterns imply that the capability to perceive various compounds during host selection was favored during adaptation to diverse host plants.

Comparison of the Chinese bamboo partridge and red Junglefowl genome sequences highlights the importance of demography in genome evolution

BackgroundRecent large-scale whole genome sequencing efforts in birds have elucidated broad patterns of avian phylogeny and genome evolution. However, despite the great interest in economically important phasianids like Gallus gallus (Red Junglefowl, the progenitor of the chicken), we know little about the genomes of closely related species. Gallus gallus is highly sexually dichromatic and polygynous, but its sister genus, Bambusicola, is smaller, sexually monomorphic, and monogamous with biparental care. We sequenced the genome of Bambusicola thoracicus (Chinese Bamboo Partridge) using a single insert library to test hypotheses about genome evolution in galliforms. Selection acting at the phenotypic level could result in more evidence of positive selection in the Gallus genome than in Bambusicola. However, the historical range size of Bambusicola was likely smaller than Gallus, and demographic effects could lead to higher rates of nonsynonymous substitution in Bambusicola than in Gallus.ResultsWe generated a genome assembly suitable for evolutionary analyses. We examined the impact of selection on coding regions by examining shifts in the average nonsynonymous to synonymous rate ratio (dN/dS) and the proportion of sites subject to episodic positive selection. We observed elevated dN/dS in Bambusicola relative to Gallus, which is consistent with our hypothesis that demographic effects may be important drivers of genome evolution in Bambusicola. We also demonstrated that alignment error can greatly inflate estimates of the number of genes that experienced episodic positive selection and heterogeneity in dN/dS. However, overall patterns of molecular evolution were robust to alignment uncertainty. Bambusicola thoracicus has higher estimates of heterozygosity than Gallus gallus, possibly due to migration events over the past 100,000 years.ConclusionsOur results emphasized the importance of demographic processes in generating the patterns of variation between Bambusicola and Gallus. We also demonstrated that genome assemblies generated using a single library can provide valuable insights into avian evolutionary history and found that it is important to account for alignment uncertainty in evolutionary inferences from draft genomes.

Examination of Signatures of Recent Positive Selection on Genes Involved in Human Sialic Acid Biology.

Sialic acids are nine carbon sugars ubiquitously found on the surfaces of vertebrate cells and are involved in various immune response-related processes. In humans, at least 58 genes spanning diverse functions, from biosynthesis and activation to recycling and degradation, are involved in sialic acid biology. Because of their role in immunity, sialic acid biology genes have been hypothesized to exhibit elevated rates of evolutionary change. Consistent with this hypothesis, several genes involved in sialic acid biology have experienced higher rates of non-synonymous substitutions in the human lineage than their counterparts in other great apes, perhaps in response to ancient pathogens that infected hominins millions of years ago (paleopathogens). To test whether sialic acid biology genes have also experienced more recent positive selection during the evolution of the modern human lineage, reflecting adaptation to contemporary cosmopolitan or geographically-restricted pathogens, we examined whether their protein-coding regions showed evidence of recent hard and soft selective sweeps. This examination involved the calculation of four measures that quantify changes in allele frequency spectra, extent of population differentiation, and haplotype homozygosity caused by recent hard and soft selective sweeps for 55 sialic acid biology genes using publicly available whole genome sequencing data from 1,668 humans from three ethnic groups. To disentangle evidence for selection from confounding demographic effects, we compared the observed patterns in sialic acid biology genes to simulated sequences of the same length under a model of neutral evolution that takes into account human demographic history. We found that the patterns of genetic variation of most sialic acid biology genes did not significantly deviate from neutral expectations and were not significantly different among genes belonging to different functional categories. Those few sialic acid biology genes that significantly deviated from neutrality either experienced soft sweeps or population-specific hard sweeps. Interestingly, while most hard sweeps occurred on genes involved in sialic acid recognition, most soft sweeps involved genes associated with recycling, degradation and activation, transport, and transfer functions. We propose that the lack of signatures of recent positive selection for the majority of the sialic acid biology genes is consistent with the view that these genes regulate immune responses against ancient rather than contemporary cosmopolitan or geographically restricted pathogens.

Selective pressure on the protein-coding genes of the pufferfish is correlated with phenotypic traits

The pufferfish accumulates neurotoxic tetrodotoxin in its body and inflates by filling its stomach with water. These traits are unique to this species, and may be a result of adaptation post-divergence of Tetraodontidae. However, evolution of the protein-coding genes in the pufferfish has not yet been well elucidated. Detection of positive selection on these genes can help us understand the mechanisms associated with functional evolution. We downloaded well-annotated gene information of two pufferfish species, Takifugu rubripes and Tetraodon nigroviridis, from the public ENSEMBL database. In order to detect selective pressure on protein-coding sequences, we performed dN/dS estimation using codeml within the PAML software package. We selected one to one orthologous genes among seven fish species (Gasterosteus aculeatus, Oryzias latipes, Poecilia formosa, Takifugu rubripes, Tetraodon nigroviridis, and Xiphophorus maculatus). Results of dN/dS analysis on orthologous genes indicate that pufferfish showed high non-synonymous substitution rate for positively selected genes, and the evolutionary rate was faster during the diversification of two pufferfishes after divergence. Additionally, a candidate mechanism for regulation of neuro-toxicity of tetrodotoxin was identified from functional annotation of positively selected genes. These results support positive selection on protein-coding genes of the pufferfish with the acquisition of specific phenotypic traits.

Characterization of major histocompatibility complex class I, and class II DRB loci of captive and wild Indian leopards (Panthera pardus fusca).

The major histocompatibility complex (MHC), in vertebrate animals, is a multi-genic protein complex that encodes various receptors. During a disease, MHC interacts with the antigen and triggers a cascade of adaptive immune responses to overcome a disease outbreak. The MHC is very important region from immunological point of view, but it is poorly characterized among Indian leopards. During this investigation, we examined genetic diversity for MHC class I (MHC-I) and MHC class II-DRB (MHC-II) among wild and captive Indian leopards. This study estimated a pool of 9 and 17 alleles for MHC-I and MHC-II, respectively. The wild group of individuals showed higher nucleotide diversity and amino acid polymorphism compared to the captive group. A phylogenetic comparison with other felids revealed a clustering in MHC-I and interspersed presence in MHC-II sequences. A test for selection also revealed a deviation from neutrality at MHC-II DRB loci and higher non-synonymous substitution rate (dN) among the individuals from wild group. Further, the wild individuals showed higher dN for both MHC I and II genes compared to the group that was bred under captive conditions. These findings suggest the role of micro-evolutionary forces, such as pathogen-mediated selection, to cause MHC variations among the two groups of Indian leopards, because the two groups have been bred in two different environments for a substantial period of time. Since, MHC diversity is often linked with the quality of immunological health; the results obtained from this study fill the gap of knowledge on disease predisposition among wild and captive Indian leopards.

Genome-wide signals of positive selection in strongylocentrotid sea urchins

BackgroundComparative genomics studies investigating the signals of positive selection among groups of closely related species are still rare and limited in taxonomic breadth. Such studies show great promise in advancing our knowledge about the proportion and the identity of genes experiencing diversifying selection. However, methodological challenges have led to high levels of false positives in past studies. Here, we use the well-annotated genome of the purple sea urchin, Strongylocentrotus purpuratus, as a reference to investigate the signals of positive selection at 6520 single-copy orthologs from nine sea urchin species belonging to the family Strongylocentrotidae paying careful attention to minimizing false positives.ResultsWe identified 1008 (15.5%) candidate positive selection genes (PSGs). Tests for positive selection along the nine terminal branches of the phylogeny identified 824 genes that showed lineage-specific adaptive diversification (1.67% of branch-sites tests performed). Positively selected codons were not enriched at exon borders or near regions containing missing data, suggesting a limited contribution of false positives caused by alignment or annotation errors. Alignments were validated at 10 loci with re-sequencing using Sanger methods. No differences were observed in the rates of synonymous substitution (dS), GC content, and codon bias between the candidate PSGs and those not showing positive selection. However, the candidate PSGs had 68% higher rates of nonsynonymous substitution (dN) and 33% lower levels of heterozygosity, consistent with selective sweeps and opposite to that expected by a relaxation of selective constraint. Although positive selection was identified at reproductive proteins and innate immunity genes, the strongest signals of adaptive diversification were observed at extracellular matrix proteins, cell adhesion molecules, membrane receptors, and ion channels. Many candidate PSGs have been widely implicated as targets of pathogen binding, inactivation, mimicry, or exploitation in other groups (notably mammals).ConclusionsOur study confirmed the widespread action of positive selection across sea urchin genomes and allowed us to reject the possibility that annotation and alignment errors (including paralogs) were responsible for creating false signals of adaptive molecular divergence. The candidate PSGs identified in our study represent promising targets for future research into the selective agents responsible for their adaptive diversification and their contribution to speciation.

Evolutionary and biogeographical implications of degraded LAGLIDADG endonuclease functionality and group I intron occurrence in stony corals (Scleractinia) and mushroom corals (Corallimorpharia).

Group I introns and homing endonuclease genes (HEGs) are mobile genetic elements, capable of invading target sequences in intron-less genomes. LAGLIDADG HEGs are the largest family of endonucleases, playing a key role in the mobility of group I introns in a process known as ‘homing’. Group I introns and HEGs are rare in metazoans, and can be mainly found inserted in the COXI gene of some sponges and cnidarians, including stony corals (Scleractinia) and mushroom corals (Corallimorpharia). Vertical and horizontal intron transfer mechanisms have been proposed as explanations for intron occurrence in cnidarians. However, the central role of LAGLIDADG motifs in intron mobility mechanisms remains poorly understood. To resolve questions regarding the evolutionary origin and distribution of group I introns and HEGs in Scleractinia and Corallimorpharia, we examined intron/HEGs sequences within a comprehensive phylogenetic framework. Analyses of LAGLIDADG motif conservation showed a high degree of degradation in complex Scleractinia and Corallimorpharia. Moreover, the two motifs lack the respective acidic residues necessary for metal-ion binding and catalysis, potentially impairing horizontal intron mobility. In contrast, both motifs are highly conserved within robust Scleractinia, indicating a fully functional endonuclease capable of promoting horizontal intron transference. A higher rate of non-synonymous substitutions (Ka) detected in the HEGs of complex Scleractinia and Corallimorpharia suggests degradation of the HEG, whereas lower Ka rates in robust Scleractinia are consistent with a scenario of purifying selection. Molecular-clock analyses and ancestral inference of intron type indicated an earlier intron insertion in complex Scleractinia and Corallimorpharia in comparison to robust Scleractinia. These findings suggest that the lack of horizontal intron transfers in the former two groups is related to an age-dependent degradation of the endonuclease activity. Moreover, they also explain the peculiar geographical patterns of introns in stony and mushroom corals.

Analysis of Plant Pan-Genomes and Transcriptomes with GET_HOMOLOGUES-EST, a Clustering Solution for Sequences of the Same Species.

The pan-genome of a species is defined as the union of all the genes and non-coding sequences found in all its individuals. However, constructing a pan-genome for plants with large genomes is daunting both in sequencing cost and the scale of the required computational analysis. A more affordable alternative is to focus on the genic repertoire by using transcriptomic data. Here, the software GET_HOMOLOGUES-EST was benchmarked with genomic and RNA-seq data of 19 Arabidopsis thaliana ecotypes and then applied to the analysis of transcripts from 16 Hordeum vulgare genotypes. The goal was to sample their pan-genomes and classify sequences as core, if detected in all accessions, or accessory, when absent in some of them. The resulting sequence clusters were used to simulate pan-genome growth, and to compile Average Nucleotide Identity matrices that summarize intra-species variation. Although transcripts were found to under-estimate pan-genome size by at least 10%, we concluded that clusters of expressed sequences can recapitulate phylogeny and reproduce two properties observed in A. thaliana gene models: accessory loci show lower expression and higher non-synonymous substitution rates than core genes. Finally, accessory sequences were observed to preferentially encode transposon components in both species, plus disease resistance genes in cultivated barleys, and a variety of protein domains from other families that appear frequently associated with presence/absence variation in the literature. These results demonstrate that pan-genome analyses are useful to explore germplasm diversity.

Evolution of the RH gene family in vertebrates revealed by brown hagfish (Eptatretus atami) genome sequences

In vertebrates, there are four major genes in the RH (Rhesus) gene family, RH, RHAG, RHBG, and RHCG. These genes are thought to have been formed by the two rounds of whole-genome duplication (2R-WGD) in the common ancestor of all vertebrates. In our previous work, where we analyzed details of the gene duplications process of this gene family, three nucleotide sequences belonging to this family were identified in Far Eastern brook lamprey (Lethenteron reissneri), and the phylogenetic positions of the genes were determined. Lampreys, along with hagfishes, are cyclostomata (jawless fishes), which is a sister group of gnathostomata (jawed vertebrates). Although those results suggested that one gene was orthologous to the gnathostome RHCG genes, we did not identify clear orthologues for other genes. In this study, therefore, we identified three novel cDNA sequences that belong to the RH gene family using de novo transcriptome analysis of another cyclostome: the brown hagfish (Eptatretus atami). We also determined the nucleotide sequences for the RHBG and RHCG genes in a red stingray (Dasyatis akajei), which belongs to the cartilaginous fishes. The phylogenetic tree showed that two brown hagfish genes, which were probably duplicated in the cyclostome lineage, formed a cluster with the gnathostome RHAG genes, whereas another brown hagfish gene formed a cluster with the gnathostome RHCG genes. We estimated that the RH genes had a higher evolutionary rate than the RHAG, RHBG, and RHCG genes. Interestingly, in the RHBG genes, only the bird lineage showed a higher rate of nonsynonymous substitutions. It is likely that this higher rate was caused by a state of relaxed functional constraints rather than positive selection nor by pseudogenization.

Genetic Characterization of Northwestern Colombian Chikungunya Virus Strains from the 2014-2015 Epidemic.

Chikungunya fever, an acute and often chronic arthralgic disease caused by the mosquito-borne alphavirus, chikungunya virus (CHIKV), spread into the Americas in late 2013. Since then it has caused epidemics in nearly all New World countries, the second largest being Colombia with over 450,000 suspected cases beginning in September, 2014, and focused in Bolivar Department in the north. We examined 32 human sera from suspected cases, including diverse age groups and both genders, and sequenced the CHIKV envelope glycoprotein genes, known determinants of vector host range. As expected for Asian lineage CHIKV strains, these isolates lacked known Aedes albopictus-adaptive mutations. All the Colombian strains were closely related to those from the Virgin Islands, Saint Lucia, Mexico, Puerto Rico, and Brazil, consistent with a single, point-source introduction from the southeast Asia/Pacific region. Two substitutions in the E2 and E1 envelope glycoprotein genes were found in the Colombian strains, especially E1-K211E involving a residue shown previously to affect epistatically the penetrance of the E1-A226V A. albopictus-adaptive substitution. We also identified two amino acid substitutions unique to all American CHIKV sequences: E2-V368A and 6K-L20M. Only one codon, 6K-47, had a high nonsynonymous substitution rate suggesting positive selection.

Large-Scale Gene Relocations following an Ancient Genome Triplication Associated with the Diversification of Core Eudicots.

Different modes of gene duplication including whole-genome duplication (WGD), and tandem, proximal and dispersed duplications are widespread in angiosperm genomes. Small-scale, stochastic gene relocations and transposed gene duplications are widely accepted to be the primary mechanisms for the creation of dispersed duplicates. However, here we show that most surviving ancient dispersed duplicates in core eudicots originated from large-scale gene relocations within a narrow window of time following a genome triplication (γ) event that occurred in the stem lineage of core eudicots. We name these surviving ancient dispersed duplicates as relocated γ duplicates. In Arabidopsis thaliana, relocated γ, WGD and single-gene duplicates have distinct features with regard to gene functions, essentiality, and protein interactions. Relative to γ duplicates, relocated γ duplicates have higher non-synonymous substitution rates, but comparable levels of expression and regulation divergence. Thus, relocated γ duplicates should be distinguished from WGD and single-gene duplicates for evolutionary investigations. Our results suggest large-scale gene relocations following the γ event were associated with the diversification of core eudicots.

Modelling the role of immunity in reversion of viral antigenic sites

Antigenic sites in viral pathogens exhibit distinctive evolutionary dynamics due to their role in evading recognition by host immunity. Antigenic selection is known to drive higher rates of non-synonymous substitution; less well understood is why differences are observed between viruses in their propensity to mutate to a novel or previously encountered amino acid. Here, we present a model to explain patterns of antigenic reversion and forward substitution in terms of the epidemiological and molecular processes of the viral population. We develop an analytical three-strain model and extend the analysis to a multi-site model to predict characteristics of observed sequence samples. Our model provides insight into how the balance between selection to escape immunity and to maintain viability is affected by the rate of mutational input. We also show that while low probabilities of reversion may be due to either a low cost of immune escape or slowly decaying host immunity, these two scenarios can be differentiated by the frequency patterns at antigenic sites. Comparison between frequency patterns of human influenza A (H3N2) and human RSV-A suggests that the increased rates of antigenic reversion in RSV-A is due to faster decaying immunity and not higher costs of escape.

Evolution of hypothalamus-pituitary growth axis among fish, amphibian, birds and mammals

Hypothalamus-pituitary growth axis (HP growth axis) regulates animal growth and development in pre-natal and post natal life governed by many factors. However, until recently, the evolutionary history of this axis among lineages is not understood. Aim of the present study was to understand the major events in evolution and evolutionary history and trend of HP growth axis. The diversity among Homo sapience, Mus musculus, Rattus norvegicus, Gallus gallus, Danio rerio and Xenopus laevis was determined for genes involved in HP growth axis in current study. Sequences of HP growth axis genes were retrieved from NCBI (http://www.ncbi.nlm.nih.gov/). Nucleotide diversity using Kimura?s two-parameter method; codon-based test of positive selection using the Nei-Gojobori; equality of evolutionary rate with Tajima's relative rate test and phylogenetic history using the RelTime method were estimated in MEGA6. Estimates of the coefficients of evolutionary differentiation based on nucleotides and amino acids substitution patterns of HP growth axis genes showed contrasting evolutionary patterns among the lineages. The results demonstrated that although these genes might have crucial functional roles in each of the species, however, their sequence divergence did not necessarily reflect similar molecular evolution among the species. Codon-based test of positive selection revealed that Human vs Mouse, Chicken vs Rat, Human vs Rat and Mouse vs Rat had similar and higher non synonymous substitutions (P > 0.05). Higher rate of non-synonymous substitutions at similar orthologs level among species indicated a similar positive selection pressure in these species. Results for relative rate test assessed with the chi-squared test showed difference on unique mutations among lineages at synonymous and non synonymous sites except Chicken vs Mouse, Human vs Mouse, Chicken vs Rat, Human vs Rat and Mouse vs Rat. This indicated that the mutagenic process that generates substitutional mutation is taking place at approximately the same rate at synonymous and non-synonymous sites these lineages. Moreover, despite of common ancestry, our results indicate a different divergent time among genes of these species. This is the first demonstration that variable rates of molecular evolution may be present within HP growth axis genes among different species. This difference could be of interest for comparative genomics analysis and physiological genes functions identification among tho comparative genomics, evolution rate, HP growth axis, positive selection se species whose HP growth axis is not explored.

Degeneration of the nonrecombining regions in the mating-type chromosomes of the anther-smut fungi.

Dimorphic mating-type chromosomes in fungi are excellent models for understanding the genomic consequences of recombination suppression. Their suppressed recombination and reduced effective population size are expected to limit the efficacy of natural selection, leading to genomic degeneration. Our aim was to identify the sequences of the mating-type chromosomes (a1 and a2) of the anther-smut fungi and to investigate degeneration in their nonrecombining regions. We used the haploid a1 Microbotryum lychnidis-dioicae reference genome sequence. The a1 and a2 mating-type chromosomes were both isolated electrophoretically and sequenced. Integration with restriction-digest optical maps identified regions of recombination and nonrecombination in the mating-type chromosomes. Genome sequence data were also obtained for 12 other Microbotryum species. We found strong evidence of degeneration across the genus in the nonrecombining regions of the mating-type chromosomes, with significantly higher rates of nonsynonymous substitution (dN/dS) than in nonmating-type chromosomes or in recombining regions of the mating-type chromosomes. The nonrecombining regions of the mating-type chromosomes also showed high transposable element content, weak gene expression, and gene losses. The levels of degeneration did not differ between the a1 and a2 mating-type chromosomes, consistent with the lack of homogametic/heterogametic asymmetry between them, and contrasting with X/Y or Z/W sex chromosomes.