Invariable Sites Research Articles

Genetic characteristics of Blastocystis sp. in cattle from Hebei Province, China

Blastocystis sp. is a protozoan parasite that infects the intestines of humans and animals, causing chronic diseases such as skin rashes, abdominal pain, and irritable bowel syndrome. A survey was conducted to determine the prevalence and genetic diversity of Blastocystis sp. infection in cattle, in Hebei Province, China. 2746 cattle fecal samples were collected from 11 cities in Hebei Province and analyzed using polymerase chain reaction targeting the Blastocystis sp. barcoding gene. MEGA, PhyloSuite, and PopART were used to analyze the subtype, sequence signature, pairwise genetic distance, and genetic diversity indices. The results showed that the Blastocystis sp. detection rate was 12.60% (346/2746). The infection rate in different herds was affected by region, age, breeding mode, and variety; that is, the infection rates in areas of southern Hebei, cattle under one year old, intensive raising, and dairy cattle were higher than the infection rates in northern Hebei, cattle over one year old, scatter feeding, and beef cattle. Seven Blastocystis subtypes were identified, namely, ST1, ST2, ST5, ST10, ST14, ST21, and ST26; ST10 was the dominant subtype, and ST14 was the second most common subtype. A total of 374 polymorphic and conserved sites were obtained, including 273 invariable (monomorphic) sites and 101 variable (polymorphic) sites, accounting for 27.01% of all nucleotides. The nucleotide diversity index (Pi) was 0.07749, and the haplotype (gene) diversity index (Hd) was 0.946. This study provides the first comprehensive information on the epidemiological situation of Blastocystis sp. infection in cattle from Hebei Province, China, and revealed rich genetic diversity of Blastocystis sp.

DNA Barcoding of Prunus Species Collection Conserved in the National Gene Bank of Egypt

Two intergenic spacers cpDNA barcoding regions were used to assess the genetic diversity and phylogenetic structure of a collection of 25 Prunus accessions. The trnH-psbA and trnL-trnF intergenic spacers were able to distinguish and identify only four Prunus species. The average aligned length was 316–352 bp and 701–756 bp for trnH-psbA and trnL-trnF, respectively. The overall evolutionary divergence was higher in trnH-psbA than trnL-trnF. The transition/transversion bias (R) recorded as 0.59 in trnL-trnF and 0.89 in trnH-psbA. The number of invariable sites, nucleotide diversity (Pi), and the average number of nucleotide differences (k) was higher in the trnH-psbA region. The trnL-trnF records was above the other region in the number of variable sites, number of singleton variable sites, and the parsimony informative sites. Phylogenetic relationships among the 25 accessions of Prunus species were investigated. Most of the different Prunus species clustered in a homogenized distribution in both regions, except for the plum (P. domestica) accession (African Rose) was assigned with the peach (P. persica) accessions. The two intergenic cpDNA trnH-psbA and trnL-trnF were able to distinguish and identify the four Prunus species accessions.

ProteinEvolverABC: coestimation of recombination and substitution rates in protein sequences by approximate Bayesian computation.

MotivationThe evolutionary processes of mutation and recombination, upon which selection operates, are fundamental to understand the observed molecular diversity. Unlike nucleotide sequences, the estimation of the recombination rate in protein sequences has been little explored, neither implemented in evolutionary frameworks, despite protein sequencing methods are largely used.ResultsIn order to accommodate this need, here I present a computational framework, called ProteinEvolverABC, to jointly estimate recombination and substitution rates from alignments of protein sequences. The framework implements the approximate Bayesian computation approach, with and without regression adjustments and includes a variety of substitution models of protein evolution, demographics and longitudinal sampling. It also implements several nuisance parameters such as heterogeneous amino acid frequencies and rate of change among sites and, proportion of invariable sites. The framework produces accurate coestimation of recombination and substitution rates under diverse evolutionary scenarios. As illustrative examples of usage, I applied it to several viral protein families, including coronaviruses, showing heterogeneous substitution and recombination rates.Availability and implementation ProteinEvolverABC is freely available from https://github.com/miguelarenas/proteinevolverabc, includes a graphical user interface for helping the specification of the input settings, extensive documentation and ready-to-use examples. Conveniently, the simulations can run in parallel on multicore machines.Supplementary information Supplementary data are available at Bioinformatics online.

Genetic Variation Analysis of Four Local Varieties of Indonesian Black Rice (Oryza sativa L.) Based on Partially rbcL cpDNA Gene Sequence

