MEGA3: Integrated software for Molecular Evolutionary Genetics Analysis and sequence alignment

With its theoretical basis firmly established in molecular evolutionary and population genetics, the comparative DNA and protein sequence analysis plays a central role in reconstructing the evolutionary histories of species and multigene families, estimating rates of molecular evolution, and inferring the nature and extent of selective forces shaping the evolution of genes and genomes. The scope of these investigations has now expanded greatly owing to the development of high-throughput sequencing techniques and novel statistical and computational methods. These methods require easy-to-use computer programs. One such effort has been to produce Molecular Evolutionary Genetics Analysis (MEGA) software, with its focus on facilitating the exploration and analysis of the DNA and protein sequence variation from an evolutionary perspective. Currently in its third major release, MEGA3 contains facilities for automatic and manual sequence alignment, web-based mining of databases, inference of the phylogenetic trees, estimation of evolutionary distances and testing evolutionary hypotheses. This paper provides an overview of the statistical methods, computational tools, and visual exploration modules for data input and the results obtainable in MEGA.     

MEGA11: Molecular Evolutionary Genetics Analysis Version 11

The Molecular Evolutionary Genetics Analysis (MEGA) software has matured to contain a large collection of methods and tools of computational molecular evolution. Here, we describe new additions that make MEGA a more comprehensive tool for building timetrees of species, pathogens, and gene families using rapid relaxed-clock methods. Methods for estimating divergence times and confidence intervals are implemented to use probability densities for calibration constraints for node-dating and sequence sampling dates for tip-dating analyses. They are supported by new options for tagging sequences with spatiotemporal sampling information, an expanded interactive Node Calibrations Editor, and an extended Tree Explorer to display timetrees. Also added is a Bayesian method for estimating neutral evolutionary probabilities of alleles in a species using multispecies sequence alignments and a machine learning method to test for the autocorrelation of evolutionary rates in phylogenies. The computer memory requirements for the maximum likelihood analysis are reduced significantly through reprogramming, and the graphical user interface has been made more responsive and interactive for very big data sets. These enhancements will improve the user experience, quality of results, and the pace of biological discovery. Natively compiled graphical user interface and command-line versions of MEGA11 are available for Microsoft Windows, Linux, and macOS from www.megasoftware.net.     

The Evolutionary Genetics of the Genes Underlying Phenotypic Associations for Loblolly Pine (Pinus taeda,Pinaceae)

A primary goal of evolutionary genetics is to discover and explain the genetic basis of fitness-related traits and how this genetic basis evolves within natural populations. Unprecedented technological advances have fueled the discovery of genetic variants associated with ecologically relevant phenotypes in many different life forms, as well as the ability to scan genomes for deviations from selectively neutral models of evolution. Theoretically, the degree of overlap between lists of genomic regions identified using each approach is related to the genetic architecture of fitness-related traits and the strength and type of natural selection molding variation at these traits within natural populations. Here we address for the first time in a plant the degree of overlap between these lists, using patterns of nucleotide diversity and divergence for >7000 unique amplicons described from the extensive expressed sequence tag libraries generated for loblolly pine (Pinus taeda L.) in combination with the >1000 published genetic associations. We show that loci associated with phenotypic traits are distinct with regard to neutral expectations. Phenotypes measured at the whole plant level (e.g., disease resistance) exhibit an approximately twofold increase in the proportion of adaptive nonsynonymous substitutions over the genome-wide average. As expected for polygenic traits, these signals were apparent only when loci were considered at the level of functional sets. The ramifications of this result are discussed in light of the continued efforts to dissect the genetic basis of quantitative traits.     

Evolutionary transitions among dioecy,androdioecy and hermaphroditism in limnadiid clam shrimp (Branchiopoda: Spinicaudata)

