Bayesian Evolutionary Analysis Research Articles

Molecular phylogenetic and estimation of evolutionary divergence and biogeography of the family Schizoparmaceae and allied families (Diaporthales, Ascomycota)

The Diaporthales includes 32 families, many of which are important plant pathogens, endophytes and saprobes, e.g., members of the families Pseudoplagiostomataceae, Pyrisporaceae and Schizoparmaceae. Nucleotide sequences derived from five genetic loci including: ITS, LSU, TEF1-α, TUB2 and RPB2 were used for Bayesian evolutionary analysis to determine divergence times and evolutionary relationships within the Schizoparmaceae. Molecular clock analyses revealed that the ancestor of Schizoparmaceae split during the Upper Cretaceous period approximately 75.7 Mya (95 % highest posterior density of 60.3–91.3 Mya). Reconstructing ancestral state in phylogenies (RASP) with using the Bayesian Binary Markov chain Monte Carlo (BBM) Method to reconstruct the historical biogeography for the family Schizoparmaceae indicated its most likely origin in Africa. Based on taxonomic and phylogenetic analyses, the Pseudoplagiostomataceae and Pyrisporaceae relationships were clarified and a total of four species described herein. For Pseudoplagiostomataceae, three new species and one known species that include, Pseudoplagiostoma fafuense sp. nov., Ps. ilicis sp. nov., Ps. sanmingense sp. nov. and Ps. bambusae are described and a key of Pseudoplagiostomataceae is provided. With respect to Pyrisporaceae, we considered Pseudoplagiostoma castaneae to be a synonym of Pyrispora castaneae. In addition, a new species of Schizoparmaceae, Coniella fujianensis sp. nov. is described and illustrated.

Quantifying the effects of the mitochondrial genome on milk production traits in dairy cows: empirical results and modelling challenges

Significant advances in livestock traits have been achieved primarily through selection strategies targeting variation in the nuclear genome, with little attention given to mitogenome variation. We analyzed the influence of the mitogenome on milk production traits of Holstein cattle in Croatia based on strategically generated next-generation sequencing data for 109 cows pedigree-linked to 7115 milk production records (milk, fat and protein yield) from 3006 cows (first 5 lactations). Since little is known about the biology of the relationship between mitogenome variation and production traits, our quantitative genetic modeling was complex. Thus, the proportion of total variance explained by mitogenome inheritance was estimated using 5 different models: (1) cytoplasmic model with maternal lineages (CYTO), (2) haplotypic model with mitogenome sequences (HAPLO), (3) amino acid model with unique amino acid sequences (AMINO), (4) evolutionary model based on a phylogenetic analysis using Bayesian Evolutionary Analysis Sampling Trees phylogenetic analysis (EVOL), and (5) mitogenome SNP model (SNPmt). The polygenic autosomal and X chromosome additive genetic effects based on pedigree were modeled, together with the effects of herd-year-season interaction, permanent environment, location, and age at first calving. The estimated proportions of phenotypic variance explained by mitogenome in 4 different models (CYTO, HAPLO, AMINO, and SNPmt) were found to be substantial given the size of mitogenome, ranging from 5% to 7% for all 3 milk traits. At the same time, a negligible proportion of the phenotypic variance was explained by mitogenome with the EVOL model. Similarly, in all models, no proportion of phenotypic variance was explained by the X chromosome. Although our results should be confirmed in other dairy cattle populations, including a large number of sequenced mitogenomes and nuclear genomes, the potential of utilizing mitogenome information in animal breeding is promising, especially as the acquisition of complete genome sequences becomes cost-effective.

Genomic surveillance reveals low-level circulation of two subtypes of genogroup C coxsackievirus A10 in Nanchang, Jiangxi Province, China, 2015-2023.

