Assessment of methods for amino acid matrix selection and their use on empirical data shows that ad hoc assumptions for choice of matrix are not justified

Chagas vectors Panstrongylus chinai (Del Ponte, 1929) and Panstrongylus howardi (Neiva, 1911): chromatic forms or true species?

Background

Chagas disease is a parasitic infection transmitted by “kissing bugs” (Hemiptera: Reduviidae: Triatominae) that has a huge economic impact in Latin American countries. The vector species with the upmost epidemiological importance in Ecuador are Rhodnius ecuadoriensis (Lent & Leon, 1958) and Triatoma dimidiata (Latreille, 1811). However, other species such as Panstrongylus howardi (Neiva, 1911) and Panstrongylus chinai (Del Ponte, 1929) act as secondary vectors due to their growing adaptation to domestic structures and their ability to transmit the parasite to humans. The latter two taxa are distributed in two different regions, they are allopatric and differ mainly by their general color. Their relative morphological similarity led some authors to suspect that P. chinai is a melanic form of P. howardi.

Methods

The present study explored this question using different approaches: antennal phenotype; geometric morphometrics of heads, wings and eggs; cytogenetics; molecular genetics; experimental crosses; and ecological niche modeling.

Results

The antennal morphology, geometric morphometrics of head and wing shape and cytogenetic analysis were unable to show distinct differences between the two taxa. However, geometric morphometrics of the eggs, molecular genetics, ecological niche modeling and experimental crosses including chromosomal analyses of the F1 hybrids, in addition to their coloration and current distribution support the hypothesis that P. chinai and P. howardi are separate species.

Conclusions

Based on the evidence provided here, P. howardi and P. chinai should not be synonymized. They represent two valid, closely related species.

Characterization of class B scavenger receptor type 1 (SRB1) in turbot (Scophthalmus maximus L.)

Class B scavenger receptor type 1 (SRB1) serves as a high-density lipoprotein (HDL) receptor essential for HDL metabolism, and plays vital roles in innate immunity. In this study, the turbot (Scophthalmus maximus) SRB1 was cloned and characterized. The gene structure consists of a coding region of 1,527 bp nucleotides dividing into 13 exons and 12 introns. Such genome structure is highly conserved among teleost fishes. The deduced SRB1 encodes 508 amino acids that mainly has a CD36 transmembrane domain. Tissue distribution of SRB1 showed the lowest expression in liver, while the highest expression was found in intestine. Significantly down-regulation pattern of SmSRB1 expression in intestine was shared after infection with Vibrio anguillarum and Streptococcus iniae. Brach and site models in CODEML program showed that SmSRB1 underwent a conservative evolutionary and three potential positive selected sites 470K, 496E, and 501Y were detected, which requires further investigation and confirmation using base-editing technologies. Subcellular localization demonstrated that turbot SRB1 was distributed in the membrane and cytoplasm. rSmSRB1 showed binding ability in vitro to bacteria, LPS, PGN, LTA and virus. Protein-protein interaction network agrees the function of SRB1 as lipoprotein receptor. Our results indicated SmSRB1 might act as co-receptors to TLRs and NLRs to modulate the immune response to pathogens. Further studies should pay attention to evaluate the specific co-receptor for SRB1 in recognition of different pathogens and selective mechanisms involved.

Transcriptomic Profiling Reveals Extraordinary Diversity of Venom Peptides in Unexplored Predatory Gastropods of the Genus Clavus

Predatory gastropods of the superfamily Conoidea number over 12,000 living species. The evolutionary success of this lineage can be explained by the ability of conoideans to produce complex venoms for hunting, defense, and competitive interactions. Whereas venoms of cone snails (family Conidae) have become increasingly well studied, the venoms of most other conoidean lineages remain largely uncharacterized. In the present study, we present the venom gland transcriptomes of two species of the genus Clavus that belong to the family Drilliidae. Venom gland transcriptomes of two specimens of Clavus canalicularis and two specimens of Clavus davidgilmouri were analyzed, leading to the identification of a total of 1,176 putative venom peptide toxins (drillipeptides). Based on the combined evidence of secretion signal sequence identity, entire precursor similarity search (BLAST), and the orthology inference, putative Clavus toxins were assigned to 158 different gene families. The majority of identified transcripts comprise signal, pro-, mature peptide, and post-regions, with a typically short (<50 amino acids) and cysteine-rich mature peptide region. Thus, drillipeptides are structurally similar to conotoxins. However, convincing homology with known groups of Conus toxins was only detected for very few toxin families. Among these are Clavus counterparts of Conus venom insulins (drillinsulins), porins (drilliporins), and highly diversified lectins (drillilectins). The short size of most drillipeptides and structural similarity to conotoxins were unexpected, given that most related conoidean gastropod families (Terebridae and Turridae) possess longer mature peptide regions. Our findings indicate that, similar to conotoxins, drillipeptides may represent a valuable resource for future pharmacological exploration.

Bat coronavirus phylogeography in the Western Indian Ocean

Bats provide key ecosystem services such as crop pest regulation, pollination, seed dispersal, and soil fertilization. Bats are also major hosts for biological agents responsible for zoonoses, such as coronaviruses (CoVs). The islands of the Western Indian Ocean are identified as a major biodiversity hotspot, with more than 50 bat species. In this study, we tested 1,013 bats belonging to 36 species from Mozambique, Madagascar, Mauritius, Mayotte, Reunion Island and Seychelles, based on molecular screening and partial sequencing of the RNA-dependent RNA polymerase gene. In total, 88 bats (8.7%) tested positive for coronaviruses, with higher prevalence in Mozambican bats (20.5% ± 4.9%) as compared to those sampled on islands (4.5% ± 1.5%). Phylogenetic analyses revealed a large diversity of α- and β-CoVs and a strong signal of co-evolution between CoVs and their bat host species, with limited evidence for host-switching, except for bat species sharing day roost sites. These results highlight that strong variation between islands does exist and is associated with the composition of the bat species community on each island. Future studies should investigate whether CoVs detected in these bats have a potential for spillover in other hosts.

Revisiting the Woolly wolf (Canis lupus chanco) phylogeny in Himalaya: Addressing taxonomy, spatial extent and distribution of an ancient lineage in Asia

Of the sub-species of Holarctic wolf, the Woolly wolf (Canis lupus chanco) is uniquely adapted to atmospheric hypoxia and widely distributed across the Himalaya, Qinghai Tibetan Plateau (QTP) and Mongolia. Taxonomic ambiguity still exists for this sub-species because of complex evolutionary history anduse of limited wild samples across its range in Himalaya. We document for the first time population genetic structure and taxonomic affinity of the wolves across western and eastern Himalayan regions from samples collected from the wild (n = 19) using mitochondrial control region (225bp). We found two haplotypes in our data, one widely distributed in the Himalaya that was shared with QTP and the other confined to Himachal Pradesh and Uttarakhand in the western Himalaya, India. After combining our data withpublished sequences (n = 83), we observed 15 haplotypes. Some of these were shared among different locations from India to QTP and a few were private to geographic locations. A phylogenetic tree indicated that Woolly wolves from India, Nepal, QTP and Mongolia are basal to other wolves with shallow divergence (K2P; 0.000-0.044) and high bootstrap values. Demographic analyses based on mismatch distribution and Bayesian skyline plots (BSP) suggested a stable population over a long time (~million years) with signs of recent declines. Regional dominance of private haplotypes across its distribution range may indicate allopatric divergence. This may be due to differences in habitat characteristics, availability of different wild prey species and differential deglaciation within the range of the Woolly wolf during historic time. Presence of basal and shallow divergence within-clade along with unique ecological requirements and adaptation to hypoxia, the Woolly wolf of Himalaya, QTP, and Mongolian regions may be considered as a distinct an Evolutionary Significant Unit (ESU). Identifying management units (MUs) is needed within its distribution range using harmonized multiple genetic data for effective conservation planning.

