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

Global In-Silico Scenario of tRNA Genes and Their Organization in Virus Genomes

Viruses are known to be highly dependent on the host translation machinery for their protein synthesis. However, tRNA genes are occasionally identified in such organisms, and in addition, few of them harbor tRNA gene clusters comprising dozens of genes. Recently, tRNA gene clusters have been shown to occur among the three domains of life. In such a scenario, the viruses could play a role in the dispersion of such structures among these organisms. Thus, in order to reveal the prevalence of tRNA genes as well as tRNA gene clusters in viruses, we performed an unbiased large-scale genome survey. Interestingly, tRNA genes were predicted in ssDNA (single-stranded DNA) and ssRNA (single-stranded RNA) viruses as well in many other dsDNA viruses of families from Caudovirales order. In the latter group, tRNA gene clusters composed of 15 to 37 tRNA genes were characterized, mainly in bacteriophages, enlarging the occurrence of such structures within viruses. These bacteriophages were from hosts that encompass five phyla and 34 genera. This in-silico study presents the current global scenario of tRNA genes and their organization in virus genomes, contributing and opening questions to be explored in further studies concerning the role of the translation apparatus in these organisms.

Niche Differentiation among Three Closely Related Competibacteraceae Clades at a Full-Scale Activated Sludge Wastewater Treatment Plant and Putative Linkages to Process Performance

Multiple clades within a microbial taxon often coexist within natural and engineered environments. Because closely related clades have similar metabolic potential, it is unclear how diversity is sustained and what factors drive niche differentiation. In this study, we retrieved three near-complete Competibacter lineage genomes from activated sludge metagenomes at a full-scale pure oxygen activated sludge wastewater treatment plant. The three genomes represent unique taxa within the Competibacteraceae A comparison of the genomes revealed differences in capacity for exopolysaccharide (EPS) biosynthesis, glucose fermentation to lactate, and motility. Using quantitative PCR (qPCR), we monitored these clades over a 2-year period. The clade possessing genes for motility and lacking genes for EPS biosynthesis (CPB_P15) was dominant during periods of suspended solids in the effluent. Further analysis of operational parameters indicate that the dominance of the CPB_P15 clade is associated with low-return activated sludge recycle rates and low wasting rates, conditions that maintain relatively high levels of biomass within the system.IMPORTANCE Members of the Competibacter lineage are relevant in biotechnology as glycogen-accumulating organisms (GAOs). Here, we document the presence of three Competibacteraceae clades in a full-scale activated sludge wastewater treatment plant and their linkage to specific operational conditions. We find evidence for niche differentiation among the three clades with temporal variability in clade dominance that correlates with operational changes at the treatment plant. Specifically, we observe episodic dominance of a likely motile clade during periods of elevated effluent turbidity, as well as episodic dominance of closely related nonmotile clades that likely enhance floc formation during periods of low effluent turbidity.

Genetic Evolution and Molecular Selection of the HE Gene of Influenza C Virus

Influenza C virus (ICV) was first identified in humans and swine, but recently also in cattle, indicating a wider host range and potential threat to both the livestock industry and public health than was originally anticipated. The ICV hemagglutinin-esterase (HE) glycoprotein has multiple functions in the viral replication cycle and is the major determinant of antigenicity. Here, we developed a comparative approach integrating genetics, molecular selection analysis, and structural biology to identify the codon usage and adaptive evolution of ICV. We show that ICV can be classified into six lineages, consistent with previous studies. The HE gene has a low codon usage bias, which may facilitate ICV replication by reducing competition during evolution. Natural selection, dinucleotide composition, and mutation pressure shape the codon usage patterns of the ICV HE gene, with natural selection being the most important factor. Codon adaptation index (CAI) and relative codon deoptimization index (RCDI) analysis revealed that the greatest adaption of ICV was to humans, followed by cattle and swine. Additionally, similarity index (SiD) analysis revealed that swine exerted a stronger evolutionary pressure on ICV than humans, which is considered the primary reservoir. Furthermore, a similar tendency was also observed in the M gene. Of note, we found HE residues 176, 194, and 198 to be under positive selection, which may be the result of escape from antibody responses. Our study provides useful information on the genetic evolution of ICV from a new perspective that can help devise prevention and control strategies.

Identification and expression analysis of strigolactone biosynthetic and signaling genes reveal strigolactones are involved in fruit development of the woodland strawberry (Fragaria vesca)

BACKGROUND

The development and ripening of fresh fruits is an important trait for agricultural production and fundamental research. Almost all plant hormones participate in this process. Strigolactones (SLs) are a new class of plant hormones that regulate plant organ development and stress tolerance, but little is known about their roles in fruit development.

RESULTS

In this study, we identified SL biosynthetic and signaling genes in woodland strawberry, a typical non-climacteric fruit, and analyzed the expression patterns of these genes in different plant tissues and developing fruits. One D27, two MAX1, and one LBO gene were identified as involved in SL biosynthesis, and one D14, one D3, and two D53 genes as related to SL signaling. The proteins encoded by these genes had similar motifs as SL biosynthetic and signaling proteins in rice and Arabidopsis. The genes had different expression levels in the root, stem, leaf, and petiole of woodland strawberry. In addition, the expression of most SL biosynthetic genes was high in developing carpel, anther, and style, while that of SL signaling genes was high in carpel and style, but low in anther, suggesting active SL biosynthesis and signaling in the developing carpel and style. Notably, the expression of SL biosynthetic and signaling genes was significantly increased in the receptacle after pollination and decreased during receptacle development. Moreover, low or no expression of these genes was detected in ripening fruits.

CONCLUSIONS

Our results suggest that SLs play a role in the early stages of woodland strawberry fruit development. Our findings provide insight into the function of SLs and will facilitate further study of the regulation by SLs of fresh fruit development.

Molecular characterization, cellular localization, and light-enhanced expression of Beta-Na+/H+ Exchanger-like in the whitish inner mantle of the giant clam, Tridacna squamosa, denote its role in light-enhanced shell formation

The fluted giant clam, Tridacna squamosa, lives in symbiosis with photosynthetic zooxanthellae, and can engage in light-enhanced growth and shell formation. Light-enhanced shell formation necessitates the elimination of excess H⁺ from the extrapallial fluid adjacent to the shell. This study aimed to clone Na⁺/H⁺Exchanger (NHE) from the whitish inner mantle adjacent to the extrapallial fluid of T. squamosa, to determine its cellular and subcellular localization, and to evaluate the effect of light exposure on its mRNA expression level and protein abundance therein. The complete coding cDNA sequence of NHE obtained was identified as a homolog of beta NHE (βNHE-like). It consisted of 2925 bp, encoding for a polypeptide of 974 amino acids and 107.1 kDa, and was expressed predominantly in the inner mantle. There, βNHE-like was localized in the apical membrane of the seawater-facing epithelium by immunofluorescence microscopy. After exposure to light for 12 h, the seawater-facing epithelium of the inner mantle displayed consistently stronger immunostaining than that of the control exposed to 12 h of darkness. Western blotting confirmed that light exposure significantly enhanced the protein abundance of βNHE-like in the inner mantle. These results denote that some of the excess H⁺ generated during light-enhanced shell formation can be excreted through the light-dependent βNHE-like of the seawater-facing epithelium to minimize the impact on the whole-body pH. Importantly, the excreted H⁺ could dehydrate exogenous HCO₃^-, and facilitate the absorption of inorganic carbon through the seawater-facing epithelium dedicated for light-enhanced shell formation due to its close proximity with the shell-facing epithelium. NUCLEOTIDE SYMBOL COMBINATIONS: Pairs: R = A/G; W = A/T; Y = C/T. Triples: D = A/G/T.

