Nucleotide Substitution Rates Research Articles

Evolutionary Dynamics of Transferred Sequences Between Organellar Genomes in Cucurbita.

Twenty-nine DNA regions of plastid origin have been previously identified in the mitochondrial genome of Cucurbita pepo (pumpkin; Cucurbitaceae). Four of these regions harbor homolog sequences of rbcL, matK, rpl20-rps12 and trnL-trnF, which are widely used as molecular markers for phylogenetic and phylogeographic studies. We extracted the mitochondrial copies of these regions based on the mitochondrial genome of C. pepo and, along with published sequences for these plastome markers from 13 Cucurbita taxa, we performed phylogenetic molecular analyses to identify inter-organellar transfer events in the Cucurbita phylogeny and changes in their nucleotide substitution rates. Phylogenetic reconstruction and tree selection tests suggest that rpl20 and rbcL mitochondrial paralogs arose before Cucurbita diversification whereas the mitochondrial matK and trnL-trnF paralogs emerged most probably later, in the mesophytic Cucurbita clade. Nucleotide substitution rates increased one order of magnitude in all the mitochondrial paralogs compared to their original plastid sequences. Additionally, mitochondrial trnL-trnF sequences obtained by PCR from nine Cucurbita taxa revealed higher nucleotide diversity in the mitochondrial than in the plastid copies, likely related to the higher nucleotide substitution rates in the mitochondrial region and loss of functional constraints in its tRNA genes.

New molecular evolutionary characteristics of H9N2 avian influenza virus in Guangdong Province, China

To understand the evolution of H9N2 avian influenza virus genotype and its molecular evolution rate, we systematically analyzed 72 H9N2 avian influenza virus sequences isolated from Guangdong province from 2014 to 2018. We found three genotypes (G57, G68, and G118) of the H9N2 avian influenza virus, of which G118 is a newly discovered genotype and G57 is the dominant genotype. The internal gene cassette of the G57 genotype H9N2 avian influenza virus is a stable combination that can easily transport internal genes to other novel avian influenza viruses, and the internal gene cassettes of the G68 and G118 are identical to those of G57.In addition, we estimated the nucleotide substitution rate of the HA and NA genes of the H9N2 influenza virus from 2014 to 2018.The nucleotide substitution rate of HA and NA genes showed an upward trend in 2015 and 2016. In the past two years, H9N2 avian influenza virus recombination has produced genotype G68, which disappeared in 2014 for one year. And very coincidentally, in 2015, there was a new genotype G118. We observed that the emergence of new genotypes was accompanied by a slight increase in overall nucleotide substitution rate. Therefore we hypothesize that the emergence of new genotypes could accelerate the molecular evolution rate of genes. Our research shows that the H9N2 avian influenza virus in Guangdong province has been undergoing intense evolution, demonstrating the need to strengthen influenza surveillance in the region.

Characteristics and evolutionary history of hepatitis B virus quasi-subgenotype B3 in Southeast Asia

To analyze the hepatitis B virus (HBV) quasi-subgenotype B3 characters and molecular evolution in Southeast Asia, 411 serum samples with HBsAg positive were collected from Xishuangbanna, China. After DNA extraction, PCR amplification and sequencing, a total of 183 HBV full-length genomes were obtained. Phylogenetic analysis showed 139 stains (76.0%) were genotype B, 41 strains were genotype C (22.4%) and 3 strains were genotype I (1.6%). Among genotype B, 34 sequences were identified as quasi-subgenotype B3. Quasi-subgenotype B3 sequences from this study and quasi-subgenotype B3 sequences from the GenBank (total of 141 complete genome sequences) were grouped into quasi-subgenotype B3 (B3, formerly B5, Chinese B6 and B7-9). Sixteen peculiar nucleotides distributed in quasi-subgenotype B3 were identified, which were differ from B1, B2, B4 and B5(formerly B6) (nt93 T, nt100C, nt355 G, nt843 T, nt861C, nt912C, nt929 T, nt930 G, nt1023 T, nt1041 T, nt2651C, nt2693 T, nt2970C, nt3054A, nt3087A and nt3171 G). Then Evolutionary dynamics analysis of HBV quasi-subgenotype B3 was performed. The mean rate of nucleotide substitution for HBV quasi-subgenotype B3 was estimated to be around 5.556–5.660 × 10−4 substitutions/site/year. Estimated time to most recent ancestor of quasi-subgenotype B3 was around the 1847–1945(95%HPD), and Yunnan strains might be the parental strains. The Bayesian sky plot showed a steady spreading of HBV quasi-genotype B3 from early of 1940s to 90 s. In summary, HBV quasi-subgenotype B3 infection is prevalent in Southeast Asia based on the current reports and still with a high prevalence rate based on the evolutionary dynamics analysis. Current vaccine and nucleotide analogues might have effective prevention and treatment for HBV quasi-subgenotype B3 based on the rare clinically relevant mutation sites included in quasi-subgenotype B3.