Mitochondrial sequence-based evolutionary analysis of riverine-swamp hybrid buffaloes of India indicates novel maternal differentiation and domestication patterns

In this study, mitochondrial D-loop sequence data on riverine, swamp and hybrid buffaloes from India have been generated and compared with other reported Indian riverine, Chinese and Bangladeshi swamp buffalo populations. Sequence analysis revealed the presence of 132 haplotypes, with a haplotype diversity of 0.9611 ± 0.0045 and a nucleotide diversity of 0.04801 ± 0.00126. For the first time, the existence of riverine-swamp hybrids among the Indian Chilika buffalo population has been recorded, having 49 chromosomes, which was also confirmed by mitochondrial haplotype sharing between Chilika and Indian swamp as well as Chinese swamp buffalo populations in the network analysis. Phylogenetic analysis documents the sharing of reported pre-domestication haplogroups 'SA1', 'SA2', 'SA3' and 'SB1' between the Chilika and swamp buffalo populations of India, China and Bangladesh, an indication of the migration of swamp buffaloes towards Bangladesh and adjoining lower parts of India and north towards Chinese domestication sites. The results have also been supplemented by multidimension scaling, grouping Indian and Chinese swamp buffaloes more closely together with Bangladeshi buffaloes, but into a separate quadrant, whereas Chilika grouped away from other riverine as well as swamp buffaloes. These findings thus confirm the previous reports that the northeast region of India, close to the Indo-China border, is the point of evolution of swamp buffaloes with multiple sites of domestication.

Evolution of endemic and sylvatic lineages of dengue virus

Recent disease outbreaks have raised awareness of tropical pathogens, especially mosquito-borne viruses. Dengue virus (DENV) is a widely studied mammalian pathogen transmitted by various species of mosquito in the genus Aedes, especially Aedes aegypti and Aedes albopictus. The prevailing view of the research community is that endemic viral lineages that cause epidemics of DENV in humans have emerged over time from sylvatic viral lineages, which persist in wild, non-human primates. These notions have been examined by researchers through phylogenetic analyses of the envelope gene (E) from the four serotypes of DENV (serotypes DENV-1 to DENV-4). In these previous reports, researchers used visual inspection of a phylogeny in order to assert that sylvatic lineages lead to endemic clades. In making this assertion, these researchers also reasserted the model of periodic sylvatic to endemic disease outbreaks. Since that study, there has been a significant increase in data both in terms of metadata (e.g., place and host of isolation) and genetic sequences of DENV. Here, we re-examine the model of sylvatic to endemic shifts in viral lineages through a phylogenetic tree search and character evolution study of metadata on the tree. We built a dataset of nucleotide sequences for 188 isolates of DENV that have metadata on sylvatic or endemic sampling along with three orthologous sequences from West Nile virus as the outgroup for the phylogenetic analysis. In contrast to previous research, we find that there are several shifts from endemic to sylvatic lineages as well as sylvatic to endemic lineages, indicating a much more dynamic model of evolution. We propose that a model that allows oscillation between sylvatic and endemic hosts better captures the dynamics of DENV transmission.

Evolution and roles of cytokinin genes in angiosperms 2: Do ancient CKXs play housekeeping roles while non-ancient CKXs play regulatory roles?

Cytokinin oxidase/dehydrogenase (CKX) is a key enzyme responsible for the degradation of endogenous cytokinins. However, the origins and roles of CKX genes in angiosperm evolution remain unclear. Based on comprehensive bioinformatic and transgenic plant analyses, we demonstrate that the CKXs of land plants most likely originated from an ancient chlamydial endosymbiont during primary endosymbiosis. We refer to the CKXs retaining evolutionarily ancient characteristics as "ancient CKXs" and those that have expanded and functionally diverged in angiosperms as "non-ancient CKXs". We show that the expression of some non-ancient CKXs is rapidly inducible within 15 min upon the dehydration of Arabidopsis, while the ancient CKX (AtCKX7) is not drought responsive. Tobacco plants overexpressing a non-ancient CKX display improved oxidative and drought tolerance and root growth. Previous mutant studies have shown that non-ancient CKXs regulate organ development, particularly that of flowers. Furthermore, ancient CKXs preferentially degrade cis-zeatin (cZ)-type cytokinins, while non-ancient CKXs preferentially target N6-(Δ2-isopentenyl) adenines (iPs) and trans-zeatins (tZs). Based on the results of this work, an accompanying study (Wang et al. 10.1038/s41438-019-0211-x) and previous studies, we hypothesize that non-ancient CKXs and their preferred substrates of iP/tZ-type cytokinins regulate angiosperm organ development and environmental stress responses, while ancient CKXs and their preferred substrates of cZs play a housekeeping role, which echoes the conclusions and hypothesis described in the accompanying report (Wang, X. et al. Evolution and roles of cytokinin genes in angiosperms 1: Doancient IPTs play housekeeping while non-ancient IPTs play regulatory roles? Hortic Res7, (2020). 10.1038/s41438-019-0211-x).

Evolution and roles of cytokinin genes in angiosperms 1: Do ancient IPTs play housekeeping while non-ancient IPTs play regulatory roles?

Isopentenyltransferase (IPT) genes, including those encoding ATP/ADP-IPTs and tRNA-IPTs, control the rate-limiting steps of the biosynthesis of N 6-(Δ2-isopentenyl)adenine (iP)-type and trans-zeatin (tZ)-type cytokinins and cis-zeatin (cZ)-type cytokinins, respectively. However, the evolution and roles of these IPTs in angiosperms are not well understood. Here, we report comprehensive analyses of the origins, evolution, expression patterns, and possible roles of ATP/ADP-IPTs and tRNA-IPTs in angiosperms. We found that Class I and II tRNA-IPTs likely coexisted in the last common ancestor of eukaryotes, while ATP/ADP-IPTs likely originated from a Class II tRNA-IPT before the divergence of angiosperms. tRNA-IPTs are conservatively retained as 2-3 copies, but ATP/ADP-IPTs exhibit considerable expansion and diversification. Additionally, tRNA-IPTs are constitutively expressed throughout the plant, whereas the expression of ATP/ADP-IPTs is tissue-specific and rapidly downregulated by abiotic stresses. Furthermore, previous studies and our present study indicate that ATP/ADP-IPTs and their products, iPs/tZs, may regulate responses to environmental stresses and organ development in angiosperms. We therefore hypothesize that tRNA-IPTs and the associated cZs play a housekeeping role, whereas ATP/ADP-IPTs and the associated iP/tZ-type cytokinins play regulatory roles in organ development and stress responses in angiosperms, which echoes the conclusions and hypothesis presented in the accompanying study by Wang, X. et al Evolution and roles of cytokinin genes in angiosperms 2: Do ancient CKXs play housekeeping roles while non-ancient CKXs play regulatory roles? Hortic Res 10.1038/s41438-020-0246-z.

The first complete mitochondrial genome of the sand dollar Sinaechinocyamus mai (Echinoidea: Clypeasteroida)

Morphologic and molecular data often lead to different hypotheses of phylogenetic relationships. Such incongruence has been found in the echinoderm class Echinoidea. In particular, the phylogenetic status of the order Clypeasteroida is not well resolved. Complete mitochondrial genomes are currently available for 29 echinoid species, but no clypeasteroid had been sequenced to date. DNA extracted from a single live individual of Sinaechinocyamus mai was sequenced with 10× Genomics technology. This first complete mitochondrial genome (mitogenome) for the order Clypeasteroida is 15,756 base pairs in length. Phylogenomic analysis based on 34 ingroup taxa belonging to nine orders of the class Echinoidea show congruence between our new genetic inference and published trees based on morphologic characters, but also includes some intriguing differences that imply the need for additional investigation.