Black rice ( Oryza sativa L.) varieties i.e. Toraja (South Sulawesi), Cempo Ireng (Yogyakarta), Wojalaka (East Nusa Tenggara), and Manggarai (East Nusa Tenggara) are four local black rice varieties in Indonesia whose character has not been widely studied, especially the character of genetic variation. Research aimed to determine the variation of the rbc L gene in the four local black rice varieties. The sample for testing the variation of the rbc L gene sequence in the form of black rice leaves six weeks after planting. Dendogram was carried out using the UPGMA method with the Kimura 2-parameter algorithmic calculation model using the MEGA5 version 5.2.2 program. The results showed that partially the rbc L gene sequence was successfully amplified on four black rice varieties with a sequence length of 487 bp. The partial rbc L sequence of black rice consisted of 26.58% tyrosine, 21.38% cytosine, 28.86% adenine, and 23.18% guanine. The value of G + C content was 0.446, with the frequency of invariable sites of 97.13%. The frequency of informative parsimony sites was 1.43% with a nucleotide diversity (Pi) value of 42-10, the number of haplotypes was 5, and the total number of mutations and polymorphic sites was 14. The ratio between transition and transversion (ts/tv ratio k) for purine bases was 1.741 and pyrimidine was 3.571, with the estimated overall ratio between transition and transversion (R) of 1.31. Based on the dendogram, the farthest genetic distance was found in Wojalaka and Manggarai varieties, which were 0.019 respectively. Keywords: black rice, genetic variation, local varieties, rbc L gen

Bioinformatic analyses of kappa casein gene in mammalian livestock species

Kappa casein (CSN3) gene is a variant of the milk protein highly conserved in mammalian species. Genetic variations in CSN3 gene of six mammalian livestock species were investigated using bioinformatics approach. A total of twenty-seven CSN3 gene sequences with corresponding amino acids belonging to the six species studied were retrieved from GenBank. Using a comparative genomic approach, we obtained 1797 bp of the CSN3 sequences from cattle, goat, horse, pig, rabbit and sheep. Alignment of twenty-seven sequences within the region of 1713 bp and containing gaps was carried out using Clustal W. The sequence alignment revealed high polymorphisms of CSN3 sequences among the six species. The percent identity and similarity between species were determined by conducting pair-wise comparisons of the sequences. Results revealed high level of identity and similarity among the ruminant animals (84-97% and 86-99%) which implied the gene is conserved in all the ruminants. However, percent identity and similarity between pig and rabbit were relatively low (45%, 61%). This revealed that CSN3 protein from them might be functionally different and divergence between them is more ancient. Higher polymorphism of amino acids was observed among the species than within species. Four non-synonymous mutations occurred in cattle, six in goat, four in horse, three in rabbit and seven in sheep. The results of DnaSP analysis indicated that the selected region (1-1797 bp) of the twenty-seven sequences from all the species had 1298 nucleotide sites excluding sites with gaps (499). There are invariable and variable sites that include 185 singleton variable sites and 56 parsimony informative sites. The nucleotide diversity (ð= 0.212) and average number of nucleotide diversity (k = 94.84) for all sequences which were lower than the highest values in rabbit (ð = 0.557, k = 229) were estimated. Divergence of the species could be inferred from high genetic diversity of the CSN3 gene. The haplotype diversity was 1.00 in goat, horse and sheep indicating abundant genetic diversity in those species. Close relationship between goat and sheep was observed in phylogenetic tree of CSN3 gene which showed that the comparability of CSN3 gene sequences was highest between the goat and sheep and they evolved from a most recent common ancestor. A more detailed study involving more divergent taxa and complete amino acid sequence will help to delineate how the mutations affect the functions of the gene.

OBAMA: OBAMA for Bayesian amino-acid model averaging.

BackgroundBayesian analyses offer many benefits for phylogenetic, and have been popular for analysis of amino acid alignments. It is necessary to specify a substitution and site model for such analyses, and often an ad hoc, or likelihood based method is employed for choosing these models that are typically of no interest to the analysis overall.MethodsWe present a method called OBAMA that averages over substitution models and site models, thus letting the data inform model choices and taking model uncertainty into account. It uses trans-dimensional Markov Chain Monte Carlo (MCMC) proposals to switch between various empirical substitution models for amino acids such as Dayhoff, WAG, and JTT. Furthermore, it switches base frequencies from these substitution models or use base frequencies estimated based on the alignment. Finally, it switches between using gamma rate heterogeneity or not, and between using a proportion of invariable sites or not.ResultsWe show that the model performs well in a simulation study. By using appropriate priors, we demonstrate both proportion of invariable sites and the shape parameter for gamma rate heterogeneity can be estimated. The OBAMA method allows taking in account model uncertainty, thus reducing bias in phylogenetic estimates. The method is implemented in the OBAMA package in BEAST 2, which is open source licensed under LGPL and allows joint tree inference under a wide range of models.

