Differences In Evolutionary Rates Research Articles

Drivers and patterns of bite force evolution in liolaemid lizards

Abstract Phenotypic variation is the result of selection on traits that are relevant in a given ecological context. Phylogenetic history, genetic drift, and any developmental or structural constraints may, however, limit variation in trait expression. It has been proposed that organismal performance traits take up a pivotal role in driving variation in morphology due to their central role in survival and reproductive success. However, how strong the links are between morphology and performance, and how the strength of this relationship impacts the rate of evolution of form and function need to be studied across a wider variety of systems to better understand the origin and evolution of biodiversity. Here we used data on the jaw system (muscle architecture and head dimensions) of liolaemid lizards to investigate the drivers of in vivo bite force variation and test for differences in evolutionary rates in morphology and performance. Our results show high rates of evolution for performance traits compared to morphological traits such as external head dimensions. Many-to-one mapping of morphology to performance, that is the possibility that different anatomical trait combinations lead to similar levels of performance, appears to be common in the jaw system of these lizards. Finally, traits showing greater mechanical sensitivity (muscle cross-sectional areas) showed higher rates of evolution compared to traits involved in other functions and that are probably subject to trade-offs (e.g. head width).

Extreme climate shifts pest dominance hierarchy through thermal evolution and transgenerational plasticity

Abstract Dominance hierarchy, the dominance ranking of members in a community, is determined by the relative population fitness of coexisting species. Climate change can shift dominance hierarchies depending on how species respond to changing climate, but ecological mechanisms underlying such shifts remain largely unknown. Specifically, dominance hierarchy shifts under climate change have rarely been linked to eco‐evolutionary responses to thermal extremes, which we consider in this study. We conducted an 8‐year field survey to link extreme high‐temperature events to the shift in dominance hierarchy of two worldwide cereal aphid species (Rhopalosiphum padi and Sitobion avenae), which may respond rapidly through evolutionary or plastic responses to thermal extremes due to their short generation duration, clonal structure and high thermal sensitivity. To further understand the mechanism involved in this change in species' dominance, we characterised aphid heat tolerance and demography of the two species sourced from relatively mild (Xinxiang) and hot (Wuhan) locations in a common garden design and compared the same measures in up‐ and down‐selected lines of the two aphid species. We found that accumulation of extreme high‐temperature events (EHTs) affected per capita growth rate of S. avenae but not R. padi, driving a shift in annual relative dominance of these two common species in wheat fields. Furthermore, evidence from common garden tests and artificial thermal selection revealed that evolutionary and plastic changes in thermal tolerance under EHTs were species‐dependent; R. padi evolved higher heat tolerance without reducing fitness while S. avenae did not evolve heat tolerance but showed detrimental transgenerational plastic changes under thermal extremes. Rhopalosiphum padi may take over the dominant position of S. avenae in a warmer future, due to a higher evolutionary potential of heat tolerance and less detrimental plasticity following selection. Field surveys and laboratory experiments together point to a shift in species dominance related to interspecific differences in species' evolutionary rates and transgenerational plasticity of thermal traits during selection events associated with extreme climatic conditions. This study highlights that eco‐evolutionary dynamics can contribute to community responses to natural thermal extremes. A free Plain Language Summary can be found within the Supporting Information of this article.

Dynamics of Nanoparticle Self-Assembly by Liquid Crystal Sorting in Two Dimensions

We study nonlinear dynamical equations for coupled conserved and non-conserved fields describing nanoparticle concentration and liquid crystal order parameter, respectively, and solve them numerically over bidimensional domains. These equations model the rapid segregation of nanoparticles away from nematic domains, which has been observed experimentally in a suspension of gold nanoparticles in 5CB below the isotropic-nematic transition temperature. We contrast the different behaviors obtained when the LC order parameter is treated as a scalar or a tensor, as well as the different rates of evolution observed with each of these. We find, after an instantaneous quench lowering the temperature below the transition one, an initial linear regime where the ordering of the nematic phase proceeds exponentially with time. Only after a lag period the nanoparticle material couples effectively to the LC order parameter and segregates to regions that are less orientationally ordered (extended domain walls for a scalar order parameter, but point disclinations for a tensor one). The lag period is followed by the onset of nonlinear dynamics and saturation of the order parameter. The choice of a scalar or tensor LC order parameter does not change this sequence but results in a clear overshooting of the nonlinear saturation level for the tensor order parameter case. These results are found to be insensitive to weak anchoring due to coupling of gradients of the conserved and non-conserved variables, for the nanoparticle concentrations and anchoring parameters studied. Our modeling approach can be extended in a straightforward manner to cases where the cooling rate is finite and to other systems where a locally conserved concentration is coupled to a orientation field, such as active Langmuir monolayers, and possibly to other examples of nonlinear dynamics in ecological or excitable media problems.

