Improved bootstrap confidence limits in large-scale phylogenies, with an example from Neo-Astragalus (Leguminosae)

Genomic analysis of two Chinese isolates of hyphantria cunea nucleopolyhedrovirus reveals a novel species of alphabaculovirus that infects hyphantria cunea drury (lepidoptera: arctiidae)

Background

Baculoviruses act as effective biological control agents against the invasive pest Hyphantria cunea Drury. In this study, two Chinese Hyphantria cunea nucleopolyhedrovirus (HycuNPV) isolates, HycuNPV-BJ and HycuNPV-HB, were deep sequenced and compared with the Japanese isolate, HycuNPV-N9, to determine whole-genome level diversity and evolutionary history.

Results

The divergence of the phylogenetic tree and the K2P distances based on 38 core-gene concatenated alignment revealed that two Chinese HycuNPV isolates were a novel species of Alphabaculovirus that infected Hyphantria cunea in China. The gene contents indicated significant differences in the HycuNPV genomes between the Chinese and Japanese isolates. The differences included gene deletions, acquisitions and structural transversions, but the main difference was the high number of single nucleotide polymorphisms (SNPs). In total, 10,393 SNPs, corresponding to approximately 8% of the entire HycuNPV-N9 genome sequence, were detected in the aligned reads. By analyzing non-synonymous variants, we found that hotspot mutation-containing genes had mainly unknown functions and most were early expressing genes. We found that the hycu78 gene which had early and late promoter was under positive selection. Biological activity assays revealed that the infectivity of HycuNPV-HB was greater than that of HycuNPV-BJ, and the killing speed of HycuNPV-HB was faster than that of HycuNPV-BJ. A comparison of molecular genetic characteristics indicated that the virulence differences between the two isolates were affected by SNP and structural variants, especially the homologous repeat regions.

Conclusions

The genomes of the two Chinese HycuNPV isolates were characterized, they belonged to a novel species of Alphabaculovirus that infected Hyphantria cunea in China. We inferred that the loss or gain of genetic material in the HycuNPV-HB and HycuNPV-BJ genomes resulted in new important adaptive capabilities to the H. cunea host. These results extend the current understanding of the genetic diversity of HycuNPV and will be useful for improving the applicability of this virus as a biological control agent.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08604-7.

Genome Size and Chromosome Number Evaluation of Astragalus L. sect. Hymenostegis Bunge (Fabaceae)

Astragalus section Hymenostegis is one of the important characteristic elements of thorn-cushion formations in the Irano-Turanian floristic region. In this paper, we examined the chromosome number of 17 species (15 new reports) and provide estimates of genome size for 62 individuals belonging to 38 taxa of A. sect. Hymenostegis, some species outside this section, plus two Oxytropis species. Based on chromosome counts 11 species were found to be diploid (2n = 16), four species tetraploid (2n = 32) and two taxa hexaploid (2n = 48). From genome size measurements on silica-gel dried material, three ploidy levels (2x, 4x and 6x) were inferred, with a majority of species being diploid. The 2C values reach from 2.07 pg in diploid Astragalus zohrabi to 7.16 pg in hexaploid A. rubrostriatus. We found indications that species might occur with different cytotypes. A phylogenetic framework using nrDNA ITS sequences was constructed to understand the evolution of ploidy changes and genome sizes. It showed that genome size values among the studied taxa differ only slightly within ploidy levels and are nearly constant within most species and groups of closely related taxa within the genus Astragalus. The results of this study show that there is a rather strong correlation between genome sizes and chromosome numbers in sect. Hymenostegis. The resolution of the ITS-based phylogenetic tree is too low to infer evolutionary or environmental correlations of genome size differences. Polyploidization seems to contribute to the high species number in Astragalus, however, in sect. Hymenostegis it is not the main driver of speciation.

Rice SIAH E3 Ligases Interact with RMD Formin and Affect Plant Morphology

Formins are actin-binding proteins that are key to maintaining the actin cytoskeleton in cells. However, molecular mechanisms controlling the stability of formin proteins in plants remain unknown. Here, we have identified six rice SIAH-type E3 ligases, named RIP1-6 (RMD Interacting Protein 1-6) respectively, with ubiquitination enzyme activity in vitro. All six proteins can form homo- and hetero-dimers with themselves, and hetero-dimers with type II formin RMD/OsFH5. In vivo assays showed that RIP1-6 proteins localize in the cytoplasm with a punctate distribution, and all of them interact with RMD to change its native diffuse cytoplasmic localization to match that of RIP1-6. To our surprise, degradation experiments revealed that RIP1, RIP5, and RIP6 decrease rather than increase the degradation rate of RMD. Genetic analyses revealed redundancy between these six genes; either single or double mutants did not show any obvious phenotypes. However, the sextuple rip1-6 mutant displayed dwarf height, wrinkled seeds and wider leaves that were similar to the previously reported rmd mutant, and defective microfilaments and increased flag leaf angles that were not reported in rmd mutant. Collectively, our study provides insights into the mechanisms determining formin protein stability in plants.

Development of a LAMP-Based Molecular Species Diagnosis Method for Four Major Agricultural Pests in the Genus Spodoptera (Lepidoptera: Noctuidae)

Simple Summary

Four major Spodoptera pests, S. exigua, S. frugiperda, S. litura, and S. littoralis, are widely distributed polyphagous pests affecting various crops. Despite different distribution areas, these four species cause serious damage to agriculture worldwide. As these species are morphologically similar at the larval stage, diagnostic methods have been developed and utilized for their identification. Here, we developed a loop-mediated isothermal amplification (LAMP) assay for rapid and effective species diagnosis, along with PCR, to identify Korean field-collected or overseas samples. The optimal conditions for the LAMP assay were 61 °C for 60 min with four LAMP primers. Additional loop primers increased the amplification efficiency in S. exigua, whereas increased non-specific amplification was found in other species. A broad range of DNA concentrations was observed in the LAMP assay, and the minimum detectable DNA concentration was 1 pg. The DNA release method for LAMP involved incubation of larval or adult tissue samples for 5 min at 95 °C, without a DNA extraction step. Considering the gradual diversification invasive pest incidence, this simple and accurate LAMP assay can be used for intensive field monitoring of invasive pests and integrated management of these species.

Abstract

Molecular-based species identification tools are helpful to identify tiny insect and lepidopteran pests that show morphological similarities in the larval stage and are essential for quarantine as well as agricultural research. Here, we focused on four major Spodoptera pests: S. exigua, S. frugiperda, S. litura, and S. littoralis. S. exigua and S. litura mitochondrial genome sequences were newly identified and species-specific sequence regions were identified in the cytochrome c oxidase subunit II and III regions. Species primers were designed and applied in loop-mediated isothermal amplification (LAMP) and PCR to identify Korean field-collected or overseas samples. The optimal incubation conditions for LAMP were 61 °C for 60 min with four LAMP primers. Additional loop primers increased the amplification efficiency for S. exigua, and the nonspecific amplification for other species. The LAMP assay could detect a wide range of DNA concentrations, with the range 1 ng–1 pg in dependence of four LAMP primers. The DNA-releasing technique, without DNA extraction, in the LAMP assay involved larval or adult tissue sample incubation at 95 °C for 5 min. The entire process takes approximately 70 min. This new molecular diagnostic method is simple and accurate, with application in the field and laboratory and for monitoring and ecological studies.

Development of a simple and accurate molecular tool for Spodoptera frugiperda species identification using LAMP

BACKGROUND