Bdellovibrio and Like Organisms in Lake Geneva: An Unseen Elephant in the Room?

When considering microbial biotic interactions, viruses as well as eukaryotic grazers are known to be important components of aquatic microbial food webs. It might be the same for bacterivorous bacteria but these groups have been comparatively less studied. This is typically the case of the Bdellovibrio and like organisms (BALOs), which are obligate bacterial predators of other bacteria. Recently, the abundance and distribution of three families of this functional group were investigated in perialpine lakes, revealing their presence and quantitative importance. Here, a more in-depth analysis is provided for Lake Geneva regarding the diversity of these bacterial predators at different seasons, sites and depths. We reveal a seasonal and spatial (vertical) pattern for BALOs. They were also found to be relatively diverse (especially Bdellovibrionaceae) and assigned to both known and unknown phylogenetic clusters. At last we found that most BALOs were positively correlated to other bacterial groups, mainly Gram-negative, in particular Myxococcales (among which many are predators of other microbes). This study is the first shedding light on this potentially important bacterial killing group in a large and deep lake.

Adaptive Radiation of the Flukes of the Family Fasciolidae Inferred from Genome-Wide Comparisons of Key Species

Liver and intestinal flukes of the family Fasciolidae cause zoonotic food-borne infections that impact both agriculture and human health throughout the world. Their evolutionary history and the genetic basis underlying their phenotypic and ecological diversity are not well understood. To close that knowledge gap, we compared the whole genomes of Fasciola hepatica, Fasciola gigantica, and Fasciolopsis buski and determined that the split between Fasciolopsis and Fasciola took place ∼90 Ma in the late Cretaceous period, and that between 65 and 50 Ma an intermediate host switch and a shift from intestinal to hepatic habitats occurred in the Fasciola lineage. The rapid climatic and ecological changes occurring during this period may have contributed to the adaptive radiation of these flukes. Expansion of cathepsins, fatty-acid-binding proteins, protein disulfide-isomerases, and molecular chaperones in the genus Fasciola highlights the significance of excretory-secretory proteins in these liver-dwelling flukes. Fasciola hepatica and Fasciola gigantica diverged ∼5 Ma near the Miocene-Pliocene boundary that coincides with reduced faunal exchange between Africa and Eurasia. Severe decrease in the effective population size ∼10 ka in Fasciola is consistent with a founder effect associated with its recent global spread through ruminant domestication. G-protein-coupled receptors may have key roles in adaptation of physiology and behavior to new ecological niches. This study has provided novel insights about the genome evolution of these important pathogens, has generated genomic resources to enable development of improved interventions and diagnosis, and has laid a solid foundation for genomic epidemiology to trace drug resistance and to aid surveillance.

The identification of tick autophagy-related genes in Ixodes scapularis responding to amino acid starvation

Ticks are obligate hematophagous arthropods and must tolerate starvation during off-host periods. Macroautophagy (hereafter autophagy) is a well-conserved self-eating mechanism of cell survival and is essential for recycling cellular contents during periods of starvation, stress, and injury in organisms. Although the genome sequence of Ixodes scapularis (Say) is available, the characteristics and functions of autophagy-related gene families remain largely unknown. To advance our understanding of autophagy in I. scapularis, we used comprehensive genomic approaches to identify Atg genes. Homologues of 14 Atg genes were identified, and their protein motif compositions were predicted. Phylogenetic analysis indicated that ATGs in I. scapularis were evolutionarily closely related to their homologues in Haemaphysalis longicornis and Rhipicephalus microplus ticks. Expression patterns of Atg genes differed across tick developmental stages. Immunofluorescence results by monodansylcadaverine (MDC) staining indicated that autophagy was activated after amino acid starvation treatments in I. scapularis embryo-derived cell lines ISE6 and IDE8. Subsequently, the expression of key Atg genes involved in autophagy pathway in both cell lines were examined. In ISE6 cells, the expression levels of three Atg genes (Atg4B, Atg6 and Atg8A) increased significantly after amino acid starvation; similarly, four Atg genes (Atg4A, Atg4B, Atg6 and Atg8B) were upregulated in IDE8 cells in response to starvation. In parallel, the MDC and lysotracker staining results indicated that autophagy was triggered after amino acid starvation treatments in R. microplus embryo-derived cell line BME26. Our observations showed that Atg family genes are highly conserved in ticks and function in autophagy pathway induced by amino acid starvation. These results also provide valuable insight for further autophagy-related research as a new strategy for blocking the transmission of tick-borne pathogens.

Comparative genomics reveal shared genomic changes in syngnathid fishes and signatures of genetic convergence with placental mammals

Syngnathids (seahorses, pipefishes and seadragons) exhibit an array of morphological innovations including loss of pelvic fins, a toothless tubular mouth and male pregnancy. They comprise two subfamilies: Syngnathinae and Nerophinae. Genomes of three Syngnathinae members have been analyzed previously. In this study, we have sequenced the genome of a Nerophinae member, the Manado pipefish (Microphis manadensis), which has a semi-enclosed brood pouch. Comparative genomic analysis revealed that the molecular evolutionary rate of the four syngnathids is higher than that of other teleosts. The loss of all but one P/Q-rich SCPP gene in the syngnathids suggests a role for the lost genes in dentin and enameloid formation in teleosts. Genome-wide comparison identified a set of 118 genes with parallel identical amino acid substitutions in syngnathids and placental mammals. Association of some of these genes with placental and embryonic development in mammals suggests a role for them in syngnathid pregnancy.

Evolutionary History of the Marchantia polymorpha Complex

The potential role of introgression in evolution has gained increased interest in recent years. Although some fascinating examples have been reported, more information is needed to generalize the importance of hybridization and introgression for adaptive divergence. As limited data exist on haploid dominant species, we analyzed genomes of three subspecies of the liverwort Marchantia polymorpha. We used available genomic data for subsp. ruderalis and carried out whole-genome (PacBio) sequencing for one individual each of subsp. montivagans and subsp. polymorpha as well as Illumina resequencing of additional genomes for all three subspecies. The three subspecies were compared against M. paleacea as outgroup. Our analyses revealed separation of the three taxa, but all three possible topologies were richly represented across the genomes, and the underlying divergence order less obvious. This uncertainty could be the result of the divergence of the three subspecies close in time, or that introgression has been frequent since divergence. In particular, we found that pseudo-chromosome 2 in subsp. montivagans was much more diverged than other parts of the genomes. This could either be explained by specific capture of chromosome 2 from an unknown related species through hybridization or by conservation of chromosome 2 despite intermittent or ongoing introgression affecting more permeable parts of the genomes. A higher degree of chromosomal rearrangements on pseudo-chromosome 2 support the second hypothesis. Species tree analyses recovered an overall topology where subsp. montivagans diverged first and subsp. ruderalis and subsp. polymorpha appeared as sister lineages. Each subspecies was associated with its own chloroplast and mitochondrial haplotype group. Our data suggest introgression but refute a previous hypothesis that subsp. ruderalis is a new stabilized hybrid between the other two subspecies.

Taxonomic diversity of Biomphalaria (Planorbidae) in São Paulo state, Brazil

Abstract: Morphological and molecular identifications were carried out for Biomphalaria occidentalis, Biomphalaria oligoza, Biomphalaria peregrina, Biomphalaria schrammi, Biomphalaria straminea and Biomphalaria tenagophila collected from 55 sites located along the upper basin of Tietê River in the Southeast Region of Brazil. Morphological analysis considered aspects of the shell, mantle, excretory organs and reproductive system. Molecular data included 122 sequences of Cytochrome C Oxidase I gene (COI). Our results showed that some shell characters, as well as other characters related to the mantle and the reproductive system, are fundamental for the identification of the six Biomphalaria species included in this study. The use of DNA barcoding together with morphological taxonomy generated more reliable results, proving to be a very useful approach, even for malacological surveillance services.

