Evolutionary genetics of pulmonary anatomical adaptations in deep-diving cetaceans.

BACKGROUND: Cetaceans, having experienced prolonged adaptation to aquatic environments, have undergone evolutionary changes in their respiratory systems. This process of evolution has resulted in the emergence of distinctive phenotypic traits, notably the abundance of elastic fibers and thickened alveolar walls in their lungs, which may facilitate alveolar collapse during diving. This structure helps selective exchange of oxygen and carbon dioxide, while minimizing nitrogen exchange, thereby reducing the risk of DCS. Nevertheless, the scientific inquiry into the mechanisms through which these unique phenotypic characteristics govern the diving behavior of marine mammals, including cetaceans, remains unresolved. RESULTS: This study entails an evolutionary analysis of 42 genes associated with pulmonary fibrosis across 45 mammalian species. Twenty-one genes in cetaceans exhibited accelerated evolution, featuring specific amino acid substitutions in 14 of them. Primarily linked to the development of the respiratory system and lung morphological construction, these genes play a crucial role. Moreover, among marine mammals, we identified eight genes undergoing positive selection, and the evolutionary rates of three genes significantly correlated with diving depth. Specifically, the SFTPC gene exhibited convergent amino acid substitutions. Through in vitro cellular experiments, we illustrated that convergent amino acid site mutations in SFTPC contribute positively to pulmonary fibrosis in marine mammals, and the presence of this phenotype can induce deep alveolar collapse during diving, thereby reducing the risk of DCS during diving. CONCLUSIONS: The study unveils pivotal genetic signals in cetaceans and other marine mammals, arising through evolution. These genetic signals may influence lung characteristics in marine mammals and have been linked to a reduced risk of developing DCS. Moreover, the research serves as a valuable reference for delving deeper into human diving physiology.

Can irisin be developed as the molecular evolutionary clock based on the origin and functions?

Irisin, a myokine identified in 2012, has garnered research interest for its capacity to induce browning of adipocytes and improve metabolic parameters. As such, the potential therapeutic applications of this exercise-induced peptide continue to be explored. Though present across diverse animal species, sequence analysis has revealed subtle variation in the irisin protein. In this review, we consider the effects of irisin on disease states in light of its molecular evolution. We summarize current evidence for irisin's influence on pathologies and discuss how sequence changes may inform development of irisin-based therapies. Furthermore, we propose that the phylogenetic variations in irisin could potentially be leveraged as a molecular clock to elucidate evolutionary relationships.

Genomic Characterization and Molecular Evolution of Sapovirus in Children under 5 Years of Age.

Sapovirus (SaV) is a type of gastroenteric virus that can cause acute gastroenteritis. It is highly contagious, particularly among children under the age of 5. In this study, a total of 712 stool samples from children under the age of 5 with acute gastroenteritis were collected. Out of these samples, 28 tested positive for SaV, resulting in a detection rate of 3.93% (28/712). Samples with Ct < 30 were collected for library construction and high-throughput sequencing, resulting in the acquisition of nine complete genomes. According to Blast, eight of them were identified as GI.1, while the remaining one was GI.6. The GI.6 strain sequence reported in our study represents the first submission of the GI.6 strain complete genome sequence from mainland China to the Genbank database, thus filling the data gap in our country. Sequence identity analysis revealed significant nucleotide variations between the two genotypes of SaV and their corresponding prototype strains. Phylogenetic and genetic evolution analyses showed no evidence of recombination events in the obtained sequences. Population dynamics analysis demonstrated potential competitive inhibition between two lineages of GI.1. Our study provides insights into the molecular epidemiological and genetic evolution characteristics of SaV prevalent in the Nantong region of China, laying the foundation for disease prevention and control, as well as pathogen tracing related to SaV in this area.

Molecular evolutionary analyses of the fusion protein gene in human respirovirus 1.