The fall armyworm, Spodoptera frugiperda is a native species of the Americas. First detected in western and central Africa in early 2016, it has become one of the most serious invasive lepidopteran pests in many African and Asian countries. S. frugiperda has spread very quickly; however, there are no molecular-based, simple and accurate diagnostic tools for identification of this species in the field. Methods to identify invasive S. frugiperda are urgently needed because farmers and agricultural managers have no prior experience with this pest.

RESULTS

Based on mitochondrial genome sequence alignment, a S. frugiperda-specific sequence region was identified in the transfer RNA-coding region between NADH dehydrogenase, ND3, and ND5. Using this unique region, species-diagnostic primers were designed and applied in a loop-mediated isothermal amplification (LAMP) assay and a conventional polymerase chain reaction to identify field-collected samples of S. frugiperda. The optimal incubation conditions for the LAMP assay were 61°C for 90 min with four LAMP primers; an additional loop primer increased the amplification efficiency. A response was obtained for a wide range of DNA concentrations in the LAMP assay and the minimum detectable DNA concentration was 10 pg.

CONCLUSIONS

We developed a new LAMP-based molecular diagnostic method that it is easy to use and accurate. The LAMP assay was used with a DNA-releasing technique for larval and adult samples, without a DNA extraction step, by incubating the tissue sample at 95°C for 5 min. This method can be applied in intensive field monitoring of S. frugiperda and its ecological studies. © 2021 Society of Chemical Industry.

Biosystematic Study on Some Egyptian Species of Astragalus L. (Fabaceae)

Astragalus L. is one of the largest angiosperm complex genera that belongs to the family Fabaceae, subfamily Papilionoideae or Faboideae under the subtribe Astragalinae of the tribe Galegeae. The current study includes the whole plant morphology, DNA barcode (ITS2), and molecular marker (SCoT). Ten taxa representing four species of Astragalus were collected from different localities in Egypt during the period from February 2018 to May 2019. Morphologically, identification and classification of collected Astragalus plants occurred by utilizing the light microscope, regarding the taxonomic revisions of the reference collected Astragalus specimens in other Egyptian Herbaria. For molecular validation, ten SCoT primers were used in this study, producing a unique banding pattern to differentiate between ten samples of Astragalus taxa which generated 212 DNA fragments with an average of 12.2 bands per 10 Astragalus samples, with 8 to 37 fragments per primer. The 212 fragments amplified were distributed as 2 monomorphic bands, 27 polymorphic without unique bands, 183 unique bands (210 Polymorphic with unique bands), and ITS2 gene sequence was showed as the optimal barcode for identifying Astragalus L. using BLAST searched on NCBI database, and afterward, analyzing the chromatogram for ITS region, 10 samples have been identified as two samples representing A. hauarensis, four samples representing A. sieberi, three samples representing A. spinosus and one sample representing A. vogelii. Based on the ITS barcode, A. hauarensis RMG1, A. hauarensis RMG2, A. sieberi RMG1, A. sieberi RMG2, A. sieberi RMG3, A. sieberi RMG4, A. spinosus RMG1, A. spinosus RMG2, A. spinosus RMG3, A. vogelii RMG were deposited into GenBank with accession # MT367587.1, MT367591.1, MT367593.1, MT367585.1, MT367586.1, MT367588.1, MT160347.1, MT367590.1, MT367589.1, MT367592.1, respectively. These results indicated the efficiency of SCoT markers and ITS2 region in identifying and determining genetic relationships between Astragalus species.

Artificial construction of the biocoenosis of deep-sea ecosystem via seeping methane

Deep-sea ecosystems, such as cold seeps and hydrothermal vents, have high biomass, even though they are located in the benthic zone, where no sunlight is present to provide energy for organism proliferation. Based on the coexistence of the reduced gases and chemoautotrophic microbes, it is inferred that the energy from the reduced gases supports the biocoenosis of deep-sea ecosystems. However, there is no direct evidence to support this deduction. Here, we developed and placed a biocoenosis generator, a device that continuously seeped methane, on the 1000-m deep-sea floor of the South China Sea to artificially construct a deep-sea ecosystem biocoenosis. The results showed that microorganisms, including bacteria and archaea, appeared in the biocoenosis generator first, followed by jellyfish and Gammaridea arthropods, indicating that a biocoenosis had been successfully constructed in the deep sea. Anaerobic methane-oxidizing archaea, which shared characteristics with the archaea of natural deep-sea cold seeps, acted as the first electron acceptors of the emitted methane; then, the energy in the electrons was transferred to downstream symbiotic archaea and bacteria and finally to animals. Nitrate-reducing bacteria served as partners to complete anaerobic oxidation of methane process. Further analysis revealed that viruses coexisted with these organisms during the origin of the deep-sea biocoenosis. Therefore, our study mimics a natural deep-sea ecosystem and provides the direct evidence to show that the chemical energy of reduced organic molecules, such as methane, supports the biocoenosis of deep-sea ecosystems.

Development of a Species Diagnostic Molecular Tool for an Invasive Pest, Mythimna loreyi, Using LAMP

The Mythimna loreyi (Duponchel) is one of the well-known invasive noctuid pests in Africa, Australia, and many Asian countries. However, it is difficult to identify the invasive and morphologically similar species, Mythimna separate, which occur at the cornfield in the larvae stage. Currently, the molecular biology method for diagnosing M. loreyi species is only using the mtCO1 universal primer (LCO1490, HCO2198), which requires a lot of time and effort, such as DNA extraction, PCR, electrophoresis, and sequencing. In this study, the LAMP assay was developed for rapid, simple, effective species identification. By analyzing the mitochondrial genome, the species-specific sequence was found at the coding region of the NADH dehydrogenase subunit 5 gene. Based on this unique sequence, four LAMP primers and two loop primers were designed. The F3 and B3 primers were able to diagnose species-specific, in general, and multiplex PCR and specifically reacted within the inner primers in LAMP assay. The optimal incubation condition of the LAMP assay was 61 °C for 60 min with four LAMP primers, though additional loop primers, BF and LF, did not significantly shorten the amplification time. The broad range of DNA concentration was workable in LAMP assay, in which the minimum detectable DNA concentration was 100 pg. DNA releasing method was applied, which took five minutes of incubation at 95 °C without the DNA extraction process. Only some pieces of tissue of larvae and adult samples were needed to extract DNA. The incidence of invasive pests is gradually diversifying. Therefore, this simple and accurate LAMP assay is possibly applied in the intensive field monitoring for invasive pests and integrated management of Mythimna loreyi.

An Evolving View of Phylogenetic Support

If all nucleotide sites evolved at the same rate within molecules and throughout the history of lineages, if all nucleotides were in equal proportion, if any nucleotide or amino acid evolved to any other with equal probability, if all taxa could be sampled, if diversification happened at well-spaced intervals, and if all gene segments had the same history, then tree building would be easy. But of course none of those conditions are true. Hence the need for evaluating the information content and accuracy of phylogenetic trees. The symposium for which this historial essay and presentation were developed focused on the importance of phylogenetic support, specifically branch support for individual clades. Here I present a timeline and review significant events in the history of systematics that set the stage for the development of the sophisticated measures of branch support and examinations of the information content of data highlighted in this symposium.

Comparisons of sampling methods for assessing intra- and inter-accession genetic diversity in three rice species using genotyping by sequencing

