The Ecology and Evolution of Cancer: The Ultra-Microevolutionary Process.

Although tumorigenesis has been accepted as an evolutionary process ( 20 , 102 ), many forces may operate differently in cancers than in organisms, as they evolve at vastly different time scales. Among such forces, natural selection, here defined as differential cellular proliferation among distinct somatic cell genotypes, is particularly interesting because its action might be thwarted in multicellular organisms ( 20 , 29 ). In this review, selection is analyzed in two stages of cancer evolution: Stage I is the evolution between tumors and normal tissues, and Stage II is the evolution within tumors. The Cancer Genome Atlas (TCGA) data show a low degree of convergent evolution in Stage I, where genetic changes are not extensively shared among cases. An equally important, albeit much less highlighted, discovery using TCGA data is that there is almost no net selection in cancer evolution. Both positive and negative selection are evident but they neatly cancel each other out, rendering total selection ineffective in the absence of recombination. The efficacy of selection is even lower in Stage II, where neutral (non-Darwinian) evolution is increasingly supported by high-density sampling studies ( 81 , 123 ). Because natural selection is not a strong deterministic force, cancers usually evolve divergently even in similar tissue environments.

Evolution under Spatially Heterogeneous Selection in Solid Tumors.

Spatial genetic and phenotypic diversity within solid tumors has been well documented. Nevertheless, how this heterogeneity affects temporal dynamics of tumorigenesis has not been rigorously examined because solid tumors do not evolve as the standard population genetic model due to the spatial constraint. We therefore, propose a neutral spatial (NS) model whereby the mutation accumulation increases toward the periphery; the genealogical relationship is spatially determined and the selection efficacy is blunted (due to kin competition). In this model, neutral mutations are accrued and spatially distributed in manners different from those of advantageous mutations. Importantly, the distinctions could be blurred in the conventional model. To test the NS model, we performed a three-dimensional multiple microsampling of two hepatocellular carcinomas. Whole-genome sequencing (WGS) revealed a 2-fold increase in mutations going from the center to the periphery. The operation of natural selection can then be tested by examining the spatially determined clonal relationships and the clonal sizes. Due to limited migration, only the expansion of highly advantageous clones can sweep through a large part of the tumor to reveal the selective advantages. Hence, even multiregional sampling can only reveal a fraction of fitness differences in solid tumors. Our results suggest that the NS patterns are crucial for testing the influence of natural selection during tumorigenesis, especially for small solid tumors.

Contrasting Patterns in the Early Stage of SARS-CoV-2 Evolution between Humans and Minks.

One of the unique features of SARS-CoV-2 is its apparent neutral evolution during the early pandemic (before February 2020). This contrasts with the preceding SARS-CoV epidemics, where viruses evolved adaptively. SARS-CoV-2 may exhibit a unique or adaptive feature which deviates from other coronaviruses. Alternatively, the virus may have been cryptically circulating in humans for a sufficient time to have acquired adaptive changes before the onset of the current pandemic. To test the scenarios above, we analyzed the SARS-CoV-2 sequences from minks (Neovision vision) and parental humans. In the early phase of the mink epidemic (April to May 2020), nonsynonymous to synonymous mutation ratio per site in the spike protein is 2.93, indicating a selection process favoring adaptive amino acid changes. Mutations in the spike protein were concentrated within its receptor-binding domain and receptor-binding motif. An excess of high-frequency derived variants produced by genetic hitchhiking was found during the middle (June to July 2020) and late phase I (August to September 2020) of the mink epidemic. In contrast, the site frequency spectra of early SARS-CoV-2 in humans only show an excess of low-frequency mutations, consistent with the recent outbreak of the virus. Strong positive selection in the mink SARS-CoV-2 implies that the virus may not be preadapted to a wide range of hosts and illustrates how a virus evolves to establish a continuous infection in a new host. Therefore, the lack of positive selection signal during the early pandemic in humans deserves further investigation.