An Investigation of the Variations in Complete Mitochondrial Genomes of Lingula anatina in the Western Pacific Region.

Simple SummaryLingula anatina has been attracted researchers because its morphological characteristics show limited changes compared to the ancestor. Even though L. anatina is a common brachiopod in the western Pacific region, a few studies were performed to investigate genetic variations of this species. To understand the mitogenomic diversity of L. anatina in the region, the present study conducted comparative analyses of mitogenome sequences from Korea, Vietnam, and Japan. The sequencing results indicated that L. anatina mitogenomes are extraordinarily longer than the typical animal mitogenome. Besides, the gene orders of L. anatina mitogenomes are variable among localities. The calculation based on protein-coding genes revealed relatively low substitution rates among mitogenome sequences from Buan (Korea), Doson (Vietnam) and Yanagawa (Japan). The phylogenetic analyses indicated the divergence of L. anatina in the western Pacific region. This investigation provides new information on the molecular systematics of L. anatina and could be helpful in exploring the diversity and evolution of mitogenomes in brachipod.Lingula anatina is a brachiopod widely distributed in the western Pacific region. Even though L. anatina has been targeted for a number of biological studies, there is still limited information on intraspecific genetic variations of L. anatina. In this study, L. anatina specimens were collected from Korea and Vietnam, and complete mitochondrial genome (mitogenome) sequences were analyzed and compared with previous records. The total mitogenomes of L. anatina were 24,875 bp and 25,305 bp in size for Korean and Vietnamese specimens, respectively. Those mitogenomes are extraordinarily longer than the typical mitogenome size for an animal but shorter than the previous record from Yanagawa (Japan) for this species. The gene orders and the sizes of the protein-coding genes are also different from those for the Japanese specimen. Furthermore, the nonsynonymous (Ka) and synonymous (Ks) substitution rates in protein-coding genes (PCGs) were calculated to test the idea of evolutionary rate differences in mitochondrial genomes. The analyses showed relatively low Ka and Ks for the complete mitogenomes from Buan (Korea), Doson (Vietnam) and Yanagawa (Japan). The Ka/Ks ratio was less than 1 in comparisons of three localities, indicating the existence of purifying selection in this species. The phylogenetic analyses showed that L. anatina diverged among localities in the western Pacific region.

Progress on the origin, function and evolutionary mechanism of RHR transcription factor family.

Transcription factors are a class of proteins that regulate gene transcription and expression by binding to gene-specific sequences and play an essential role in regulating the biological activities of cells. The RHR (Rel-homology region) transcription factor family is the primary member of the IF (immunoglobulin fold) transcription factor superfamily, whose members contain the conserved Rel domain and IPT (immunoglobulin-like fold) domain. As an ancient transcription factor family, the RHR family continues differentiation on gene gain and loss through gene duplication, mutation, and silencing, accompanied with the evolution of diverse species. Natural selection has led to different rates of evolution among members of the family, and some domains of the protein family have shown unique mechanisms of evolution. However, the current reviews about the origin and differentiation of RHR family are rare. In this review, we summarize the research results on the distribution, classification, function, and evolution of the members of the RHR family in order to provide a reference and new idea for studying the evolution mechanism of the whole transcription factor family and the evolutionary relationship among species.

English