To minimize the cost of sample preparation and genotyping, most genebank genomics studies in self-pollinating species are conducted on a single individual to represent an accession, which may be heterogeneous with larger than expected intra-accession genetic variation. Here, we compared various population genetics parameters among six DNA (leaf) sampling methods on 90 accessions representing a wild species (O. barthii), cultivated and landraces (O. glaberrima, O. sativa), and improved varieties derived through interspecific hybridizations. A total of 1,527 DNA samples were genotyped with 46,818 polymorphic single nucleotide polymorphisms (SNPs) using DArTseq. Various statistical analyses were performed on eleven datasets corresponding to 5 plants per accession individually and in a bulk (two sets), 10 plants individually and in a bulk (two sets), all 15 plants individually (one set), and a randomly sampled individual repeated six times (six sets). Overall, we arrived at broadly similar conclusions across 11 datasets in terms of SNP polymorphism, heterozygosity/heterogeneity, diversity indices, concordance among genetic dissimilarity matrices, population structure, and genetic differentiation; there were, however, a few discrepancies between some pairs of datasets. Detailed results of each sampling method, the concordance in their outputs, and the technical and cost implications of each method were discussed.

Identification of the prion-related protein gene (PRNT) sequences in various species of the Cervidae family

Chronic wasting disease (CWD) is caused by abnormal deleterious prion protein (PrP^Sc), and transmissible spongiform encephalopathy occurs in the Cervidae family. In recent studies, the susceptibility of prion disease has been affected by polymorphisms of the prion gene family. However, the study of the prion-related protein gene (PRNT) is rare, and the DNA sequence of this gene was not fully reported in all Cervidae families. In the present study, we amplified and first identified PRNT DNA sequences in the Cervidae family, including red deer, elk, sika deer and Korean water deer, using polymerase chain reaction (PCR). We aligned nucleotide sequences of the PRNT gene and the amino acid sequences of prion-related protein (Prt) protein among several species. In addition, we performed phylogenetic analysis to measure the evolutionary relationships of the PRNT gene in the Cervidae family. Furthermore, we performed homology modeling of the Prt protein using SWISS-MODEL and compared the structure of Prt protein between sheep and the Cervidae family using the Swiss-PdbViewer program. We obtained much longer PRNT sequences of red deer compared to the PRNT gene sequence registered in GenBank. Korean water deer denoted more close evolutionary distances with goats and cattle than the Cervidae family. We found 6 Cervidae family-specific amino acids by the alignment of Prt amino acid sequences. There are significantly different distributions of hydrogen bonds and the atomic distance of the N-terminal tail and C-terminal tail between sheep and the Cervidae family. We also detected the mRNA expression of PRNT gene in 3 tissues investigated. To our knowledge, this report is the first genetic study of the PRNT gene in the Cervidae family.

Genome Analysis of a Novel Clade II.b Alphabaculovirus Obtained from Artaxa digramma

Artaxa digramma is a lepidopteran pest distributed throughout southern China, Myanmar, Indonesia, and India. Artaxa digramma nucleopolyhedrovirus (ArdiNPV) is a specific viral pathogen of A. digramma and deemed as a promising biocontrol agent against the pest. In this study, the complete genome sequence of ArdiNPV was determined by deep sequencing. The genome of ArdiNPV contains a double-stranded DNA (dsDNA) of 161,734 bp in length and 39.1% G+C content. Further, 149 hypothetical open reading frames (ORFs) were predicted to encode proteins >50 amino acids in length, covering 83% of the whole genome. Among these ORFs, 38 were baculovirus core genes, 22 were lepidopteran baculovirus conserved genes, and seven were unique to ArdiNPV, respectively. No typical baculoviral homologous regions (hrs) were identified in the genome. ArdiNPV had five multi-copy genes including baculovirus repeated ORFs (bros), calcium/sodium antiporter B (chaB), DNA binding protein (dbp), inhibitor of apoptosis protein (iap), and p26. Interestingly, phylogenetic analyses showed that ArdiNPV belonged to Clade II.b of Group II Alphabaculoviruses, which all contain a second copy of dbp. The genome of ArdiNPV was the closest to Euproctis pseudoconspersa nucleopolyhedrovirus, with 57.4% whole-genome similarity. Therefore, these results suggest that ArdiNPV is a novel baculovirus belonging to a newly identified cluster of Clade II.b Alphabaculoviruses.

Genome sequence analysis and organization of the Hyphantria cunea granulovirus (HycuGV-Hc1) from Turkey

The fall webworm (Hyphantria cunea) impacts a wide variety of crops and cultivated broadleaf plant species. The pest is native to North America, was introduced to Europe and has since spread further as far as central Asia. Despite several attempts to control its distribution, the pest continues to spread causing damage all over the world. A naturally occurring baculovirus, Hyphantria cunea granulovirus (HycuGV-Hc1), isolated from the larvae of H. cunea in Turkey appears to have a potential as microbial control agent against this pest. In this report we describe the complete genome sequence and organization of the granulovirus isolate (HycuGV-Hc1) that infects the larval stages and compare it to other baculovirus genomes. The HycuGV-Hc1 genome is a circular double-stranded DNA of 114,825 bp in size with a nucleotide distribution of 39.3% G + C. Bioinformatics analysis predicted 132 putative open reading frames of (ORFs) ≥ 150 nucleotides. There are 24 ORFs with unknown function. Seven homologous repeated regions (hrs) and two bro genes (bro-1 and bro-2) were identified in the genome. Comparison to other baculovirus genomes, HycuGV-Hc1 revealed some differences in gene content and organization. Gene parity plots and phylogenetics confirmed that HycuGV-Hc1 is a Betabaculovirus and is closely related to Plutella xylostella granulovirus. This study expands our knowledge on the genetic variation of HycuGV isolates and provides further novel knowledge on the nature of granuloviruses.

Genome analysis of a novel Group I alphabaculovirus obtained from Oxyplax ochracea

Oxyplax ochracea (Moore) is a pest that causes severe damage to a wide range of crops, forests and fruit trees. The complete genome sequence of Oxyplax ochracea nucleopolyhedrovirus (OxocNPV) was determined using a Roche 454 pyrosequencing system. OxocNPV has a double-stranded DNA (dsDNA) genome of 113,971 bp with a G+C content of 31.1%. One hundred and twenty-four putative open reading frames (ORFs) encoding proteins of >50 amino acids in length and with minimal overlapping were predicted, which covered 92% of the whole genome. Six baculoviral typical homologous regions (hrs) were identified. Phylogenetic analysis and gene parity plot analysis showed that OxocNPV belongs to clade “a” of Group I alphabaculoviruses, and it seems to be close to the most recent common ancestor of Group I alphabaculoviruses. Three unique ORFs (with no homologs in the National Center for Biotechnology Information database) were identified. Interestingly, OxocNPV lacks three auxiliary genes (lef7, ie-2 and pcna) related to viral DNA replication and RNA transcription. In addition, OxocNPV has significantly different sequences for several genes (including ie1 and odv-e66) in comparison with those of other baculoviruses. However, three dimensional structure prediction showed that OxocNPV ODV-E66 contain the conserved catalytic residues, implying that it might possess polysaccharide lyase activity as AcMNPV ODV-E66. All these unique features suggest that OxocNPV represents a novel species of the Group I alphabaculovirus lineage.

Why Do Phylogenomic Data Sets Yield Conflicting Trees? Data Type Influences the Avian Tree of Life more than Taxon Sampling