Mitogenomic evidence of close relationships between New Zealand's extinct giant raptors and small-sized Australian sister-taxa

Prior to human arrival in the 13th century, two large birds of prey were the top predators in New Zealand. In the absence of non-volant mammals, the extinct Haast's eagle (Hieraaetus moorei), the largest eagle in the world, and the extinct Eyles' harrier (Circus teauteensis) the largest harrier in the world, had filled ecological niches that are on other landmasses occupied by animals such as large cats or canines. The evolutionary and biogeographic history of these island giants has long been a mystery. Here we reconstruct the origin and evolution of New Zealand's giant raptors using complete mitochondrial genome data. We show that both Eyles' harrier and Haast's eagle diverged from much smaller, open land adapted Australasian relatives in the late Pliocene to early Pleistocene. These events coincided with the development of open habitat in the previously densely forested islands of New Zealand. Our study provides evidence of rapid evolution of island gigantism in New Zealand's extinct birds of prey. Early Pleistocene climate and environmental changes were likely to have triggered the establishment of Australian raptors into New Zealand. Our results shed light on the evolution of two of the most impressive cases of island gigantism in the world.

A novel polyketide synthase gene cluster in the plant pathogenic fungus Pseudocercospora fijiensis

Pseudocercospora fijiensis, causal agent of black Sigatoka of banana, produces polyketide synthase (PKS) pathways shown to be important in disease development by related Dothideomycete fungi. Genome analysis of the P. fijiensis PKS8-1 gene identified it as part of a gene cluster including genes encoding two transcription factors, a regulatory protein, a glyoxylase/beta-lactamase-like protein, an MFS transporter, a cytochrome P450, two aldo/keto reductases, a dehydrogenase, and a decarboxylase. Genome analysis of the related pathogens Pseudocercospora musae, Pseudocercospora eumusae, and Pseudocercospora pini-densiflorae, identified orthologous clusters containing a nearly identical combination of genes. Phylogenetic analysis of PKS8-1 identified homology to PKS proteins in the monodictyphenone and cladofulvin pathways in Aspergillus nidulans and Cladosporium fulvum, respectively. Analysis of clustered genes showed that the PKS8-1 cluster shares genes for enzymes involved in the production of the emodin intermediate in the monodictyphenone and cladofulvin pathways, but differs in many genes, suggesting production of a different metabolic product. Time course analysis of gene expression in infected banana showed up-regulation of PKS8-1 and four of eight clustered genes as early as 2 weeks post-inoculation and remaining high through 9 weeks. Overexpression of the pathway through constitutive expression of an aflR-like transcription factor gene in the cluster resulted in increased expression in culture of PKS8-1 as well as the four clustered genes that are up-regulated in infected plants. No differences were seen in timing or severity of disease symptoms with the overexpression strains relative to controls, however gene expression analysis showed no difference in expression in planta by an overexpression strain relative to controls. Thus constitutive expression of the aflR-like gene is not sufficient to upregulate the pathway above normal expression in planta. Pathway expression during all phases of disease development and conservation of the pathway in related Pseudocercospora species support a role for this pathway in disease.

Evolution of the modular, disordered stress proteins known as dehydrins

Dehydrins, plant proteins that are upregulated during dehydration stress conditions, have modular sequences that can contain three conserved motifs (the Y-, S-, and K-segments). The presence and order of these motifs are used to classify dehydrins into one of five architectures: Kn, SKn, KnS, YnKn, and YnSKn, where the subscript n describes the number of copies of that motif. In this study, an architectural and phylogenetic analysis was performed on 426 dehydrin sequences that were identified in 53 angiosperm and 3 gymnosperm genomes. It was found that angiosperms contained all five architectures, while gymnosperms only contained Kn and SKn dehydrins. This suggests that the ancestral dehydrin in spermatophytes was either Kn or SKn, and the Y-segment containing dehydrins first arose in angiosperms. A high-level split between the YnSKn dehydrins from either the Kn or SKn dehydrins could not be confidently identified, however, two lower level architectural divisions appear to have occurred after different duplication events. The first likely occurred after a whole genome duplication, resulting in the duplication of a Y3SK2 dehydrin; the duplicate subsequently lost an S- and K- segment to become a Y3K1 dehydrin. The second split occurred after a tandem duplication of a Y1SK2 dehydrin, where the duplicate lost both the Y- and S- segment and gained four K-segments, resulting in a K6 dehydrin. We suggest that the newly arisen Y3K1 dehydrin is possibly on its way to pseudogenization, while the newly arisen K6 dehydrin developed a novel function in cold protection.

Genomes and seroprevalence of severe fever with thrombocytopenia syndrome virus and Nairobi sheep disease virus in Haemaphysalis longicornis ticks and goats in Hubei, China

Ticks are medically-important arthropods that maintain and transmit numerous emerging viruses. China suffers severely from tick-borne viral diseases such as tick-borne encephalitis and severe fever with thrombocytopenia syndrome (SFTS), but the background of tick-borne viruses is very limited. Here we report the virome profiling of ticks and goat sera from SFTS-epidemic areas, and serological investigation of SFTS virus (SFTSV) and Nairobi sheep disease virus (NSDV). Results revealed divergent viruses in ticks and goat sera, including SFTSV and NSDV. Sequence and phylogenetic analyses showed that the SFTSV identified here was most closely related to human SFTSV in sampling and surrounding areas, and the NSDV to the previously identified NSDV from northeast China. Serological investigation of SFTSV infection in goats revealed intensive activity in those areas. Surprisingly, two different methods of NSDV serological investigation showed no sera positive for this virus.

Phylogenetic relationship of prophages is affected by CRISPR selection in Group A Streptococcus

BACKGROUND

Group A Streptococcus (GAS) is a major human pathogen, which is associated with a wide spectrum of invasive diseases, such as pharyngitis, scarlet fever, rheumatic fever, and streptococcal toxic shock syndrome (STSS). It is hypothesized that differences in GAS pathogenicity are related to the acquisition of diverse bacteriophages (phages). Nevertheless, the GAS genome also harbors clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated (cas) genes, which play an important role in eliminating foreign DNA, including those of phages. However, the structure of prophages in GAS strains is mosaic, and the phylogenetic relationship between prophages and CRISPR is not clear. In this study, we analyzed CRISPR and prophage structure using 118 complete genome sequences of GAS strains to elucidate the relationship between two genomic elements. Additionally, phylogenetic and M-type analyses were performed.

RESULTS