Current study was conducted on family Aeshnidae from Hazara region of Pakistan. During the survey a total of 125 members were collected and identified into 2 species under single genus. The present study focuses on molecular characterizations and phylogenetics of family Aeshnidae. Phylogenies of the analyzed taxa were elaborated with maximum likelihood, maximum parsimony and Bayesian analysis. We sequenced both mitochondrial genes i.e., COI and 16S rRNA, separate and combined CO1+16S data sets revealed evolutionary relationship within Aeshnidae at the species and genera level. Mean Pairwise Distances (MPD) of each species were ranged from 0.00 to 84.60%. Evolutionary rate differences among two categories Gamma distribution and Invariant were recorded as 0.07 and 1.20 substitutions per site. DNA based identification using CO1, 16S and combined CO1+16S data set, for all Aeshnidae species shared genetic similarities having bootstrap values MLB=70–100%, MPB= 52–100% and BPP=0.75–1% respectively. The analysis of the combined (COI+16S) data set produced trees with complete stronger bootstrap support than analyses of either gene alone. These findings had shown that the taxonomic position of Aeshnidae species based on morphological characters could be verified, further improved and confirmed by the use of modern molecular biological tools which involve the nucleotide sequences of genes used in phylogenetic investigations. © 2021 Friends Science Publishers

Mutational Pressure Drives Differential Genome Conservation in Two Bacterial Endosymbionts of Sap-Feeding Insects

Compared with free-living bacteria, endosymbionts of sap-feeding insects have tiny and rapidly evolving genomes. Increased genetic drift, high mutation rates, and relaxed selection associated with host control of key cellular functions all likely contribute to genome decay. Phylogenetic comparisons have revealed massive variation in endosymbiont evolutionary rate, but such methods make it difficult to partition the effects of mutation versus selection. For example, the ancestor of Auchenorrhynchan insects contained two obligate endosymbionts, Sulcia and a betaproteobacterium (BetaSymb; called Nasuia in leafhoppers) that exhibit divergent rates of sequence evolution and different propensities for loss and replacement in the ensuing ∼300 Ma. Here, we use the auchenorrhynchan leafhopper Macrosteles sp. nr. severini, which retains both of the ancestral endosymbionts, to test the hypothesis that differences in evolutionary rate are driven by differential mutagenesis. We used a high-fidelity technique known as duplex sequencing to measure and compare low-frequency variants in each endosymbiont. Our direct detection of de novode novo mutations reveals that the rapidly evolving endosymbiont (Nasuia) has a much higher frequency of single-nucleotide variants than the more stable endosymbiont (Sulcia) and a mutation spectrum that is potentially even more AT-biased than implied by the 83.1% AT content of its genome. We show that indels are common in both endosymbionts but differ substantially in length and distribution around repetitive regions. Our results suggest that differences in long-term rates of sequence evolution in Sulcia versus BetaSymb, and perhaps the contrasting degrees of stability of their relationships with the host, are driven by differences in mutagenesis.

The differential diversification of Mongolic

Abstract This paper examines the Mongolic family of languages from the point of view of their different paths and rates of evolution, and with a view on the general problem concerning the speed of language change. All extant Mongolic languages descend from a relatively recent ancestral form of speech, Proto-Mongolic, spoken by the historical Mongols in the twelfth to fourteenth centuries and documented in a number of written sources conventionally known as Middle Mongol. A comparison of the modern Mongolic languages with Proto-Mongolic and Middle Mongol reveals considerable differences in their rates of evolution, with some languages being highly innovative, while others are conspicuously conservative. These differences are evident at all levels of linguistic structure and substance, including phonology, grammar, and lexicon. The reasons for the different rates of evolution can be sought in a variety of linguistic and extralinguistic factors, including not only the linguistic environment, but also the political and geographical context of the speakers, as well as their demographic structure, economic and cultural profile, and degree of mobility.

Male genitalia's evolutionary rate is higher than those of body traits: the case of two Liolaemus lizards' group