Factors associated with viral rebound among COVID-19 patients receiving oral antivirals.

BACKGROUND: COVID-19 rebound is usually reported among patients experiencing concurrent symptomatic and viral rebound. But longitudinal viral RT-PCR results from early stage to rebound of COVID-19 was less characterized. Further, identifying the factors associated with viral rebound after nirmatrelvir-ritonavir (NMV/r) and molnupiravir may expand understanding of COVID-19 rebound. METHODS: We retrospectively analyzed clinical data and sequential viral RT-PCR results from COVID-19 patients receiving oral antivirals between April and May, 2022. Viral rebound was defined by the degree of viral load increase (ΔCt ≥ 5 units). RESULTS: A total of 58 and 27 COVID-19 patients taking NMV/r and molnupiravir, respectively, were enrolled. Patients receiving NMV/r were younger, had fewer risk factors for disease progression and faster viral clearance rate compared to those receiving molnupiravr (All P < 0.05). The overall proportion of viral rebound (n = 11) was 12.9%, which was more common among patients receiving NMV/r (10 [17.2%] vs. 1 [3.7%], P = 0.16). Of them, 5 patients experienced symptomatic rebound, suggesting the proportion of COVID-19 rebound was 5.9%. The median interval to viral rebound was 5.0 (interquartile range, 2.0-8.0) days after completion of antivirals. Initial lymphopenia (<0.8 × 109/L) was associated with viral rebound among overall population (adjusted odds ratio [aOR], 5.34; 95% confidence interval [CI], 1.33-21.71), and remained significant (aOR, 4.50; 95% CI, 1.05-19.25) even when patients receiving NMV/r were considered. CONCLUSION: Our data suggest viral rebound after oral antivirals may be more commonly observed among lymphopenic individuals in the context of SARS-CoV-2 Omicron BA.2 variant.

The origin and past demography of murine astrovirus 1 in laboratory mice.

Astroviruses are non-enveloped, positive-sense, ssRNA viruses and often associated with gastrointestinal diseases. Murine astrovirus (MuAstV) was first confirmed in a laboratory mouse colony in 2011. Although infected mice do not present significant clinical symptoms, the virus might interfere with research results. A recent surveillance has shown that MuAstV is highly prevalent in laboratory mice. The aims of the present study were to identify and characterize MuAstV strains as well as to investigate the prevalence rate of viral RNA in laboratory mice in Taiwan, and to estimate the origin and past population demography of MuAstVs. Based on molecular surveillance, MuAstV RNA was detected in 45.7 % of laboratory mice (48/105) from seven of nine colonies. Three fully sequenced MuAstV strains, MuAstV TW1, TW2 and TW3, exhibited 89.1-94.4 % and 89.1-90.0 % nucleotide identities with the reference strains MuAstV STL1 and STL2, respectively. Phylogenetic analyses of the partial regions of the RNA-dependent RNA polymerase (RdRp) and capsid protein (CP) genes of 18 Taiwan strains along with other astroviruses revealed that there are three distinct lineages of mouse astrovirus, MuAstV1, MuAstV2 and mouse astrovirus JF755422. The mutation rates of MuAstV1 were 2.6×10-4 and 6.2×10-4 substitutions/site/year for the RdRp and CP regions, respectively. Based on the above molecular clock, the colonization of MuAstV1 in laboratory mice was between 1897 and 1912, in good agreement with the establishment of 'modern' laboratory mouse facilities. Since its initial infection, the population size of MuAstV1 has increased 15-60-fold, probably consistent with the increased use of laboratory mice. In conclusion, MuAstV1 has been associated with modern laboratory mice since the beginning, and its influence on research results may require further investigation.

Genetic Load and Potential Mutational Meltdown in Cancer Cell Populations.