Rapid evolution of Mexican H7N3 highly pathogenic avian influenza viruses in poultry.

Highly pathogenic avian influenza (HPAI) virus subtype H7N3 has been circulating in poultry in Mexico since 2012 and vaccination has been used to control the disease. In this study, eight Mexican H7N3 HPAI viruses from 2015–2017 were isolated and fully sequenced. No evidence of reassortment was detected with other avian influenza (AI) viruses, but phylogenetic analyses show divergence of all eight gene segments into three genetic clusters by 2015, with 94.94 to 98.78 percent nucleotide homology of the HA genes when compared to the index virus from 2012. The HA protein of viruses from each cluster showed a different number of basic amino acids (n = 5–7) in the cleavage site, and six different patterns at the predicted N-glycosylation sites. Comparison of the sequences of the Mexican lineage H7N3 HPAI viruses and American ancestral wild bird AI viruses to characterize the virus evolutionary dynamics showed that the nucleotide substitution rates in PB2, PB1, PA, HA, NP, and NS genes greatly increased once the virus was introduced into poultry. The global nonsynonymous and synonymous ratios imply strong purifying selection driving the evolution of the virus. Forty-nine positively selected sites out of 171 nonsynonymous mutations were identified in the Mexican H7N3 HPAI viruses, including 7 amino acid changes observed in higher proportion in North American poultry origin AI viruses isolates than in wild bird-origin viruses. Continuous monitoring and molecular characterization of the H7N3 HPAI virus is important for better understanding of the virus evolutionary dynamics and further improving control measures including vaccination.

Genome-wide analysis of cotton C2H2-zinc finger transcription factor family and their expression analysis during fiber development

BackgroundC2H2-zinc finger protein family is commonly found in the plant, and it is known as the key actors in the regulation of transcription and vital component of chromatin structure. A large number of the C2H2-zinc finger gene members have not been well characterized based on their functions and structure in cotton. However, in other plants, only a few C2H2-zinc finger genes have been studied.ResultsIn this work, we performed a comprehensive analysis and identified 386, 196 and 195 C2H2-zinc finger genes in Gossypium hirsutum (upland cotton), Gossypium arboreum and Gossypium raimondii, respectively. Phylogenetic tree analysis of the C2H2-zinc finger proteins encoding the C2H2-zinc finger genes were classified into seven (7) subgroups. Moreover, the C2H2-zinc finger gene members were distributed in all cotton chromosomes though with asymmetrical distribution patterns. All the orthologous genes were detected between tetraploid and the diploid cotton, with 154 orthologous genes pair detected between upland cotton and Gossypium arboreum while 165 orthologous genes were found between upland cotton and Gossypium raimondii. Synonymous (Ks) and non-synonymous (Ka) nucleotide substitution rates (Ka/Ks) analysis indicated that the cotton C2H2-zinc finger genes were highly influenced mainly by negative selection, which maintained their protein levels after the duplication events. RNA-seq data and RT-qPCR validation of the RNA seq result revealed differential expression pattern of some the C2H2-zinc finger genes at different stages of cotton fiber development, an indication that the C2H2-zinc finger genes play an important role in initiating and regulating fiber development in cotton.ConclusionsThis study provides a strong foundation for future practical genome research on C2H2-zinc finger genes in upland cotton. The expression levels of C2H2-zinc finger genes family is a pointer of their involvement in various biochemical and physiological functions which are directly related to cotton fiber development during initiation and elongation stages. This work not only provides a basis for determining the nominal role of the C2H2-zinc finger genes in fiber development but also provide valuable information for characterization of potential candidate genes involved in regulation of cotton fiber development.