AbstractDifferent studies suggest that reproductive characters evolve faster than non‐reproductive characters. Males in the order Squamata have paired copulatory organs called hemipenes, with high morphological diversity, including differences in size, shape, and ornamentation. Some studies in the species‐rich lizard genus, Anolis suggest that genital traits evolve faster than the rest of the body. However, these studies were made considering only a few traits, across a wide phylogeny, without considering species relatedness, which may inflate differences in evolutionary rates. Here, we study two phylogenetic distantly related lizard groups, which differ in the number of species, but have similar divergence times. We evaluate as follows: (1) evolutionary rate, models of evolution and phylogenetic signal among the different genital and non‐genital traits; (2) which kind of traits (genital and non‐genital) are divergent across sister species and (3) whether the species‐rich group shows a faster rate of trait change. We studied 24 Liolaemus lizard species, belonging to two monophyletic groups that differ in species number: L. elongatus' clade, which has more species than L. lineomaculatus' clade. We studied 20 different traits (9 genital and 11 non‐genital) and calculated their phylogenetic signal, evolutionary rate of change and models that best explain the evolutionary change. Our results show that: (1) in general, genital traits evolve faster than non‐genital ones in both groups, and both phylogenetic signal and best evolutionary model vary depending on the trait. (2) Genital traits diverged more among sister species within the L. lineomaculatus group, but within the L. elongatus group, both sets of traits show similar degrees of divergence. Finally, (3) the species‐rich group (L. elongatus), has the highest genital evolutionary rate but also the highest non‐genital evolutionary rate.

Early Diverging and Core Bromelioideae (Bromeliaceae) Reveal Contrasting Patterns of Genome Size Evolution and Polyploidy.

The subfamily Bromelioideae is one of the most diverse groups among the neotropical Bromeliaceae. Previously, key innovations have been identified which account for the extraordinary radiation and species richness of this subfamily, especially in the so-called core Bromelioideae. However, in order to extend our understanding of the evolutionary mechanisms, the genomic mechanisms (e.g. polyploidy, dysploidy) that potentially underlie this accelerated speciation also need to be tested. Here, using PI and DAPI staining and flow cytometry we estimated genome size and GC content of 231 plants covering 30 genera and 165 species and combined it with published data. The evolutionary and ecological significance of all three genomic characters was tested within a previously generated dated phylogenetic framework using ancestral state reconstructions, comparative phylogenetic methods, and multiple regressions with climatic variables. The absolute genome size (2C) of Bromelioideae varied between 0.59 and 4.11 pg, and the GC content ranged between 36.73 and 41.43%. The monoploid genome sizes (Cx) differed significantly between core and early diverging lineages. The occurrence of dysploidy and polyploidy was, with few exceptions, limited to the phylogenetically isolated early diverging tank-less lineages. For Cx and GC content Ornstein–Uhlenbeck models outperformed the Brownian motion models suggesting adaptive potential linked to the temperature conditions. 2C-values revealed different rates of evolution in core and early diverging lineages also related to climatic conditions. Our results suggest that polyploidy is not associated with higher net diversification and fast radiation in core bromelioids. On the other hand, although coupled with higher extinction rates, dysploidy, polyploidy, and resulting genomic reorganizations might have played a role in the survival of the early diverging bromelioids in hot and arid environments.

Differential Genome Size and Repetitive DNA Evolution in Diploid Species of Melampodium sect. Melampodium (Asteraceae).

Plant genomes vary greatly in composition and size mainly due to the diversity of repetitive DNAs and the inherent propensity for their amplification and removal from the host genome. Most studies addressing repeatome dynamics focus on model organisms, whereas few provide comprehensive investigations across the genomes of related taxa. Herein, we analyze the evolution of repeats of the 13 species in Melampodium sect. Melampodium, representing all but two of its diploid taxa, in a phylogenetic context. The investigated genomes range in size from 0.49 to 2.27 pg/1C (ca. 4.5-fold variation), despite having the same base chromosome number (x = 10) and very strong phylogenetic affinities. Phylogenetic analysis performed in BEAST and ancestral genome size reconstruction revealed mixed patterns of genome size increases and decreases across the group. High-throughput genome skimming and the RepeatExplorer pipeline were utilized to determine the repeat families responsible for the differences in observed genome sizes. Patterns of repeat evolution were found to be highly correlated with phylogenetic position, namely taxonomic series circumscription. Major differences found were in the abundances of the SIRE (Ty1-copia), Athila (Ty3-gypsy), and CACTA (DNA transposon) lineages. Additionally, several satellite DNA families were found to be highly group-specific, although their overall contribution to genome size variation was relatively small. Evolutionary changes in repetitive DNA composition and genome size were complex, with independent patterns of genome up- and downsizing throughout the evolution of the analyzed diploids. A model-based analysis of genome size and repetitive DNA composition revealed evidence for strong phylogenetic signal and differential evolutionary rates of major lineages of repeats in the diploid genomes.