Examinations of breeding system transitions have primarily concentrated on the transition from hermaphroditism to dioecy, likely because of the preponderance of this transition within flowering plants. Fewer studies have considered the reverse transition: dioecy to hermaphroditism. A fruitful approach to studying this latter transition can be sought by studying clades in which transitions between dioecy and hermaphroditism have occurred multiple times. Freshwater crustaceans in the family Limnadiidae comprise dioecious, hermaphroditic and androdioecious (males + hermaphrodites) species, and thus this family represents an excellent model system for the assessment of the evolutionary transitions between these related breeding systems. Herein we report a phylogenetic assessment of breeding system transitions within the family using a total evidence comparative approach. We find that dioecy is the ancestral breeding system for the Limnadiidae and that a minimum of two independent transitions from dioecy to hermaphroditism occurred within this family, leading to (1) a Holarctic, all‐hermaphrodite species, Limnadia lenticularis and (2) mixtures of hermaphrodites and males in the genus Eulimnadia. Both hermaphroditic derivatives are essentially females with only a small amount of energy allocated to male function. Within Eulimnadia, we find several all‐hermaphrodite populations/species that have been independently derived at least twice from androdioecious progenitors within this genus. We discuss two adaptive (based on the notion of ‘reproductive assurance’) and one nonadaptive explanations for the derivation of all‐hermaphroditism from androdioecy. We propose that L. lenticularis likely represents an all‐hermaphrodite species that was derived from an androdioecious ancestor, much like the all‐hermaphrodite populations derived from androdioecy currently observed within the Eulimnadia. Finally, we note that the proposed hypotheses for the dioecy to hermaphroditism transition are unable to explain the derivation of a fully functional, outcrossing hermaphroditic species from a dioecious progenitor.     

鼠疫耶尔森氏菌的进化研究:从系统发育学到系统发育基因组学

20世纪90年代末起,基因组学在细菌研究中应用越来越广泛,尤其在进化领域,取得了一系列革命性的发现.本文以鼠疫耶尔森氏菌进化研究为例,介绍了从利用基因组中少数特定片段(等位基因)多态性进行分析的传统系统发育学,到基于大量菌株全基因组序列进行系统发育基因组学的研究发展历程,回顾讨论了基因组学技术的进步为鼠疫菌进化研究领域带来的成果.     

Spatial analysis method (sam): a software tool combining molecular and environmental data to identify candidate loci for selection

sam is a Windows program designed to detect candidate loci for selection in whole-genome scans. It also gives valuable clues as regards the ecological factors at stake in the selection process. The method used is based on multiple univariate logistic regression models to test for association between allelic frequencies at marker loci and environmental variables. The software reads matrices constituted of presence/absence of molecular markers, and of the corresponding environmental parameters at sampling locations. It provides dynamic analysis tables to process the results. The tool is freely available for download at http://www.econogene.eu/software/sam/.     

A chromosome‐level genome assembly of the woolly apple aphid,Eriosoma lanigerum Hausmann (Hemiptera: Aphididae)

Woolly apple aphid (WAA, Eriosoma lanigerum Hausmann) (Hemiptera: Aphididae) is a major pest of apple trees (Malus domestica, order Rosales) and is critical to the economics of the apple industry in most parts of the world. Here, we generated a chromosome‐level genome assembly of WAA—representing the first genome sequence from the aphid subfamily Eriosomatinae—using a combination of 10X Genomics linked‐reads and in vivo Hi‐C data. The final genome assembly is 327 Mb, with 91% of the assembled sequences anchored into six chromosomes. The contig and scaffold N50 values are 158 kb and 71 Mb, respectively, and we predicted a total of 28,186 protein‐coding genes. The assembly is highly complete, including 97% of conserved arthropod single‐copy orthologues based on Benchmarking Universal Single‐Copy Orthologs (busco ) analysis. Phylogenomic analysis of WAA and nine previously published aphid genomes, spanning four aphid tribes and three subfamilies, reveals that the tribe Eriosomatini (represented by WAA) is recovered as a sister group to Aphidini + Macrosiphini (subfamily Aphidinae). We identified syntenic blocks of genes between our WAA assembly and the genomes of other aphid species and find that two WAA chromosomes (El5 and El6) map to the conserved Macrosiphini and Aphidini X chromosome. Our high‐quality WAA genome assembly and annotation provides a valuable resource for research in a broad range of areas such as comparative and population genomics, insect–plant interactions and pest resistance management.     