Of the 118 GAS strains, 80 harbored type I-C and/or II-A CRISPR/cas loci. A total of 553 spacer sequences were identified from CRISPR/cas loci and sorted into 229 patterns. We identified and classified 373 prophages into 14 groups. Some prophage groups shared a common integration site, and were related to M-type. We further investigated the correlation between spacer sequences and prophages. Of the 229 spacer sequence patterns, 203 were similar to that of other GAS prophages. No spacer showed similarity with that of a specific prophage group with mutL integration site. Moreover, the average number of prophages in strains with type II-A CRISPR was significantly less than that in type I-C CRISPR and non-CRISPR strains. However, there was no statistical difference between the average number of prophages in type I-C strains and that in non-CRISPR strains.

CONCLUSIONS

Our results indicated that type II-A CRISPR may play an important role in eliminating phages and that the prophage integration site may be an important criterion for the acceptance of foreign DNA by GAS. M type, spacer sequence, and prophage group data were correlated with the phylogenetic relationships of GAS. Therefore, we hypothesize that genetic characteristics and/or phylogenetic relationships of GAS may be estimated by analyzing its spacer sequences.

Gut Virome Analysis of Cameroonians Reveals High Diversity of Enteric Viruses, Including Potential Interspecies Transmitted Viruses

Diarrhea remains one of the most common causes of deaths in children. A limited number of studies have investigated the prevalence of enteric pathogens in Cameroon, and as in many other African countries, the cause of many diarrheal episodes remains unexplained. A proportion of these unknown cases of diarrhea are likely caused by yet-unidentified viral agents, some of which could be the result of (recent) interspecies transmission from animal reservoirs, like bats. Using viral metagenomics, we screened fecal samples of 221 humans (almost all with gastroenteritis symptoms) between 0 and 89 years of age with different degrees of bat contact. We identified viruses belonging to families that are known to cause gastroenteritis such as Adenoviridae, Astroviridae, Caliciviridae, Picornaviridae, and Reoviridae Interestingly, a mammalian orthoreovirus, picobirnaviruses, a smacovirus, and a pecovirus were also found. Although there was no evidence of interspecies transmission of the most common human gastroenteritis-related viruses (Astroviridae, Caliciviridae, and Reoviridae), the phylogenies of the identified orthoreovirus, picobirnavirus, and smacovirus indicate a genetic relatedness of these viruses identified in stools of humans and those of bats and/or other animals. These findings points out the possibility of interspecies transmission or simply a shared host of these viruses (bacterial, fungal, parasitic, …) present in both animals (bats) and humans. Further screening of bat viruses in humans or vice versa will elucidate the epidemiological potential threats of animal viruses to human health. Furthermore, this study showed a huge diversity of highly divergent novel phages, thereby expanding the existing phageome considerably.IMPORTANCE Despite the availability of diagnostic tools for different enteric viral pathogens, a large fraction of human cases of gastroenteritis remains unexplained. This could be due to pathogens not tested for or novel divergent viruses of potential animal origin. Fecal virome analyses of Cameroonians showed a very diverse group of viruses, some of which are genetically related to those identified in animals. This is the first attempt to describe the gut virome of humans from Cameroon. Therefore, the data represent a baseline for future studies on enteric viral pathogens in this area and contribute to our knowledge of the world's virome. The studies also highlight the fact that more viruses may be associated with diarrhea than the typical known ones. Hence, it provides meaningful epidemiological information on diarrhea-related viruses in this area.

Phylogeny and species delimitation of near Eastern Neurergus newts (Salamandridae) based on genome-wide RADseq data analysis

We reconstruct the molecular phylogeny of Near Eastern mountain brook newts of the genus Neurergus (family Salamandridae) based on newly determined RADseq data, and compare the outcomes of concatenation-based phylogenetic reconstruction with species-tree inference. Furthermore, we test the current taxonomy of Neurergus (with four species: Neurergus strauchii, N. crocatus, N. kaiseri, and N. derjugini) against coalescent-based species-delimitation approaches of our genome-wide genetic data set. While the position of N. strauchii as sister species to all other Neurergus species was consistent in all of our analyses, the phylogenetic relationships between the three remaining species changed depending on the applied method. The concatenation approach, as well as quartet-based species-tree inference, supported a topology with N. kaiseri as the closest relative to N. derjugini, while full-coalescent species-tree inference approaches supported N. crocatus as sister species of N. derjugini. Investigating the individual signal of gene trees highlighted an extensive variation among gene histories, most likely resulting from incomplete lineage sorting. Coalescent-based species-delimitation models suggest that the current taxonomy might underestimate the species richness within Neurergus and supports seven species. Based on the current sampling, our analysis suggests that N. strauchii, N. derjugini and N. kaiseri might each be subdivided into further species. However, as amphibian species are known to be composed of deep conspecific lineages that do not always warrant species status, these results need to be cautiously interpreted in an integrative taxonomic framework. We hypothesize that the rather shallow divergences detected within N. kaiseri and N. derjugini likely reflect an ongoing speciation process and thus require further investigation. On the contrary, the much deeper genetic divergence found between the two morphologically and geographically differentiated subspecies of N. strauchii leads us to propose that N. s. barani should be considered a distinct species, Neurergus barani Öz, 1994.

Functional and structural profiles of GST gene family from three Populus species reveal the sequence-function decoupling of orthologous genes

A common assumption in comparative genomics is that orthologous genes are functionally more similar than paralogous genes. However, the validity of this assumption needs to be assessed using robust experimental data. We conducted tissue-specific gene expression and protein function analyses of orthologous groups within the glutathione S-transferase (GST) gene family in three closely related Populus species: Populus trichocarpa, Populus euphratica and Populus yatungensis. This study identified 21 GST orthologous groups in the three Populus species. Although the sequences of the GST orthologous groups were highly conserved, the divergence in enzymatic functions was prevalent. Through site-directed mutagenesis of orthologous proteins, this study revealed that nonsynonymous substitutions at key amino acid sites played an important role in the divergence of enzymatic functions. In particular, a single amino acid mutation (Arg39→Trp39) contributed to P. euphratica PeGSTU30 possessing high enzymatic activity via increasing the hydrophobicity of the active cavity. This study provided experimental evidence showing that orthologues belonging to the gene family have functional divergences. The nonsynonymous substitutions at a few amino acid sites resulted in functional divergence of the orthologous genes. Our findings provide new insights into the evolution of orthologous genes in closely related species.

Diversity of Biomphalaria spp. freshwater snails and associated mollusks in areas with schistosomiasis risk, using molecular and spatial analysis tools