Few evolutionary studies of the human respiratory virus (HRV) have been conducted, but most of them have focused on HRV3. In this study, the full-length fusion (F) genes in HRV1 strains collected from various countries were subjected to time-scaled phylogenetic, genome population size, and selective pressure analyses. Antigenicity analysis was performed on the F protein. The time-scaled phylogenetic tree using the Bayesian Markov Chain Monte Carlo method estimated that the common ancestor of the HRV1 F gene diverged in 1957 and eventually formed three lineages. Phylodynamic analyses showed that the genome population size of the F gene has doubled over approximately 80 years. Phylogenetic distances between the strains were short (< 0.02). No positive selection sites were detected for the F protein, whereas many negative selection sites were identified. Almost all conformational epitopes of the F protein, except one in each monomer, did not correspond to the neutralising antibody (NT-Ab) binding sites. These results suggest that the HRV1 F gene has constantly evolved over many years, infecting humans, while the gene may be relatively conserved. Mismatches between computationally predicted epitopes and NT-Ab binding sites may be partially responsible for HRV1 reinfection and other viruses such as HRV3 and respiratory syncytial virus.

Evaluation of codon usage patterns and molecular evolution dynamics in Japanese encephalitis virus: An integrated bioinformatics approach.

In the recent survey, Japanese encephalitis (JE) is one of the most common mosquito-borne diseases, accounting for ∼30% of fatalities. The outbreaks of the JE virus (JEV) suggests that exhaustive study is essential for the prevention and management of the disease. The disease mainly spreads in humans and pigs by the vector: mosquito; as this is a major concern, this study had employed various bioinformatics tools to investigate the codon usage bias, evolutionary inference and selection pressure analysis of the Japanese encephalitis virus disease. The results indicated that the JE virus was biased and natural selection was the main factor shaping the codon usage that was determined and confirmed with the Nc, neutrality, PR2 plots and correlation analysis. The evolutionary analysis revealed that the virus had a substitution rate of 1.54 × 10-4 substitution/site/year and the tMRCA was found to be in 1723. The transmission of the virus in the map found transmissions mostly from China and transmitted across Asia and Africa. The selection pressure analysis employed three methods which had 969th codon site as diversifying site and had many purifying sites that shows the virus had evolved rapidly. The inferences from this study would aid people to employ this methodology on various diseases and also perform insilico studies in the field of vaccinology and immunoinformatics.

Molecular Evolutionary Growth of Ultralong Semiconducting Double-Walled Carbon Nanotubes.

The self-assembling preparation accompanied with template auto-catalysis loop and the ability to gather energy, induces the appearance of chirality and entropy reduction in biotic systems. However, an abiotic system with biotic characteristics is of great significance but still missing. Here, it is demonstrated that the molecular evolution is characteristic of ultralong carbon nanotube preparation, revealing the advantage of chiral assembly through template auto-catalysis growth, stepwise-enriched chirality distribution with decreasing entropy, and environmental effects on the evolutionary growth. Specifically, the defective and metallic nanotubes perform inferiority to semiconducting counterparts, among of which the ones with double walls and specific chirality (n, m) are more predominant due to molecular coevolution. An explicit evolutionary trend for tailoring certain layer chirality is presented toward perfect near-(2n, n)-containing semiconducting double-walled nanotubes. These findings extend our conceptual understanding for the template auto-catalysis assembly of abiotic carbon nanotubes, and provide an inspiration for preparing chiral materials with kinetic stability by evolutionary growth.

Molecular Evolutionary Rate Predicts Intraspecific Genetic Polymorphism and Species-Specific Selection.