Phylogenomics, the use of large-scale data matrices in phylogenetic analyses, has been viewed as the ultimate solution to the problem of resolving difficult nodes in the tree of life. However, it has become clear that analyses of these large genomic data sets can also result in conflicting estimates of phylogeny. Here, we use the early divergences in Neoaves, the largest clade of extant birds, as a "model system" to understand the basis for incongruence among phylogenomic trees. We were motivated by the observation that trees from two recent avian phylogenomic studies exhibit conflicts. Those studies used different strategies: 1) collecting many characters [$\sim$ 42 mega base pairs (Mbp) of sequence data] from 48 birds, sometimes including only one taxon for each major clade; and 2) collecting fewer characters ($\sim$ 0.4 Mbp) from 198 birds, selected to subdivide long branches. However, the studies also used different data types: the taxon-poor data matrix comprised 68% non-coding sequences whereas coding exons dominated the taxon-rich data matrix. This difference raises the question of whether the primary reason for incongruence is the number of sites, the number of taxa, or the data type. To test among these alternative hypotheses we assembled a novel, large-scale data matrix comprising 90% non-coding sequences from 235 bird species. Although increased taxon sampling appeared to have a positive impact on phylogenetic analyses the most important variable was data type. Indeed, by analyzing different subsets of the taxa in our data matrix we found that increased taxon sampling actually resulted in increased congruence with the tree from the previous taxon-poor study (which had a majority of non-coding data) instead of the taxon-rich study (which largely used coding data). We suggest that the observed differences in the estimates of topology for these studies reflect data-type effects due to violations of the models used in phylogenetic analyses, some of which may be difficult to detect. If incongruence among trees estimated using phylogenomic methods largely reflects problems with model fit developing more "biologically-realistic" models is likely to be critical for efforts to reconstruct the tree of life. [Birds; coding exons; GTR model; model fit; Neoaves; non-coding DNA; phylogenomics; taxon sampling.].

Complete Pathway Elucidation and Heterologous Reconstitution of Rhodiola Salidroside Biosynthesis

Salidroside is a bioactive tyrosine-derived phenolic natural product found in medicinal plants under the Rhodiola genus. In addition to their anti-fatigue and anti-anoxia roles in traditional medicine, Rhodiola total extract and salidroside have also displayed medicinal properties as anti-cardiovascular diseases and anti-cancer agents. The resulting surge in global demand of Rhodiola plants and salidroside has driven some species close to extinction. Here, we report the full elucidation of the Rhodiola salidroside biosynthetic pathway utilizing the first comprehensive transcriptomics and metabolomics datasets for Rhodiola rosea. Unlike the previously proposed pathway involving separate decarboxylation and deamination enzymatic steps from tyrosine to the key intermediate 4-hydroxyphenylacetaldehyde (4-HPAA), Rhodiola contains a pyridoxal phosphate-dependent 4-HPAA synthase that directly converts tyrosine to 4-HPAA. We further identified genes encoding the subsequent 4-HPAA reductase and tyrosol:UDP-glucose 8-O-glucosyltransferase, respectively, to complete salidroside biosynthesis in Rhodiola. We show that heterologous production of salidroside can be achieved in the yeast Saccharomyces cerevisiae as well as the plant Nicotiana benthamiana through transgenic expression of Rhodiola salidroside biosynthetic genes. This study provides new tools for engineering sustainable production of salidroside in heterologous hosts.

Using Minimum Bootstrap support for Splits to Construct Confidence Regions for Trees

Many of the estimated topologies in phylogenetic studies are presented with the bootstrap support for each of the splits in the topology indicated. If phylogenetic estimation is unbiased, high bootstrap support for a split suggests that there is a good deal of certainty that the split actually is present in the tree and low bootstrap support suggests that one or more of the taxa on one side of the estimated split might in reality be located with taxa on the other side. In the latter case the follow-up questions about how many and which of the taxa could reasonably be incorrectly placed as well as where they might alternatively be placed are not addressed through the presented bootstrap support. We present here an algorithm that finds the set of all trees with minimum bootstrap support for their splits greater than some given value. The output is a ranked list of trees, ranked according to the minimum bootstrap supports for splits in the trees. The number of such trees and their topologies provides useful supplementary information in bootstrap analyses about the reasons for low bootstrap support for splits. We also present ways of quantifying low bootstrap support by considering the set of all topologies with minimum bootstrap greater than some quantity as providing a confidence region of topologies. Using a double bootstrap we are able to choose a cutoff so that the set of topologies with minimum bootstrap support for a split greater than that cutoff gives an approximate 95% confidence region. As with bootstrap support one advantage of the methods is that they are generally applicable to the wide variety of phylogenetic estimation methods.

Genomic adaptation to polyphagy and insecticides in a major East Asian noctuid pest

The tobacco cutworm, Spodoptera litura, is among the most widespread and destructive agricultural pests, feeding on over 100 crops throughout tropical and subtropical Asia. By genome sequencing, physical mapping and transcriptome analysis, we found that the gene families encoding receptors for bitter or toxic substances and detoxification enzymes, such as cytochrome P450, carboxylesterase and glutathione-S-transferase, were massively expanded in this polyphagous species, enabling its extraordinary ability to detect and detoxify many plant secondary compounds. Larval exposure to insecticidal toxins induced expression of detoxification genes, and knockdown of representative genes using short interfering RNA (siRNA) reduced larval survival, consistent with their contribution to the insect's natural pesticide tolerance. A population genetics study indicated that this species expanded throughout southeast Asia by migrating along a South India-South China-Japan axis, adapting to wide-ranging ecological conditions with diverse host plants and insecticides, surviving and adapting with the aid of its expanded detoxification systems. The findings of this study will enable the development of new pest management strategies for the control of major agricultural pests such as S. litura.

Genetic basis for soma is present in undifferentiated volvocine green algae

Somatic cellular differentiation plays a critical role in the transition from unicellular to multicellular life, but the evolution of its genetic basis remains poorly understood. By definition, somatic cells do not reproduce to pass on genes and so constitute an extreme form of altruistic behaviour. The volvocine green algae provide an excellent model system to study the evolution of multicellularity and somatic differentiation. In Volvox carteri, somatic cell differentiation is controlled by the regA gene, which is part of a tandem duplication of genes known as the reg cluster. Although previous work found the reg cluster in divergent Volvox species, its origin and distribution in the broader group of volvocine algae has not been known. Here, we show that the reg cluster is present in many species without somatic cells and determine that the genetic basis for soma arose before the phenotype at the origin of the family Volvocaceae approximately 200 million years ago. We hypothesize that the ancestral function was involved in regulating reproduction in response to stress and that this function was later co-opted to produce soma. Determining that the reg cluster was co-opted to control somatic cell development provides insight into how cellular differentiation, and with it greater levels of complexity and individuality, evolves.

Deciphering MCR-2 Colistin Resistance

Antibiotic resistance is a prevalent problem in public health worldwide. In general, the carbapenem β-lactam antibiotics are considered a final resort against lethal infections by multidrug-resistant bacteria. Colistin is a cationic polypeptide antibiotic and acts as the last line of defense for treatment of carbapenem-resistant bacteria. Very recently, a new plasmid-borne colistin resistance gene, mcr-2, was revealed soon after the discovery of the paradigm gene mcr-1, which has disseminated globally. However, the molecular mechanisms for MCR-2 colistin resistance are poorly understood. Here we show a unique transposon unit that facilitates the acquisition and transfer of mcr-2 Evolutionary analyses suggested that both MCR-2 and MCR-1 might be traced to their cousin phosphoethanolamine (PEA) lipid A transferase from a known polymyxin producer, Paenibacillus Transcriptional analyses showed that the level of mcr-2 transcripts is relatively higher than that of mcr-1 Genetic deletions revealed that the transmembrane regions (TM1 and TM2) of both MCR-1 and MCR-2 are critical for their location and function in bacterial periplasm, and domain swapping indicated that the TM2 is more efficient than TM1. Matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS) confirmed that all four MCR proteins (MCR-1, MCR-2, and two chimeric versions [TM1-MCR-2 and TM2-MCR-1]) can catalyze chemical modification of lipid A moiety anchored on lipopolysaccharide (LPS) with the addition of phosphoethanolamine to the phosphate group at the 4' position of the sugar. Structure-guided site-directed mutagenesis defined an essential 6-residue-requiring zinc-binding/catalytic motif for MCR-2 colistin resistance. The results further our mechanistic understanding of transferable colistin resistance, providing clues to improve clinical therapeutics targeting severe infections by MCR-2-containing pathogens.IMPORTANCE Carbapenem and colistin are the last line of refuge in fighting multidrug-resistant Gram-negative pathogens. MCR-2 is a newly emerging variant of the mobilized colistin resistance protein MCR-1, posing a potential challenge to public health. Here we report transfer of the mcr-2 gene by a unique transposal event and its possible origin. Distribution of MCR-2 in bacterial periplasm is proposed to be a prerequisite for its role in the context of biochemistry and the colistin resistance. We also define the genetic requirement of a zinc-binding/catalytic motif for MCR-2 colistin resistance. This represents a glimpse of transferable colistin resistance by MCR-2.