Evolutionary rates of and selective constraints on the mitochondrial genomes of Orthoptera insects with different wing types

Orthoptera is the most diverse order of polyneopterans, and the forewing and hindwing of its members exhibit extremely variability from full length to complete loss in many groups; thus, this order provides a good model for studying the effects of insect flight ability on the evolutionary constraints on and evolutionary rate of the mitochondrial genome. Based on a data set of mitochondrial genomes from 171 species, including 43 newly determined, we reconstructed Orthoptera phylogenetic relationships and estimated the divergence times of this group. The results supported Caelifera and Ensifera as two monophyletic groups, and revealed that Orthoptera originated in the Carboniferous (298.997 Mya). The date of divergence between the suborders Caelifera and Ensifera was 255.705 Mya, in the late Permian. The major lineages of Acrididae seemed to have radiated in the Cenozoic, and the six patterns of rearrangement of 171 Orthoptera mitogenomes mostly occurred in the Cretaceous and Cenozoic. Based on phylogenetic relationships and ancestral state reconstruction, we analysed the evolutionary selection pressure on and evolutionary rate of mitochondrial protein-coding genes (mPCGs). The results indicated that during approximately 300 Mya of evolution, these genes experienced purifying selection to maintain their function. Flightless orthopteran insects accumulated more non-synonymous mutations than flying species and experienced more relaxed evolutionary constraints. The different wing types had different evolutionary rates, and the mean evolutionary rate of Orthoptera mitochondrial mPCGs was 13.554 × 10−9 subs/s/y. The differences in selection pressures and evolutionary rates observed between the mitochondrial genomes suggested that functional constraints due to locomotion play an important role in the evolution of mitochondrial DNA in orthopteran insects with different wing types.

Childhood immune imprinting to influenza A shapes birth year-specific risk during seasonal H1N1 and H3N2 epidemics.

Across decades of co-circulation in humans, influenza A subtypes H1N1 and H3N2 have caused seasonal epidemics characterized by different age distributions of cases and mortality. H3N2 causes the majority of severe, clinically attended cases in high-risk elderly cohorts, and the majority of overall deaths, whereas H1N1 causes fewer deaths overall, and cases shifted towards young and middle-aged adults. These contrasting age profiles may result from differences in childhood imprinting to H1N1 and H3N2 or from differences in evolutionary rate between subtypes. Here we analyze a large epidemiological surveillance dataset to test whether childhood immune imprinting shapes seasonal influenza epidemiology, and if so, whether it acts primarily via homosubtypic immune memory or via broader, heterosubtypic memory. We also test the impact of evolutionary differences between influenza subtypes on age distributions of cases. Likelihood-based model comparison shows that narrow, within-subtype imprinting shapes seasonal influenza risk alongside age-specific risk factors. The data do not support a strong effect of evolutionary rate, or of broadly protective imprinting that acts across subtypes. Our findings emphasize that childhood exposures can imprint a lifelong immunological bias toward particular influenza subtypes, and that these cohort-specific biases shape epidemic age distributions. As a consequence, newer and less "senior" antibody responses acquired later in life do not provide the same strength of protection as responses imprinted in childhood. Finally, we project that the relatively low mortality burden of H1N1 may increase in the coming decades, as cohorts that lack H1N1-specific imprinting eventually reach old age.

Postembryonic Development of Cladocera Males of the Family Podonidae Mordukhai-Boltovskoi, 1968 (Cladocera, Onychopoda)