It is unknown what determines genetic diversity and how genetic diversity is associated with various biological traits. In this work, we provide insight into these issues. By comparing genetic variation of 14,671 mammalian gene trees with thousands of individual human, chimpanzee, gorilla, mouse, and dog/wolf genomes, we found that intraspecific genetic diversity can be predicted by long-term molecular evolutionary rates rather than de novo mutation rates. This relationship was established during the early stage of mammalian evolution. Moreover, we developed a method to detect fluctuations of species-specific selection on genes based on the deviations of intraspecific genetic diversity predicted from long-term rates. We showed that the evolution of epithelial cells, rather than connective tissue, mainly contributed to morphological evolution of different species. For humans, evolution of the immune system and selective sweeps caused by infectious diseases are the most representative examples of adaptive evolution.

Epidemiological characteristics and molecular evolution analysis of HIV-1 CRF07_BC in Nanjing from 2015 to 2019 / 中华微生物学和免疫学杂志

Objective:To investigate the epidemiology, evolution and transmission characteristics of HIV-1 CRF07_BC in Nanjing between 2015 and 2019 to provide scientific basis for precise prevention and control of HIV-1 transmission.Methods:Pol gene sequences of 319 patients with HIV-1 CRF07_BC infection in Nanjing were amplified and sequenced and a maximum likelihood phylogenetic tree was then constructed. Markov Chain Monte Carlo sampling method was used to generate Maximum Clade Credibility Tree. Bayesian Skyline Plot was used to reconstruct the changing trend of the effective population size. Pairwise gene distance method was used to construct molecular network to investigate the transmission features. Results:Among the 319 patients, 303 (95.0%) were male; 264 (82.8%) had multiple sexual partners; only 14 (4.4%) had been using condoms. Most of the infections were occurred in men who have sex with men (MSM), accounting for 77.4%. The maximum likelihood phylogenetic analysis of HIV-1 CRF07_BC revealed two clusters: Cluster1 and Cluster2. Cluster1 mainly contained the strains isolated from MSM and Cluster2 mainly consisted of the strains isolated from heterosexual patients. The most recent common ancestor was 2002.47(1999.91, 2005.43) year for Cluster1 and 1996.38(1992.55, 1999.76) year for Cluster2. The evolutionary rates (95% highest posterior density, 95%HPD) of Cluster1 and Cluster2 were 1.73×10 -3 (1.36×10 -3-2.16×10 -3) substitutions·site -1·year -1 and 2.09×10 -3 (1.50×10 -3 -2.79×10 -3 ) substitutions·site -1·year -1, respectively. The effective population sizes of Cluster1 and Cluster2 tended to be stable after 2002 and 2003, respectively. In addition, Cluster1 and Cluster2 formed eleven and eight unique branches, respectively, suggesting the possibility of divergent epidemics of this genotype. A total of 35 propagation clusters were formed in the molecular propagation network, including 92 Nanjing sequences with an average degree of 4.3. Males, MSM and people with multiple sexual partners were more likely to be connected to the network. Students and young patients were more likely to be connected to the network. Conclusions:HIV-1 CRF07_BC infection was characterized by low age, multiple sexual partners, unprotected behaviors and transmission among MSM in Nanjing from 2015 to 2019. It was recommended to pay more attention to students and young people, to formulate more effective prevention and control measures for high-risk sexual behaviors, and to carry out continuous molecular monitoring of CRF07_BC strain, so as to provide scientific basis for the prevention and control of HIV CRF07_BC.

Territorywide Study of Early Coronavirus Disease Outbreak, Hong Kong, China.

Initial cases of coronavirus disease in Hong Kong were imported from mainland China. A dramatic increase in case numbers was seen in February 2020. Most case-patients had no recent travel history, suggesting the presence of transmission chains in the local community. We collected demographic, clinical, and epidemiologic data from 50 patients, who accounted for 53.8% of total reported case-patients as of February 28, 2020. We performed whole-genome sequencing to determine phylogenetic relationship and transmission dynamics of severe acute respiratory syndrome coronavirus 2 infections. By using phylogenetic analysis, we attributed the community outbreak to 2 lineages; 1 harbored a common mutation, Orf3a-G251V, and accounted for 88.0% of the cases in our study. The estimated time to the most recent common ancestor of local coronavirus disease outbreak was December 24, 2019, with an evolutionary rate of 3.04 × 10-3 substitutions/site/year. The reproduction number was 1.84, indicating ongoing community spread.