Large genomes with elevated mutation rates are prone to accumulating deleterious mutations more rapidly than natural selection can purge (Muller's ratchet). As a consequence, it may lead to the extinction of small populations. Relative to most unicellular organisms, cancer cells, with large and nonrecombining genome and high mutation rate, could be particularly susceptible to such "mutational meltdown." However, the most common type of mutation in organismal evolution, namely, deleterious mutation, has received relatively little attention in the cancer biology literature. Here, by monitoring single-cell clones from HeLa cell lines, we characterize deleterious mutations that retard the rate of cell proliferation. The main mutation events are copy number variations (CNVs), which, estimated from fitness data, happen at a rate of 0.29 event per cell division on average. The mean fitness reduction, estimated reaching 18% per mutation, is very high. HeLa cell populations therefore have very substantial genetic load and, at this level, natural population would likely face mutational meltdown. We suspect that HeLa cell populations may avoid extinction only after the population size becomes large enough. Because CNVs are common in most cell lines and tumor tissues, the observations hint at cancer cells' vulnerability, which could be exploited by therapeutic strategies.

Lymphocyte Antigen 6 Complex, Locus C+ Monocytes and Kupffer Cells Orchestrate Liver Immune Responses Against Hepatitis B Virus in Mice.

To understand the mechanism(s) of age-dependent outcomes of hepatitis B virus (HBV) infection in humans, we previously established an age-related HBV mouse model in which 6-week-old (N6W) C3H/HeN mice exhibited virus tolerance whereas 12-week-old (N12W) counterparts presented virus clearance. By investigating the hepatic myeloid cell dynamics in mice of these two ages, we aim to identify factors associated with HBV clearance. C3H/HeN mice were transfected with an HBV plasmid by hydrodynamic injection. Serum HBV markers were monitored weekly. Hepatic leucocyte populations and their cytokine/chemokine productions were examined at baseline, day 3 (D3), day 7 (D7), and day 14 after injection. C-C chemokine receptor type 2 (CCR2) antagonist and clodronate (CLD) were respectively administered to N12W and N6W mice to study the roles of lymphocyte antigen 6 complex, locus C (Ly6C)+ monocytes and Kupffer cells (KCs) in viral clearance. N12W mice had a significantly higher number of TNF-α-secreting Ly6C+ monocytes and fewer IL-10-secreting KCs at D3 in the liver than their younger N6W counterparts after HBV transfection. In addition, the elevated number of interferon-γ+ TNF-α+ CD8+ T cells at D7 was only seen in the older cohort. The enhanced Ly6C+ monocyte induction in N12W mice resulted from elevated C-C motif chemokine ligand 2 (CCL2) secretion by hepatocytes. CCR2 antagonist administration hampered Ly6C+ monocyte recruitment and degree of KC reduction and delayed HBV clearance in N12W animals. Depletion of KCs by CLD liposomes enhanced Ly6C+ monocyte recruitment and accelerated HBV clearance in N6W mice. Conclusions: Ly6C+ monocytes and KCs may, respectively, represent the resistance and tolerance arms of host defenses. These two cell types play an essential role in determining HBV clearance/tolerance. Manipulation of these cells is a promising avenue for immunotherapy of HBV-related liver diseases.

The origin and underlying driving forces of the SARS-CoV-2 outbreak.

BACKGROUND: SARS-CoV-2 began spreading in December 2019 and has since become a pandemic that has impacted many aspects of human society. Several issues concerning the origin, time of introduction to humans, evolutionary patterns, and underlying force driving the SARS-CoV-2 outbreak remain unclear. METHOD: Genetic variation in 137 SARS-CoV-2 genomes and related coronaviruses as of 2/23/2020 was analyzed. RESULT: After correcting for mutational bias, the excess of low frequency mutations on both synonymous and nonsynonymous sites was revealed which is consistent with the recent outbreak of the virus. In contrast to adaptive evolution previously reported for SARS-CoV during its brief epidemic in 2003, our analysis of SARS-CoV-2 genomes shows signs of relaxation. The sequence similarity in the spike receptor binding domain between SARS-CoV-2 and a sequence from pangolin is probably due to an ancient intergenomic introgression that occurred approximately 40 years ago. The current outbreak of SARS-CoV-2 was estimated to have originated on 12/11/2019 (95% HPD 11/13/2019-12/23/2019). The effective population size of the virus showed an approximately 20-fold increase from the onset of the outbreak to the lockdown of Wuhan (1/23/2020) and ceased to increase afterwards, demonstrating the effectiveness of social distancing in preventing its spread. Two mutations, 84S in orf8 protein and 251 V in orf3 protein, occurred coincidentally with human intervention. The former first appeared on 1/5/2020 and plateaued around 1/23/2020. The latter rapidly increased in frequency after 1/23/2020. Thus, the roles of these mutations on infectivity need to be elucidated. Genetic diversity of SARS-CoV-2 collected from China is two times higher than those derived from the rest of the world. A network analysis found that haplotypes collected from Wuhan were interior and had more mutational connections, both of which are consistent with the observation that the SARS-CoV-2 outbreak originated in China. CONCLUSION: SARS-CoV-2 might have cryptically circulated within humans for years before being discovered. Data from the early outbreak and hospital archives are needed to trace its evolutionary path and determine the critical steps required for effective spreading.

Clonality of Fluconazole-Nonsusceptible Candida tropicalis in Bloodstream Infections, Taiwan, 2011-2017.

Candida tropicalis is the leading cause of non-C. albicans candidemia in tropical Asia and Latin America. We evaluated isolates from 344 patients with an initial episode of C. tropicalis candidemia. We found that 58 (16.9%) patients were infected by fluconazole-nonsusceptible (FNS) C. tropicalis with cross resistance to itraconazole, voriconazole, and posaconazole; 55.2% (32/58) of patients were azole-naive. Multilocus sequence typing analysis revealed FNS isolates were genetically closely related, but we did not see time- or place-clustering. Among the diploid sequence types (DSTs), we noted DST225, which has been reported from fruit in Taiwan and hospitals in Beijing, China, as well as DST376 and DST505-7, which also were reported from hospitals in Shanghai, China. Our findings suggest cross-boundary expansion of FNS C. tropicalis and highlight the importance of active surveillance of clinical isolates to detect dissemination of this pathogen and explore potential sources in the community.

Effects of indoor residual spraying and outdoor larval control on Anopheles coluzzii from São Tomé and Príncipe, two islands with pre-eliminated malaria.

BACKGROUND: Vector control is a key component of malaria prevention. Two major vector control strategies have been implemented in São Tomé and Príncipe (STP), indoor residual spraying (IRS) and outdoor larval control using Bacillus thuringiensis israelensis (Bti). This study evaluated post-intervention effects of control strategies on vector population density, composition, and knockdown resistance mutation, and their implications for malaria epidemiology in STP. METHODS: Mosquitoes were collected by indoor and outdoor human landing catches and mosquito light traps in seven districts. Mosquito density was calculated by numbers of captured adult mosquitoes/house/working hour. Mitochondrial cytochrome c oxidase subunit I (COI) was PCR amplified and sequenced to understand the spatial-temporal population composition of malaria vector in STP. Knockdown resistance L1014F mutation was detected using allele-specific PCR. To estimate the malaria transmission risks in STP, a negative binomial regression model was constructed. The response variable was monthly incidence, and the explanatory variables were area, rainfall, entomological inoculation rate (EIR), and kdr mutation frequency. RESULTS: Malaria vector in STP is exophilic Anopheles coluzzii with significant population differentiation between Príncipe and São Tomé (mean FST = 0.16, p < 0.001). Both vector genetic diversity and knockdown resistance mutation were relatively low in Príncipe (mean of kdr frequency = 15.82%) compared to São Tomé (mean of kdr frequency = 44.77%). Annual malaria incidence rate in STP had been rapidly controlled from 37 to 2.1% by three rounds of country-wide IRS from 2004 to 2007. Long-term application of Bti since 2007 kept the mosquito density under 10 mosquitoes/house/hr/month, and malaria incidence rate under 5% after 2008, except for a rising that occurred in 2012 (incidence rate = 6.9%). Risk factors of area (São Tomé compared to Príncipe), rainfall, outdoor EIR, and kdr mutation frequency could significantly increase malaria incidence by 9.33-11.50, 1.25, 1.07, and 1.06 fold, respectively. CONCLUSIONS: Indoor residual spraying could rapidly decrease Anopheles density and malaria incidence in STP. Outdoor larval control using Bti is a sustainable approach for controlling local vector with exophilic feature and insecticide resistance problem. Vector control interventions should be intensified especially at the north-eastern part of São Tomé to minimize impacts of outbreaks.

Age-related immune clearance of hepatitis B virus infection requires the establishment of gut microbiota.

A unique feature of hepatitis B virus (HBV) infection in humans is that viral clearance heavily depends on the age of exposure. However, the reason for this remains unclear. Here we show that gut microbiota contribute to the age dependence of HBV immunity in a hydrodynamic transfection mouse model. Although adult (12-wk-old) C3H/HeN mice cleared HBV within 6 wk postinjection (wpi), their young (6-wk-old) counterparts remained HBV-positive at 26 wpi. Sterilization of gut microbiota from 6 to 12 wk of age using antibiotics prevented adult mice from rapidly clearing HBV. Young mice with the Toll-like-receptor (TLR) 4 mutation (C3H/HeJ) exhibited rapid HBV clearance. The results suggest that an immuno-tolerating pathway to HBV prevailed in young mice, before the establishment of gut bacteria, through a TLR4-dependent pathway and that the maturation of gut microbiota in adult mice stimulated liver immunity, resulting in rapid HBV clearance.

De novo assembly of highly polymorphic metagenomic data using in situ generated reference sequences and a novel BLAST-based assembly pipeline.

BACKGROUND: The accuracy of metagenomic assembly is usually compromised by high levels of polymorphism due to divergent reads from the same genomic region recognized as different loci when sequenced and assembled together. A viral quasispecies is a group of abundant and diversified genetically related viruses found in a single carrier. Current mainstream assembly methods, such as Velvet and SOAPdenovo, were not originally intended for the assembly of such metagenomics data, and therefore demands for new methods to provide accurate and informative assembly results for metagenomic data. RESULTS: In this study, we present a hybrid method for assembling highly polymorphic data combining the partial de novo-reference assembly (PDR) strategy and the BLAST-based assembly pipeline (BBAP). The PDR strategy generates in situ reference sequences through de novo assembly of a randomly extracted partial data set which is subsequently used for the reference assembly for the full data set. BBAP employs a greedy algorithm to assemble polymorphic reads. We used 12 hepatitis B virus quasispecies NGS data sets from a previous study to assess and compare the performance of both PDR and BBAP. Analyses suggest the high polymorphism of a full metagenomic data set leads to fragmentized de novo assembly results, whereas the biased or limited representation of external reference sequences included fewer reads into the assembly with lower assembly accuracy and variation sensitivity. In comparison, the PDR generated in situ reference sequence incorporated more reads into the final PDR assembly of the full metagenomics data set along with greater accuracy and higher variation sensitivity. BBAP assembly results also suggest higher assembly efficiency and accuracy compared to other assembly methods. Additionally, BBAP assembly recovered HBV structural variants that were not observed amongst assembly results of other methods. Together, PDR/BBAP assembly results were significantly better than other compared methods. CONCLUSIONS: Both PDR and BBAP independently increased the assembly efficiency and accuracy of highly polymorphic data, and assembly performances were further improved when used together. BBAP also provides nucleotide frequency information. Together, PDR and BBAP provide powerful tools for metagenomic data studies.

New insights into the evolutionary rate of hepatitis B virus at different biological scales.

UNLABELLED: The evolutionary rates of hepatitis B virus (HBV) estimated using contemporary sequences are 10(2) to 10(4) times higher than those derived from archaeological and genetic evidence. This discrepancy makes the origin of HBV and the time scale of its spread, both of which are critical for studying the burden of HBV pathogenicity, largely unresolved. To evaluate whether the dual demands (i.e., adaptation within hosts and colonization between hosts) of the viral life cycle affect this conundrum, the HBV quasispecies dynamics within and among hosts from a family consisting of a grandmother, her 5 children, and her 2 granddaughters, all of whom presumably acquired chronic HBV through mother-to-infant transmission, were examined by PCR cloning and next-generation sequencing methods. We found that the evolutionary rate of HBV between hosts was considerably lower than that within hosts. Moreover, the between-host substitution rates of HBV decreased as transmission numbers between individuals increased. Both observations were due primarily to changes at nonsynonymous rather than synonymous sites. There were significantly more multiple substitutions than expected for random mutation processes, and 97% of substitutions were changed from common to rare amino acid residues in the database. Continual switching between colonization and adaptation resulted in a rapid accumulation of mutations at a limited number of positions, which quickly became saturated, whereas substitutions at the remaining regions occurred at a much lower rate. Our study may help to explain the time-dependent HBV substitution rates reported in the literature and provide new insights into the origin of the virus. IMPORTANCE: It is known that the estimated hepatitis B virus (HBV) substitution rate is time dependent, but the reason behind this observation is still elusive. We hypothesize that owing to the small genome size of HBV, transmission between hosts and adaptation within hosts must exhibit high levels of fitness trade-offs for the virus. By studying the HBV quasispecies dynamics for a chain of sequentially infected transmissions within a family, we found the HBV substitution rate between patients to be negatively correlated with the number of transmissions. Continual switching between hosts resulted in a rapid accumulation of mutations at a limited number of genomic sites, which quickly became saturated in the short term. Nevertheless, substitutions at the remaining regions occurred at a much lower rate. Therefore, the HBV substitution rate decreased as the divergence time increased.

Micro-evolution of the hepatitis B virus genome in hepatitis B e-antigen-positive carriers: comparison of genotypes B and C at various immune stages.

BACKGROUND AND AIM: Patients with hepatitis B virus (HBV) genotype B infection experience hepatitis B e-antigen (HBeAg) seroconversion at an earlier stage than do patients with genotype C infection. Therefore, this study investigated whether the differential phenotypes are related to HBV genomic evolution. METHODS: Thirty-three HBeAg-positive patients with a mean follow-up of 3.1 years were enrolled: 16 at the immune tolerance stage (group I) and 17 at the immune clearance stage (group II). The evolution rates of paired viral genomes at enrollment and at the final follow-up in the full-length genome (µf), nonoverlapping regions (synonymous [µs] and nonsynonymous [µa]), and overlapping regions (µ) were calculated. The evolution rates were then compared according to serum alanine aminotransferase (ALT) levels and HBV genotype. RESULTS: The overall µf evolution rate was lower in group I than in group II (1.4 × 10(-5) ± 3.3 × 10(-5) vs 1.2 × 10(-3) ± 1.2 × 10(-3) nucleotide substitution/site/year, P < 0.001). We observed similar results for the µs, µa, and µ evolution rates. All evolution parameters were comparable between genotypes B and C. We determined a positive correlation between µa/y and the area under the average ALT time curve in genotype B (R(2) = 0.6935, P < 0.0001), but not in genotype C (R(2) = 0.1606, P = 0.124). CONCLUSION: The evolution rate of the HBV genome is higher at the immune clearance stage than at the immune tolerance stage. Host immune selection might play a role in triggering evolution of genotype B.

Molecular characterization of cryptically circulating rabies virus from ferret badgers, Taiwan.

After the last reported cases of rabies in a human in 1959 and a nonhuman animal in 1961, Taiwan was considered free from rabies. However, during 2012-2013, an outbreak occurred among ferret badgers in Taiwan. To examine the origin of this virus strain, we sequenced 3 complete genomes and acquired multiple rabies virus (RABV) nucleoprotein and glycoprotein sequences. Phylogeographic analyses demonstrated that the RABV affecting the Taiwan ferret badgers (RABV-TWFB) is a distinct lineage within the group of lineages from Asia and that it has been differentiated from its closest lineages, China I (including isolates from Chinese ferret badgers) and the Philippines, 158-210 years ago. The most recent common ancestor of RABV-TWFB originated 91-113 years ago. Our findings indicate that RABV could be cryptically circulating in the environment. An understanding of the underlying mechanism might shed light on the complex interaction between RABV and its host.

Influence of gene flow on divergence dating - implications for the speciation history of Takydromus grass lizards.

Dating the time of divergence and understanding speciation processes are central to the study of the evolutionary history of organisms but are notoriously difficult. The difficulty is largely rooted in variations in the ancestral population size or in the genealogy variation across loci. To depict the speciation processes and divergence histories of three monophyletic Takydromus species endemic to Taiwan, we sequenced 20 nuclear loci and combined with one mitochondrial locus published in GenBank. They were analysed by a multispecies coalescent approach within a Bayesian framework. Divergence dating based on the gene tree approach showed high variation among loci, and the divergence was estimated at an earlier date than when derived by the species-tree approach. To test whether variations in the ancestral population size accounted for the majority of this variation, we conducted computer inferences using isolation-with-migration (IM) and approximate Bayesian computation (ABC) frameworks. The results revealed that gene flow during the early stage of speciation was strongly favoured over the isolation model, and the initiation of the speciation process was far earlier than the dates estimated by gene- and species-based divergence dating. Due to their limited dispersal ability, it is suggested that geographical isolation may have played a major role in the divergence of these Takydromus species. Nevertheless, this study reveals a more complex situation and demonstrates that gene flow during the speciation process cannot be overlooked and may have a great impact on divergence dating. By using multilocus data and incorporating Bayesian coalescence approaches, we provide a more biologically realistic framework for delineating the divergence history of Takydromus.

Multilocus approach to clarify species status and the divergence history of the Bemisia tabaci (Hemiptera: Aleyrodidae) species complex.

The sweet potato whitefly, Bemisia tabaci, is a highly differentiated species complex. Despite consisting of several morphologically indistinguishable entities and frequent invasions on all continents with important associated economic losses, the phylogenetic relationships, species status, and evolutionary history of this species complex is still debated. We sequenced and analyzed one mitochondrial and three single-copy nuclear genes from 9 of the 12 genetic groups of B. tabaci and 5 closely related species. Bayesian species delimitation was applied to investigate the speciation events of B. tabaci. The species statuses of the different genetic groups were strongly supported under different prior settings and phylogenetic scenarios. Divergence histories were estimated by a multispecies coalescence approach implemented in (*)BEAST. Based on mitochondrial locus, B. tabaci was originated 6.47 million years ago (MYA). Nevertheless, the time was 1.25MYA based on nuclear loci. According to the method of approximate Bayesian computation, this difference is probably due to different degrees of migration among loci; i.e., although the mitochondrial locus had differentiated, gene flow at nuclear loci was still possible, a scenario similar to parapatric mode of speciation. This is the first study in whiteflies using multilocus data and incorporating Bayesian coalescence approaches, both of which provide a more biologically realistic framework for delimiting species status and delineating the divergence history of B. tabaci. Our study illustrates that gene flow during species divergence should not be overlooked and has a great impact on divergence time estimation.

Rapid growth of a hepatocellular carcinoma and the driving mutations revealed by cell-population genetic analysis of whole-genome data.

We present the analysis of the evolution of tumors in a case of hepatocellular carcinoma. This case is particularly informative about cancer growth dynamics and the underlying driving mutations. We sampled nine different sections from three tumors and seven more sections from the adjacent nontumor tissues. Selected sections were subjected to exon as well as whole-genome sequencing. Putative somatic mutations were then individually validated across all 9 tumor and 7 nontumor sections. Among the mutations validated, 24 were amino acid changes; in addition, 22 large indels/copy number variants (>1 Mb) were detected. These somatic mutations define four evolutionary lineages among tumor cells. Separate evolution and expansion of these lineages were recent and rapid, each apparently having only one lineage-specific protein-coding mutation. Hence, by using a cell-population genetic definition, this approach identified three coding changes (CCNG1, P62, and an indel/fusion gene) as tumor driver mutations. These three mutations, affecting cell cycle control and apoptosis, are functionally distinct from mutations that accumulated earlier, many of which are involved in inflammation/immunity or cell anchoring. These distinct functions of mutations at different stages may reflect the genetic interactions underlying tumor growth.

Genetic and physiological data suggest demographic and adaptive responses in complex interactions between populations of figs (Ficus pumila) and their pollinating wasps (Wiebesia pumilae).

To study interactions between host figs and their pollinating wasps and the influence of climatic change on their genetic structures, we sequenced cytoplasmic and nuclear genes and genotyped nuclear microsatellite loci from two varieties of Ficus pumila, the widespread creeping fig and endemic jelly fig, and from their pollinating wasps, Wiebesia pumilae, found in Taiwan and on nearby offshore islands. Great divergence in the mitochondrial cytochrome c oxidase subunit I (mtCOI) with no genetic admixture in nuclear markers indicated that creeping- and jelly-fig wasps are genetically distinct. Compared with creeping-fig wasps, jelly-fig wasps also showed better resistance under cold (20 °C) than warm (25 and 30 °C) conditions in a survival test, indicating their adaptation to a cold environment, which may have facilitated population expansion during the ice age as shown by a nuclear intron and 10 microsatellite loci. An excess of amino acid divergence and a pattern of too many rare mtCOI variants of jelly-fig wasps as revealed by computer simulations and neutrality tests implied the effect of positive selection, which we hypothesize was associated with the cold-adaptation process. Chloroplast DNA of the two fig plants was completely segregated, with signs of genetic admixture in nuclear markers. As creeping- and jelly-fig wasps can pollinate creeping figs, occasional gene flow between the two figs is thus possible. Therefore, it is suggested that pollinating wasps may be playing an active role in driving introgression between different types of host fig.

Clinical features and phylogenetic analysis of Coxsackievirus A9 in Northern Taiwan in 2011.

BACKGROUND: Coxsackievirus A9 (CA9) was one of the most prevalent serotype of enteroviral infections in Taiwan in 2011. After several patient series were reported in the 1960s and 1970s, few studies have focused on the clinical manifestations of CA9 infections. Our study explores and deepens the current understanding of CA9. METHODS: We analyzed the clinical presentations of 100 culture-proven CA9-infected patients in 2011 by reviewing their medical records and depicted the CA9 phylogenetic tree. RESULTS: Of the 100 patients with culture-proven CA9 infections, the mean (SD) age was 4.6 (3.4) years and the male to female ratio was 1.9. For clinical manifestations, 96 patients (96%) had fever and the mean (SD) duration of fever was 5.9 (3.4) days. Sixty one patients (61%) developed a skin rash, and the predominant pattern was a generalized non-itchy maculopapular rash without vesicular changes. While most patients showed injected throat, oral ulcers were found in only 19 cases (19%), among whom, 6 were diagnosed as herpangina. Complicated cases included: aseptic meningitis (n=8), bronchopneumonia (n=6), acute cerebellitis (n=1), and polio-like syndrome (n=1). Phylogenetic analysis for current CA9 strains is closest to the CA9 isolate 27-YN-2008 from the border area of mainland China and Myanmar. CONCLUSIONS: The most common feature of CA9 during the 2011 epidemic in Taiwan is generalized febrile exanthema rather than herpangina or hand, foot, and mouth disease. Given that prolonged fever and some complications are possible, caution should be advised in assessing patients as well as in predicting the clinical course.