In recent years, coxsackievirus (CV) A10 has been associated with increasing sporadic hand, foot, and mouth disease (HFMD) cases and outbreaks globally. In addition to mild symptoms such as pharyngitis and herpangina, CVA10-related complications or even fatality can occur. Currently, systematic phylogenetic studies of CVA10 are limited. In this study, we first explored the epidemiological and genetic characteristics of CVA10 in Nanchang, an inland southeastern city of China, based on the HFMD surveillance network from 2015-2023. Among 3429 enterovirus-positive cases, 110 (3.04%) were associated with CVA10, with a male-to-female ratio of 1.62. The median age of the CVA10 patients was 2.3 years (interquartile range, IQR 1.0-4.0), with 94.55% (104/110) of the patients aged less than 5 years. Phylogenetic analyses using the full-length VP1, 5'UTR, P1, P2, P3 sequences and near full-length genomes indicated that CVA10 strains (n = 57) isolated in Nanchang belonged to genogroup C; two strains identified in 2017 belonged to C1 subtypes clustered with strains from Vietnam, Madagascar, France and Spain; and the others belonged to C2 subtypes interdigitating with CVA10 isolates from mainland China, the United States and Australia. Through extensive analysis, we identified a rare F168Y mutation in epitope 4 of VP1 in a Madagascar strain of genogroup F and a Chinese strain of genogroup C. Based on Bayesian evolutionary analyses, the average nucleotide substitution rate for the VP1 gene of CV10 strains was 3.07×10-3 substitutions/site/year. The most recent common ancestor (tMRCA) of genogroup C was dated 1990.84, and the tMRCA of CVA10 strains from Nanchang was dated approximately 2003.16, similar to strains circulating in other regions of China, suggesting that the viruses were likely introduced and cryptically circulated in China before the establishment of the HFMD surveillance network. Recombination analysis indicated intertypic recombination of the Nanchang strain with the genogroup G strain in the 3D region. Given the shifting dominance of viral genotypes and frequent recombination events, the existing surveillance system needs to be regulated to enhance genomic surveillance efforts on a more diverse spectrum of genotypes in the future.

Deciphering the evolutionary landscape of severe fever with thrombocytopenia syndrome virus across East Asia.

Severe fever with thrombocytopenia syndrome virus (SFTSV) poses a significant public health challenge in East Asia, necessitating a deeper understanding of its evolutionary dynamics to effectively manage its spread and pathogenicity. This study provides a comprehensive analysis of the genetic diversity, recombination patterns, and selection pressures across the SFTSV genome, utilizing an extensive dataset of 2041 sequences from various hosts and regions up to November 2023. Employing maximum likelihood and Bayesian evolutionary analysis by sampling trees (BEAST), we elucidated the phylogenetic relationships among nine distinct SFTSV genotypes (A, B1, B2, B3, B4, C, D, E, and F), revealing intricate patterns of viral evolution and genotype distribution across China, South Korea, and Japan. Furthermore, our analysis identified 34 potential reassortments, underscoring a dynamic genetic interplay among SFTSV strains. Genetic recombination was observed most frequently in the large segment and least in the small segment, with notable recombination hotspots characterized by stem-loop hairpin structures, indicative of a structural propensity for genetic recombination. Additionally, selection pressure analysis on critical viral genes indicated a predominant trend of negative selection, with specific sites within the RNA-dependent RNA polymerase and glycoprotein genes showing positive selection. These sites suggest evolutionary adaptations to host immune responses and environmental pressures. This study sheds light on the intricate evolutionary mechanisms shaping SFTSV, offering insights into its adaptive strategies and potential implications for vaccine development and therapeutic interventions.

Molecular evolutionary analysis of novel NSP4 mono-reassortant G1P[8]-E2 rotavirus strains that caused a discontinuous epidemic in Japan in 2015 and 2018.