The Effects of Ecological Traits on the Rate of Molecular Evolution in Ray-Finned Fishes: A Multivariable Approach.

Myriad environmental and biological traits have been investigated for their roles in influencing the rate of molecular evolution across various taxonomic groups. However, most studies have focused on a single trait, while controlling for additional factors in an informal way, generally by excluding taxa. This study utilized a dataset of cytochrome c oxidase subunit I (COI) barcode sequences from over 7000 ray-finned fish species to test the effects of 27 traits on molecular evolutionary rates. Environmental traits such as temperature were considered, as were traits associated with effective population size including body size and age at maturity. It was hypothesized that these traits would demonstrate significant correlations with substitution rate in a multivariable analysis due to their associations with mutation and fixation rates, respectively. A bioinformatics pipeline was developed to assemble and analyze sequence data retrieved from the Barcode of Life Data System (BOLD) and trait data obtained from FishBase. For use in phylogenetic regression analyses, a maximum likelihood tree was constructed from the COI sequence data using a multi-gene backbone constraint tree covering 71% of the species. A variable selection method that included both single- and multivariable analyses was used to identify traits that contribute to rate heterogeneity estimated from different codon positions. Our analyses revealed that molecular rates associated most significantly with latitude, body size, and habitat type. Overall, this study presents a novel and systematic approach for integrative data assembly and variable selection methodology in a phylogenetic framework.

Computational analysis of SARS-CoV-2, SARS-CoV, and MERS-CoV genome using MEGA.

The novel coronavirus pandemic that has originated from China and spread throughout the world in three months. Genome of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) predecessor, severe acute respiratory syndrome coronavirus (SARS-CoV) and Middle East respiratory syndrome coronavirus (MERS-CoV) play an important role in understanding the concept of genetic variation. In this paper, the genomic data accessed from National Center for Biotechnology Information (NCBI) through Molecular Evolutionary Genetic Analysis (MEGA) for statistical analysis. Firstly, the Bayesian information criterion (BIC) and Akaike information criterion (AICc) are used to evaluate the best substitution pattern. Secondly, the maximum likelihood method used to estimate of transition/transversions (R) through Kimura-2, Tamura-3, Hasegawa-Kishino-Yano, and Tamura-Nei nucleotide substitutions model. Thirdly and finally nucleotide frequencies computed based on genomic data of NCBI. The results indicate that general times reversible model has the lowest BIC and AICc score 347,394 and 347,287, respectively. The transition/transversions bias for nucleotide substitutions models varies from 0.56 to 0.59 in MEGA output. The average nitrogenous bases frequency of U, C, A, and G are 31.74, 19.48, 28.04, and 20.74, respectively in percentages. Overall the genomic data analysis of SARS-CoV-2, SARS-CoV, and MERS-CoV highlights the close genetic relationship.

Role of BgaA as a Pneumococcal Virulence Factor Elucidated by Molecular Evolutionary Analysis.