Heterologous biosynthesis and manipulation of crocetin in Saccharomyces cerevisiae

BACKGROUND

Due to excellent performance in antitumor, antioxidation, antihypertension, antiatherosclerotic and antidepressant activities, crocetin, naturally exists in Crocus sativus L., has great potential applications in medical and food fields. Microbial production of crocetin has received increasing concern in recent years. However, only a patent from EVOVA Inc. and a report from Lou et al. have illustrated the feasibility of microbial biosynthesis of crocetin, but there was no specific titer data reported so far. Saccharomyces cerevisiae is generally regarded as food safety and productive host, and manipulation of key enzymes is critical to balance metabolic flux, consequently improve output. Therefore, to promote crocetin production in S. cerevisiae, all the key enzymes, such as CrtZ, CCD and ALD should be engineered combinatorially.

RESULTS

By introduction of heterologous CrtZ and CCD in existing β-carotene producing strain, crocetin biosynthesis was achieved successfully in S. cerevisiae. Compared to culturing at 30 °C, the crocetin production was improved to 223 μg/L at 20 °C. Moreover, an optimal CrtZ/CCD combination and a titer of 351 μg/L crocetin were obtained by combinatorial screening of CrtZs from nine species and four CCDs from Crocus. Then through screening of heterologous ALDs from Bixa orellana (Bix_ALD) and Synechocystis sp. PCC6803 (Syn_ALD) as well as endogenous ALD6, the crocetin titer was further enhanced by 1.8-folds after incorporating Syn_ALD. Finally a highest reported titer of 1219 μg/L at shake flask level was achieved by overexpression of CCD2 and Syn_ALD. Eventually, through fed-batch fermentation, the production of crocetin in 5-L bioreactor reached to 6278 μg/L, which is the highest crocetin titer reported in eukaryotic cell.

CONCLUSIONS

Saccharomyces cerevisiae was engineered to achieve crocetin production in this study. Through combinatorial manipulation of three key enzymes CrtZ, CCD and ALD in terms of screening enzymes sources and regulating protein expression level (reaction temperature and copy number), crocetin titer was stepwise improved by 129.4-fold (from 9.42 to 1219 μg/L) as compared to the starting strain. The highest crocetin titer (6278 μg/L) reported in microbes was achieved in 5-L bioreactors. This study provides a good insight into key enzyme manipulation involved in serial reactions for microbial overproduction of desired compounds with complex structure.

Expression map of a complete set of gustatory receptor genes in chemosensory organs of Bombyx mori

Most lepidopteran species are herbivores, and interaction with host plants affects their gene expression and behavior as well as their genome evolution. Gustatory receptors (Grs) are expected to mediate host plant selection, feeding, oviposition and courtship behavior. However, due to their high diversity, sequence divergence and extremely low level of expression it has been difficult to identify precisely a complete set of Grs in Lepidoptera. By manual annotation and BAC sequencing, we improved annotation of 43 gene sequences compared with previously reported Grs in the most studied lepidopteran model, the silkworm, Bombyx mori, and identified 7 new tandem copies of BmGr30 on chromosome 7, bringing the total number of BmGrs to 76. Among these, we mapped 68 genes to chromosomes in a newly constructed chromosome distribution map and 8 genes to scaffolds; we also found new evidence for large clusters of BmGrs, especially from the bitter receptor family. RNA-seq analysis of diverse BmGr expression patterns in chemosensory organs of larvae and adults enabled us to draw a precise organ specific map of BmGr expression. Interestingly, most of the clustered genes were expressed in the same tissues and more than half of the genes were expressed in larval maxillae, larval thoracic legs and adult legs. For example, BmGr63 showed high expression levels in all organs in both larval and adult stages. By contrast, some genes showed expression limited to specific developmental stages or organs and tissues. BmGr19 was highly expressed in larval chemosensory organs (especially antennae and thoracic legs), the single exon genes BmGr53 and BmGr67 were expressed exclusively in larval tissues, the BmGr27-BmGr31 gene cluster on chr7 displayed a high expression level limited to adult legs and the candidate CO₂ receptor BmGr2 was highly expressed in adult antennae, where few other Grs were expressed. Transcriptional analysis of the Grs in B. mori provides a valuable new reference for finding genes involved in plant-insect interactions in Lepidoptera and establishing correlations between these genes and vital insect behaviors like host plant selection and courtship for mating.

Molecular epidemiology of HIV-1 among the HIV infected people of Manipur, Northeastern India: Emergence of unique recombinant forms

According to the Joint National Programme on HIV/AIDS (UNAIDS), the northeastern region of India has the highest HIV prevalence in the country. This study was conducted to determine the current HIV-1 molecular epidemiology of Manipur, a state in northeast India. Blood samples from HIV-1 seropositive subjects were collected between June 2011 and February 2014. The partial regions of HIV-1 genes; pol and tat-vpu-env were independently amplified, sequenced, analyzed, and genotyped. Based on all sequences generated from 110 samples using pol and/or tat-vpu-env gene, the overall HIV-1 genotypes distribution of Manipur was as follows: 65.45% (72/110) subtype C, 32.73% (36/110) unique recombinant forms (URFs), and 1.82% (2/110) subtype B. The distribution of HIV-1 genotypes among the risk groups was: heterosexual: 58.33% (35/60) subtype C, 38.33% (23/60) URFs, and 3.34% (2/60) subtype B; intravenous drug users (IDUs): 85.36% (35/41) subtype C, 9.76% (4/41) URFs, and 4.88% (2/41) subtype B; mother to child (MTC): 50% (3/6) URFs and 50% (3/6) subtype C and blood transfusion: 100% (3/3) subtype C. The findings for the first time revealed the emergence of URFs of HIV-1 in Manipur which is predominant among the sexual and MTC risk groups as compared to IDUs. Taking together, this study illustrated that Manipur is the "recombinant hotspot of HIV" of India. The results will provide the clinical importance for continuous monitoring of HIV-infections in order to design appropriate prevention measures to limit the spread of new HIV infections.

The global spread of HIV-1 subtype B epidemic

Human immunodeficiency virus type 1 (HIV-1) was discovered in the early 1980s when the virus had already established a pandemic. For at least three decades the epidemic in the Western World has been dominated by subtype B infections, as part of a sub-epidemic that traveled from Africa through Haiti to United States. However, the pattern of the subsequent spread still remains poorly understood. Here we analyze a large dataset of globally representative HIV-1 subtype B strains to map their spread around the world over the last 50years and describe significant spread patterns. We show that subtype B travelled from North America to Western Europe in different occasions, while Central/Eastern Europe remained isolated for the most part of the early epidemic. Looking with more detail in European countries we see that the United Kingdom, France and Switzerland exchanged viral isolates with non-European countries than with European ones. The observed pattern is likely to mirror geopolitical landmarks in the post-World War II era, namely the rise and the fall of the Iron Curtain and the European colonialism. In conclusion, HIV-1 spread through specific migration routes which are consistent with geopolitical factors that affected human activities during the last 50years, such as migration, tourism and trade. Our findings support the argument that epidemic control policies should be global and incorporate political and socioeconomic factors.