Abstract: The Middle Paranapanema River region of São Paulo, Brazil is home to significant diversity of Biomphalaria species and is very vulnerable to health and environmental impacts such as schistosomiasis. This study updates freshwater malacological surveys for ecosystems in one portion of the Middle Paranapanema River Basin, with emphasis on the genus Biomphalaria. Snails were collected from 114 distinct bodies of water between 2015 and 2018. Biomphalaria specimens were identified according to morphological and molecular characteristics, while animals in other genera (Drepanotrema, Lymnaea, Melanoides, Physa and Pomacea) were identified solely according to shell characteristics. A geographic information system was used to update intermediate host colonization sites and consequently assist in identifying probable hotspots for intermediate hosts of schistosomiasis. The sequences of the COI gene relating to the DNA barcode stretch were tested for similarity against sequences found in GenBank, for monophyly through Maximum Likelihood phylogenetic inference, and analyzed in ABDG, bPTP and GMYC for the delimitation of putative species. Of the 10,722 snails collected, 86.7% were in the Planorbidae family (75.5% Biomphalaria and 11.2% Drepanotrema) and 13.3% were other non-Planorbidae species (Lymnaea, Melanoides, Physa and Pomacea). The taxonomic COI reference sequences in the NCBI nucleotide database used for DNA sequence comparison, and phylogenetic analysis used to test the monophyly of the groups, resulted in more reliable taxonomic units than delimitation of the COI sequences in MOTUs using statistical taxonomic models. Analysis of the species distribution shows that B. glabrata and B. tenagophila are heterogeneously distributed in the study area. B. glabrata colonizes only five water bodies, in the study area, most of them in Ourinhos, while B. tenagophila predominates in water bodies in Ipaussu. Contrasting with this, B. straminea, B. occidentalis and B. peregrina are evenly distributed throughout the study area.

Assignment of a dubious gene cluster to melanin biosynthesis in the tomato fungal pathogen Cladosporium fulvum

Pigments and phytotoxins are crucial for the survival and spread of plant pathogenic fungi. The genome of the tomato biotrophic fungal pathogen Cladosporium fulvum contains a predicted gene cluster (CfPKS1, CfPRF1, CfRDT1 and CfTSF1) that is syntenic with the characterized elsinochrome toxin gene cluster in the citrus pathogen Elsinoë fawcettii. However, a previous phylogenetic analysis suggested that CfPks1 might instead be involved in pigment production. Here, we report the characterization of the CfPKS1 gene cluster to resolve this ambiguity. Activation of the regulator CfTSF1 specifically induced the expression of CfPKS1 and CfRDT1, but not of CfPRF1. These co-regulated genes that define the CfPKS1 gene cluster are orthologous to genes involved in 1,3-dihydroxynaphthalene (DHN) melanin biosynthesis in other fungi. Heterologous expression of CfPKS1 in Aspergillus oryzae yielded 1,3,6,8-tetrahydroxynaphthalene, a typical precursor of DHN melanin. Δcfpks1 deletion mutants showed similar altered pigmentation to wild type treated with DHN melanin inhibitors. These mutants remained virulent on tomato, showing this gene cluster is not involved in pathogenicity. Altogether, our results showed that the CfPKS1 gene cluster is involved in the production of DHN melanin and suggests that elsinochrome production in E. fawcettii likely involves another gene cluster.

Two more Posterior Hox genes and Hox cluster dispersal in echinoderms

BACKGROUND

Hox genes are key elements in patterning animal development. They are renowned for their, often, clustered organisation in the genome, with supposed mechanistic links between the organisation of the genes and their expression. The widespread distribution and comparable functions of Hox genes across the animals has led to them being a major study system for comparing the molecular bases for construction and divergence of animal morphologies. Echinoderms (including sea urchins, sea stars, sea cucumbers, feather stars and brittle stars) possess one of the most unusual body plans in the animal kingdom with pronounced pentameral symmetry in the adults. Consequently, much interest has focused on their development, evolution and the role of the Hox genes in these processes. In this context, the organisation of echinoderm Hox gene clusters is distinctive. Within the classificatory system of Duboule, echinoderms constitute one of the clearest examples of Disorganized (D) clusters (i.e. intact clusters but with a gene order or orientation rearranged relative to the ancestral state).

RESULTS

Here we describe two Hox genes (Hox11/13d and e) that have been overlooked in most previous work and have not been considered in reconstructions of echinoderm Hox complements and cluster organisation. The two genes are related to Posterior Hox genes and are present in all classes of echinoderm. Importantly, they do not reside in the Hox cluster of any species for which genomic linkage data is available.

CONCLUSION

Incorporating the two neglected Posterior Hox genes into assessments of echinoderm Hox gene complements and organisation shows that these animals in fact have Split (S) Hox clusters rather than simply Disorganized (D) clusters within the Duboule classification scheme. This then has implications for how these genes are likely regulated, with them no longer covered by any potential long-range Hox cluster-wide, or multigenic sub-cluster, regulatory mechanisms.

Regeneration of the acorn worm pygochord with the implication for its convergent evolution with the notochord

The origin of the notochord is a central issue in chordate evolution. This study examined the development of the acorn worm pygochord, a putative homologue of the notochord. Because the pygochord differentiates only after metamorphosis, the developmental was followed process by inducing regeneration after artificial amputation in Ptychodera flava. It was found that although the regeneration of the posterior part of the body did not proceed via formation of an obvious regeneration bud, pygochord regeneration was observed within a few weeks, possibly via trans-differentiation of endoderm cells. The expression of the fibrillary collagen gene (Fcol) and elav in the pygochord during regeneration was detected. This indicates that pygochord cells are not part of gut epithelial cells, but that they differentiated as a distinct cell type. Our gene expression analyses do not provide supporting evidence for the homology between the pygochord and notochord, but rather favored the convergent evolution between them.

Identification and partial characterization of a novel serpin from Eudiplozoon nipponicum (Monogenea, Polyopisthocotylea)

Background: Serpins are a superfamily of serine peptidase inhibitors that participate in the regulation of many physiological and cell peptidase-mediated processes in all organisms (e.g. in blood clotting, complement activation, fibrinolysis, inflammation, and programmed cell death). It was postulated that in the blood-feeding members of the monogenean family Diplozoidae, serpins could play an important role in the prevention of thrombus formation, activation of complement, inflammation in the host, and/or in the endogenous regulation of protein degradation.

Results: In silico analysis showed that the DNA and primary protein structures of serpin from Eudiplozoon nipponicum (EnSerp1) are similar to other members of the serpin superfamily. The inhibitory potential of EnSerp1 on four physiologically-relevant serine peptidases (trypsin, factor Xa, kallikrein, and plasmin) was demonstrated and its presence in the worm’s excretory-secretory products (ESPs) was confirmed.

Conclusion: EnSerp1 influences the activity of peptidases that play a role in blood coagulation, fibrinolysis, and complement activation. This inhibitory potential, together with the serpin’s presence in ESPs, suggests that it is likely involved in host-parasite interactions and could be one of the molecules involved in the control of feeding and prevention of inflammatory responses.

Polychromophilus melanipherus and haemoplasma infections not associated with clinical signs in southern bent-winged bats (Miniopterus orianae bassanii) and eastern bent-winged bats (Miniopterus orianae oceanensis)

While bats are often viewed as carriers of infectious disease agents, little research has been conducted on the effects these pathogens may have on the bat populations themselves. The southern bent-winged bat (Miniopterus orianae bassanii) is a critically endangered subspecies endemic to south-eastern Australia. Population numbers of this bat have declined over the past 50 years, but the reasons for this are unclear. As part of a larger study to determine if disease could be a contributing factor to this decline, southern bent-winged bats from several locations in Victoria and South Australia were captured and examined for the presence of the blood parasite, Polychromophilus melanipherus, and haemoplasmas (Mycoplasma sp.). Results were compared with those obtained from populations of the more common, partially sympatric, eastern bent-winged bat (Miniopterus orianae oceanensis) from three different locations in Victoria. Both organisms were found in both subspecies (prevalence of P. melanipherus 60% by PCR for southern bent-winged bats compared with 46% for eastern bent-winged bats; prevalence of haemoplasmas 10% for southern bent-winged bats compared with 8% for eastern bent-winged bats), with no association between the probability of infection, body weight, abnormal blood parameters or any other indicators of ill health. However, Victorian southern bent-winged bats had heavier burdens of P. melanipherus than both the South Australian southern bent-winged bats and eastern bent-winged bats. Further investigations are required to determine if these differences are impacting population health.