Sequencing the exon 1 of MSTN in Mongolian horse (Equus caballus)

Myostatin (MSTN) is a protein responsible for muscular tissue differentiation, development, and growth in mammals. Studies on cattle, mice, pig, and dog have provided sample demonstrations of mutations on myostatin, which strongly affected animal Ankole93200phenotype, particularly muscle development. In horse, myostatin gene is usually studied with respect to racing performance. Most Mongolian native horses do race for long distance and such racing is a very popular sport in Mongolia. On the other hand, the MSTN gene also affects meat yield and quality. We partially sequenced of MSTN gene to explore its polymorphic features in Mongolian horse. A 238 bp long segment of exon 1 was sequenced, using 3130 Genetic Analyser and MEGAX. The sequence was registered with accession number LC216412 in the Gene bank database of NCBI and aligned to four reference sequences (GQ183900.1, KC708233.1, AY840554.2, NM 001081817) from this database. The nucleotide diversity was 0.0084, Tajima test statistic was -1.12. so that there were 233 invariable (monomorphic) sites and 5 variable (polymorphic) sites, in the group of animals compared. The Mongolian horse differed from the reference sequence at 3 sites. There was an amino acid difference in the partial sequence, in which isoleucine in the reference sequence was converted into valine. The results of the study indicate that the Mongolian horse seems to differ substantially from the reference genome of equus caballus. Sequencing of the full length of exon 1 or even the whole MSTN gene, including introns, is strongly advised, to get more insights with regard to the identification and characterization of variants at the myostatin locus in Mongolian horses.

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.

Complex Models of Sequence Evolution Require Accurate Estimators as Exemplified with the Invariable Site Plus Gamma Model.

The invariable site plus n}{}Gamma model (In}{}+Gamma) is widely used to model rate heterogeneity among alignment sites in maximum likelihood and Bayesian phylogenetic analyses. The proof that the In}{}+ continuous n}{}Gamma model is identifiable (model parameters can be inferred correctly given enough data) has increased the creditability of its application to phylogeny reconstruction. However, most phylogenetic software implement the In}{}+ discrete n}{}Gamma model, whose identifiability is likely but unproven. How well the parameters of the In}{}+ discrete n}{}Gamma model are estimated is still disputed. Especially the correlation between the fraction of invariable sites and the fractions of sites with a slow evolutionary rate is discussed as being problematic. We show that optimization heuristics as implemented in frequently used phylogenetic software (PhyML, RAxML, IQ-TREE, and MrBayes) cannot always reliably estimate the shape parameter, the proportion of invariable sites, and the tree length. Here, we propose an improved optimization heuristic that accurately estimates the three parameters. While research efforts mainly focus on tree search methods, our results signify the equal importance of verifying and developing effective estimation methods for complex models of sequence evolution.

CoMuS: simulating coalescent histories and polymorphic data from multiple species.

The simultaneous analysis of intra- and interspecies variation is challenging mainly because our knowledge about patterns of polymorphisms where both intra- and interspecies samples coexist is limited. In this study, we present CoMuS (Coalescent of Multiple Species), a multispecies coalescent software that can simulate intra- and interspecies polymorphisms. CoMuS supports a variety of speciation models and demographic scenarios related to the history of each species. In CoMuS, speciation can be accompanied by either instant or gradual isolation between sister species. Sampling may also occur in the past, and thus, we can study simultaneously extinct and extant species. Our software supports both the infinite- and the finite-site model, with substitution rate heterogeneity among sites and a user-defined proportion of invariable sites. We demonstrate the usage of CoMuS in various applications: species delimitation, software testing, model selection and parameter inference involving present-day and ancestral samples, comparison between gradual and instantaneous isolation models, estimation of speciation time between human and chimpanzee using both intra- and interspecies variation. We expect that CoMuS will be particularly useful for studies where species have been separated recently from their common ancestor and phenomena such as incomplete lineage sorting or introgression still occur.