Demographic reconstruction from ancient DNA supports rapid extinction of the great auk

The great auk was once abundant and distributed across the North Atlantic. It is now extinct, having been heavily exploited for its eggs, meat, and feathers. We investigated the impact of human hunting on its demise by integrating genetic data, GPS-based ocean current data, and analyses of population viability. We sequenced complete mitochondrial genomes of 41 individuals from across the species' geographic range and reconstructed population structure and population dynamics throughout the Holocene. Taken together, our data do not provide any evidence that great auks were at risk of extinction prior to the onset of intensive human hunting in the early 16th century. In addition, our population viability analyses reveal that even if the great auk had not been under threat by environmental change, human hunting alone could have been sufficient to cause its extinction. Our results emphasise the vulnerability of even abundant and widespread species to intense and localised exploitation.

Characterization of the mitochondrial genome of the pathogenic fungus Scytalidium auriculariicola (Leotiomycetes) and insights into its phylogenetics

Scytalidium auriculariicola is the causative pathogen of slippery scar disease in the cultivated cloud ear fungus, Auricularia polytricha. In the present study, the mitogenome of S. auriculariicola was sequenced and assembled by next-generation sequencing technology. The circular mitogenome is 96,857 bp long and contains 56 protein-coding genes, 2 ribosomal RNA genes, and 30 transfer RNA genes (tRNAs). The high frequency of A and T used in codons contributed to the high AT content (73.70%) of the S. auriculariicola mitogenome. Comparative analysis indicated that the base composition and the number of introns and protein-coding genes in the S. auriculariicola mitogenome varied from that of other Leotiomycetes mitogenomes, including a uniquely positive AT skew. Five distinct groups were found in the gene arrangements of Leotiomycetes. Phylogenetic analyses based on combined gene datasets (15 protein-coding genes) yielded well-supported (BPP = 1) topologies. A single-gene phylogenetic tree indicated that the nad4 gene may be useful as a molecular marker to analyze the phylogenetic relationships of Leotiomycetes species. This study is the first report on the mitochondrial genome of the genus Scytalidium, and it will contribute to our understanding of the population genetics and evolution of S. auriculariicola and related species.

Description of a human Bocavirus recombinant strain in the Americas

Human bocaviruses (HBoV) are mainly associated with respiratory and gastroenteric infections. These viruses belong to the family Parvoviridae, genus Bocaparvovirus and are classified in four subtypes (HBoV1-4). Recombination and point mutation have been described as basis of parvovirus evolution. In this study three viral sequences were obtained from positives HBoV sewage samples collected in two Uruguayan cities and were characterised by different methods as recombinant strains. This recombination event was localised in the 5’ end of VP1 gene and the parental strains belonged to subtypes 3 and 4. These three Uruguayan strains are identical at the nucleotide sequences in the analysed genome region of the virus. As far as we known, this study represents the first detection of HBoV recombinants strains in the Americas.

The global diversity of Haemonchus contortus is shaped by human intervention and climate

Haemonchus contortus is a haematophagous parasitic nematode of veterinary interest. We have performed a survey of its genome-wide diversity using single-worm whole genome sequencing of 223 individuals sampled from 19 isolates spanning five continents. We find an African origin for the species, together with evidence for parasites spreading during the transatlantic slave trade and colonisation of Australia. Strong selective sweeps surrounding the β-tubulin locus, a target of benzimidazole anthelmintic drug, are identified in independent populations. These sweeps are further supported by signals of diversifying selection enriched in genes involved in response to drugs and other anthelmintic-associated biological functions. We also identify some candidate genes that may play a role in ivermectin resistance. Finally, genetic signatures of climate-driven adaptation are described, revealing a gene acting as an epigenetic regulator and components of the dauer pathway. These results begin to define genetic adaptation to climate in a parasitic nematode.

Based on single worm whole genome sequencing, the authors here characterise the global evolution of the gastrointestinal parasite Haemonchus contortus and identify genes that play a role in drug resistance as well as climate-driven adaptations involving an epigenetic regulator.

An evolutionary model motivated by physicochemical properties of amino acids reveals variation among proteins

Motivation

The relative rates of amino acid interchanges over evolutionary time are likely to vary among proteins. Variation in those rates has the potential to reveal information about constraints on proteins. However, the most straightforward model that could be used to estimate relative rates of amino acid substitution is parameter-rich and it is therefore impractical to use for this purpose.

Results

A six-parameter model of amino acid substitution that incorporates information about the physicochemical properties of amino acids was developed. It showed that amino acid side chain volume, polarity and aromaticity have major impacts on protein evolution. It also revealed variation among proteins in the relative importance of those properties. The same general approach can be used to improve the fit of empirical models such as the commonly used PAM and LG models.

Availability and implementation

Perl code and test data are available from https://github.com/ebraun68/sixparam.

Supplementary information

Supplementary data are available at Bioinformatics online.

Sympatric western lowland gorillas, central chimpanzees and humans are infected with different trichomonads

We investigated intestinal trichomonads in western lowland gorillas, central chimpanzees and humans cohabiting the forest ecosystem of Dzanga-Sangha Protected Area in Central African Republic, using the internal transcribed spacer (ITS) region and SSU rRNA gene sequences. Trichomonads belonging to the genus Tetratrichomonas were detected in 23% of the faecal samples and in all host species. Different hosts were infected with different genotypes of Tetratrichomonas. In chimpanzees, we detected tetratrichomonads from 'novel lineage 2', which was previously reported mostly in captive and wild chimpanzees. In gorillas, we found two different genotypes of Tetratrichomonas. The ITS region sequences of the more frequent genotype were identical to the sequence found in a faecal sample of a wild western lowland gorilla from Cameroon. Sequences of the second genotype from gorillas were almost identical to sequences previously obtained from an anorexic French woman. We provide the first report of the presence of intestinal tetratrichomonads in asymptomatic, apparently healthy humans. Human tetratrichomonads belonged to the lineage 7, which was previously reported in domestic and wild pigs and a domestic horse. Our findings suggest that the ecology and spatial overlap among hominids in the tropical forest ecosystem has not resulted in exchange of intestinal trichomonads among these hosts.

Next-Generation Genome Sequencing of Sedum plumbizincicola Sheds Light on the Structural Evolution of Plastid rRNA Operon and Phylogenetic Implications within Saxifragales

The genus Sedum, with about 470 recognized species, is classified in the family Crassulaceae of the order Saxifragales. Phylogenetic relationships within the Saxifragales are still unresolved and controversial. In this study, the plastome of S. plumbizincicola was firstly presented, with a focus on the structural analysis of rrn operon and phylogenetic implications within the order Saxifragaceae. The assembled complete plastome of S. plumbizincicola is 149,397 bp in size, with a typical circular, double-stranded, and quadripartite structure of angiosperms. It contains 133 genes, including 85 protein-coding genes (PCGs), 36 tRNA genes, 8 rRNA genes, and four pseudogenes (one ycf1, one rps19, and two ycf15). The predicted secondary structure of S. plumbizincicola 16S rRNA includes three main domains organized in 74 helices. Further, our results confirm that 4.5S rRNA of higher plants is associated with fragmentation of 23S rRNA progenitor. Notably, we also found the sequence of putative rrn5 promoter has some evolutionary implications within the order Saxifragales. Moreover, our phylogenetic analyses suggested that S. plumbizincicola had a closer relationship with S. sarmentosum than S. oryzifolium, and supported the taxonomic revision of Phedimus. Our findings of the present study will be useful for further investigation of the evolution of plastid rRNA operon and phylogenetic relationships within Saxifragales.