Streptococcus pneumoniae is a major cause of pneumonia, sepsis, and meningitis. Previously, we identified a novel virulence factor by investigating evolutionary selective pressure exerted on pneumococcal choline-binding cell surface proteins. Herein, we focus on another pneumococcal cell surface protein. Cell wall-anchoring proteins containing the LPXTG motif are conserved in Gram-positive bacteria. Our evolutionary analysis showed that among the examined genes, nanA and bgaA had high proportions of codon that were under significant negative selection. Both nanA and bgaA encode a multi-functional glycosidase that aids nutrient acquisition in a glucose-poor environment, pneumococcal adherence to host cells, and evasion from host immunity. However, several studies have shown that the role of BgaA is limited in a mouse pneumonia model, and it remains unclear if BgaA affects pneumococcal pathogenesis in a mouse sepsis model. To evaluate the distribution and pathogenicity of bgaA, we performed phylogenetic analysis and intravenous infection assay. In both Bayesian and maximum likelihood phylogenetic trees, the genetic distances between pneumococcal bgaA was small, and the cluster of pneumococcal bgaA did not contain other bacterial orthologs except for a Streptococcus gwangjuense gene. Evolutionary analysis and BgaA structure indicated BgaA active site was not allowed to change. The mouse infection assay showed that the deletion of bgaA significantly reduced host mortality. These results indicated that both nanA and bgaA encode evolutionally conserved pneumococcal virulence factors and that molecular evolutionary analysis could be a useful alternative strategy for identification of virulence factors.

Emerging Frontiers in the Study of Molecular Evolution.

A collection of the editors of Journal of Molecular Evolution have gotten together to pose a set of key challenges and future directions for the field of molecular evolution. Topics include challenges and new directions in prebiotic chemistry and the RNA world, reconstruction of early cellular genomes and proteins, macromolecular and functional evolution, evolutionary cell biology, genome evolution, molecular evolutionary ecology, viral phylodynamics, theoretical population genomics, somatic cell molecular evolution, and directed evolution. While our effort is not meant to be exhaustive, it reflects research questions and problems in the field of molecular evolution that are exciting to our editors.

Hyper-production of taxol from Aspergillus fumigatus, an endophytic fungus isolated from Taxus sp. of the Northern Himalayan region.

Taxol® (generic name Paclitaxel) is a chemotherapeutic drug, effective against head, neck, breast, lung, bladder, ovary, and cervix cancers. Rising demands in chemotherapy and limited supply of natural taxol have ultimately increased the cost of the drug. Semi synthesis using taxol precursors is not able to meet the global supply and has intensified the need to find alternative ways of taxol production. In the present study, 34 different endophytes were isolated from Taxus sp. collected from Shimla, Himachal Pradesh (India). Primary screening of taxol-producing fungi was carried out based on the presence of dbat gene, essential for the taxol biosynthetic pathway. A fungal isolate TPF-06 was screened to be a taxol-producing strain based on the PCR amplification results. It was characterized and identified as Aspergillus fumigatus by 18S rRNA (Accession No. KU-837249). Multiple sequence alignment (MSA) of nuclear ribosomal internal transcribed spacer (ITS) region and phylogenetic analysis confirmed that strain belonged to A. fumigatus clade (Accession No. MF-374798) and is endophytic in nature. Presence of taxol was detected and quantified by High-Performance Liquid Chromatography (HPLC) and characterized by using Thin Layer Chromatography (TLC), Ultraviolet (UV) spectroscopy, Mass spectrometry (MS), Fourier-Transform Infrared Spectroscopy (FTIR) and Nuclear Magnetic Resonance (NMR) spectroscopy. Microbial fermentation in the S7 medium yielded 1.60 g/L of taxol, which to the best of our knowledge is the highest taxol production from an endophytic fungus. Findings of the present study suggest that the A. fumigatus is an excellent alternate source for taxol supply, and it may become a highly potent strain on a commercial scale. The involvement of dbat gene in A. fumigatus KU-837249 strain further suggested a way of increasing taxol yield in fungi by medium engineering and recombinant DNA technology in the future.

Genetic Analysis of Avian Coronavirus Infectious Bronchitis Virus in Yellow Chickens in Southern China over the Past Decade: Revealing the Changes of Genetic Diversity, Dominant Genotypes, and Selection Pressure.