The first taxonomic revaluation of the Iranian water frogs of the genus Pelophylax (Anura: Ranidae) using sequences of the mitochondrial genome

The Eurasian water frog species and their geographic ranges have undergone considerable changes in the last four decades, but the Iranian populations have largely remained unknown. All the Iranian populations of water frogs, despite their vast distribution range have attributed to a single species: Rana ridibunda. In order to understand the phylogenetic relationships and taxonomic status of water frogs of Iran, we collected samples from many populations across the country and used the mitochondrial DNA sequence variation. A data set with a final sequence length of 616 nucleotides was generated for Cyt b from 70 individuals of Pelophylax in which there are 422 invariable sites, 174 variable sites of which 123 were parsimony informative. In total, 43 haplotypes were found (Hd: 0.9752). The result demonstrated that, two major clades with strong support can be identified within the Iranian water frogs. One of these clades that include north western and southwestern populations forms a monophyletic group along with P. bedriagae samples from Turkey. The second clade consists of water frog populations of north and northeastern parts of Iran which in turn is subdivided into two subclades. Inclusion of water frog samples from adjacent areas showed that the second clade of our study is, most likely, a distinct taxonomic entity at species rank with its two subclades indicating two diagnosable subspecies for the clade. In conclusion, we suggest that two distinct species, P. bedriagae and Pelophylax sp., with its two subspecies, should be identified as water frogs of Iran. In Addition, another traditionally reported water frog of Iran, P.ridibundus, most likely should be excluded from the Iranian water frog’s checklist.

First Report of the Spiral Nematode Rotylenchus incultus (Nematoda: Hoplolaimidae) from Cultivated Olive in Tunisia, with Additional Molecular Data on Rotylenchus eximius.

Spiral nematode species of the genus Rotylenchus have been reported on olive (Olea europaea L.) in several Mediterranean countries (Castillo et al., 2010; Ali et al., 2014). Nematological surveys for plant-parasitic nematodes on olive trees were carried out in Tunisia between 2013 and 2014, and two nematode species of Rotylenchus were collected from the rhizosphere of olive cv. Chemlali in several localities of Tunisia (Tables 1,2 [Table: see text] [Table: see text] ). Twenty-two soil samples of 3 to 4 kg were collected with a shovel from the upper 50 cm of soil from arbitrarily chosen olive trees. Nematodes were extracted from 500 cm3 of soil by centrifugal flotation method (Coolen, 1979). Specimens were heat killed by adding hot 4% formaldehyde solution and processed to pure glycerin using the De Grisse's (1969) method. Measurements were done using a drawing tube attached to a Zeiss III compound microscope. Nematode DNA was extracted from single individuals and PCR assays were conducted as described by Castillo et al. (2003). Moderate-to-low soil populations of these spiral nematodes were detected (5.5-11.5, 1.5-5.0 individuals/500 cm3 of soil, respectively). This prompted us to undertake a detailed morphological and molecular comparative study with previous reported data. Morphological and molecular analyses of females identified these species as Rotylenchus eximius Siddiqi, 1964, and Rotylenchus incultus Sher, 1965. The morphology of R. eximius females (five specimens studied) was characterized by having a hemispherical lip region clearly off set, with four to five annuli, body without longitudinal striations, lateral fields areolated in the pharyngeal region only, stylet 32 to 36 μm long, and broadly rounded tail. The morphology of R. incultus females (51 females and 16 males; Table 2) was characterized by a hemispherical lip region with the basal annulus subdivided by irregular longitudinal striations, with three, rarely four annuli; stylet 21.5 to 27.5 μm long, female tail hemispherical with terminus regularly annulated; phasmids anterior to anus level (3-6 annuli above). The morphology of the isolated nematodes agreed with previous descriptions of R. eximius (Siddiqi, 1964; Castillo and Vovlas, 2005) and R. incultus (Sher, 1965; Castillo and Vovlas, 2005; Vovlas et al., 2008), respectively. A single individual was used for DNA extraction. Primers and PCR conditions used in this research were specified in Cantalapiedra-Navarrete et al. (2013), and a single amplicon of 800, 1,100, and 450 bp was obtained and sequenced for D2 to D3, ITS1, and cytochrome c oxidase subunit 1 (coxI), respectively. Sequence alignments for D2 to D3 (KX669231-KX669233), ITS1 (KX669238-KX669240), and coxI (KX669244-KX669245) from R. eximius, showed 99% to 97%, 98% to 94%, 93% similarity to other sequences of R. eximius deposited in GenBank (EU280794-DQ328741, EU373663-EU373664, JX015401-JX015402, respectively). Similarly, D2 to D3 (KX669234-KX669237), ITS1 (KX669241-KX669243), and coxI (KX669246-KX669249) sequence alignments from R. incultus, showed 99%, 99% to 95%, 99% to 90% similarity, respectively, to other sequences of R. incultus deposited in GenBank (EU280797, EU373672-EU373673, JX015403, respectively). The best fitted model of DNA evolution was obtained using jModelTest v. 2.1.7 (Darriba et al. 2012) with the Akaike information criterion. BI analyses were performed under the general time reversible (GTR) with invariable sites and a gamma-shaped distribution of substitution rates (GTR + I + G) model for ITS1 and coxI. Phylogenetic analyses of ITS1 and coxI using Bayesian inference (BI) placed R. eximius and R. incultus from Tunisia in subclades that included all R. eximius and R. incultus sequences deposited in GenBank (Fig. 1Fig. 1Phylogenetic relationships within Rotylenchus species found in Tunisia and other species from GenBank. Bayesian 50% majority rule consensus trees as inferred from ITS1 and coxI sequences alignments under the GTR + I + G model. Posterior probabilities more than 0.70 are given for appropriate clades. Newly obtained sequences in this study are in bold. Scale bar = expected changes per site.), which agrees with previous results (Cantalapiedra-Navarrete et al., 2013). Morphology, morphometry, and molecular and phylogenetic data obtained from these samples were consistent with R. eximius and R. incultus identification. To our knowledge, this is the first report of R. incultus in Tunisia. Consequently, all these data suggest that spiral nematode species of the genus Rotylenchus are predominant in olive as previously reported in other Mediterranean areas (Ali et al., 2014).