Evolutionary distances corrected for purifying selection and ancestral polymorphisms

Evolutionary distance formulas that take into account effects due to ancestral polymorphisms and purifying selection are obtained on the basis of the full solution of Jukes–Cantor and Kimura DNA substitution models. In the case of purifying selection two different methods are developed. It is shown that avoiding the dimensional reduction implicitly carried out in the conventional model solving is instrumental to incorporate the quoted effects into the formalism. The problem of estimating the numerical values of the model parameters, as well as those of the correction terms, is not addressed.

Many human RNA viruses show extraordinarily stringent selective constraints on protein evolution

How negative selection, positive selection, and population size contribute to the large variation in nucleotide substitution rates among RNA viruses remains unclear. Here, we studied the ratios of nonsynonymous-to-synonymous substitution rates (dN/dS) in protein-coding genes of human RNA and DNA viruses and mammals. Among the 21 RNA viruses studied, 18 showed a genome-average dN/dS from 0.01 to 0.10, indicating that over 90% of nonsynonymous mutations are eliminated by negative selection. Only HIV-1 showed a dN/dS (0.31) higher than that (0.22) in mammalian genes. By comparing the dN/dS values among genes in the same genome and among species or strains, we found that both positive selection and population size play significant roles in the dN/dS variation among genes and species. Indeed, even in flaviviruses and picornaviruses, which showed the lowest ratios among the 21 species studied, positive selection appears to have contributed significantly to dN/dS We found the view that positive selection occurs much more frequently in influenza A subtype H3N2 than subtype H1N1 holds only for the hemagglutinin and neuraminidase genes, but not for other genes. Moreover, we found no support for the view that vector-borne RNA viruses have lower dN/dS ratios than non-vector-borne viruses. In addition, we found a correlation between dN and dS, implying a correlation between dN and the mutation rate. Interestingly, only 2 of the 8 DNA viruses studied showed a dN/dS < 0.10, while 4 showed a dN/dS > 0.22. These observations increase our understanding of the mechanisms of RNA virus evolution.

The HPV16 Genome Is Stable in Women Who Progress to In Situ or Invasive Cervical Cancer: A Prospective Population-Based Study.

The human papillomavirus (HPV) rate of evolution is essential for cancer-preventive strategies targeting HPV. We analyzed variability over time in a prospective, population-based nested case-control study of in situ (CIS) and invasive squamous cervical cancer (SCC). Among 757,690 women who participated in cervical screening in Sweden during 1969 to 2002, we identified 94 women who had HPV16 persistence in two serial cervical screening samples (median 24 months apart, range 0.5-178 months) and later were diagnosed with CIS (n = 59), SCC (n = 32), or remained healthy (n = 3). Whole-HPV16-genome sequencing and comparison of sequences in the serial samples revealed that all women had the same HPV16 lineage, particularly lineage A, in both serial smears. Fifty-six percent of women had an identical 7,906 base pair HPV16 sequence in both samples, and no woman had more than 15 nucleotide substitutions. The median substitution rate was 0 substitutions/site/year (95% confidence interval, 0-0.00008), with no variation between quartiles of follow-up. We concluded that in most women with HPV16 persistence preceding disease, the nucleotide substitution rate was not measurable within up to 15-years follow-up. This slow rate of evolution has important implications for both HPV-based screening and HPV vaccination. SIGNIFICANCE: These findings show there is no genomic variation over time in HPV16 infections progressing to cervical cancer, which could influence risk stratification of women when screening for cervical cancer and inform HPV vaccination strategies.

Computational Analysis of Distance based Phylogenetic Tree for Azotobacter Species