The high mutation rates of infectious bronchitis virus (IBV) pose economic threats to the poultry industry. In order to track the genetic evolutionary of IBV isolates circulating in yellow chickens, we continued to conduct the genetic analyses of the structural genes S1, E, M, and N from 64 IBV isolates in southern China during 2009-2017. The results showed that the dominant genotypes based on the four genes had changed when compared with those during 1985-2008. Based on the S1 gene phylogenetic tree, LX4-type (GI-19) was the most dominant genotype, which was different from that during 1985-2008. The second most dominant genotype was LDT3-A-type, but this genotype disappeared after 2012. New-type 1 (GVI-1) isolates showed increasing tendency and there were four aa (QKEP) located in the hypervariable region (HVR) III and one aa (S) insertion in all the New-type 1 isolates. Both the analyses of amino acid entropy and molecular evolutionary rate revealed that the variations from large to small were S1, E, M, and N. Purifying selection was detected in the S1, E, M, and N gene proteins, which was different from the positive selection during 1985-2008. Six isolates were confirmed to be recombinants, possibly generated from a vaccine virus of the 4/91-type or LDT3-A-type and a circulating virus. The estimated times for the most recent common ancestors based on the S1, E, M, and N genes were the years of 1744, 1893, 1940, and 1945, respectively. Bayesian skyline analysis revealed a sharp decrease in genetic diversity of all the four structural genes after 2010 and since late 2015, the viral population rapidly rose. In conclusion, the IBVs circulating in southern China over the past decade have experienced a remarkable change in genetic diversity, dominant genotypes, and selection pressure, indicating the importance of permanent monitoring of circulating strains and the urgency for developing new vaccines to counteract the emerging LX4-type and New-type IBVs.

Epidemiological and molecular analysis of avian influenza A(H7N9) virus in Shanghai, China, 2013-2017.

BACKGROUND: Human infections with a novel avian influenza A virus (H7N9) were reported in Shanghai municipality, China, at the beginning of 2013. High-pathogenic avian influenza (HPAI) H7N9 virus emerged in late February 2017 along with existing low-pathogenic avian influenza (LPAI) H7N9 virus, and this has the potential to develop into a pandemic that could be harmful to humans. METHODS: To elucidate the epidemiological characteristics of H7N9-infected cases from 2013 to 2017 in Shanghai, data on the 59 laboratory-confirmed human cases and 26 bird and environmental contamination cases were collected from the WHO website and Food and Agriculture Organization Emergency Prevention System for Animal Health (FAO EMPRES-AH). Full-length sequences of H7N9 viruses that emerged in Shanghai were collected from the Global Initiative on Sharing Avian Influenza Data to analyze the evolutionary and genetic features. RESULTS: We found that genetically different strains emerged in every epidemic in Shanghai, and most of the circulating H7N9 strains had affinity to human-type receptors, with the characteristics of high-virulence and low-pathogenic influenza viruses. Furthermore, our findings suggest that the Shanghai chicken strains are closely related to the HPAI H7N9 virus A/Guangdong/17SF003/2016, indicating that this viral strain is of avian origin and generated from the LPAI H7N9 viruses in Shanghai. The gradual decrease in H7N9 human infection in Shanghai was probably due to the control measures taken by the Shanghai government and the enhanced public awareness leading to a reduced risk of H7N9 virus infection. However, LPAI H7N9 viruses from poultry and environmental samples were continually detected in Shanghai across the epidemics, increasing the risk of new emerging H7N9 outbreaks. CONCLUSION: It is important to consistently obtain sufficient surveillance data and implement prevention measures against H7N9 viruses in Shanghai municipality.

Evolutional conservation of molecular structure and antiviral function of a type I interferon, IFN-kappa, in poultry.