Enrichment and Genomic Characterization of a N2O-Reducing Chemolithoautotroph From a Deep-Sea Hydrothermal Vent

Nitrous oxide (N₂O) is a greenhouse gas and also leads to stratospheric ozone depletion. In natural environments, only a single N₂O sink process is the microbial reduction of N₂O to N₂, which is mediated by nitrous oxide reductase (NosZ) encoded by nosZ gene. The nosZ phylogeny has two distinct clades, clade I and formerly overlooked clade II. In deep-sea hydrothermal environments, several members of the class Campylobacteria are shown to harbor clade II nosZ gene and perform the complete denitrification of nitrate to N₂; however, little is known about their ability to grow on exogenous N₂O as the sole electron acceptor. Here, we obtained an enrichment culture from a deep-sea hydrothermal vent in the Southern Mariana Trough, which showed a respiratory N₂O reduction with H₂ as an electron donor. The single amplicon sequence variant (ASV) presenting 90% similarity to Hydrogenimonas species within the class Campylobacteria was predominant throughout the cultivation period. Metagenomic analyses using a combination of short-read and long-read sequence data succeeded in reconstructing a complete genome of the dominant ASV, which encoded clade II nosZ gene. This study represents the first cultivation analysis that shows the occurrence of N₂O-respiring microorganisms in a deep-sea hydrothermal vent and provides the opportunity to assess their capability to reduce N₂O emission from the environments.

Identification of insect-specific flaviviruses in areas of Brazil and Paraguay experiencing endemic arbovirus transmission and the description of a novel flavivirus infecting Sabethes belisarioi

Viral infection was examined with pan-flavivirus and pan-alphavirus sets of primers in mosquitoes collected in four South American regions with confirmed pathogenic arbovirus circulation. Positive pools for flavivirus infection were sequenced and screened for specific arboviruses, which were not detected. However, NS5 gene sequencing showed that most sequences corresponded to the insect-specific Culex flavivirus. One sequence retrieved from an Aedes albopictus pool grouped with the insect-specific Aedes flavivirus and two Sabethes belisarioi pools were infected by a previously unknown flavivirus, tentatively named Sabethes flavivirus (SbFV). Phylogenetic inference placed SbFV as ancestral to a clade formed by Culiseta flavivirus, Mercadeo, and Calbertado. SbFV polyprotein showed an average aminoacidic identity of 51% in comparison to these flaviviruses. In vitro studies suggest that SbFV infects insect cells, but not vertebrate cells, therefore, we propose it as a new insect-specific flavivirus. These results highlight the wide distribution of insect-specific flaviviruses concomitant with the circulation of emergent arboviruses.

Deciphering the evolutionary signatures of pinnipeds using novel genome sequences: The first genomes of Phoca largha, Callorhinus ursinus, and Eumetopias jubatus

The pinnipeds, which comprise seals, sea lions, and walruses, are a remarkable group of marine animals with unique adaptations to semi-aquatic life. However, their genomes are poorly characterized. In this study, we sequenced and characterized the genomes of three pinnipeds (Phoca largha, Callorhinus ursinus, and Eumetopias jubatus), focusing on site-wise sequence changes. We detected rapidly evolving genes in pinniped lineages and substitutions unique to pinnipeds associated with amphibious sound perception. Phenotypic convergence-related sequence convergences are not common in marine mammals. For example, FASN, KCNA5, and IL17RA contain substitutions specific to pinnipeds, yet are potential candidates of phenotypic convergence (blubber, response to hypoxia, and immunity to pathogens) in all marine mammals. The outcomes of this study will provide insight into targets for future studies of convergent evolution or gene function.

Delineation of the Genera Haemoproteus and Plasmodium Using RNA-Seq and Multi-gene Phylogenetics

Members of the order Haemosporida are protist parasites that infect mammals, reptiles and birds. This group includes the causal agents of malaria, Plasmodium parasites, the genera Leucocytozoon and Fallisia, as well as the species rich genus Haemoproteus with its two subgenera Haemoproteus and Parahaemoproteus. Some species of Haemoproteus cause severe disease in avian hosts, and these parasites display high levels of diversity worldwide. This diversity emphasizes the need for accurate evolutionary information. Most molecular studies of wildlife haemosporidians use a bar coding approach by sequencing a fragment of the mitochondrial cytochrome b gene. This method is efficient at differentiating parasite lineages but insufficient for accurate phylogenetic inferences in highly diverse taxa such as haemosporidians. Recent studies have utilized multiple mitochondrial genes (cyt b, cox1 and cox3), sometimes combined with a few apicoplast and nuclear genes. These studies have been highly successful with one notable exception: the evolutionary relationships of the genus Haemoproteus remain unresolved. Here we describe the transcriptome of Haemoproteus columbae and investigate its phylogenetic position recovered from a multi-gene dataset (600 genes). This genomic approach restricts the taxon sampling to 18 species of apicomplexan parasites. We employed Bayesian inference and maximum likelihood methods of phylogenetic analyses and found H. columbae and a representative from the subgenus Parahaemoproteus to be sister taxa. This result strengthens the hypothesis of genus Haemoproteus being monophyletic; however, resolving this question will require sequences of orthologs from, in particular, representatives of Leucocytozoon species.

Molecular phylogeny of Caudofoveata (Mollusca) challenges traditional views

The shell-less, worm-shaped Caudofoveata (=Chaetodermomorpha) is one of the least known groups of molluscs. The taxon consists of 141 recognized species found from intertidal environments to the deep-sea where they live burrowing in sediment. Evolutionary relationships of the group have been debated, but few studies based on morphological or molecular data have investigated the phylogeny of the group. Here we use molecular phylogenetics to resolve relationships among and within families of Caudofoveata. Phylogenetic analyses were performed using selected mitochondrial and nuclear genes from species from all recognized families of Caudofoveata. In resulting trees and contrary to traditional views, Prochaetodermatidae forms the sister clade to a clade containing the other two currently recognized families, Chaetodermatidae and Limifossoridae. The monophyly of Prochaetodermatidae is highly supported, but Limifossoridae and Chaetodermatidae are not recovered as monophyletic. Most of the caudofoveate genera are also not recovered as monophyletic in our analyses. Thus results from our molecular data suggest that the current classification of Caudofoveata is in need of revision, and indicate evolutionary scenarios that differ from previously proposed hypotheses based on morphology.

Hidden diversity and evolution of viruses in market fish