Genome Sequencing and Analysis of Catopsilia pomona nucleopolyhedrovirus: A Distinct Species in Group I Alphabaculovirus

The genome sequence of Catopsilia pomona nucleopolyhedrovirus (CapoNPV) was determined by the Roche 454 sequencing system. The genome consisted of 128,058 bp and had an overall G+C content of 40%. There were 130 hypothetical open reading frames (ORFs) potentially encoding proteins of more than 50 amino acids and covering 92% of the genome. Among all the hypothetical ORFs, 37 baculovirus core genes, 23 lepidopteran baculovirus conserved genes and 10 genes conserved in Group I alphabaculoviruses were identified. In addition, the genome included regions of 8 typical baculoviral homologous repeat sequences (hrs). Phylogenic analysis showed that CapoNPV was in a distinct branch of clade “a” in Group I alphabaculoviruses. Gene parity plot analysis and overall similarity of ORFs indicated that CapoNPV is more closely related to the Group I alphabaculoviruses than to other baculoviruses. Interesting, CapoNPV lacks the genes encoding the fibroblast growth factor (fgf) and ac30, which are conserved in most lepidopteran and Group I baculoviruses, respectively. Sequence analysis of the F-like protein of CapoNPV showed that some amino acids were inserted into the fusion peptide region and the pre-transmembrane region of the protein. All these unique features imply that CapoNPV represents a member of a new baculovirus species.

Inhibition of dengue virus replication by diisopropyl chrysin-7-yl phosphate

Dengue fever is a tropical disease and caused by dengue virus (DENV), which is transmitted by mosquitoes and infects about 400 million people annually. With the development of international trade and travel, China is facing a growing threat. Over 40 thousands of people were infected during the 2014 DENV outbreak in Guangdong. Neither licensed vaccine nor therapeutic drug has been available. In this report, we isolated two clinical DENV strains. The full-length genome was sequenced and characterized. We also applied a flavonoid, CPI, into an anti-DENV assay. Replication of viral RNA and expression of viral protein was all strongly inhibited. These results indicated that CPI may serve as potential protective agents in the treatment of patients with chronic DENV infection.

Genetic diversity, population structure and association analysis in cut chrysanthemum (Chrysanthemum morifolium Ramat.)

Characterizing the genetic diversity present in a working set of plant germplasm can contribute to its effective management and genetic improvement. The cut flower chrysanthemum (Chrysanthemum morifolium Ramat.) is an economically important ornamental species. With the repeated germplasm exchange and intensive breeding activities, it remains a major task in genetic research. The purpose of the present study was to characterize the genetic diversity and the population structure of a worldwide collection of 159 varieties, and to apply an association mapping approach to identify DNA-based markers linked to five plant architecture traits and six inflorescence traits. The genotyping demonstrated that there was no lack of genetic diversity in the collection and that pair-wise kinship values were relatively low. The clustering based on a Bayesian model of population structure did not reflect known variation in either provenance or inflorescence type. A principal coordinate analysis was, however, able to discriminate most of the varieties according to both of these criteria. About 1 in 100 marker pairs exhibited a degree of linkage disequilibrium. The association analysis identified a number of markers putatively linked to one or more of the traits. Some of these associations were robust over two seasons. The findings provide an in-depth understanding of genetic diversity and population structure present in cut flower chrysanthemum varieties, and an insight into the genetic control of plant architecture and inflorescence-related traits.

Genome sequencing and analysis of a granulovirus isolated from the Asiatic rice leafroller, Cnaphalocrocis medinalis

The complete genome of Cnaphalocrocis medinalis granulovirus (CnmeGV) from a serious migratory rice pest, Cnaphalocrocis medinalis (Lepidoptera: Pyralidae), was sequenced using the Roche 454 Genome Sequencer FLX system (GS FLX) with shotgun strategy and assembled by Roche GS De Novo assembler software. Its circular double-stranded genome is 111,246 bp in size with a high A+T content of 64.8% and codes for 118 putative open reading frames (ORFs). It contains 37 conserved baculovirus core ORFs, 13 unique ORFs, 26 ORFs that were found in all Lepidoptera baculoviruses and 42 common ORFs. The analysis of nucleotide sequence repeats revealed that the CnmeGV genome differs from the rest of sequenced GVs by a 23 kb and a 17kb gene block inversions, and does not contain any typical homologous region (hr) except for a region of non-hr-like sequence. Chitinase and cathepsin genes, which are reported to have major roles in the liquefaction of the hosts, were not found in the CnmeGV genome, which explains why CnmeGV infected insects do not show the phenotype of typical liquefaction. Phylogenetic analysis, based on the 37 core baculovirus genes, indicates that CnmeGV is closely related to Adoxophyes orana granulovirus. The genome analysis would contribute to the functional research of CnmeGV, and would benefit to the utilization of CnmeGV as pest control reagent for rice production.

Epigenetic Mutation of RAV6 Affects Leaf Angle and Seed Size in Rice

Heritable epigenetic variants of genes, termed epialleles, can broaden genetic and phenotypic diversity in eukaryotes. Epialleles may also provide a new source of beneficial traits for crop breeding, but very few epialleles related to agricultural traits have been identified in crops. Here, we identified Epi-rav6, a gain-of-function epiallele of rice (Oryza sativa) RELATED TO ABSCISIC ACID INSENSITIVE3 (ABI3)/VIVIPAROUS1 (VP1) 6 (RAV6), which encodes a B3 DNA-binding domain-containing protein. The Epi-rav6 plants show larger lamina inclination and smaller grain size; these agronomically important phenotypes are inherited in a semidominant manner. We did not find nucleotide sequence variation of RAV6. Instead, we found hypomethylation in the promoter region of RAV6, which caused ectopic expression of RAV6 in Epi-rav6 plants. Bisulfite analysis revealed that cytosine methylation of four CG and two CNG loci within a continuous 96-bp region plays essential roles in regulating RAV6 expression; this region contains a conserved miniature inverted repeat transposable element transposon insertion in cultivated rice genomes. Overexpression of RAV6 in the wild type phenocopied the Epi-rav6 phenotype. The brassinosteroid (BR) receptor BR INSENSITIVE1 and BR biosynthetic genes EBISU DWARF, DWARF11, and BR-DEFICIENT DWARF1 were ectopically expressed in Epi-rav6 plants. Also, treatment with a BR biosynthesis inhibitor restored the leaf angle defects of Epi-rav6 plants. This indicates that RAV6 affects rice leaf angle by modulating BR homeostasis and demonstrates an essential regulatory role of epigenetic modification on a key gene controlling important agricultural traits. Thus, our work identifies a unique rice epiallele, which may represent a common phenomenon in complex crop genomes.

Quantitative assessment of gene expression network module-validation methods