Detecting recent selective sweeps while controlling for mutation rate and background selection.

A composite likelihood ratio test implemented in the program sweepfinder is a commonly used method for scanning a genome for recent selective sweeps. sweepfinder uses information on the spatial pattern (along the chromosome) of the site frequency spectrum around the selected locus. To avoid confounding effects of background selection and variation in the mutation process along the genome, the method is typically applied only to sites that are variable within species. However, the power to detect and localize selective sweeps can be greatly improved if invariable sites are also included in the analysis. In the spirit of a Hudson–Kreitman–Aguadé test, we suggest adding fixed differences relative to an out‐group to account for variation in mutation rate, thereby facilitating more robust and powerful analyses. We also develop a method for including background selection, modelled as a local reduction in the effective population size. Using simulations, we show that these advances lead to a gain in power while maintaining robustness to mutation rate variation. Furthermore, the new method also provides more precise localization of the causative mutation than methods using the spatial pattern of segregating sites alone.

The application of fluorescence in situ hybridization in different ploidy levels cross-breeding of lily.

21 crossing were conducted between Asiatic Lily with different ploidy levels, the results showed that the interploidy hybridization between diploid and tetraploid lilies was not as successful as intraploidy hybridization. Regardless of male sterility, triploid lilies could be used as female parents in the hybridization which the progenies were aneuploidy. 3x×4x crosses could be cultured more successfully than 3x×2x crosses. 45S rDNA was mapped on the chromosomes of seven Lilium species and their progenies using fluorescence in situ hybridization (FISH). FISH revealed six to sixteen 45S rDNA gene loci, and normally the sites were not in pairs. The asymmetry indexes of LA (Longiflorum hybrids × Asiatic hybrids) hybrids was higher than Asiatic hybrids, the evolution degree was LA hybrids > Asiatic hybrids. 45S rDNA distributed variably on chromosome 1-10 and 12 among Asiatic hybrids. Chromosome 1 had invariable sites of 45S rDNA in all Asiatic hybrids, which could be considered as the characteristic of Asiatic hybrids. LA hybrid ‘Freya’ had two sites of 45S rDNA on one homologous chromosome 5, and also it could be found in the progenies. The karyotype and fluorescence in situ hybridization with 45S rDNA as probe were applied to identify the different genotypes of 9 hybrids. Typical chromosomes with parental signal sites could be observed in all the genotypes of hybrids, it was confirmed that all the hybrids were true.