The postembryonic development of males of marine Cladocera (Crustacea) of the family Podonidae was studied for the first time, based on two common species from the White Sea, Evadne nordmanni Loven 1836 and Pleopis polyphemoides (Leuckart 1859). The development of podonid males includes two juvenile stages which differ in the size of the testes, the shape and size of the penis, and the morphology of the copulatory hook on the distal segment of the endopodite of thoracic limb I. The females of Podonidae are pedogenetic and lack juvenile stages. Males of Podonidae retain the plesiomorphic postembryonal development, suggesting different rates of evolution in both sexes in this family.

Post K-Pg diversification of the mammalian order Eulipotyphla as suggested by phylogenomic analyses of ultra-conserved elements

The origin of the mammalian order Eulipotyphla has been debated intensively with arguments around whether they began diversifying before or after the Cretaceous-Palaeogene (K-Pg) boundary at 66 Ma. Here, we used an in-solution nucleotide capture method and next generation DNA sequencing to determine the sequence of hundreds of ultra-conserved elements (UCEs), and conducted phylogenomic and molecular dating analyses for the four extant eulipotyphlan lineages—Erinaceidae, Solenodontidae, Soricidae, and Talpidae. Concatenated maximum-likelihood analyses with single or partitioned models and a coalescent species-tree analysis showed that divergences among the four major eulipotyphlan lineages occurred within a short period of evolutionary time, but did not resolve the interrelationships among them. Alternative suboptimal phylogenetic hypotheses received consistently the same amount of support from different UCE loci, and were not significantly different from the maximum likelihood tree topology, suggesting the prevalence of stochastic lineage sorting. Molecular dating analyses that incorporated among-lineage evolutionary rate differences supported a scenario where the four eulipotyphlan families diversified between 57.8 and 63.2 Ma. Given short branch lengths with low support values, traces of rampant genome-wide stochastic lineage sorting, and post K-Pg diversification, we concluded that the crown eulipotyphlan lineages arose through a rapid diversification after the K-Pg boundary when novel niches were created by the mass extinction of species.

Speciation of cultivated temperate and tropical Pleurotus species –an in silico prediction using conserved sequences

Conserved sequences are molecular fossil record of any organism and their evolution can be traced using the sequences. Genus Pleurotus is one of the most cultivated mushroom species for food purposes. To trace the evolution in Pleurotus species, two conserved sequences (ITS and LSU) of twenty three species were analysed and phylogenetic relationship was studied to see evolution pattern linked with domestication. The sequence were analysed using the Molecular Evolutionary Genetics Analysis (MEGA) and Bayesian phylogenetic analysis (Mr. Bayes). Analysis of LSU and ITS categorized the Pleurotus species into two groups. The differences in evolutionary rate among various species were ascertained. Grouping of sequences confirmed that cultivation requirements of Pleurotus species might be due to allopatric variation occurred during the evolution. Our results suggested that summer and winter species of Pleurotus might be evolved at different evolutionary rate.

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.

Comprehensive Genomic Survey, Characterization and Expression Analysis of the HECT Gene Family in Brassica rapa L. and Brassica oleracea L.

The HECT-domain protein family is one of the most important classes of E3 ligases. While the roles of this family in human diseases have been intensively studied, the information for plant HECTs is limited. In the present study, we performed the identification of HECT genes in Brassica rapa and Brassica oleracea, followed by analysis of phylogeny, gene structure, additional domains, putative cis-regulatory elements, chromosomal location, synteny, and expression. Ten and 13 HECT genes were respectively identified in B. rapa and B. oleracea and then resolved into seven groups along with their Arabidopsis orthologs by phylogenetic analysis. This classification is well supported by analyses of gene structure, motif composition within the HECT domain and additional protein domains. Ka/Ks ratio analysis showed that these HECT genes primarily underwent purifying selection with varied selection pressures resulting in different rates of evolution. RNA-Seq data analysis showed that the overwhelming majority of them were constitutively expressed in all tested tissues. qRT-PCR based expression analysis of the 10 B. rapa HECT genes under salt and drought stress conditions showed that all of them were responsive to the two stress treatments, which was consistent with their promoter sequence analysis revealing the presence of an important number of phytohormone-responsive and stress-related cis-regulatory elements. Our study provides useful information and lays the foundation for further functional determination of each HECT gene across the Brassica species.