Validation of pluripotent modules in diverse networks holds enormous potential for systems biology and network pharmacology. An arising challenge is how to assess the accuracy of discovering all potential modules from multi-omic networks and validating their architectural characteristics based on innovative computational methods beyond function enrichment and biological validation. To display the framework progress in this domain, we systematically divided the existing Computational Validation Approaches based on Modular Architecture (CVAMA) into topology-based approaches (TBA) and statistics-based approaches (SBA). We compared the available module validation methods based on 11 gene expression datasets, and partially consistent results in the form of homogeneous models were obtained with each individual approach, whereas discrepant contradictory results were found between TBA and SBA. The TBA of the Zsummary value had a higher Validation Success Ratio (VSR) (51%) and a higher Fluctuation Ratio (FR) (80.92%), whereas the SBA of the approximately unbiased (AU) p-value had a lower VSR (12.3%) and a lower FR (45.84%). The Gray area simulated study revealed a consistent result for these two models and indicated a lower Variation Ratio (VR) (8.10%) of TBA at 6 simulated levels. Despite facing many novel challenges and evidence limitations, CVAMA may offer novel insights into modular networks.

Integrative genomic signatures of hepatocellular carcinoma derived from nonalcoholic Fatty liver disease

Nonalcoholic fatty liver disease (NAFLD) is a risk factor for Hepatocellular carcinoma (HCC), but he transition from NAFLD to HCC is poorly understood. Feature selection algorithms in human and genetically modified mice NAFLD and HCC microarray data were applied to generate signatures of NAFLD progression and HCC differential survival. These signatures were used to study the pathogenesis of NAFLD derived HCC and explore which subtypes of cancers that can be investigated using mouse models. Our findings show that: (I) HNF4 is a common potential transcription factor mediating the transcription of NAFLD progression genes (II) mice HCC derived from NAFLD co-cluster with a less aggressive human HCC subtype of differential prognosis and mixed etiology (III) the HCC survival signature is able to correctly classify 95% of the samples and gives Fgf20 and Tgfb1i1 as the most robust genes for prediction (IV) the expression values of genes composing the signature in an independent human HCC dataset revealed different HCC subtypes showing differences in survival time by a Logrank test. In summary, we present marker signatures for NAFLD derived HCC molecular pathogenesis both at the gene and pathway level.

STBase: one million species trees for comparative biology

Comprehensively sampled phylogenetic trees provide the most compelling foundations for strong inferences in comparative evolutionary biology. Mismatches are common, however, between the taxa for which comparative data are available and the taxa sampled by published phylogenetic analyses. Moreover, many published phylogenies are gene trees, which cannot always be adapted immediately for species level comparisons because of discordance, gene duplication, and other confounding biological processes. A new database, STBase, lets comparative biologists quickly retrieve species level phylogenetic hypotheses in response to a query list of species names. The database consists of 1 million single- and multi-locus data sets, each with a confidence set of 1000 putative species trees, computed from GenBank sequence data for 413,000 eukaryotic taxa. Two bodies of theoretical work are leveraged to aid in the assembly of multi-locus concatenated data sets for species tree construction. First, multiply labeled gene trees are pruned to conflict-free singly-labeled species-level trees that can be combined between loci. Second, impacts of missing data in multi-locus data sets are ameliorated by assembling only decisive data sets. Data sets overlapping with the user's query are ranked using a scheme that depends on user-provided weights for tree quality and for taxonomic overlap of the tree with the query. Retrieval times are independent of the size of the database, typically a few seconds. Tree quality is assessed by a real-time evaluation of bootstrap support on just the overlapping subtree. Associated sequence alignments, tree files and metadata can be downloaded for subsequent analysis. STBase provides a tool for comparative biologists interested in exploiting the most relevant sequence data available for the taxa of interest. It may also serve as a prototype for future species tree oriented databases and as a resource for assembly of larger species phylogenies from precomputed trees.

Molecular epidemiology of HIV-1 in Jilin province, northeastern China: emergence of a new CRF07_BC transmission cluster and intersubtype recombinants

Objective

To investigate the HIV-1 molecular epidemiology among newly diagnosed HIV-1 infected persons living in the Jilin province of northeastern China.

Methods

Plasma samples from 189 newly diagnosed HIV-1 infected patients were collected between June 2010 and August 2011 from all nine cities of Jilin province. HIV-1 nucleotide sequences of gag P17–P24 and env C2–C4 gene regions were amplified using a multiplex RT-PCR method and sequenced. Phylogenetic and recombination analyses were used to determine the HIV-1 genotypes.

Results

Based on all sequences generated, the subtype/CFR distribution was as follows: CRF01_AE (58.1%), CRF07_BC (13.2%), subtype B’ (13.2%), recombinant viruses (8.1%), subtype B (3.7%), CRF02_AG (2.9%), subtype C (0.7%). In addition to finding CRF01_AE strains from previously reported transmission clusters 1, 4 and 5, a new transmission cluster was described within the CRF07_BC radiation. Among 11 different recombinants identified, 10 contained portions of gene regions from the CRF01_AE lineage. CRF02_AG was found to form a transmission cluster of 4 in local Jilin residents.

Conclusions

Our study presents a molecular epidemiologic investigation describing the complex structure of HIV-1 strains co-circulating in Jilin province. The results highlight the critical importance of continuous monitoring of HIV-infections, along with detailed socio-demographic data, in order to design appropriate prevention measures to limit the spread of new HIV infections.

Why some clades have low bootstrap frequencies and high Bayesian posterior probabilities

Bayesian posterior probabilities are wrongly considered by many systematists as indicative of character support, and equivalent to non-parametric bootstrap frequencies. Here I argue against this view. Non-parametric bootstrap is indicative of the amount of evidence in a data matrix supporting each clade in the tree, while Bayesian posterior probabilities are not intended to represent that property. Clades with high posterior probability may not have a large amount of characters favouring them, and their frequencies are the result of the particular sampling procedure of the Bayesian Markov chain Monte Carlo method, which tends to sample very similar topologies according to their posterior probabilities. Both metrics may relate to the notion of confidence, but depict different properties.

The phylogeography of Eurasian Fraxinus species reveals ancient transcontinental reticulation

To investigate the biogeographical history of ashes species of the Eurasian section Fraxinus and to test the hypothesis of ancient reticulations, we sequenced nuclear DNA (nETS and nITS, 1075 bp) for 533 samples and scored AFLP for 63 samples of Eurasian ashes within the section Fraxinus. The nITS phylogeny retrieved the classical view of the evolution of the section, whereas nETS phylogeny indicated an unexpected separation of F. angustifolia in two paraphyletic groups, respectively found in southeastern Europe and in the other parts of the Mediterranean basin. In the nETS phylogeny, the former group was closely related to F. excelsior, whereas the later was closely related to F. mandshurica, a species which is restricted nowadays to northeastern Asia. This topological incongruence between the two loci indicated the occurrence of an ancient reticulation between European and Asian ash species. Several other ancient reticulation events between the two European species and the other species of the section were supported by the posterior predictive checking method. Some of these reticulation events would have occurred during the Miocene, when climatic variations may have lead these species to expand their distribution range and come into contact. The recurrent reticulations observed among Eurasian ash species indicate that they should be considered as conspecific taxa, with subspecific status for some groups. Altogether, the results of the present study provide a rare documented evidence for the occurrence of multiple ancient reticulations within a group of temperate tree taxa with modern disjunct distributions in Eurasia. These ancient reticulation events indicate that the speciation process is slow in ashes, necessitating long periods of geographical isolation. The implications for speciation processes in temperate trees with similar life history and reproductive biology are discussed.

Does more sequence data improve estimates of galliform phylogeny? Analyses of a rapid radiation using a complete data matrix