Molecular and Evolutionary Analysis of NEAr-Iron Transporter (NEAT) Domains

Iron is essential for bacterial survival, being required for numerous biological processes. NEAr-iron Transporter (NEAT) domains have been studied in pathogenic Gram-positive bacteria to understand how their proteins obtain heme as an iron source during infection. While a 2002 study initially discovered and annotated the NEAT domain encoded by the genomes of several Gram-positive bacteria, there remains a scarcity of information regarding the conservation and distribution of NEAT domains throughout the bacterial kingdom, and whether these domains are restricted to pathogenic bacteria. This study aims to expand upon initial bioinformatics analysis of predicted NEAT domains, by exploring their evolution and conserved function. This information was used to identify new candidate domains in both pathogenic and nonpathogenic organisms. We also searched metagenomic datasets, specifically sequence from the Human Microbiome Project. Here, we report a comprehensive phylogenetic analysis of 343 NEAT domains, encoded by Gram-positive bacteria, mostly within the phylum Firmicutes, with the exception of Eggerthella sp. (Actinobacteria) and an unclassified Mollicutes bacterium (Tenericutes). No new NEAT sequences were identified in the HMP dataset. We detected specific groups of NEAT domains based on phylogeny of protein sequences, including a cluster of novel clostridial NEAT domains. We also identified environmental and soil organisms that encode putative NEAT proteins. Biochemical analysis of heme binding by a NEAT domain from a protein encoded by the soil-dwelling organism Paenibacillus polymyxa demonstrated that the domain is homologous in function to NEAT domains encoded by pathogenic bacteria. Together, this study provides the first global bioinformatics analysis and phylogenetic evidence that NEAT domains have a strong conservation of function, despite group-specific differences at the amino acid level. These findings will provide information useful for future projects concerning the structure and function of NEAT domains, particularly in pathogens where they have yet to be studied.     

Ultraconserved element (UCE) probe set design: Base genome and initial design parameters critical for optimization

Targeted capture and enrichment approaches have proven effective for phylogenetic study. Ultraconserved elements (UCEs) in particular have exhibited great utility for phylogenomic analyses, with the software package phyluce being among the most utilized pipelines for UCE phylogenomics, including probe design. Despite the success of UCEs, it is becoming increasing apparent that diverse lineages require probe sets tailored to focal taxa in order to improve locus recovery. However, factors affecting probe design and methods for optimizing probe sets to focal taxa remain underexplored. Here, we use newly available beetle (Coleoptera) genomic resources to investigate factors affecting UCE probe set design using phyluce . In particular, we explore the effects of stringency during initial design steps, as well as base genome choice on resulting probe sets and locus recovery. We found that both base genome choice and initial bait design stringency parameters greatly alter the number of resultant probes included in final probe sets and strongly affect the number of loci detected and recovered during in silico testing of these probe sets. In addition, we identify attributes of base genomes that correlated with high performance in probe design. Ultimately, we provide a recommended workflow for using Phyluce to design an optimized UCE probe set that will work across a targeted lineage, and use our findings to develop a new, open‐source UCE probe set for beetles of the suborder Adephaga.     

Molecular genetic identification of southern hemisphere beaked whales (Cetacea: Ziphiidae)