Aquaculture is the fastest growing industry worldwide. Aquatic diseases have had enormous economic and environmental impacts in the recent past and the emergence of new aquatic pathogens, particularly viruses, poses a continuous threat. Nevertheless, little is known about the diversity, abundance and evolution of fish viruses. We used a meta-transcriptomic approach to help determine the virome of seemingly healthy fish sold at a market in Sydney, Australia. Specifically, by identifying and quantifying virus transcripts we aimed to determine (i) the abundance of viruses in market fish, (ii) test a key component of epidemiological theory that large and dense host populations harbour a greater number of viruses compared to their more solitary counterparts and (iii) reveal the relative roles of virus–host co-divergence and cross-species transmission in the evolution of fish viruses. The species studied comprised both shoaling fish—eastern sea garfish (Hyporhamphus australis) and Australasian snapper (Chrysophrys auratus)—and more solitary fish—eastern red scorpionfish (Scorpaena jacksoniensis) and largetooth flounder (Pseudorhombus arsius). Our analysis identified twelve potentially novel viruses, eight of which were likely vertebrate-associated across four viral families and that exhibited frequent cross-species transmission. Notably, the most solitary of the fish species studied, the largetooth flounder, harboured the least number of viruses while eastern sea garfish, a densely shoaling fish, had the highest number of viruses. These results support the emerging view that fish harbour a large and largely uncharacterised virome.

First evidence of polychaete intermediate hosts for Neospirorchis spp. marine turtle blood flukes (Trematoda: Spirorchiidae)

Life cycles of spirorchiids that infect the vascular system of turtles are poorly understood. Few life cycles of these blood flukes have been elucidated and all intermediate hosts reported are gastropods (Mollusca), regardless of whether the definitive host is a freshwater or a marine turtle. During a recent survey of blood fluke larvae in polychaetes on the coast of South Carolina, USA, spirorchiid-like cercariae were found to infect the polychaetes Amphitrite ornata (Terebellidae) and Enoplobranchus sanguineus (Polycirridae). Cercariae were large, furcate, with a ventral acetabulum, but no eyespots were observed. Partial sequences of D1-D2 domains of the large ribosomal subunit, the internal transcribed spacer 2, and the mitochondrial cytochrome oxidase 1 genes allowed the identification of sporocysts and cercariae as belonging to two unidentified Neospirorchis species reported from the green turtle, Chelonia mydas, in Florida: Neospirorchis sp. (Neogen 13) in A. ornata and Neospirorchis sp. (Neogen 14) in E. sanguineus. Phylogenetic analysis suggests that infection of annelids by blood flukes evolved separately in aporocotylids and spirorchiids. Our results support the contention that the Spirorchiidae is not a valid family and suggest that Neospirorchis is a monophyletic clade within the paraphyletic Spirorchiidae. Since specificity of spirorchiids for their intermediate hosts is broader than it was thus far assumed, surveys of annelids in turtle habitats are necessary to further our understanding of the life history of these pathogenic parasites.

Complex virome in feces from Amerindian children in isolated Amazonian villages

The number of viruses circulating in small isolated human populations may be reduced by viral extinctions and rare introductions. Here we used viral metagenomics to characterize the eukaryotic virome in feces from healthy children from a large urban center and from three Amerindian villages with minimal outside contact. Numerous human enteric viruses, mainly from the Picornaviridae and Caliciviridae families, were sequenced from each of the sites. Multiple children from the same villages shed closely related viruses reflecting frequent transmission clusters. Feces of isolated villagers also contained multiple viral genomes of unknown cellular origin from the Picornavirales order and CRESS-DNA group and higher levels of nematode and protozoan DNA. Despite cultural and geographic isolation, the diversity of enteric human viruses was therefore not reduced in these Amazonian villages. Frequent viral introductions and/or increased susceptibility to enteric infections may account for the complex fecal virome of Amerindian children in isolated villages.

The first complete mitochondrial genome from the family Hygrophoraceae (Hygrophorus russula) by next-generation sequencing and phylogenetic implications

Hygrophorus russula (Schaeff.) Kauffman is an edible ectomycorrhizal fungus that is widely distributed in the world. In this study, the mitogenome of H. russula was sequenced and assembled. The mitogenome of H. russula is composed of circular DNA molecules, with a total size of 55,769 bp. Further analysis indicated that the frequent use of A and T in codons contributes to the high AT content (80.87%) in the H. russula mitogenome. Comparative analysis indicated that the length and base composition of the core protein-encoding genes, and the number of tRNA genes in the H. russula mitogenome varied from that of other Agaricales mitogenomes. Gene arrangement analysis revealed a novel gene order in the H. russula mitogenome. In addition, the expansion of the mitogenome in Agaricales was found to be closely related to the increase in the number of introns. Phylogenetic analysis of the combined mitochondrial gene set showed strong support for tree topologies, and H. russula was determined to be relatively distant from other Agaricales species. This study is the first report on the mitogenome of a member of genus Hygrophorus as well as family Hygrophoraceae, which improves our understanding of mitochondrial differentiation and evolution in the important ectomycorrhizal fungi Hygrophorus species.

BioWorkbench: a high-performance framework for managing and analyzing bioinformatics experiments

Advances in sequencing techniques have led to exponential growth in biological data, demanding the development of large-scale bioinformatics experiments. Because these experiments are computation- and data-intensive, they require high-performance computing techniques and can benefit from specialized technologies such as Scientific Workflow Management Systems and databases. In this work, we present BioWorkbench, a framework for managing and analyzing bioinformatics experiments. This framework automatically collects provenance data, including both performance data from workflow execution and data from the scientific domain of the workflow application. Provenance data can be analyzed through a web application that abstracts a set of queries to the provenance database, simplifying access to provenance information. We evaluate BioWorkbench using three case studies: SwiftPhylo, a phylogenetic tree assembly workflow; SwiftGECKO, a comparative genomics workflow; and RASflow, a RASopathy analysis workflow. We analyze each workflow from both computational and scientific domain perspectives, by using queries to a provenance and annotation database. Some of these queries are available as a pre-built feature of the BioWorkbench web application. Through the provenance data, we show that the framework is scalable and achieves high-performance, reducing up to 98% of the case studies execution time. We also show how the application of machine learning techniques can enrich the analysis process.

PIWI genes and piRNAs are ubiquitously expressed in mollusks and show patterns of lineage-specific adaptation

PIWI proteins and PIWI-interacting RNAs (piRNAs) suppress transposon activity in animals, thus protecting their genomes from detrimental insertion mutagenesis. Here, we reveal that PIWI genes and piRNAs are ubiquitously expressed in mollusks, similar to the situation in arthropods. We describe lineage-specific adaptations of transposon composition in piRNA clusters in the great pond snail and the pacific oyster, likely reflecting differential transposon activity in gastropods and bivalves. We further show that different piRNA clusters with unique transposon composition are dynamically expressed during oyster development. Finally, bioinformatics analyses suggest that different populations of piRNAs presumably bound to different PIWI paralogs participate in homotypic and heterotypic ping-pong amplification loops in a tissue- and sex-specific manner. Together with recent findings from other animal species, our results support the idea that somatic piRNA expression represents the ancestral state in metazoans.

A Revised Spiralian Homeobox Gene Classification Incorporating New Polychaete Transcriptomes Reveals a Diverse TALE Class and a Divergent Hox Gene