Phylogenetic tree is a pictorial representation of evolutionary relationships between organisms. It is important method to analyze the biological data. Phylogenetic trees are based on two methods : Distance based and Character based. Phylogenetic tree are used comparative analysis of any organism like human Beings, Animals, Bacteria, Viruses and Fungi’s etc. In this paper we compare 12 different nucleotide sequences of Azotobacter species having linear DNA of 999 BP as maximum size using substitution model and phylogenetic model. In this study two different models name P-Distance and Jukes cantor model are used and helped in finding UPGMA or Neighbour joining method efficiency in evaluating the similarity and dissimilarity of bacterial species. This paper gives influence in reconciliation of Azotobacter species to produce phylogram with informative branch lengths. This further leads to analyze and understand various expressive characters of Azotobacter in agriculture field.

Estimation of AT and GC content distributions of nucleotide substitution rates in bacterial core genomes

BackgroundGenomic GC content varies both within and, substantially, between microbial genomes. While some of this variation can be explained by evolutionary divergence and environmental factors, a notable portion is not understood. To investigate further, we explore a non-linear mathematical model (gcMOD) of single-nucleotide polymorphism (SNP) GC content (sbGC, the GC content of substituted bases) as a function of core genome GC content (cgGC). We estimate the model’s parameters using Bayesian inference on empirical genetic data from the microbial core genomes of 35 bacterial species, each of which contains at least 10 representative strains. We utilize 716 bacterial genomes in total. We also explore some possible implications that result from the mathematical properties of gcMOD.ResultsWe find that the median GC → AT substitution rates (β) are almost always considerably higher than the corresponding AT → GC substitution rates (α) for all 35 core genomes. The distribution of β is also noticeably more concentrated (i.e. thinner) than the corresponding distribution of α for almost all species, excepting the bacteria with the most GC-rich genomes. We also demonstrate that at the singularity point of gcMOD (where α = β), the model is reduced to a linear equation. By analyzing the linear model, we show that due to the constraints on gcMOD, the mutation rates can have profound influence on both cgGC as well as sbGC. Moreover, by examining the mathematical properties of gcMOD’s inverse function, we find that change in cgGC, and hence in genomic GC content, can potentially occur quite rapidly.ConclusionsExamining the distributions of the GC → AT and AT → GC substitution rates for 35 bacterial species, we demonstrate that the former (β) are remarkably similar for all species examined. In addition, GC → AT substitution rate distributions were considerably more concentrated for all species, with the mode consistently peaking at higher rates than for AT → GC substitution rates.

A63 Quantifying the dynamics of evolutionary rates through time

The availability of evolutionary rate estimates in recent years led to the observation that they may depend on the time-scale on which they are measured. Specifically, RNA virus evolutionary rates are frequently estimated to be low towards the past and high towards the present. This time-dependent rate phenomenon (TDRP) has important implications for evolutionary studies as it could severely bias divergence time estimates. While recent studies are providing insights into the relationship between viral evolutionary rate and time, formal probabilistic models to draw inference under TDRP scenarios remain lacking. Here, we adopt epoch-modelling to develop a Bayesian model of discrete rate changes through time in an unknown evolutionary history and combine this with a log-linear parameterization of rates as a function of times in the past. We provide an implementation for nucleotide substitution rates as well as for nonsynonymous rates change in a codon substitution model. Using a foamy virus dataset for which internal node calibrations can be applied based on host-virus co-divergence, we estimate a significant decline in evolutionary rates as a function of time into the past for nucleotide substitutions as well as for non-synonymous substitutions in a codon model. We also estimate a deep evolutionary history for primate Lentiviruses by combining an HIV-1 group M node calibration and a biogeographic calibration for viruses in drill monkeys in the TDRP model. Our analyses lead to the conclusion that studies of evolutionary timescales require a reconsideration of substitution rates, in either codon and nucleotide substitution model, as a dynamic feature of molecular evolution.

Dynamic virulence-related regions of the plant pathogenic fungus Verticillium dahliae display enhanced sequence conservation.