To assist in the species-level identification of stranded and hunted beaked whales, we compiled a database of 'reference' sequences from the mitochondrial DNA control region for 15 of the 20 described ziphiid species. Reference samples for eight species were obtained from stranded animals in New Zealand and South Australia. Sequences for a further seven species were obtained from a previously published report. This database was used to identify 20 'test' samples obtained from incompletely documented strandings around New Zealand. Analyses showed that four of these 'test' specimens (20%) had initially been misidentified. These included two animals of particular interest: (i) a Blainville's beaked whale ( Mesoplodon densirostris) , the first record of this species in New Zealand waters; and, (ii) a juvenile Andrews' beaked whale ( Mesoplodon bowdoini ), a species known from just over 20 strandings worldwide. A published sequence from a beaked whale product purchased in the Republic of Korea was identified as a Cuvier's beaked whale ( Ziphius cavirostris ). Levels of intra- and interspecific variation were compared to determine the potential for misidentification when the database or taxonomy is incomplete. Intraspecific variation was generally <2%, and interspecific divergence was generally >4.7%. Exceptions were within-species variation in Hyperoodon planifrons , southern bottlenosed whale (4.12%), which exceeded the variation between the two species of Berardius (3.78%), and variation between the two specimens assigned to M. hectori , Hector's beaked whale (7.14%). The latter case appears to be an error in species identification, and could represent the discovery of a new species of beaked whale.     

Molecular techniques in phytoplankton research: from allozyme electrophoresis to genomics

Molecular techniques have become a valuable tool in phytoplankton studies over the past decades. The appropriate choice of a technique from an increasing array of methods can be rather complex, because different techniques are suitable for different questions or problems in ecology and evolution. Each technique has its particular strengths and weaknesses and is based upon different (theoretical) assumptions. Our aim is to give a better insight in the (correct) use of various molecular techniques in phytoplankton research, with special emphasis on the fields of strain identification, differentiation of populations and the establishment of phylogenetic relationships. The basic steps in the development of molecular techniques like allozyme electrophoresis, RFLP, DGGE, SSCP, RAPD, AFLP and microsatellites, and the application of these techniques in phytoplankton research, are discussed. Furthermore, recent developments in molecular biology, that have so far only found limited application in phytoplankton studies, such as single-cell PCR, PCR assays combined with molecular probes (Heteroduplex Mobility Assays or DNA arrays), Real-time PCR, complete genome sequencing, multi-gene expression studies using microarrays, and Single Nucleotide Polymorphism (SNPs), are discussed. We emphasise the relevance of fundamental and applied molecular studies on phytoplankton for a wider community of ecologists and evolutionary biologists.     

杜鹃花科白珠树属分子系统学和生物地理学研究进展

白珠树属(Gaultheria)在杜鹃花科(Ericaceae)系统演化中占有十分重要的地位,其系统位置和演化关系一直备受争议.最近的分子系统学研究认为,白珠树属已经不再属于传统上的越橘亚科(Vaccinioideae)的綟木族(Andromedeae),而是与一些相关属组成了白珠树族(Gaultherieae).对白珠树属产于美洲的类群和相关类群的分子系统学的初步研究则表明,该属与Diplycosia、Tepuia和Pernettya等属(均为"常绿类群")关系密切,可能应将这几个属并入到白珠树属中,但其属下分类系统关系还需要对产于亚洲的类群进行深入的研究后才能确定.白珠树属与其近缘属的进化历史和生物地理学关系较为复杂,与杜鹃花科其他大多数属不同,白珠树属为典型的环太平洋分布.关于白珠树属的起源问题存在两种不同的推测:一种观点认为该属起源于南半球的冈瓦纳古陆;另一种观点则认为其起源于北半球的劳亚古大陆.本文概述了近年来白珠树属的分子系统学和生物地理学研究进展,并对该属尚存的一些问题进行了讨论.     

Molecular phylogenetics of Triculine snails (Gastropoda: Pomatiopsidae) from southern China