The diversity of mechanisms and capacity for regeneration across the Metazoa present an intriguing challenge in evolutionary biology, impacting on the burgeoning field of regenerative medicine. Broad taxonomic sampling is essential to improve our understanding of regeneration, and studies outside of the traditional model organisms have proved extremely informative. Within the historically understudied Spiralia, the Annelida have an impressive variety of tractable regenerative systems. The biomeralizing, blastema-less regeneration of the head appendage (operculum) of the serpulid polychaete keelworm Spirobranchus (formerly Pomatoceros) lamarcki is one such system. To profile potential regulatory mechanisms, we classified the homeobox gene content of opercular regeneration transcriptomes. As a result of retrieving several difficult-to-classify homeobox sequences, we performed an extensive search and phylogenetic analysis of the TALE and PRD-class homeobox gene content of a broad selection of lophotrochozoan genomes. These analyses contribute to our increasing understanding of the diversity, taxonomic extent, rapid evolution, and radical flexibility of these recently discovered homeobox gene radiations. Our expansion and integration of previous nomenclature systems helps to clarify their cryptic orthology. We also describe an unusual divergent S. lamarcki Antp gene, a previously unclassified lophotrochozoan orphan gene family (Lopx), and a number of novel Nk class orphan genes. The expression and potential involvement of many of these lineage- and clade-restricted homeobox genes in S. lamarcki operculum regeneration provides an example of diversity in regenerative mechanisms, as well as significantly improving our understanding of homeobox gene evolution.

Relative Evolutionary Rates in Proteins Are Largely Insensitive to the Substitution Model

The relative evolutionary rates at individual sites in proteins are informative measures of conservation or adaptation. Often used as evolutionarily aware conservation scores, relative rates reveal key functional or strongly selected residues. Estimating rates in a phylogenetic context requires specifying a protein substitution model, which is typically a phenomenological model trained on a large empirical data set. A strong emphasis has traditionally been placed on selecting the "best-fit" model, with the implicit understanding that suboptimal or otherwise ill-fitting models might bias inferences. However, the pervasiveness and degree of such bias has not been systematically examined. We investigated how model choice impacts site-wise relative rates in a large set of empirical protein alignments. We compared models designed for use on any general protein, models designed for specific domains of life, and the simple equal-rates Jukes Cantor-style model (JC). As expected, information theoretic measures showed overwhelming evidence that some models fit the data decidedly better than others. By contrast, estimates of site-specific evolutionary rates were impressively insensitive to the substitution model used, revealing an unexpected degree of robustness to potential model misspecification. A deeper examination of the fewer than 5% of sites for which model inferences differed in a meaningful way showed that the JC model could uniquely identify rapidly evolving sites that models with empirically derived exchangeabilities failed to detect. We conclude that relative protein rates appear robust to the applied substitution model, and any sensible model of protein evolution, regardless of its fit to the data, should produce broadly consistent evolutionary rates.

Pretreatment Hepatitis C Virus NS5A/NS5B Resistance-Associated Substitutions in Genotype 1 Uruguayan Infected Patients

Hepatitis C Virus (HCV) infection treatment has dramatically changed with the advent of direct-acting antiviral agents (DAAs). However, the efficacy of DAAs can be attenuated by the presence of resistance-associated substitutions (RASs) before and after treatment. Indeed, RASs detected in DAA treatment-naïve HCV-infected patients could be useful for clinical management and outcome prediction. Although the frequency of naturally occurring HCV NS5A and NS5B RASs has been addressed in many countries, there are only a few reports on their prevalence in the South American region. The aim of this study was to investigate the presence of RASs to NS5A and NS5B inhibitors in a DAA treatment naïve cohort of Uruguayan patients infected with chronic hepatitis C and compare them with reports from other South American countries. Here, we found that naturally occurring substitutions conferring resistance to NS5A and NS5B inhibitors were present in 8% and 19.2%, respectively, of treatment-naïve HCV genotype 1 infected patients. Importantly, the baseline substitutions in NS5A and NS5B herein identified differ from the studies previously reported in Brazil. Furthermore, Uruguayan strains subtype 1a clustered within all major world clades, showing that HCV variants currently circulating in this country are characterized by a remarkable genetic diversity.

Biodiversity estimates and ecological interpretations of meiofaunal communities are biased by the taxonomic approach

Accurate assessments of biodiversity are crucial to advising ecosystem-monitoring programs and understanding ecosystem function. Nevertheless, a standard operating procedure to assess biodiversity accurately and consistently has not been established. This is especially true for meiofauna, a diverse community (>20 phyla) of small benthic invertebrates that have fundamental ecological roles. Recent studies show that metabarcoding is a cost-effective and time-effective method to estimate meiofauna biodiversity, in contrast to morphological-based taxonomy. Here, we compare biodiversity assessments of a diverse meiofaunal community derived by applying multiple taxonomic methods based on comparative morphology, molecular phylogenetic analysis, DNA barcoding of individual specimens, and metabarcoding of environmental DNA. We show that biodiversity estimates are strongly biased across taxonomic methods and phyla. Such biases affect understanding of community structures and ecological interpretations. This study supports the urgency of improving aspects of environmental high-throughput sequencing and the value of taxonomists in correctly understanding biodiversity estimates.

Ancient goat genomes reveal mosaic domestication in the Fertile Crescent

Current genetic data are equivocal as to whether goat domestication occurred multiple times or was a singular process. We generated genomic data from 83 ancient goats (51 with genome-wide coverage) from Paleolithic to Medieval contexts throughout the Near East. Our findings demonstrate that multiple divergent ancient wild goat sources were domesticated in a dispersed process that resulted in genetically and geographically distinct Neolithic goat populations, echoing contemporaneous human divergence across the region. These early goat populations contributed differently to modern goats in Asia, Africa, and Europe. We also detect early selection for pigmentation, stature, reproduction, milking, and response to dietary change, providing 8000-year-old evidence for human agency in molding genome variation within a partner species.

Characterization and comparative mitogenomic analysis of six newly sequenced mitochondrial genomes from ectomycorrhizal fungi (Russula) and phylogenetic analysis of the Agaricomycetes

In this study, the mitochondrial genomes of six Russula species were sequenced using next generation sequencing. The six mitogenomes were all composed of circular DNA molecules, with lengths ranging from 40,961 bp to 69,423 bp. The length and number of protein coding genes (PCGs), GC content, AT skew, and GC skew varied among the six mitogenomes. The increased number and total size of introns likely contributed to the size expansion of mitogenomes in some Russula species. Gene synteny analysis revealed some gene rearrangements among the six mitochondrial genomes. The nad4L gene had the lowest K2P genetic distance of the 15 core PCGs among the six Russula species, indicating that this gene was highly conserved. The Ka/Ks values for all 15 core PCGs were <1, suggesting that they were all subject to purifying selection. Phylogenetic analyses based on two gene datasets (15 core PCGs, and 15 core PCGs + rnl + rns) recovered identical and well-supported trees. In addition, cox1 was identified as a potential single-gene molecular marker for the phylogenetic analysis of relationships among Agaricomycetes species. This study provides the first report of mitogenomes from the Russulaceae family and facilitates the investigation of population genetics and evolution of other ectomycorrhizal fungi.