Identifiability of the unrooted species tree topology under the coalescent model with time-reversible substitution processes, site-specific rate variation, and invariable sites.

The inference of the evolutionary history of a collection of organisms is a problem of fundamental importance in evolutionary biology. The abundance of DNA sequence data arising from genome sequencing projects has led to significant challenges in the inference of these phylogenetic relationships. Among these challenges is the inference of the evolutionary history of a collection of species based on sequence information from several distinct genes sampled throughout the genome. It is widely accepted that each individual gene has its own phylogeny, which may not agree with the species tree. Many possible causes of this gene tree incongruence are known. The best studied is the incomplete lineage sorting, which is commonly modeled by the coalescent process. Numerous methods based on the coalescent process have been proposed for the estimation of the phylogenetic species tree given DNA sequence data. However, use of these methods assumes that the phylogenetic species tree can be identified from DNA sequence data at the leaves of the tree, although this has not been formally established. We prove that the unrooted topology of the n-leaf phylogenetic species tree is generically identifiable given observed data at the leaves of the tree that are assumed to have arisen from the coalescent process under a time-reversible substitution process with the possibility of site-specific rate variation modeled by the discrete gamma distribution and a proportion of invariable sites.

Mixture models of nucleotide sequence evolution that account for heterogeneity in the substitution process across sites and across lineages.