The resolution of rapid evolutionary radiations or "bushes" in the tree of life has been one of the most difficult and interesting problems in phylogenetics. The avian order Galliformes appears to have undergone several rapid radiations that have limited the resolution of prior studies and obscured the position of taxa important both agriculturally and as model systems (chicken, turkey, Japanese quail). Here we present analyses of a multi-locus data matrix comprising over 15,000 sites, primarily from nuclear introns but also including three mitochondrial regions, from 46 galliform taxa with all gene regions sampled for all taxa. The increased sampling of unlinked nuclear genes provided strong bootstrap support for all but a small number of relationships. Coalescent-based methods to combine individual gene trees and analyses of datasets that are independent of published data indicated that this well-supported topology is likely to reflect the galliform species tree. The inclusion or exclusion of mitochondrial data had a limited impact upon analyses upon analyses using either concatenated data or multispecies coalescent methods. Some of the key phylogenetic findings include support for a second major clade within the core phasianids that includes the chicken and Japanese quail and clarification of the phylogenetic relationships of turkey. Jackknifed datasets suggested that there is an advantage to sampling many independent regions across the genome rather than obtaining long sequences for a small number of loci, possibly reflecting the differences among gene trees that differ due to incomplete lineage sorting. Despite the novel insights we obtained using this increased sampling of gene regions, some nodes remain unresolved, likely due to periods of rapid diversification. Resolving these remaining groups will likely require sequencing a very large number of gene regions, but our analyses now appear to support a robust backbone for this order.

Genome sequence and analysis of Buzura suppressaria nucleopolyhedrovirus: a group II Alphabaculovirus

The genome of Buzura suppressaria nucleopolyhedrovirus (BusuNPV) was sequenced by 454 pyrosequencing technology. The size of the genome is 120,420 bp with 36.8% G+C content. It contains 127 hypothetical open reading frames (ORFs) covering 90.7% of the genome and includes the 37 conserved baculovirus core genes, 84 genes found in other baculoviruses, and 6 unique ORFs. No typical baculoviral homologous repeats (hrs) were present but the genome contained a region of repeated sequences. Gene Parity Plots revealed a 28.8 kb region conserved among the alpha- and beta-baculoviruses. Overall comparisons of BusuNPV to other baculoviruses point to a distinct species in group II Alphabaculovirus.

Functional characterization of protein 4.1 homolog in amphioxus: defining a cryptic spectrin-actin-binding site

Vertebrate 4.1 proteins have a spectrin-actin-binding (SAB) domain, which is lacking in all the invertebrate 4.1 proteins indentified so far, and it was therefore proposed that the SAB domain emerged with the advent of vertebrates during evolution. Here we demonstrated for the first time that amphioxus (an invertebrate chordate) protein 4.1, though lacking a recognizable SAB, was able to bind both spectrin and actin, with a binding capacity comparable to that of human protein 4.1. Detailed structure-activity analyses revealed that the unique domain U2/3 was a newly identified SAB-like domain capable of interacting with spectrin and actin, suggesting the presence of a “cryptic” SAB domain in amphioxus 4.1 protein. We also showed that amphioxus 4.1 protein gene was the common ancestor of vertebrate 4.1 protein genes, from which 4.1R, 4.1N, 4.1G, and 4.1B genes originated. This work will encourage further study on the structure-activity of invertebrate 4.1 protein and its interacting proteins.

Assessing statistical reliability of phylogenetic trees via a speedy double bootstrap method

Evaluating the reliability of estimated phylogenetic trees is of critical importance in the field of molecular phylogenetics, and for other endeavors that depend on accurate phylogenetic reconstruction. The bootstrap method is a well-known computational approach to phylogenetic tree assessment, and more generally for assessing the reliability of statistical models. However, it is known to be biased under certain circumstances, calling into question the accuracy of the method. Several advanced bootstrap methods have been developed to achieve higher accuracy, one of which is the double bootstrap approach, but the computational burden of this method has precluded its application to practical problems of phylogenetic tree selection. We address this issue by proposing a simple method called the speedy double bootstrap, which circumvents the second-tier resampling step in the regular double bootstrap approach. We also develop an implementation of the regular double bootstrap for comparison with our speedy method. The speedy double bootstrap suffers no significant loss of accuracy compared with the regular double bootstrap, while performing calculations significantly more rapidly (at minimum around 371 times faster, based on analysis of mammalian mitochondrial amino acid sequences and 12S and 16S rRNA genes). Our method thus enables, for the first time, the practical application of the double bootstrap technique in the context of molecular phylogenetics. The approach can also be used more generally for model selection problems wherever the maximum likelihood criterion is used.

A time lag insensitive approach for estimating HIV-1 transmission direction

OBJECTIVES

Identifying the direction of transmission in transmission pairs is important both for forensic investigations and for the monitoring of HIV epidemics, however, reliable methods are not yet available due to the long time lag between infection and sampling in most real cases.

DESIGNS

Based on bottleneck effect and coreceptor switching, we aimed at identifying an estimator from sequences of viral gp120 proteins to determine transmission direction between transmission pairs. The estimator should be changed with HIV transmission but was independent of disease progression in an individual.

METHODS

Here, we present a novel and reliable approach for identifying transmission direction. We derived a set of conserved patterns, called common patterns, from the sequences of viruses, which differed in their coreceptor usage. The number of unique common patterns in viral sequences decreased with transmission but remained almost constant with the progress of disease in an individual. We used this number as an estimator to determine transmission direction in 73 transmission pairs for which the transmission direction was already known.

RESULTS

Our method predicted transmission direction with an accuracy of up to 94.5%. Of greater importance, our approach was not influenced by time lags between infection and sampling, and even transmission direction for transmission pairs with long time lags ranging from 2 years to more than 18 years were correctly determined.

CONCLUSION

Our approach for accurately determining transmission direction between transmission pairs is irrespective of the time lag between infection and sampling, which means a promising applications prospect.

Radical instability and spurious branch support by likelihood when applied to matrices with non-random distributions of missing data

Non-random distributions of missing data are a general problem for likelihood-based statistical analyses, including those in a phylogenetic context. Extensive non-randomly distributed missing data are particularly problematic in supermatrix analyses that include many terminals and/or loci. It has been widely reported that missing data can lead to loss of resolution, but only very rarely create misleading or otherwise unsupported results in a parsimony context. Yet this does not hold for all parametric-based analyses because of their assumption of homogeneity across characters and lineages, which can lead to both long-branch attraction and long-branch repulsion. Contrived examples were used to demonstrate that non-random distributions of missing data, even without rate heterogeneity among characters and a well fitting model, can provide misleading likelihood-based topologies and branch-support values that are radically unstable based on slight modifications to character sampling. The same can occur despite complete absence of parsimony-informative characters. Otherwise unsupported resolution and high branch support for these clades were found to occur frequently in 22 empirical examples derived from a published supermatrix. Partitioning characters based on the distribution of missing data helped to decrease, but did not eliminate, these artifacts. These artifacts were exacerbated by low quality tree searches, particularly when holding only a single optimal tree that must be fully resolved.

Towards resolving the complete fern tree of life

In the past two decades, molecular systematic studies have revolutionized our understanding of the evolutionary history of ferns. The availability of large molecular data sets together with efficient computer algorithms, now enables us to reconstruct evolutionary histories with previously unseen completeness. Here, the most comprehensive fern phylogeny to date, representing over one-fifth of the extant global fern diversity, is inferred based on four plastid genes. Parsimony and maximum-likelihood analyses provided a mostly congruent results and in general supported the prevailing view on the higher-level fern systematics. At a deep phylogenetic level, the position of horsetails depended on the optimality criteria chosen, with horsetails positioned as the sister group either of Marattiopsida-Polypodiopsida clade or of the Polypodiopsida. The analyses demonstrate the power of using a 'supermatrix' approach to resolve large-scale phylogenies and reveal questionable taxonomies. These results provide a valuable background for future research on fern systematics, ecology, biogeography and other evolutionary studies.