A simple method for data partitioning based on relative evolutionary rates.

BackgroundMultiple studies have demonstrated that partitioning of molecular datasets is important in model-based phylogenetic analyses. Commonly, partitioning is done a priori based on some known properties of sequence evolution, e.g. differences in rate of evolution among codon positions of a protein-coding gene. Here we propose a new method for data partitioning based on relative evolutionary rates of the sites in the alignment of the dataset being analysed. The rates are inferred using the previously published Tree Independent Generation of Evolutionary Rates (TIGER), and the partitioning is conducted using our novel python script RatePartitions. We conducted simulations to assess the performance of our new method, and we applied it to eight published multi-locus phylogenetic datasets, representing different taxonomic ranks within the insect order Lepidoptera (butterflies and moths) and one phylogenomic dataset, which included ultra-conserved elements as well as introns.MethodsWe used TIGER-rates to generate relative evolutionary rates for all sites in the alignments. Then, using RatePartitions, we partitioned the data into partitions based on their relative evolutionary rate. RatePartitions applies a simple formula that ensures a distribution of sites into partitions following the distribution of rates of the characters from the full dataset. This ensures that the invariable sites are placed in a partition with slowly evolving sites, avoiding the pitfalls of previously used methods, such as k-means. Different partitioning strategies were evaluated using BIC scores as calculated by PartitionFinder.ResultsSimulations did not highlight any misbehaviour of our partitioning approach, even under difficult parameter conditions or missing data. In all eight phylogenetic datasets, partitioning using TIGER-rates and RatePartitions was significantly better as measured by the BIC scores than other partitioning strategies, such as the commonly used partitioning by gene and codon position. We compared the resulting topologies and node support for these eight datasets as well as for the phylogenomic dataset.DiscussionWe developed a new method of partitioning phylogenetic datasets without using any prior knowledge (e.g. DNA sequence evolution). This method is entirely based on the properties of the data being analysed and can be applied to DNA sequences (protein-coding, introns, ultra-conserved elements), protein sequences, as well as morphological characters. A likely explanation for why our method performs better than other tested partitioning strategies is that it accounts for the heterogeneity in the data to a much greater extent than when data are simply subdivided based on prior knowledge.

Remarkable geographical variations between India and Europe in carriage of the staphylococcal surface protein-encoding sasX/sesI and in the population structure of methicillin-resistant Staphylococcus aureus belonging to clonal complex 8

ObjectivessasX is a colonization-virulence factor that potentially underlies the success of methicillin-resistant Staphylococcus aureus (MRSA) sequence type (ST) 239 in Asia. We aimed to study the spread of sasX and the population structure of MRSA in two geographically distinct regions, Europe and India. MethodsMRSA (n = 128) from screening and clinical samples from tertiary care patients in 12 European countries (n = 119), and from India (n = 9) were multilocus-sequence-typed and screened for sasX and its carrier φSPβ-like prophage by PCR. Whole genome sequencing was performed on sasX-harbouring strains from India (n = 5) and Europe (n = 2) and on a selection non-harbouring sasX (n = 36) (2 × 150 bp, Miseq, Illumina). Reads were mapped to the ST239 reference strain, TW20. ResultssasX and sesI, a sasX homologue native to Staphylococcus epidermidis, were detected in five of the nine Indian MRSA belonging to ST239 and to other sequence types of CC8. In contrast, sasX was restricted to two ST239 strains in Europe. The intact sasX and sesI carrier φSPβ-like prophages were ∼80 kb and ∼118 kb, and integrated in the yeeE gene. We identified ‘novel’ ST239 clades in India and Serbia that showed significant differences in base substitution frequencies (0.130 and 0.007, respectively, Tamura–Nei model) (p <0.05). ConclusionsOur data highlight dissemination of sasX to non-ST239 sequence types of CC8. Detection of the S. epidermidis-associated sesI in MRSA provided unquestionable evidence of transfer between the two species. Stark differences in evolutionary rates between the novel Indian and Serbian ST239 clades identified here might be due to inherent clade characteristics or influenced by other environmental differences such as antibiotic use.