Molecular phylogenetic studies of homologous sequences of nucleotides often assume that the underlying evolutionary process was globally stationary, reversible, and homogeneous (SRH), and that a model of evolution with one or more site-specific and time-reversible rate matrices (e.g., the GTR rate matrix) is enough to accurately model the evolution of data over the whole tree. However, an increasing body of data suggests that evolution under these conditions is an exception, rather than the norm. To address this issue, several non-SRH models of molecular evolution have been proposed, but they either ignore heterogeneity in the substitution process across sites (HAS) or assume it can be modeled accurately using the distribution. As an alternative to these models of evolution, we introduce a family of mixture models that approximate HAS without the assumption of an underlying predefined statistical distribution. This family of mixture models is combined with non-SRH models of evolution that account for heterogeneity in the substitution process across lineages (HAL). We also present two algorithms for searching model space and identifying an optimal model of evolution that is less likely to over- or underparameterize the data. The performance of the two new algorithms was evaluated using alignments of nucleotides with 10 000 sites simulated under complex non-SRH conditions on a 25-tipped tree. The algorithms were found to be very successful, identifying the correct HAL model with a 75% success rate (the average success rate for assigning rate matrices to the tree's 48 edges was 99.25%) and, for the correct HAL model, identifying the correct HAS model with a 98% success rate. Finally, parameter estimates obtained under the correct HAL-HAS model were found to be accurate and precise. The merits of our new algorithms were illustrated with an analysis of 42 337 second codon sites extracted from a concatenation of 106 alignments of orthologous genes encoded by the nuclear genomes of Saccharomyces cerevisiae, S. paradoxus, S. mikatae, S. kudriavzevii, S. castellii, S. kluyveri, S. bayanus, and Candida albicans. Our results show that second codon sites in the ancestral genome of these species contained 49.1% invariable sites, 39.6% variable sites belonging to one rate category (V1), and 11.3% variable sites belonging to a second rate category (V2). The ancestral nucleotide content was found to differ markedly across these three sets of sites, and the evolutionary processes operating at the variable sites were found to be non-SRH and best modeled by a combination of eight edge-specific rate matrices (four for V1 and four for V2). The number of substitutions per site at the variable sites also differed markedly, with sites belonging to V1 evolving slower than those belonging to V2 along the lineages separating the seven species of Saccharomyces. Finally, sites belonging to V1 appeared to have ceased evolving along the lineages separating S. cerevisiae, S. paradoxus, S. mikatae, S. kudriavzevii, and S. bayanus, implying that they might have become so selectively constrained that they could be considered invariable sites in these species.

The impact of modelling rate heterogeneity among sites on phylogenetic estimates of intraspecific evolutionary rates and timescales.

Phylogenetic analyses of DNA sequence data can provide estimates of evolutionary rates and timescales. Nearly all phylogenetic methods rely on accurate models of nucleotide substitution. A key feature of molecular evolution is the heterogeneity of substitution rates among sites, which is often modelled using a discrete gamma distribution. A widely used derivative of this is the gamma-invariable mixture model, which assumes that a proportion of sites in the sequence are completely resistant to change, while substitution rates at the remaining sites are gamma-distributed. For data sampled at the intraspecific level, however, biological assumptions involved in the invariable-sites model are commonly violated. We examined the use of these models in analyses of five intraspecific data sets. We show that using 6–10 rate categories for the discrete gamma distribution of rates among sites is sufficient to provide a good approximation of the marginal likelihood. Increasing the number of gamma rate categories did not have a substantial effect on estimates of the substitution rate or coalescence time, unless rates varied strongly among sites in a non-gamma-distributed manner. The assumption of a proportion of invariable sites provided a better approximation of the asymptotic marginal likelihood when the number of gamma categories was small, but had minimal impact on estimates of rates and coalescence times. However, the estimated proportion of invariable sites was highly susceptible to changes in the number of gamma rate categories. The concurrent use of gamma and invariable-site models for intraspecific data is not biologically meaningful and has been challenged on statistical grounds; here we have found that the assumption of a proportion of invariable sites has no obvious impact on Bayesian estimates of rates and timescales from intraspecific data.

Ab-initio tiling and atomic structure for decagonal ZnMgY quasicrystal

We discover the detailed atomic structure of -MgZnY, a stable decagonal quasicrystal alloy of the layered Frank–Kasper type, and related phases, using the ‘tiling and decoration’ approach. The atoms have invariable sites in the rectangle and triangle tiles of a 10-fold-symmetric planar tiling. To discover the lowest-energy structures, we combine the methods of density functional theory (DFT) total energy calculations, empirical oscillating pair potentials (fitted to DFT), fitting effective Hamiltonians for tilings, and discovering optimum tiling structures using a nonlocal tile-reshuffling algorithm. We find a family of almost stable binary compounds with varying composition, including the decagonal quasicrystal and the known MgZn; in addition, we show the ternary versions can be stable: -MgZnY at high temperature, and -MgZnSc even at .

Assessing Molecular Signature for Some Potential Date (Phoenix dactylifera L.) Cultivars from Saudi Arabia, Based on Chloroplast DNA Sequences rpoB and psbA-trnH

Phoenix dactylifera L. (date palm), being economically very important, is widely cultivated in the Middle East and North Africa, having about 400 different cultivars. Assessment of date cultivars under trading and farming is a widely accepted problem owing to lack of a unique molecular signature for specific date cultivars. In the present study, eight different cultivars of dates viz., Khodry, Khalas, Ruthana, Sukkari, Sefri, Segae, Ajwa and Hilali were sequenced for rpoB and psbA-trnH genes and analyzed using bioinformatics tools to establish a cultivar-specific molecular signature. The combined aligned data matrix was of 1147 characters, of which invariable and variable sites were found to be 958 and 173, respectively. The analysis clearly reveals three major groups of these cultivars: (i) Khodary, Sefri, Ajwa, Ruthana and Hilali (58% BS); (ii) Sukkari and Khalas (64% BS); and (iii) Segae. The economically most important cultivar Ajwa showed similarity with Khodary and Sefri (67% BS).The sequences of the date cultivars generated in the present study showed bootstrap values between 38% and 70% so these sequences could be carefully used as molecular signature for potential date cultivars under trading and selection of genuine cultivars at the seedling stage for farming.

Two Stationary Nonhomogeneous Markov Models of Nucleotide Sequence Evolution

The general Markov model (GMM) of nucleotide substitution does not assume the evolutionary process to be stationary, reversible, or homogeneous. The GMM can be simplified by assuming the evolutionary process to be stationary. A stationary GMM is appropriate for analyses of phylogenetic data sets that are compositionally homogeneous; a data set is considered to be compositionally homogeneous if a statistical test does not detect significant differences in the marginal distributions of the sequences. Though the general time-reversible (GTR) model assumes stationarity, it also assumes reversibility and homogeneity. We propose two new stationary and nonhomogeneous models--one constrains the GMM to be reversible, whereas the other does not. The two models, coupled with the GTR model, comprise a set of nested models that can be used to test the assumptions of reversibility and homogeneity for stationary processes. The two models are extended to incorporate invariable sites and used to analyze a seven-taxon hominoid data set that displays compositional homogeneity. We show that within the class of stationary models, a nonhomogeneous model fits the hominoid data better than the GTR model. We note that if one considers a wider set of models that are not constrained to be stationary, then an even better fit can be obtained for the hominoid data. However, the methods for reducing model complexity from an extremely large set of nonstationary models are yet to be developed.