Different maternal lineages revealed by ancient mitochondrial genome of Camelus bactrianus from China

Domestic Bactrian camel (Camelus bactrianus) used to be one of the most important livestock species in Chinese history, as well as the major transport carrier on the ancient Silk Road. However, archeological studies on Chinese C. bactrianus are still limited, and molecular biology research on this species is mainly focused on modern specimens. In this study, we retrieved the complete mitochondrial genome from a C. bactrianus specimen, which was excavated from northwestern China and dated at 1290-1180 cal. years before present (yBP). Phylogenetic analyses using 18 mitochondrial genomes indicated that the C. bactrianus clade was divided into two maternal lineages. The majority of samples originating from Iran to Japan and Mongolia belong to subclade A1, while our sample together with two Mongolian individuals formed the much smaller subclade A2. Furthermore, the divergence time of these two maternal lineages was estimated as 165 Kya (95% credibility interval 117-222 Kya), this might indicate that several different evolutionary lineages were incorporated into the domestic gene pool during the initial domestication process. Bayesian skyline plot (BSP) analysis suggest a slow increase in female effective population size of C. bactrianus from 5000 years ago, which corresponds to the beginning of domestication of C. bactrianus. The present study also revealed that there were extensive exchanges of genetic information among C. bactrianus populations in regions along the Silk Road.

Characterization and Expression Analysis of ERF Genes in Fragaria vesca Suggest Different Divergences of Tandem ERF Duplicates

Ethylene-responsive factors (ERFs) play important roles in plant growth and development and in responses to abiotic stresses. However, little information was available about the ERF genes in woodland strawberry (Fragaria vesca), a genetic model plant for the Fragaria genus and Rosaceae family. In this study, 91 FveERF genes were identified, including 35 arrayed in tandem, indicating that tandem duplication is a major mechanism for the expansion of the FveERF family. According to their phylogenetic relationships with AtERFs from Arabidopsis thaliana, the tandem FveERF genes could be grouped into ancestral and lineage-specific tandem ones. The ancestral tandem FveERFs are likely derived from tandem duplications that occurred in the common ancestor of F. vesca and A. thaliana, whereas the lineage-specific ones are specifically present in the F. vesca lineage. The lineage-specific tandem FveERF duplicates are more conserved than the ancestral ones in sequence and structure. However, their expression in flowers and fruits is similarly diversified, indicating that tandem FveERFs have diverged rapidly after duplication in this respect. The lineage-specific tandem FveERFs display the same response patterns with only one exception under drought or cold, whereas the ancestral tandem ones are largely differentially expressed, suggesting that divergence of tandem FveERF expression under stress may have occurred later in the reproductive development. Our results provide evidence that the retention of tandem FveERF duplicates soon after their duplication may be related to their divergence in the regulation of reproductive development. In contrast, their further divergence in expression pattern likely contributes to plant response to abiotic stress.

Strigolactone synthesis is ancestral in land plants, but canonical strigolactone signalling is a flowering plant innovation

Background

Strigolactones (SLs) are an important class of carotenoid-derived signalling molecule in plants, which function both as exogenous signals in the rhizosphere and as endogenous plant hormones. In flowering plants, SLs are synthesized by a core pathway of four enzymes and are perceived by the DWARF14 (D14) receptor, leading to degradation of SMAX1-LIKE7 (SMXL7) target proteins in a manner dependent on the SCF^MAX2 ubiquitin ligase. The evolutionary history of SLs is poorly understood, and it is not clear whether SL synthesis and signalling are present in all land plant lineages, nor when these traits evolved.

Results

We have utilized recently-generated genomic and transcriptomic sequences from across the land plant clade to resolve the origin of each known component of SL synthesis and signalling. We show that all enzymes in the core SL synthesis pathway originated at or before the base of land plants, consistent with the previously observed distribution of SLs themselves in land plant lineages. We also show that the late-acting enzyme LATERAL BRANCHING OXIDOREDUCTASE (LBO) may be considerably more ancient than previously thought. We perform a detailed phylogenetic analysis of SMXL proteins and show that specific SL target proteins only arose in flowering plants. We also assess diversity and protein structure in the SMXL family, identifying several previously unknown clades.

Conclusions

Overall, our results suggest that SL synthesis is much more ancient than canonical SL signalling, consistent with the idea that SLs first evolved as rhizosphere signals and were only recruited much later as hormonal signals.

Electronic supplementary material

The online version of this article (10.1186/s12915-019-0689-6) contains supplementary material, which is available to authorized users.

The Bacterial Symbionts of Closely Related Hydrothermal Vent Snails With Distinct Geochemical Habitats Show Broad Similarity in Chemoautotrophic Gene Content

Symbiosis has evolved between a diversity of invertebrate taxa and chemosynthetic bacterial lineages. At the broadest level, these symbioses share primary function: the bacterial symbionts use the energy harnessed from the oxidation of reduced chemicals to power the fixation of inorganic carbon and/or other nutrients, providing the bulk of host nutrition. However, it is unclear to what extent the ecological niche of the host species is influenced by differences in symbiont traits, particularly those involved in chemoautotrophic function and interaction with the geochemical environment. Hydrothermal vents in the Lau Basin (Tonga) are home to four morphologically and physiologically similar snail species from the sister genera Alviniconcha and Ifremeria. Here, we assembled nearly complete genomes from their symbionts to determine whether differences in chemoautotrophic capacity exist among these symbionts that could explain the observed distribution of these snail species into distinct geochemical habitats. Phylogenomic analyses confirmed that the symbionts have evolved from four distinct lineages in the classes γ-proteobacteria or Campylobacteria. The genomes differed with respect to genes related to motility, adhesion, secretion, and amino acid uptake or excretion, though were quite similar in chemoautotrophic function, with all four containing genes for carbon fixation, sulfur and hydrogen oxidation, and oxygen and nitrate respiration. This indicates that differences in the presence or absence of symbiont chemoautotrophic functions does not likely explain the observed geochemical habitat partitioning. Rather, differences in gene expression and regulation, biochemical differences among these chemoautotrophic pathways, and/or differences in host physiology could all influence the observed patterns of habitat partitioning.

Molecular evolution of mammalian genes with epistatic interactions in fertilization

BACKGROUND

Genes that encode proteins associated with sperm competition, fertilization, and sexual conflicts of interest are often among the most rapidly evolving parts of animal genomes. One family of sperm-expressed genes (Zp3r, C4bpa) in the mammalian gene cluster called the regulator of complement activation (RCA) encodes proteins that bind eggs and mediate reproductive success, and are therefore expected to show high relative rates of nonsynonymous nucleotide substitution in response to sexual selection in comparison to other genes not involved in gamete binding at fertilization. We tested that working hypothesis by using phylogenetic models of codon evolution to identify episodes of diversifying positive selection. We used a comparative approach to quantify the evidence for episodic diversifying selection acting on RCA genes with known functions in fertilization (and sensitivity to sexual selection), and contrast them with other RCA genes in the same gene family that function in innate immunity (and are not sensitive to sexual selection).

RESULTS

We expected but did not find evidence for more episodes of positive selection on Zp3r in Glires (the rodents and lagomorphs) or on C4BPA in Primates, in comparison to other paralogous RCA genes in the same taxon, or in comparison to the same orthologous RCA gene in the other taxon. That result was not unique to RCA genes: we also found little evidence for more episodes of diversifying selection on genes that encode selective sperm-binding molecules in the egg coat or zona pellucida (Zp2, Zp3) in comparison to members of the same gene family that encode structural elements of the egg coat (Zp1, Zp4). Similarly, we found little evidence for episodic diversifying selection acting on two other recently discovered genes (Juno, Izumo1) that encode essential molecules for sperm-egg fusion.

CONCLUSIONS