Plant pathogens continuously evolve to evade host immune responses. During host colonization, many fungal pathogens secrete effectors to perturb such responses, but these in turn may become recognized by host immune receptors. To facilitate the evolution of effector repertoires, such as the elimination of recognized effectors, effector genes often reside in genomic regions that display increased plasticity, a phenomenon that is captured in the two‐speed genome hypothesis. The genome of the vascular wilt fungus Verticillium dahliae displays regions with extensive presence/absence polymorphisms, so‐called lineage‐specific regions, that are enriched in in planta‐induced putative effector genes. As expected, comparative genomics reveals differential degrees of sequence divergence between lineage‐specific regions and the core genome. Unanticipated, lineage‐specific regions display markedly higher sequence conservation in coding as well as noncoding regions than the core genome. We provide evidence that disqualifies horizontal transfer to explain the observed sequence conservation and conclude that sequence divergence occurs at a slower pace in lineage‐specific regions of the V. dahliae genome. We hypothesize that differences in chromatin organisation may explain lower nucleotide substitution rates in the plastic, lineage‐specific regions of V. dahliae.

Molecular evolution of mammalian genes with epistatic interactions in fertilization

BackgroundGenes that encode proteins associated with sperm competition, fertilization, and sexual conflicts of interest are often among the most rapidly evolving parts of animal genomes. One family of sperm-expressed genes (Zp3r, C4bpa) in the mammalian gene cluster called the regulator of complement activation (RCA) encodes proteins that bind eggs and mediate reproductive success, and are therefore expected to show high relative rates of nonsynonymous nucleotide substitution in response to sexual selection in comparison to other genes not involved in gamete binding at fertilization. We tested that working hypothesis by using phylogenetic models of codon evolution to identify episodes of diversifying positive selection. We used a comparative approach to quantify the evidence for episodic diversifying selection acting on RCA genes with known functions in fertilization (and sensitivity to sexual selection), and contrast them with other RCA genes in the same gene family that function in innate immunity (and are not sensitive to sexual selection).ResultsWe expected but did not find evidence for more episodes of positive selection on Zp3r in Glires (the rodents and lagomorphs) or on C4BPA in Primates, in comparison to other paralogous RCA genes in the same taxon, or in comparison to the same orthologous RCA gene in the other taxon. That result was not unique to RCA genes: we also found little evidence for more episodes of diversifying selection on genes that encode selective sperm-binding molecules in the egg coat or zona pellucida (Zp2, Zp3) in comparison to members of the same gene family that encode structural elements of the egg coat (Zp1, Zp4). Similarly, we found little evidence for episodic diversifying selection acting on two other recently discovered genes (Juno, Izumo1) that encode essential molecules for sperm–egg fusion.ConclusionsThese negative results help to illustrate the importance of a comparative context for this type of codon model analysis. The results may also point to other phylogenetic contexts in which the effects of selection acting on these fertilization proteins might be more readily discovered and documented in mammals and other taxa.

Phylogeographic and genetic characterization of porcine circovirus type 2 in Taiwan from 2001\u20132017

Porcine circovirus type 2 (PCV2) is an important pathogen that causes significant economic losses in the swine industry worldwide. Five major PCV2 genotypes have been identified, including PCV2a, PCV2b, PCV2c, PCV2d, and PCV2e. To investigate the prevalence and phylodynamics of the different PCV2 genotypes in Taiwan, 214 PCV2 ORF2 sequences from Taiwan and other countries were analyzed. Genotypic differences were observed among PCV2a, 2b, and 2d at amino acid position 89 in ORF2, with isoleucine (I), arginine (R), and leucine (L), respectively. Similar to other countries, a genotypic shift was also observed in Taiwan, where the predominant genotype shifted from PCV2b to 2d after 2010. The estimated nucleotide substitution rate of Taiwanese strains in the ORF2 region was 8.467 × 10−4 substitutions per site per year. This rapid evolution rate of PCV2 may lead to the genotypic shift observed in Taiwan. The times to the most recent common ancestor (TMRCA) for PCV2a, -2b, and -2d-2 was dated to 1970, 1992 and 2004, respectively. Thus, the PCV2a, -2b, and -2d genotypes were already present in Taiwan before the introduction of the PCV2 vaccine.