IFN-kappa (IFN-κ) is a type I IFN expressed by keratinocytes, monocytes and dendritic cells with important roles during the innate immune response period. This research was conducted to elaborate the evolution and characteristics of IFN-κ in poultry. Chicken IFN-κ is located on the sex-determining Z chromosome, which is greatly different from mammals. Poultry IFN-κ cluster together in a species-specific manner through positive selection pressure and share only 19-33% homology with mammalian IFN-κ and poultry other type I IFN. Both chicken and duck IFN-κ was constitutively expressed in spleen, skin, lung, and peripheral blood mononuclear cells (PBMC), as well as being significantly induced after treatment with virus in PBMC. Biologically, poultry IFN-κ has antiviral activity against VSV in chicken embryonic fibroblasts and duck embryonic fibroblasts (CEF and DEF) cells, and induces the expression of IFN stimulated genes (ISGs). After treatment with JAK1 inhibitor, the ISGs expression can be down-regulated. Overall, our research on poultry IFN-κ not only enriches the knowledge about IFN-κ but also facilitates further research on the role of type I IFNs in antiviral defense responses in poultry.

Functional redundancy in Echinocandin B in-cluster transcription factor ecdB of Emericella rugulosa NRRL 11440.

Echinocandin B is a potent antifungal against the majority of fungal pathogens and its biosynthesis occurred by ecd and hty gene clusters in Emericella rugulosa NRRL 11440. We elucidated the functional necessity of in-clustered transcription factor; ecdB in the production of echinocandin B. We deleted the ecdB gene and found that ΔecdB mutant has no significant effect on echinocandin B production. The expression level of most of the ecd and hty cluster genes was not significantly altered except few of them up-regulated in knockout strain. The complete abrogation in ecdB gene expression was observed in ΔecdB strain. However, the interactions of purified EcdB protein with DNA sequence of ecdA, ecdH, ecdK and ecdI promoter was confirmed in-vitro. Our results conclude that EcdB protein in-vitro binds to the ecdA, ecdH, ecdK and ecdI promoter but in-vivo, it could not significantly affect the gene expression and echinocandin B production in Emericella rugulosa.

Phylogenetic estimates of diversification rate are affected by molecular rate variation.

Molecular phylogenies are increasingly being used to investigate the patterns and mechanisms of macroevolution. In particular, node heights in a phylogeny can be used to detect changes in rates of diversification over time. Such analyses rest on the assumption that node heights in a phylogeny represent the timing of diversification events, which in turn rests on the assumption that evolutionary time can be accurately predicted from DNA sequence divergence. But there are many influences on the rate of molecular evolution, which might also influence node heights in molecular phylogenies, and thus affect estimates of diversification rate. In particular, a growing number of studies have revealed an association between the net diversification rate estimated from phylogenies and the rate of molecular evolution. Such an association might, by influencing the relative position of node heights, systematically bias estimates of diversification time. We simulated the evolution of DNA sequences under several scenarios where rates of diversification and molecular evolution vary through time, including models where diversification and molecular evolutionary rates are linked. We show that commonly used methods, including metric-based, likelihood and Bayesian approaches, can have a low power to identify changes in diversification rate when molecular substitution rates vary. Furthermore, the association between the rates of speciation and molecular evolution rate can cause the signature of a slowdown or speedup in speciation rates to be lost or misidentified. These results suggest that the multiple sources of variation in molecular evolutionary rates need to be considered when inferring macroevolutionary processes from phylogenies.

Genome sequence characterization of pigeon circoviruses in China.

Pigeon circovirus (PiCV) was detected by PCR in pigeons from China. Altogether, 48 out of 244 pigeons tested positive for PiCV (positive rate, 19.67%), suggesting that the virus was prevalent in China. From the 48 PiCV-positive samples, about 2040bp complete genome fragments were obtained by full length genome amplification and sequenced with a next-generation sequencing platform. Characteristics of the ORFs from different PiCV strains tested in this study were analyzed. Several insertion, deletion or substitutions were discovered during the analysis of the nucleotide sequence compared with sequences reported previously. In phylogenetic tree analysis, 48 sequences isolated in this study could be further divided into five clades (A, B, C, D, and F), clade E includes reference sequences only. Two major groups were found in the six clades, distinguished by ATA and ATG initiation codons. Most of the viruses isolated in the study were in the ATG group, with fewer in the ATA branch.