These negative results help to illustrate the importance of a comparative context for this type of codon model analysis. The results may also point to other phylogenetic contexts in which the effects of selection acting on these fertilization proteins might be more readily discovered and documented in mammals and other taxa.

Selection, periodicity and potential function for Highly Iterative Palindrome-1 (HIP1) in cyanobacterial genomes

Highly Iterated Palindrome 1 (HIP1, GCGATCGC) is hyper-abundant in most cyanobacterial genomes. In some cyanobacteria, average HIP1 abundance exceeds one motif per gene. Such high abundance suggests a significant role in cyanobacterial biology. However, 20 years of study have not revealed whether HIP1 has a function, much less what that function might be. We show that HIP1 is 15- to 300-fold over-represented in genomes analyzed. More importantly, HIP1 sites are conserved both within and between open reading frames, suggesting that their overabundance is maintained by selection rather than by continual replenishment by neutral processes, such as biased DNA repair. This evidence for selection suggests a functional role for HIP1. No evidence was found to support a functional role as a peptide or RNA motif or a role in the regulation of gene expression. Rather, we demonstrate that the distribution of HIP1 along cyanobacterial chromosomes is significantly periodic, with periods ranging from 10 to 90 kb, consistent in scale with periodicities reported for co-regulated, co-expressed and evolutionarily correlated genes. The periodicity we observe is also comparable in scale to chromosomal interaction domains previously described in other bacteria. In this context, our findings imply HIP1 functions associated with chromosome and nucleoid structure.

Inadvertent Paralog Inclusion Drives Artifactual Topologies and Timetree Estimates in Phylogenomics

Increasingly, large phylogenomic data sets include transcriptomic data from nonmodel organisms. This not only has allowed controversial and unexplored evolutionary relationships in the tree of life to be addressed but also increases the risk of inadvertent inclusion of paralogs in the analysis. Although this may be expected to result in decreased phylogenetic support, it is not clear if it could also drive highly supported artifactual relationships. Many groups, including the hyperdiverse Lissamphibia, are especially susceptible to these issues due to ancient gene duplication events and small numbers of sequenced genomes and because transcriptomes are increasingly applied to resolve historically conflicting taxonomic hypotheses. We tested the potential impact of paralog inclusion on the topologies and timetree estimates of the Lissamphibia using published and de novo sequencing data including 18 amphibian species, from which 2,656 single-copy gene families were identified. A novel paralog filtering approach resulted in four differently curated data sets, which were used for phylogenetic reconstructions using Bayesian inference, maximum likelihood, and quartet-based supertrees. We found that paralogs drive strongly supported conflicting hypotheses within the Lissamphibia (Batrachia and Procera) and older divergence time estimates even within groups where no variation in topology was observed. All investigated methods, except Bayesian inference with the CAT-GTR model, were found to be sensitive to paralogs, but with filtering convergence to the same answer (Batrachia) was observed. This is the first large-scale study to address the impact of orthology selection using transcriptomic data and emphasizes the importance of quality over quantity particularly for understanding relationships of poorly sampled taxa.

Paternally Expressed Imprinted Genes under Positive Darwinian Selection in Arabidopsis thaliana

Genomic imprinting is an epigenetic phenomenon where autosomal genes display uniparental expression depending on whether they are maternally or paternally inherited. Genomic imprinting can arise from parental conflicts over resource allocation to the offspring, which could drive imprinted loci to evolve by positive selection. We investigate whether positive selection is associated with genomic imprinting in the inbreeding species Arabidopsis thaliana. Our analysis of 140 genes regulated by genomic imprinting in the A. thaliana seed endosperm demonstrates they are evolving more rapidly than expected. To investigate whether positive selection drives this evolutionary acceleration, we identified orthologs of each imprinted gene across 34 plant species and elucidated their evolutionary trajectories. Increased positive selection was sought by comparing its incidence among imprinted genes with nonimprinted controls. Strikingly, we find a statistically significant enrichment of imprinted paternally expressed genes (iPEGs) evolving under positive selection, 50.6% of the total, but no such enrichment for positive selection among imprinted maternally expressed genes (iMEGs). This suggests that maternally- and paternally expressed imprinted genes are subject to different selective pressures. Almost all positively selected amino acids were fixed across 80 sequenced A. thaliana accessions, suggestive of selective sweeps in the A. thaliana lineage. The imprinted genes under positive selection are involved in processes important for seed development including auxin biosynthesis and epigenetic regulation. Our findings support a genomic imprinting model for plants where positive selection can affect paternally expressed genes due to continued conflict with maternal sporophyte tissues, even when parental conflict is reduced in predominantly inbreeding species.

Characterization of the complete mitochondrial genome of Brentisentisyangtzensis Yu & Wu, 1989 (Acanthocephala, Illiosentidae)

The mitogenome of Brentisentisyangtzensis is 13,864 bp in length and has the circular structure typical of metazoans. It contains 36 genes: 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs) and 12 protein-encoding genes (PCGs). All genes are transcribed from the same strand. Thirteen overlapping regions were found in the mitochondrial genome. The overall A+T content of B.yangtzensis is 68.3% versus 31.7% of G+C content (A = 27.8%, T = 40.5%, C = 9.0%, G = 22.7%). B.yangtzenensis (Illiosentidae) and Leptorhynchoidesthecatus (Rhadinorhynchidae) form a sister clade, showing the relatively close relationship between the Illiosentidae and the Rhadinorhynchidae. The mitochondrial gene arrangements of acanthocephalan species are relatively conserved, with only a few translocations of tRNAs (trnS1, trnS2, trnV, and trnK) detected. An identical gene order was found both in a sister clade (Centrorhynchusaluconis and Plagiorhynchustransversus) and across different classes (B.yangtzensis (Palaeacanthocephala), Acanthosentischeni (Eoacanthocephala) and Macracanthorhynchushirudinaceus (Archiacanthocephala), Oncicolaluehei and L.thecatus (Palaeacanthocephala)). More studies and more sequences of acanthocephalan species are needed to gain a clear understanding of the phylogenetic relationships.

Characterization and phylogenomics of the complete mitochondrial genome of the polyzoic cestode Gangesia oligonchis (Platyhelminthes: Onchoproteocephalidea)

The order Onchoproteocephalidea (Eucestoda) was recently erected to accommodate the hook-bearing tetraphyllideans and the proteocephalideans, which are characterized by internal proglottization and a tetra-acetabulate scolex. The recognized subfamilies in the Proteocephalidae appeared to be non-monophyletic based on 28S recombinant DNA (rDNA) sequence data. Other molecular markers with higher phylogenetic resolution, such as large mitochondrial DNA fragments and multiple genes, are obviously needed. Thus the mitochondrial genome of Gangesia oligonchis, belonging to the putative earliest diverging group of the Proteocephalidae, was sequenced. The circular mitogenome of G. oligonchis was 13,958 bp in size, and contained the standard 36 genes: 22 transfer RNA genes, two rRNA genes and 12 protein-coding genes, as well as two major non-coding regions. A short NCR and a large NCR (lNCR) region were 216 bp and 419 bp in size, respectively. Highly repetitive regions in the lNCR region were detected with that of 11 repeat units. The mitogenome of G. oligonchis shared 71.1% nucleotide identity with Testudotaenia sp. WL-2016. Phylogenetic analyses of the complete mitochondrial genomes with Bayesian inference and maximum likelihood methods indicated that G. oligonchis formed a sister clade with Testudotaenia sp. WL-2016 with maximum support. The ordinal topology is (Caryophyllidea, (Diphyllobothriidea, (Bothriocephalidea, (Onchoproteocephalidea, Cyclophyllidea)))). The mitogenomic gene arrangement of G. oligonchis was identical to that of Testudotaenia sp. WL-2016. Both mitogenomic and nuclear sequence data for many more taxa are required to effectively explore the inter-relationships among the Onchoproteocephalidea.