Comparative mitogenomics of the Decapoda reveals evolutionary heterogeneity in architecture and composition

The emergence of cost-effective and rapid sequencing approaches has resulted in an exponential rise in the number of mitogenomes on public databases in recent years, providing greater opportunity for undertaking large-scale comparative genomic and systematic research. Nonetheless, current datasets predominately come from small and disconnected studies on a limited number of related species, introducing sampling biases and impeding research of broad taxonomic relevance. This study contributes 21 crustacean mitogenomes from several under-represented decapod infraorders including Polychelida and Stenopodidea, which are used in combination with 225 mitogenomes available on NCBI to investigate decapod mitogenome diversity and phylogeny. An overview of mitochondrial gene orders (MGOs) reveals a high level of genomic variability within the Decapoda, with a large number of MGOs deviating from the ancestral arthropod ground pattern and unevenly distributed among infraorders. Despite the substantial morphological and ecological variation among decapods, there was limited evidence for correlations between gene rearrangement events and species ecology or lineage specific nucleotide substitution rates. Within a phylogenetic context, predicted scenarios of rearrangements show some MGOs to be informative synapomorphies for some taxonomic groups providing strong independent support for phylogenetic relationships. Additional comparisons for a range of mitogenomic features including nucleotide composition, strand asymmetry, unassigned regions and codon usage indicate several clade-specific trends that are of evolutionary and ecological interest.

Application of a sequence-based taxonomic classification method to uncultured and unclassified marine single-stranded RNA viruses in the order Picornavirales.

Metagenomics has altered our understanding of microbial diversity and ecology. This includes its applications to viruses in marine environments that have demonstrated their enormous diversity. Within these are RNA viruses, many of which share genetic features with members of the order Picornavirales; yet, very few of these have been taxonomically classified. The only recognized family of marine RNA viruses is the Marnaviridae, which was founded based on discovery and characterization of the species Heterosigma akashiwo RNA virus. Two additional genera of marine RNA viruses, Labyrnavirus (one species) and Bacillarnavirus (three species), were subsequently defined within the order Picornavirales but not assigned to a family. We have defined a sequence-based framework for taxonomic classification of twenty marine RNA viruses into the family Marnaviridae. Using RNA-dependent RNA polymerase (RdRp) phylogeny and distance-based analyses, we assigned the genera Labyrnavirus and Bacillarnavirus to the family Marnaviridae and created four additional genera in the family: Locarnavirus (four species), Kusarnavirus (one species), Salisharnavirus (four species) and Sogarnavirus (six species). We used pairwise capsid protein comparisons to delineate species within families, with 75 per cent identity as the species demarcation threshold. The family displays high sequence diversities and Jukes–Cantor distances for both the RdRp and capsid genes, suggesting that the classified viruses are not representative of all of the virus diversity within the family and that there are many more extant taxa. Our proposed taxonomic framework provides a sound classification system for this group of viruses that will have broadly applicable principles for other viral groups. It is based on sequence data alone and provides a robust taxonomic framework to include viruses discovered via metagenomic studies, thereby greatly expanding the realm of viruses subject to taxonomic classification.

Do angiosperms with highly divergent mitochondrial genomes have altered mitochondrial function?

Angiosperm mitochondrial (mt) genes are generally slow-evolving, but multiple lineages have undergone dramatic accelerations in rates of nucleotide substitution and extreme changes in mt genome structure. While molecular evolution in these lineages has been investigated, very little is known about their mt function. Some studies have suggested altered respiration in individual taxa, although there are several reasons why mt variation might be neutral in others. Here, we develop a new protocol to characterize respiration in isolated plant mitochondria and apply it to species of Silene with mt genomes that are rapidly evolving, highly fragmented, and exceptionally large (~11 Mbp). This protocol, complemented with traditional measures of plant fitness, cytochrome c oxidase activity assays, and fluorescence microscopy, was also used to characterize inter- and intraspecific variation in mt function. Contributions of the individual “classic” OXPHOS complexes, the alternative oxidase, and external NADH dehydrogenases to overall mt respiratory flux were found to be similar to previously studied angiosperms with more typical mt genomes. Some differences in mt function could be explained by inter- and intraspecific variation. This study suggests that Silene species with peculiar mt genomes still show relatively normal mt respiration. This may be due to strong purifying selection on mt variants, coevolutionary responses in the nucleus, or a combination of both. Future experiments should explore such questions using a comparative framework and investigating other lineages with unusual mitogenomes.