In the 2010s, several unusual rotavirus strains emerged, causing epidemics worldwide. This study reports a comprehensive molecular epidemiological study of rotaviruses in Japan based on full-genome analysis. From 2014 to 2019, a total of 489 rotavirus-positive stool specimens were identified, and the associated viral genomes were analyzed by next-generation sequencing. The genotype constellations of those strains were classified into nine patterns (G1P[8] (Wa), G1P[8]-E2, G1P[8] (DS-1), G2P[4] (DS-1), G3P[8] (Wa), G3P[8] (DS-1), G8P[8] (DS-1), G9P[8] (Wa), and G9P[8]-E2). The major prevalent genotype differed by year, comprising G8P[8] (DS-1) (37% of that year's isolates) in 2014, G1P[8] (DS-1) (65%) in 2015, G9P[8] (Wa) (72%) in 2016, G3P[8] (DS-1) (66%) in 2017, G1P[8]-E2 (53%) in 2018, and G9P[8] (Wa) (26%) in 2019. The G1P[8]-E2 strains (G1-P[8]-I1-R1-C1-M1-A1-N1-T1-E2-H1) isolated from a total of 42 specimens in discontinuous years (2015 and 2018), which were the newly-emerged NSP4 mono-reassortant strains. Based on the results of the Bayesian evolutionary analyses, G1P[8]-E2 and G9P[8]-E2 were hypothesized to have been generated from distinct independent inter-genogroup reassortment events. The G1 strains detected in this study were classified into multiple clusters, depending on the year of detection. A comparison of the predicted amino acid sequences of the VP7 epitopes revealed that the G1 strains detected in different years encoded VP7 epitopes harboring distinct mutations. These mutations may be responsible for immune escape and annual changes in the prevalent strains.

The maternal genetic origin and diversity of the extant populations of the Ladakh region in India

Ladakh lies at a strategic location between the Indus River valley and the Hindu Khush Mountains, which makes the “Land of high passes” one of the major routes of movement. Through the years the region has faced multi-layered cultural movements, genetic assimilation and demographic changes. The initial settlement in the years goes back to the early Neolithic age and still continues despite its harsh, unhospitable and cold climate. Previous studies mostly covered the patrilineal markers of the region and an in-depth study lacked to represent the matrilineal ancestry and possible genetic inflow in the region. Hence, our current study first time generated complete mitogenomes of 108 unrelated individuals from Ladakh belonging to three population groups namely, Changpa (n = 38), Brokpa (n = 32) and Monpa (n = 38). In the in-depth analysis, we found that the mitogenome of the three Ladakhi groups are highly diverse in terms of maternal haplogroup distribution carrying lineages specific to East Asia (M9a), Tibbet (A21) and South Asia (M3, M30, U2). In our analysis we found that Changpa and Monpa probably have shared maternal ancestry compared to Brokpa, which is very distinct and also later suffered possible historical Bottleneck. Bayesian evolutionary and Network analysis indicates more ancient maternal lineage of Changpa and Monpa in terms of M9a haplotypes, but they also share some genetic history with Tibeto-Burman speakers in past. These findings conclusively indicate possible matrilineal genetic inflow in Ladakh from three directions, primarily from East Asia or South East Asia during post-glacial, West Eurasia and also from South Asia.

Bridging Evolutionary History and Conservation of New World Vultures.

The New World Vultures (Cathartidae) include seven species of obligate scavengers that, despite their ecological relevance, present critical information gaps around their evolutionary history and conservation. Insights into their phylogenetic relationships in recent years has enabled the addressing of such information gaps through approaches based on phylogeny. We reconstructed the ancestral area in America of the current species using two regionalization schemes and methods: Biogeography with Bayesian Evolutionary Analysis (BioGeoBears) and Bayesian Binary Model-Monte Carlo Markov Chains (BBM-MCMC). Then, we identified the priority species and areas for conservation by means of the Evolutionary Distinctiveness index (ED), as a proxy of the uniqueness of species according to phylogeny, and the Global Endangerment index (GE), mapping phylogenetic diversity. We found that the ancestral area of New World Vultures in America corresponds to South America, with dispersal processes that led to a recolonization of North America by Coragyps atratus, Gymnogyps californianus and Cathartes aura. We identified the Black Vulture, G. californianus and Vultur gryphus as priority species based on ED and "Evolutionary Distinct Globally Endangered" (EDGE) indexes, and the lowlands of Amazon River basin and the Orinoco basin and some tributaries areas of the Guiana Shield were identified as the priority areas when mapping the phylogenetic diversity. This study highlights the importance of filling knowledge gaps of species of conservation concern through the integration of evolutionary and ecological information and tools and, thus, developing adequate strategies to enhance the preservation of these species in the face of the current loss of biodiversity.

Molecular Identification and Evolutionary Divergence of the Sri Lankan Sambar Deer, Rusa unicolor (Kerr 1792).

The Sambar is one of the largest deer species distributed mainly in Asia, and it has been listed as a vulnerable species. Taxonomy based on morphological characterization has been the gold standard method used to identify the Sambar deer species. Yet, morphological identification is challenging and requires expertise. To conduct species identification and taxonomic decisions, we performed the molecular identification of R. unicolor found in Sri Lanka using DNA barcodes, COI, and Cyt b to compare the Sri Lankan R. unicolor with the Indian R. unicolor and other R. unicolor subspecies. We obtained mitochondrial DNA sequences from COI and Cyt b from blood samples collected from the wet zone in Sri Lanka. A phylogenetic tree was constructed based on the Bayesian analyses using MrBayes 3.2.7. Molecular dating was implemented in Bayesian Evolutionary Analysis Sampling Trees (BEAST v1.8.2) on the concatenated sequence using a log-normal relaxed clock and Yule species tree prior, with four categories. The results showed that the Sri Lankan R. unicolor is genetically different from the Indian R. unicolor and other R. unicolor subspecies. The divergence occurred approximately 1.1 MYA (million years ago) in the Pleistocene era. The results are essential for designing new conservation platforms for these Sambar deer species.

Deep Sequencing to Reveal Phylo-Geographic Relationships of Juquitiba Virus in Paraguay.

Several hantaviruses result in zoonotic infections of significant public health concern, causing hemorrhagic fever with renal syndrome (HFRS) or hantavirus cardiopulmonary syndrome (HCPS) in the Old and New World, respectively. Given a 35% case fatality rate, disease-causing New World hantaviruses require a greater understanding of their biology, genetic diversity, and geographical distribution. Juquitiba hantaviruses have been identified in Oligoryzomys nigripes in Brazil, Paraguay, and Uruguay. Brazil has reported the most HCPS cases associated with this virus. We used a multiplexed, amplicon-based PCR strategy to screen and deep-sequence the virus harbored within lung tissues collected from Oligoryzomys species during rodent field collections in southern (Itapúa) and western (Boquerón) Paraguay. No Juquitiba-like hantaviruses were identified in Boquerón. Herein, we report the full-length S and M segments of the Juquitiba hantaviruses identified in Paraguay from O. nigripes. We also report the phylogenetic relationships of the Juquitiba hantaviruses in rodents collected from Itapúa with those previously collected in Canindeyú. We showed, using the TN93 nucleotide substitution model, the coalescent (constant-size) population tree model, and Bayesian inference implemented in the Bayesian evolutionary analysis by sampling trees (BEAST) framework, that the Juquitiba virus lineage in Itapúa is distinct from that in Canindeyú. Our spatiotemporal analysis showed significantly different time to the most recent ancestor (TMRA) estimates between the M and S segments, but a common geographic origin. Our estimates suggest the additional geographic diversity of the Juquitiba virus within the Interior Atlantic Forest and highlight the need for more extensive sampling across this biome.

Phylogenomic and population genetics analyses of extant tomato yellow leaf curl virus strains on a global scale

Tomato yellow leaf curl virus (TYLCV) is a monopartite DNA virus with a genome size of ~ 2,800 base pairs. The virus belongs to the genus Begomovirus within the family Geminiviridae. Extant TYLCV strains are differentiated based on an established threshold of 94% genome-wide pairwise nucleotide identity. The phylogenetic relationships, diversification mechanisms, including recombination, and extent of spread within and from the center of origin for TYLCV have been reported in previous studies. However, the evolutionary relationships among strains, strains’ distribution and genomic diversification, and genetic mechanisms shaping TYLCV strains’ evolution have not been re-evaluated to consider globally representative genome sequences in publicly available sequence database, including herein newly sequenced genomes from the U.S. and Middle East, respectively. In this study, full-length genome sequences for the extant strains and isolates of TYLCV (n=818) were downloaded from the GenBank database. All previously published genome sequences, and newly sequenced TYLCV genomes of TYLCV isolates from Kuwait and USA, determined herein (n=834), were subjected to recombination analysis. To remove the ‘phylogenetic noise’ imparted by interspecific recombination, the recombinant genomes were removed from the data set, and the remaining non-recombinant genome sequences (n=423) were subjected to population genetics and Bayesian analyses. Results of the phylogeographical analysis indicated that the type strain, TYLCV-Israel, and TYLCV-Mild strain, were globally distributed, spanning Africa, America, Asia, Australia/Oceania, Europe, and New Caledonia, while the other TYLCV strains were prevalent only throughout the Middle East. The results of Bayesian evolutionary (ancestral) analysis predicted that TYLCV-Israel represents the oldest, most recent common ancestor (MRCA) (41,795 years), followed by TYLCV-Mild at 39,808 years. These were closely followed by two Iranian strains viz., TYLCV-Kerman and TYLCV-Iran at 37,529 and 36,420 years, respectively. In contrast, the most recently evolving strains were TYLCV-Kuwait and TYLCV-Kahnooj at 12,445 and 298 years, respectively. Results of the neutrality test indicated that TYLCV-Israel and TYLCV-Mild populations are undergoing purifying selection and/or population expansion, although statistically significant selection was documented for only TYLCV-Israel, based on positive selection acting on five codons.

Genomic Surveillance of Recent Dengue Outbreaks in Colombo, Sri Lanka.

All four serotypes of the dengue virus (DENV1-4) cause a phenotypically similar illness, but serial infections from different serotypes increase the risk of severe disease. Thus, genomic surveillance of circulating viruses is important to detect serotype switches that precede community outbreaks of disproportionate magnitude. A phylogenetic analysis was conducted on near full length DENV genomes sequenced from serum collected from a prospective cohort study from the Colombo district, Sri Lanka during a 28-month period using Oxford nanopore technology, and the consensus sequences were analyzed using maximum likelihood and Bayesian evolutionary analysis. From 523 patients, 328 DENV sequences were successfully generated (DENV1: 43, DENV2: 219, DENV3:66). Most circulating sequences originated from a common ancestor that was estimated to have existed from around 2010 for DENV2 and around 2015/2016 for DENV1 and DENV3. Four distinct outbreaks coinciding with monsoon rain seasons were identified during the observation period mostly driven by DENV2 cosmopolitan genotype, except for a large outbreak in 2019 contributed by DENV3 genotype I. This serotype switch did not result in a more clinically severe illness. Phylogeographic analyses showed that all outbreaks started within Colombo city and then spread to the rest of the district. In 2019, DENV3 genotype I, previously, rarely reported in Sri Lanka, is likely to have contributed to a disease outbreak. However, this did not result in more severe disease in those infected, probably due to pre-existing DENV3 immunity in the community. Targeted vector control within Colombo city before anticipated seasonal outbreaks may help to limit the geographic spread of outbreaks.

Tracing the evolutionary history of hepatitis B virus genotype H endemic to Mexico.

Hepatitis B virus (HBV) spreads efficiently among all human populations worldwide. HBV is classified into ten genotypes (A to J) with their geographic distribution and clinical features. In Mexico, HBV genotype H is the leading cause of hepatitis B and has been detected in indigenous populations, suggesting that HBV genotype H may be native to Mexico. However, little is known about the evolutionary history of HBV genotype H. Thus, we aimed to determine the age of HBV genotype H in Mexico using molecular dating techniques. Ninety-two HBV sequences of the reverse transcriptase (RT) domain of the polymerase gene (~1,251 bp) were analyzed; 48 were genotype H, 43 were genotype F, and the oldest HBV sequence from America was included as the root. All sequences were aligned, and the most recent common ancestor (TMRCA) time was calculated using the Bayesian Skyline Evolutionary Analysis. Our results estimate a TMRCA for the genotype H in Mexico of 2070.9 (667.5-4489.2) years before the present (YBP). We identified four major diversification events in genotype H, named H1, H2, H3, and H4. The TMRCA of H1 was 1213.0 (253.3-2638.3) YBP, followed by H2 1175.5 (557.5-2424.2) YBP, H3 949.6 (279.3-2105.0) YBP, and H4 1230.5 (336.3, 2756.7) YBP. We estimated that genotype H diverged from its sister genotype F around 8140.8 (1867.5-18012.8) YBP. In conclusion, this study found that genotype H in Mexico has an estimated age of 2070.9 (667.5-4489.2) YBP and has experienced at least four major diversification events since then.

Genomic characterization and phylogenetic evolution of the canine parvoviruses in Iranian dogs, a nationwide study: CPV evolutionary analysis in Iran

Canine Parvo Virus 2 (CPV-2) culminated in lots of fatalities in domestic dogs since its emergence in 1978. Mainly, it is responsible for severe hemorrhagic diarrhea, vomiting, and dehydration. CPV-2 has three main variants known as 2a, 2b, and 2c. Due to the necessity of monitoring the evolutionary parameters of the virus, and also the lack of comprehensive study of CPV2 in Iran, this study is done for the first time in this country not only to characterize Iranian CPV genomes but also to study the evolutionary parameters and phylodynamics of CPV. The phylogenetic trees were constructed using the Maximum Likelihood (ML) method. By the use of the Bayesian Monte Carlo Markov Chain (BMCMC) method, evolutionary analysis and phylodynamics of the virus were investigated. Phylogenetic results showed that all Iranian isolates were classified in the CPV-2a variant. The central part of Iran was suggested to be the origin of the virus, especially the Alborz province. Before its prevalence throughout the country, the virus circulated in the central part, in Thran, Karaj, and Qom. Mutational analysis showed a positive selection pressure of CPV-2a. Investigating the evolutionary parameters of the virus proposed 1970 to be the date of birth of the virus, with a 95% credible interval between 1953 and 1987. The effective number of infections increased dramatically from 2012 to 2015, then faced a slightly decreasing trend from 2015 to 2019. A considerable up warding pattern was witnessed from the middle of 2019, which can be taken as a concern about the risk of vaccination failure.

Sunda–Sahul floristic exchange and pathways into the Southwest Pacific: New insights from wet tropical forest trees

AbstractAimRecent investigations on the floristic exchange between Southeast Asia and Australia have shown a clear dispersal directionality bias (West to East) of wet‐adapted plant taxa. However, dispersal routes and directions of wet forest taxa into the South Pacific remain insufficiently known. We here aimed to establish the most likely routes and directions of plant dispersal into the Southwest Pacific islands.LocationSoutheast Asia, East Asia, Australia, Southwest Pacific.TaxonDysoxylum s.l. (Meliaceae). This includes Dysoxylum s.s., Didymocheton, Epicharis, Goniocheton, Pseudocarapa, and Prasoxylon.MethodWe sampled 75% of the species diversity in Dysoxylum s.l., covering the entire distribution range, all genera and major lineages. Phylogenetic relationships of 149 accessions were reconstructed using Bayesian Evolutionary Analysis and two internal constraints. The dispersal–extinction–cladogenesis variant, founder‐event speciation (DEC+J), was used for reconstructing the biogeographic history, and 100 Biogeographical Stochastic Mappings were simulated.ResultsDysoxylum s.l. originated and firstly diversified in the western part of its current distribution range (including Indochina) during the Miocene to Pliocene, followed by an overall eastern range expansion towards Malesia, Australia and the Southwest Pacific in the Pliocene.Main ConclusionsThe south‐eastward expansion of lineages into Wallacea and Australia is in temporal agreement with the convergence of the Asian and Australian tectonic plates since the Miocene. Long‐distance dispersal is the main mechanism that led to the current distribution. Two dispersal pathways into the Southwest Pacific are identified, (1) through New Guinea and the Solomon Islands to Fiji, and (2) from New Zealand to Fiji. For both routes, Fiji was an important secondary source area for dispersal into the Southwest Pacific.

Detection and genetic characterization of canine kobuvirus from stray dogs in Shanghai, China.

In this study, rectal samples collected from 60 stray dogs in dog shelters were screened for canine kobuvirus and other enteroviruses by quantitative real-time reverse transcription polymerase chain reaction. Canine kobuvirus was detected in 25% (15/60) of the samples. In the 15 positive samples, the coinfection rates of canine distemper virus, canine coronavirus, canine astrovirus, canine norovirus, and canine rotavirus were 26.67%, 20.00%, 73.33%, 0%, and 20.00%, respectively. Phylogenetic analysis based on partial VP1 sequences identified a novel canine kobuvirus that was a recombinant of canine and feline kobuvirus. Bayesian evolutionary analysis revealed that the rate of evolution of the VP1 gene of canine kobuvirus was 1.36×10-4 substitutions per site per year (95% highest posterior density interval, 6.28×10-7-4.30×10-4 substitutions per site per year). Finally, the divergence time of VP1 was around 19.44years ago (95% highest posterior density interval, 12.96-27.57years).

Climacocystaceae fam. nov. and Gloeoporellaceae fam. nov., two new families of Polyporales (Basidiomycota).

Polyporales is a diverse group of Agaricomycetes including more than 2,500 species belonging to 255 genera and 18 families. Recently, many studies focused on the classification of Polyporales, but the familial placements of some taxa remain uncertain. In this study, two new families, Climacocystaceae and Gloeoporellaceae of Polyporales, are proposed based on morphological characters and molecular data. Phylogenetic analyses of the two new families are inferred from the DNA sequences of the internal transcribed spacer regions (ITS), the large subunit of nuclear ribosomal RNA gene (nLSU), the largest subunit of RNA polymerase II gene (RPB1), the second largest subunit of RNA polymerase II gene (RPB2), and the translation elongation factor 1-α gene (TEF1). Furthermore, the divergence time of Polyporales was estimated as an additional taxonomic criterion based on the conserved regions of five DNA fragments (5.8S, nLSU, RPB1, RPB2, and TEF1). Bayesian evolutionary analysis revealed that the ancestor of Polyporales splits with a mean stem age of 136.53 Mya with a 95% highest posterior density (HPD) of 118.08-158.06 Mya. The mean stem ages of the families within Polyporales originated between 66.02 and 119.22 Mya, of which Climacocystaceae occurred in a mean stem age of 77.49 Mya with a 95% HPD of 61.45-93.16 Mya, and Gloeoporellaceae occurred in a mean stem age of 88.06 Mya with a 95% HPD of 67.15-107.76 Mya.

Out of the Qinghai-Tibetan plateau: Origin, evolution and historical biogeography of Morchella (both Elata and Esculenta clades).

Morchella has become a research hotspot because of its wide distribution, delicious taste, and phenotypic plasticity. The Qinghai-Tibet Plateau subkingdoms (QTPs) are known as the cradle of Ice age biodiversity. However, the diversity of Morchella in the QTPs has been poorly investigated, especially in phylogenetic diversity, origin, and biogeography. The genealogical concordance phylogenetic species recognition (GCPSR, based on Bayesian evolutionary analysis using sequences from the internal transcribed spacer (ITS), nuclear large subunit rDNA (nrLSU), translation elongation factor 1-α (EF1-α), and the largest and second largest subunits of RNA polymerase II (RPB1 and RPB2)), differentiation time estimation, and ancestral region reconstruction were used to infer Morchella's phylogenetic relationships and historical biogeography in the QTPs. Firstly, a total of 18 Morchella phylogenetic species are recognized in the QTPs, including 10 Elata clades and 8 Esculenta clades of 216 individuals Secondly, the divergences of the 18 phylogenetic species were 50.24-4.20 Mya (Eocene-Pliocene), which was closely related to the geological activities in the QTPs. Furthermore, the ancestor of Morchella probably originated in the Northern regions (Qilian Shan, Elata cade) and southwestern regions (Shangri-La, Esculenta clade) of QTPs and might have migrated from North America (Rufobrunnea clade) via Beringian Land Bridge (BLB) and Long-Distance Dispersal (LDD) expansions during the Late Cretaceous. Moreover, as the cradle of species origin and diversity, the fungi species in the QTPs have spread out and diffused to Eurasia and South Africa starting in the Paleogene Period. This is the first report that Esculenta and Elata clade of Morchella originated from the QTPs because of orogenic, and rapid differentiation of fungi is strongly linked to geological uplift movement and refuge in marginal areas of the QTPs. Our findings contribute to increasing the diversity of Morchella and offer more evidence for the origin theory of the QTPs.

A pioneering method to identify bovine horn trophy: A combined morphometric and DNA-based approach in wildlife forensics

We describe a protocol for identifying wild versus domestic bovid horns using multidisciplinary approaches for implementing applicable wildlife laws. We used DNA and morphometry methods to identify four wild and two domestic bovid species from horn trophies. It helped in testing a horn trophy suspected from a wild buffalo that was seized by authorities and sent to us for species identification. Two measurements, circumference at the base (CAB) and length on the front curve (LOFC), were measured for the seized and different bovids horns. The three-dimensional (3-D) scatter plot, generated by the values of CAB, LOFC and CAB/LOFC, differentiated the different bovid species into distinct clusters and placed the seized horn in the proximity of domestic water buffalos. The Bayesian evolutionary analysis of the partial D-loop gene (521 bp) placed the seized horn in a clade with swamp buffalos, which are considered domestic. Hence, the morphometric and the DNA analyses showed consensus in concluding that the seized horn originated from a domestic buffalo. The current protocol may differentiate between wild water buffalo, wild yak, guar, takin and domestic water buffalos (river and swamp). It may also help address the illegal wildlife trade of different bovid species at national and international levels.

Phylodynamics of Highly Pathogenic Avian Influenza A(H5N1) Virus Circulating in Indonesian Poultry.

After its first detection in 1996, the highly pathogenic avian influenza A(H5Nx) virus has spread extensively worldwide. HPAIv A(H5N1) was first detected in Indonesia in 2003 and has been endemic in poultry in this country ever since. However, Indonesia has limited information related to the phylodynamics of HPAIv A(H5N1) in poultry. The present study aimed to increase the understanding of the evolution and temporal dynamics of HPAIv H5N1 in Indonesian poultry between 2003 and 2016. To this end, HPAIv A(H5N1) hemagglutinin sequences of viruses collected from 2003 to 2016 were analyzed using Bayesian evolutionary analysis sampling trees. Results indicated that the common ancestor of Indonesian poultry HPAIv H5N1 arose approximately five years after the common ancestor worldwide of HPAI A(H5Nx). In addition, this study indicated that only two introductions of HPAIv A(H5N1) occurred, after which these viruses continued to evolve due to extensive spread among poultry. Furthermore, this study revealed the divergence of H5N1 clade 2.3.2.1c from H5N1 clade 2.3.2.1b. Both clades 2.3.2.1c and 2.3.2.1b share a common ancestor, clade 1, suggesting that clade 2.3.2.1 originated and diverged from China and other Asian countries. Since there was limited sequence and surveillance data for the HPAIv A(H5N1) from wild birds in Indonesia, the exact role of wild birds in the spread of HPAIv in Indonesia is currently unknown. The evolutionary dynamics of the Indonesian HPAIv A(H5N1) highlight the importance of continuing and improved genomic surveillance and adequate control measures in the different regions of both the poultry and wild birds. Spatial genomic surveillance is useful to take adequate control measures. Therefore, it will help to prevent the future evolution of HPAI A(H5N1) and pandemic threats.

Genomic and evolutionary characteristics of G9P[8], the dominant group a rotavirus in China (2016-2018).

G9P[8] became the predominant rotavirus A (RVA) genotype in China in 2012. To evaluate its genetic composition at the whole-genome level, 115 G9P[8] RVA strains isolated from children under 5 years old were sequenced and characterized. All 13 strains in 2016 and 2017 and an additional 54 strains in 2018 were genotyped as G9-P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1. The other 48 strains in 2018 were all genotyped as G9-P[8]-I1-R1-C1-M1-A1-N1-T1-E2-H1, with the NSP4 gene characterized as a DS-1-like genotype. The time of the most recent common ancestor (tMRCA) and evolution rates of the VP7, VP4, and NSP4 (E1 and E2) genes of these strains were estimated by Bayesian evolutionary dynamics analysis. We estimated the evolution rates (nt substitutions per site per year) as 1.38 × 10–3 [the 95% highest posterior density (HPD) was 1.09–1.72 × 10–3] for VP7, 0.87 × 10–3 (95% HPD: 0.75–1.00 × 10–3) for VP4, 0.56 × 10–3 (95% HPD: 0.41–0.73 × 10–3) for NSP4-E1, and 1.35 × 10–3 (95% HPD: 0.92–1.86 × 10–3) for NSP4-E2. The tMRCA was estimated to be 1935.4 (95% HPD: 1892.4–1961.3) for VP7, 1894.3 (95% HPD: 1850.5–1937.8) for VP4, 1929.4 (95% HPD: 1892.4–1961.3) for NSP4-E1, and 1969.2 (95% HPD: 1942.2–1985.3) for NSP4-E2. The baseline genetic information in this study is expected to improve our understanding of the genomic and evolutionary characteristics of the rotavirus genome. Furthermore, it will provide a basis for the development of next-generation rotavirus vaccines for humans.