Co-infection of Newcastle disease virus genotype XIII with low pathogenic avian influenza exacerbates clinical outcome of Newcastle disease in vaccinated layer poultry flocks

Newcastle disease (ND) and avian influenza (AI) are economically important infectious diseases of poultry. Sometime, concomitant secondary viral/or bacterial infections significantly alters the pathobiology of ND and AI in poultry. As of now, the disease patterns and dynamics of co-infections caused by ND virus (NDV, genotype XIII) and Low Pathogenic AI viruses (LPAI, H9N2) are explicitly elusive. Thus, we examined the clinicopathological disease conditions due to these two economically important viruses to understand the complex disease outcomes by virus-virus interactions in vaccinated flocks. The findings of clinicopathological and molecular investigations carried on 37 commercial ND vaccinated poultry flocks revealed simultaneous circulation of NDV and AIV in same flock/bird. Further, molecular characterization of hemagglutinin (HA) and neuraminidase (NA) genes confirmed that all the identified AIVs were of low pathogenicity H9N2 subtype and fusion (F) gene analysis of detected NDVs belong to NDV class II, genotype XIII, a virulent type. The NDV and H9N2 alone or co-infected flocks (NDV + LPAI) exhibit clinical signs and lesions similar to that of virulent NDV except the degree of severity, which was higher in H9N2-NDV co-infected flocks. Additionally, avian pathogenic E. coli and mycoplasma infections were detected in majority of the ailing/dead birds from the co-infected flocks during progression of the clinical disease. Overall, the findings highlight the multi-factorial disease complexity in commercial poultry and suggest the importance of NDV genotype XIII in intensifying the clinical disease in vaccinated birds.

Physicochemical Amino Acid Properties Better Describe Substitution Rates in Large Populations

Substitutions between chemically distant amino acids are known to occur less frequently than those between more similar amino acids. This knowledge, however, is not reflected in most codon substitution models, which treat all nonsynonymous changes as if they were equivalent in terms of impact on the protein. A variety of methods for integrating chemical distances into models have been proposed, with a common approach being to divide substitutions into radical or conservative categories. Nevertheless, it remains unclear whether the resulting models describe sequence evolution better than their simpler counterparts. We propose a parametric codon model that distinguishes between radical and conservative substitutions, allowing us to assess if radical substitutions are preferentially removed by selection. Applying our new model to a range of phylogenomic data, we find differentiating between radical and conservative substitutions provides significantly better fit for large populations, but see no equivalent improvement for smaller populations. Comparing codon and amino acid models using these same data shows that alignments from large populations tend to select phylogenetic models containing information about amino acid exchangeabilities, whereas the structure of the genetic code is more important for smaller populations. Our results suggest selection against radical substitutions is, on average, more pronounced in large populations than smaller ones. The reduced observable effect of selection in smaller populations may be due to stronger genetic drift making it more challenging to detect preferences. Our results imply an important connection between the life history of a phylogenetic group and the model that best describes its evolution.

Low genetic diversity in a critically endangered primate: shallow evolutionary history or recent population bottleneck?

BACKGROUND

Current patterns of population genetic variation may have been shaped by long-term evolutionary history and contemporary demographic processes. Understanding the underlying mechanisms that yield those patterns is crucial for informed conservation of endangered species. The critically endangered white-headed langur, Trachypithecus leucocephalus, is endemic to a narrow range in southwest China. This species shows very low genetic diversity in its 2 main relict populations, Fusui and Chongzuo. Whether this has been caused by a short evolutionary history or recent population declines is unknown. Therefore, we investigated the contributions of historical and recent population demographic changes to population genetic diversity by using 15 nuclear microsatellite markers and mitochondrial DNA (mtDNA) control region sequences.

RESULTS

Using genetic data from 214 individuals we found a total of 9 mtDNA haplotypes in the Fusui population but only 1 haplotype in the Chongzuo population, and we found an overall low genetic diversity (haplotype and nucleotide diversities: h = 0.486 ± 0.036; π = 0.0028 ± 0.0003). The demographic history inferred from mtDNA and microsatellite markers revealed no evidence for historical population size fluctuations or recent population bottlenecks. Simulations of possible population divergence histories inferred by DIYABC analysis supported a recent divergence of the Chongzuo population from the Fusui population and no population bottlenecks.

CONCLUSIONS

Despite severe population declines caused by anthropogenic activities in the last century, the low genetic diversity of the extant white-headed langur populations is most likely primarily due to the species' shallow evolutionary history and to a recent, local population founder event.

Genetic Analysis and Evolutionary Changes of the Torque teno sus Virus

The torque teno sus virus (TTSuV) is an emerging virus threating the Suidae species of unclear pathogenicity, although it was previously reported as a worsening factor of other porcine diseases, in particular, porcine circovirus associated disease (PCVAD). Here, a comprehensive codon usage analysis of the open reading frame 1 (ORF1), which encodes the viral capsid protein, was undertaken for the first time to reveal its evolutionary history. We revealed independent phylogenetic processes for the two genera during TTSuV evolution, which was confirmed by principal component analysis (PCA). A low codon usage bias was observed in different genera and different species, with Kappatorquevirus a (TTSuVk2a) displaying the highest, which was mainly driven by mutation pressure and natural selection, especially natural selection. Overall, ATs were more abundant than GCs, along with more A-ended synonymous codons in relative synonymous codon usage (RSCU) analysis. To further confirm the role of natural selection and TTSuV adaptation to the Suidae species, codon adaptation index (CAI), relative codon deoptimization index (RCDI), and similarity index (SiD) analyses were performed, which showed different adaptations for different TTSuVs. Importantly, we identified a more dominant role of Sus scrofa in the evolution of Iotatorquevirus (TTSuV1), with the highest CAI values and lowest RCDI values compared to Sus scrofa domestica. However, in TTSuVk2, the roles of Sus scrofa and Sus scrofa domestica were the same, regarding codon usage, with similar CAI and RCDI values. Our study provides a new perspective of the evolution of TTSuV and valuable information to develop control measures against TTSuV.

Tracking Five Millennia of Horse Management with Extensive Ancient Genome Time Series

Horse domestication revolutionized warfare and accelerated travel, trade, and the geographic expansion of languages. Here, we present the largest DNA time series for a non-human organism to date, including genome-scale data from 149 ancient animals and 129 ancient genomes (≥1-fold coverage), 87 of which are new. This extensive dataset allows us to assess the modern legacy of past equestrian civilizations. We find that two extinct horse lineages existed during early domestication, one at the far western (Iberia) and the other at the far eastern range (Siberia) of Eurasia. None of these contributed significantly to modern diversity. We show that the influence of Persian-related horse lineages increased following the Islamic conquests in Europe and Asia. Multiple alleles associated with elite-racing, including at the MSTN "speed gene," only rose in popularity within the last millennium. Finally, the development of modern breeding impacted genetic diversity more dramatically than the previous millennia of human management.

Detection and characterization of a novel marine birnavirus isolated from Asian seabass in Singapore

Background

Lates calcarifer, known as seabass in Asia and barramundi in Australia, is a widely farmed species internationally and in Southeast Asia and any disease outbreak will have a great economic impact on the aquaculture industry. Through disease investigation of Asian seabass from a coastal fish farm in 2015 in Singapore, a novel birnavirus named Lates calcarifer Birnavirus (LCBV) was detected and we sought to isolate and characterize the virus through molecular and biochemical methods.

Methods

In order to propagate the novel birnavirus LCBV, the virus was inoculated into the Bluegill Fry (BF-2) cell line and similar clinical signs of disease were reproduced in an experimental fish challenge study using the virus isolate. Virus morphology was visualized using transmission electron microscopy (TEM). Biochemical analysis using chloroform and 5-Bromo-2′-deoxyuridine (BUDR) sensitivity assays were employed to characterize the virus. Next-Generation Sequencing (NGS) was also used to obtain the virus genome for genetic and phylogenetic analyses.