Triculine (Gastropoda: Rissooidea: Pomatiopsidae) snails are involved in the transmission of schistosomiasis and paragonimiasis; their distributions are mainly across southeastern Asia and southern China. In the present investigation, partial sequences of COI, 16S, and 28S were examined to infer the phylogenetic relationships among the species rich and poorly understood gastropod. Samples were collected from 12 geographic locations in six provinces of southern China. Several methods such as maximum parsimony, maximum likelihood and distance analysis were used in phylogenetic analyses among these taxa. The resultant phylogenetic trees showed a similar topology irrespective of the phylogenetic methods used. The taxa fell into two clades, with those from Fujian, Guangxi, and Zhejiang Provinces in one clade and those from Hunan, Sichuan and Hubei in the other. Among the taxa in Hubei Province, five formed a monophyletic clade, but Tricula sp. H-SHY fell into a sister clade of Tricula hortensis of Sichuan, whilst Tricula hongshanensis formed a single clade. Sister taxa Tricula pingi and Tricula hsiangi formed well-supported clade within almost all the trees. These results, while preliminary, represent the first attempt to reconstruct a phylogeny for Triculinae across China.     

A Molecular Phylogeny of Costaceae (Zingiberales)

The phylogenetic relationships of Costaceae, a tropical monocotyledonous family sister to the gingers (Zingiberaceae), were investigated with a combination of two chloroplast loci (the trnL-F locus, including the trnL intron, the 3'trnL exon, and the trnL-F intergenic spacer, and the trnK locus, including the trnK intron and the matK coding region) and one nuclear locus (ITS1-5.8s-ITS2). The resulting parsimony analysis of selected taxa that demonstrate the range of floral morphological variation in the family shows that the Cadalvena-type [corrected] floral morphology is ancestral to the group and that both Tapeinochilos species and a Monocostus + Dimerocostus clade represent recent divergences. The genus Costus is broadly paraphyletic but Costus subgenus Eucostus K. Schum. represents a large monophyletic radiation that is poorly resolved. Within this clade, secondary analyses suggest that pollination syndrome, traditionally used for taxonomic and classification purposes within the genus Costus, is a relatively plastic trait of limited phylogenetic utility. This represents the first detailed investigation into intrageneric and interspecific evolutionary relationships within the family Costaceae and presents some novel evolutionary trends with respect to floral morphology and biogeography.     

Inferring the Evolutionary History of Your Favorite Protein: A Guide for Molecular Biologists

Comparative genomics has proven a fruitful approach to acquire many functional and evolutionary insights into core cellular processes. Here it is argued that in order to perform accurate and interesting comparative genomics, one first and foremost has to be able to recognize, postulate, and revise different evolutionary scenarios. After all, these studies lack a simple protocol, due to different proteins having different evolutionary dynamics and demanding different approaches. The authors here discuss this challenge from a practical (what are the observations?) and conceptual (how do these indicate a specific evolutionary scenario?) viewpoint, with the aim to guide investigators who want to analyze the evolution of their protein(s) of interest. By sharing how the authors draft, test, and update such a scenario and how it directs their investigations, the authors hope to illuminate how to execute molecular evolution studies and how to interpret them. Also see the video abstract here https://youtu.be/VCt3l2pbdbQ .     

一个实用的群体遗传学分析软件包——GENEPOP 3.1版

ＧＥＮＥＰＯＰ是一个非常实用的群体遗传学分析软件包,适用于对大量的群体遗传学数据进行分析。它主要有以下３个方面的用途：１）进行正合检验,如对哈迪－温伯格平衡、种群差异和位点间的连锁不平衡进行检验;２）估算经典的群体遗传学参数,如Ｆｓｔ和其它相关指数及基因频率等;３）可把ＧＥＮＥＰＯＰ的文件转换为常用的群体遗传学分析软件包（如ＢＩＳＹＳ、ＦＳＴＡＴ和ＬＩＮＫＤＯＳ）所要求的输入文件格式。与软件ＢＩＯ