Highly accelerated rates of genomic rearrangements and nucleotide substitutions in plastid genomes of Passiflora subgenus Decaloba

Plastid genomes (plastomes) of photosynthetic angiosperms are for the most part highly conserved in their organization, mode of inheritance and rates of nucleotide substitution. A small number of distantly related lineages share a syndrome of features that deviate from this general pattern, including extensive genomic rearrangements, accelerated rates of nucleotide substitution, biparental inheritance and plastome-genome incompatibility. Previous studies of plastomes in Passiflora with limited taxon sampling suggested that the genus exhibits this syndrome. To examine this phenomenon further, 15 new plastomes from Passiflora were sequenced and combined with previously published data to examine the phylogenetic relationships, genome organization and evolutionary rates across all five subgenera and the sister genus Adenia. Phylogenomic analyses using 68 protein-coding genes shared by Passiflora generated a fully resolved and strongly supported tree that is congruent with previous phylogenies based on a few plastid and nuclear loci. This phylogeny was used to examine the distribution of plastome rearrangements across Passiflora. Multiple gene and intron losses and inversions were identified in Passiflora with some occurring in parallel and others that extended across the Passifloraceae. Furthermore, extensive expansions and contractions of the inverted repeat (IR) were uncovered and in some cases this resulted in exclusion of all ribosomal RNA genes from the IR. The most highly rearranged lineage was subgenus Decaloba, which experienced extensive IR expansion that incorporated up to 25 protein-coding genes usually located in large single copy region. Nucleotide substitution rate analyses of 68 protein-coding genes across the genus showed lineage- and locus-specific acceleration. Significant increase in dS, dN and dN/dS was detected for clpP across the genus and for ycf4 in certain lineages. Significant increases in dN and dN/dS for ribosomal subunits and plastid-encoded RNA polymerase genes were detected in the branch leading to the expanded IR-clade in subgenus Decaloba. This subgenus displays the syndrome of unusual features, making it an ideal system to investigate the dynamic evolution of angiosperm plastomes.

Number of Microstates of Dark DNAs in Extra Dimensions for Normal and Cancerous Cells

Recently, Hargreaves (New Scientist, Volume 237, Issue 3168, March 2018, Pages 29-31) has argued that some animal genomes seem to be missing certain genes, ones that appear in other similar species and must be present to keep the animals alive. He called these missing genes as “dark DNA”. On the other hand, Sepehri and his collaborations (Open Physics, 16(1), pp. 463-475) has discussed that some biological events like DNA teleportation and water memory may be due to existence of some extra genes in extra dimensions. From these results, we can conclude that origin of some cancers may be evolutions of dark DNA in extra dimension. To show this, we propose a model for calculating number of microstates of a DNA for a chick embryo in extra dimension and compare with experimental data. We show that number of microstates of a DNA in four dimensions for a normal chick embryo is more than number of microstates for a cancerous chick embryo. This is because that some extra microstates are transformed to four dimensions.

Modeling of Mutational Events in the Evolution of Viruses.

Diverse studies of viral evolution have led to the recognition that the evolutionary rates of viral taxa observed are dependent on the time scale being investigated—with short-term studies giving fast substitution rates, and orders of magnitude lower rates for deep calibrations. Although each of these factors may contribute to this time dependent rate phenomenon, a more fundamental cause should be considered. We sought to test computationally whether the basic phenomena of virus evolution (mutation, replication, and selection) can explain the relationships between the evolutionary and phylogenetic distances. We tested, by computational inference, the hypothesis that the phylogenetic distances between the pairs of sequences are functions of the evolutionary path lengths between them. A Basic simulation revealed that the relationship between simulated genetic and mutational distances is non-linear, and can be consistent with different rates of nucleotide substitution at different depths of branches in phylogenetic trees.