Resolving deep-sea pelagic saccopharyngiform eel mysteries: Identification of Neocyema and Monognathidae leptocephali and establishment of a new fish family "Neocyematidae" based on larvae, adults and mitogenomic gene orders

Deep-sea midwater "saccopharyngiform" eels of the families Cyematidae, Monognathidae, Eurypharyngidae and Saccopharyngidae (order Anguilliformes) are extraordinary fishes having major skeletal reductions and modifications compared to the general anguilliform body structure. Little is known about most aspects of the systematics, phylogeny, and ecology of these families, and few of the approximately 30 species described from adult specimens have been matched with their leptotocephalus larvae. Based on mitogenomic sequence data from rare new specimens, we show that the long-speculated-about larval form referred to as "Leptocephalus holti", which was thought to possibly be the larva of the rare orange-colored eels of Neocyema (5 known specimens; speculated to belong to the Cyematidae) are actually the larvae of the one-jaw eels of the family Monognathidae. One of the 5 types of L. holti larvae that were collected in the Pacific is genetically matched with Monognathus jesperseni, but multiple species exist based on larval sequence data and the morphology of adult specimens. A rare leptocephalus from the Sargasso Sea, with unique morphological characteristics including many small orange spots on the gut, was found to be the larva of Neocyema, which is presently only known from the Atlantic Ocean. We demonstrate that Neocyema constitutes a separate family being most closely related to Eurypharyngidae and Saccopharyngidae based on mitogenomic DNA sequences and unique mitochondrial gene orders.

Bat Astrovirus in Mozambique

Astroviruses (AstVs) are responsible for infection of a large diversity of mammalian and avian species, including bats, aquatic birds, livestock and humans. We investigated AstVs circulation in bats in Mozambique and Mayotte, a small island in the Comoros Archipelago located between east Africa and Madagascar. Biological material was collected from 338 bats and tested for the presence of the AstV RNA-dependent RNA-polymerase gene with a pan-AstV semi-nested polymerase chain reaction assay. None of the 79 samples obtained from Mayotte bats (Pteropus seychellensis comorensis and Chaerephon pusillus) tested positive; however, 20.1% of bats sampled in Mozambique shed AstVs at the time of sampling and significant interspecific variation in the proportion of positive bats was detected. Many AstVs sequences obtained from a given bat species clustered in different phylogenetic lineages, while others seem to reflect some level of host-virus association, but also with AstVs previously reported from Malagasy bats. Our findings support active circulation of a large diversity of AstVs in bats in the western Indian Ocean islands, including the southeastern African coast, and highlight the need for more detailed assessment of its risk of zoonotic transmission to human populations.

Pax3/7 duplicated and diverged independently in amphioxus, the basal chordate lineage

The Pax3/7 transcription factor family is integral to developmental gene networks contributing to important innovations in vertebrate evolution, including the neural crest. The basal chordate lineage of amphioxus is ideally placed to understand the dynamics of the gene regulatory network evolution that produced these novelties. We report here the discovery that the cephalochordate lineage possesses two Pax3/7 genes, Pax3/7a and Pax3/7b. The tandem duplication is ancestral to all extant amphioxus, occurring in both Asymmetron and Branchiostoma, but originated after the split from the lineage leading to vertebrates. The two paralogues are differentially expressed during embryonic development, particularly in neural and somitic tissues, suggesting distinct regulation. Our results have implications for the study of amphioxus regeneration, neural plate and crest evolution, and differential tandem paralogue evolution.

Lack of host specificity of copepod crustaceans associated with mushroom corals in the Red Sea

The radiation of symbiotic copepods (Crustacea: Copepoda) living in association with stony corals (Cnidaria: Scleractinia) is considered host-specific and linked to the phylogenetic diversification of their hosts. However, symbiotic copepods are poorly investigated, occurrence records are mostly anecdotal, and no explicit analysis exists regarding their relationship with the hosts. Here, we analysed the occurrence of symbiotic copepods on different co-occurring and phylogenetically closely related scleractinian corals. We used an innovative approach of DNA extraction from single microscopic specimens that preserves the shape of the organisms for integrative morphological studies. The rationale of the study involved: (i) sampling of mushroom corals (Fungiidae) belonging to 13 species and eight genera on different reefs along the Saudi coastline in the Red Sea, (ii) extraction of all the associated copepods, (iii) morphological screening and identification of copepod species, (iv) use of DNA taxonomy on mitochondrial and nuclear markers to determine species boundaries for morphologically unknown copepod species, (v) reconstruction of phylogenies to understand their evolutionary relationships, and (vi) analysis of the ecological drivers of the occurrence, diversity and host specificity of the copepods. The seven species of coral-associated copepods, all new to science, did not show any statistically significant evidence of host-specificity or other pattern of ecological association. We thus suggest that, contrary to most assumptions and previous anecdotal evidence on this coral-copepod host-symbiont system, the association between copepods and their host corals is rather labile, not strict, and not phylogenetically constrained, changing our perception on evolutionary patterns and processes in symbiotic copepods.

Whole-genome resequencing reveals world-wide ancestry and adaptive introgression events of domesticated cattle in East Asia

Cattle domestication and the complex histories of East Asian cattle breeds warrant further investigation. Through analysing the genomes of 49 modern breeds and eight East Asian ancient samples, worldwide cattle are consistently classified into five continental groups based on Y-chromosome haplotypes and autosomal variants. We find that East Asian cattle populations are mainly composed of three distinct ancestries, including an earlier East Asian taurine ancestry that reached China at least ~3.9 kya, a later introduced Eurasian taurine ancestry, and a novel Chinese indicine ancestry that diverged from Indian indicine approximately 36.6–49.6 kya. We also report historic introgression events that helped domestic cattle from southern China and the Tibetan Plateau achieve rapid adaptation by acquiring ~2.93% and ~1.22% of their genomes from banteng and yak, respectively. Our findings provide new insights into the evolutionary history of cattle and the importance of introgression in adaptation of cattle to new environmental challenges in East Asia.

There are various indigenous cattle breeds in East Asia which have a complex history. Here, the authors analyse the genomes of 49 modern breeds and eight ancient samples and identify three distinct ancestries and multiple adaptive introgressions from other bovine species.

A new minute ectosymbiotic harpacticoid copepod living on the sea cucumberEupentacta fraudatrixin the East/Japan Sea

The ectosymbiotic copepods,Vostoklaophonte eupentagen. & sp. nov. associated with the sea cucumberEupentacta fraudatrix, was found in the subtidal zone of Peter the Great Bay, East/Japan Sea. The new genus,Vostoklaophonte, is similar toMicrocheloniain the flattened body form, reduced mandible, maxillule and maxilla, but with well-developed prehensile maxilliped, and in the reduced segmentation and setation of legs 1–5. Most appendages of the new genus are more primitive than those ofMicrochelonia. The inclusion of the symbiotic generaMicrocheloniaandVostoklaophontegen. nov. in Laophontidae, as well as their close phylogenetic relationships, are supported by morphological observations and molecular data. This is the third record of laophontid harpacticoid copepods living in symbiosis with sea cucumbers recorded from the Korean and Californian coasts.