Results

The LCBV-infected BF-2 cell line showed cytopathic effects such as rounding and granulation of cells, localized cell death and detachment of cells observed at 3 to 5 days’ post-infection. The propagated virus, when injected intra-peritoneally into naïve Asian seabass under experimental conditions, induced lesions similar to fish naturally infected with LCBV. Morphology of LCBV, visualized under TEM, revealed icosahedral particles around 50 nm in diameter. Chloroform and BUDR sensitivity assays confirmed the virus to be a non-enveloped RNA virus. Further genome analysis using NGS identified the virus to be a birnavirus with two genome segments. Phylogenetic analyses revealed that LCBV is more closely related to the Blosnavirus genus than to the Aquabirnavirus genus within the Birnaviridae family.

Conclusions

These findings revealed the presence of a novel birnavirus that could be linked to the disease observed in the Asian seabass from the coastal fish farms in Singapore. This calls for more studies on disease transmission and enhanced surveillance programs to be carried out to understand pathogenicity and epidemiology of this novel virus. The gene sequences data obtained from the study can also pave way to the development of PCR-based diagnostic test methods that will enable quick and specific identification of the virus in future disease investigations.

Composite Analysis of the Virome and Bacteriome of HIV/HPV Co-Infected Women Reveals Proxies for Immunodeficiency

The human cervical microbiome is complex, and its role in health and disease has just begun to be elucidated. In this study, 57 cervical swab samples from 19 HIV/HPV co-infected women were analyzed for both virome and bacteriome composition. Virome analysis focused on circular DNA viruses through rolling circle amplification followed by next-generation sequencing (NGS). Data were assigned to virus families and genera, and HPV types were identified. NGS data of bacterial 16S from a subset of 24 samples were assigned to operational taxonomic units and classified according to vaginal microbiome community state types (CSTs). Four viral families were found: Papillomaviridae, Anelloviridae, Genomoviridae, and Herpesviridae. Papillomavirus reads were more abundant in women with premalignant cervical lesions, which were also strongly associated with multiple (≥3) high-risk HPV infection. Anellovirus read abundance was negatively correlated with host CD4+ T-cell counts. The bacteriome revealed the presence of CST III and CST IV, and women with ≥1% frequency of genomovirus or herpesvirus reads displayed an increased risk of carrying CST IV. By characterizing the composition of the cervical circular DNA viruses and the bacteriome of HIV/HPV co-infected women, we identified putative interactions between these two microorganism communities and their associations with patients’ clinical characteristics, notably immunodeficiency status.

Estimating HIV-1 Genetic Diversity in Brazil Through Next-Generation Sequencing

Approximately 36.7 million people were living with the human immunodeficiency virus (HIV) at the end of 2016 according to UNAIDS, representing a global prevalence rate of 0.8%. In Brazil, an HIV prevalence of 0.24% has been estimated, which represents approximately 830,000 individuals living with the virus. As a touristic and commercial hub in Latin America, Brazil harbors an elevated HIV genetic variability, further contributed by the selective pressure exerted by the host immune system and by antiretroviral treatment. Through the progress of the next-generation sequencing (NGS) techniques, it has been possible to expand the study of HIV genetic diversity, evolutionary, and epidemic processes, allowing the generation of HIV complete or near full-length genomes (NFLG) and improving the characterization of intra- and interhost diversity of viral populations. Greater sensitivity in the detection of viral recombinant forms represents one of the major improvements associated with this development. It is possible to identify unique or circulating recombinant forms using the near full-length viral genomes with increasing accuracy. It also permits the characterization of multiple viral infections within individual hosts. Previous Brazilian studies using NGS to analyze HIV diversity were able to identify several distinct unique and circulating recombinant forms and evidenced dual infections. These data unveiled unprecedented high rates of viral recombination and highlighted that novel recombinants are continually arising in the Brazilian epidemic. In the pooled analysis depicted in this report, HIV subtypes have been determined from HIV-positive patients in five states of Brazil with some of the highest HIV prevalence, three in the Southeast (Rio de Janeiro, São Paulo, and Minas Gerais), one in the Northeast (Pernambuco) and one in the South (Rio Grande do Sul). Combined data analysis showed a significant prevalence of recombinant forms (29%; 101/350), and a similar 26% when only NFLGs were considered. Moreover, the analysis was able to evidence the occurrence of multiple infections in some individuals. Our data highlight the great HIV genetic diversity found in Brazil and unveils a more accurate scenario of the HIV evolutionary dynamics in the region.

Nuclear protein phylogenies support the monophyly of the three bryophyte groups (Bryophyta Schimp.)

Unraveling the phylogenetic relationships between the four major lineages of terrestrial plants (mosses, liverworts, hornworts, and vascular plants) is essential for an understanding of the evolution of traits specific to land plants, such as their complex life cycles, and the evolutionary development of stomata and vascular tissue. Well supported phylogenetic hypotheses resulting from different data and methods are often incongruent due to processes of nucleotide evolution that are difficult to model, for example substitutional saturation and composition heterogeneity. We reanalysed a large published dataset of nuclear data and modelled these processes using degenerate-codon recoding and tree-heterogeneous composition substitution models. Our analyses resolved bryophytes as a monophyletic group and showed that the nonnonmonophyly of the clade that is supported by the analysis of nuclear nucleotide data is due solely to fast-evolving synonymous substitutions. The current congruence among phylogenies of both nuclear and chloroplast analyses lent considerable support to the conclusion that the bryophytes are a monophyletic group. An initial split between bryophytes and vascular plants implies that the bryophyte life cycle (with a dominant gametophyte nurturing an unbranched sporophyte) may not be ancestral to all land plants and that stomata are likely to be a symplesiomorphy among embryophytes.

Genome-wide identification, expression and functional analysis of the phosphofructokinase gene family in Chinese white pear (Pyrus bretschneideri)

Phosphofructokinase plays an essential role in sugar metabolism in plants. Plants possess two types of phosphofructokinase proteins for phosphorylation of fructose-6-phosphate, the pyrophosphate-dependent fructose-6-phosphate phosphotransferase (PFP), and the ATP-dependent phosphofructokinase (PFK). Until now, the gene evolution, expression patterns, and functions of phosphofructokinase proteins were unknown in pear. In this report, 14 phosphofructokinase genes were identified in pear. The phylogenetic tree indicated that the phosphofructokinase gene family could be grouped into two subfamilies, with 10 genes belonging to the PbPFK subfamily, and 4 genes belonging to the PbPFP subfamily. Conserved motifs and exon numbers of the phosphofructokinase were found in pear and other six species. The evolution analysis indicated that WGD/Segmental and dispersed duplications were the main duplication models for the phosphofructokinase genes expansion in pear and other six species. Analysis of cis-regulatory element sequences of all phosphofructokinase genes identified light regulation and the MYB binding site in the promoter of all pear phosphofructokinase genes, suggesting that phosphofructokinase might could be regulated by light and MYB transcription factors (TFs). Gene expression patterns revealed that PbPFP1 showed similar pattern with sorbitol contents, suggesting important contributions to sugar accumulation during fruit development. Further functional analysis indicated that the phosphofructokinase gene PbPFP1 was localized on plasma membrane compartment, indicating that PbPFP1 had function in plasma membrane. Transient transformation of PbPFP1 in pear fruits led to significant increases of fructose and sorbitol compared to controls. Overall, our study provides important insights into the gene expression patterns and important potential functions of phosphofructokinase for sugar accumulation in pear fruits, which will help to enrich understanding of sugar-related bio-pathways and lay the molecular basis for fruit quality improvement.