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

Characterization and phylogenetic analysis of the complete mitochondrial genome of the medicinal fungus Laetiporus sulphureus

The medicinal fungus Laetiporus sulphureus is widely distributed worldwide. To screen for molecular markers potentially useful for phylogenetic analyses of this species and related species, the mitochondrial genome of L. sulphureus was sequenced and assembled. The complete circular mitochondrial genome was 101,111 bp long, and contained 38 protein-coding genes (PCGs), 2 rRNA genes, and 25 tRNA genes. Our BLAST search aligned about 6.1 kb between the mitochondrial and nuclear genomes of L. sulphureus, indicative of possible gene transfer events. Both the GC and AT skews in the L. sulphureus mitogenome were negative, in contrast to the other seven Polyporales species tested. Of the 15 PCGs conserved across the seven species of Polyporales, the lengths of 11 were unique in the L. sulphureus mitogenome. The Ka/Ks of these 15 PCGs were all less than 1, indicating that PCGs were subject to purifying selection. Our phylogenetic analysis showed that three single genes (cox1, cob, and rnl) were potentially useful as molecular markers. This study is the first publication of a mitochondrial genome in the family Laetiporaceae, and will facilitate the study of population genetics and evolution in L. sulphureus and other species in this family.

Virological Sampling of Inaccessible Wildlife with Drones

There is growing interest in characterizing the viromes of diverse mammalian species, particularly in the context of disease emergence. However, little is known about virome diversity in aquatic mammals, in part due to difficulties in sampling. We characterized the virome of the exhaled breath (or blow) of the Eastern Australian humpback whale (Megaptera novaeangliae). To achieve an unbiased survey of virome diversity, a meta-transcriptomic analysis was performed on 19 pooled whale blow samples collected via a purpose-built Unmanned Aerial Vehicle (UAV, or drone) approximately 3 km off the coast of Sydney, Australia during the 2017 winter annual northward migration from Antarctica to northern Australia. To our knowledge, this is the first time that UAVs have been used to sample viruses. Despite the relatively small number of animals surveyed in this initial study, we identified six novel virus species from five viral families. This work demonstrates the potential of UAVs in studies of virus disease, diversity, and evolution.

Genomic Characterization and Seroprevalence Studies on Alphaviruses in Uruguay

Alphaviruses (Togaviridae) are arboviruses frequently associated with emerging infectious diseases. In this study, we aimed to investigate the presence of alphaviruses in Uruguay by detecting the viral genome in mosquitoes and neutralizing antibodies in equines. A total of 3,575 mosquitoes were analyzed for alphavirus genome detection. Serologic studies were performed on 425 horse sera by plaque reduction neutralization test (PRNT80) against Venezuelan equine encephalitis virus (VEEV) subtype IAB, Pixuna virus (PIXV), Rio Negro virus (RNV), western equine encephalitis virus (WEEV), and Madariaga virus (MADV). Mosquitoes belonging to six genera were captured and 82.9% were identified as Culex pipiens. Two Cx. pipiens pools collected in Fray Bentos and Las Toscas localities were alphavirus positive, and phylogenetic analyses showed that the sequences grouped into two different clusters: the lineage I of eastern equine encephalitis virus and RNV (VEEV complex), respectively. Plaque reduction neutralization test assays showed antibodies against strains of the VEEV complex, MADV, and WEEV. Rio Negro virus was the most geographically widespread virus, showing higher seroprevalences (up to 20%). Seroprevalences against VEEV IAB ranged between 4.6% and 13%; antibodies against PIXV, WEEV, and MADV were less frequent (3-4%). In conclusion, RNV exhibited the highest seroprevalence in horses, a wide geographical distribution, and viral genome was detected in Cx. pipiens mosquitoes. Madariaga virus had a low seroprevalence in equines, but an epizootic lineage typical of North America was detected in Cx. pipiens mosquitoes. Taken together, our results show that alphaviruses are present in Uruguay with variable occurrence and geographical distribution being a potential threat for human and equine health.

Comparative mitogenomics supports synonymy of the genera Ligula and Digramma (Cestoda: Diphyllobothriidae)

Background

After observing differences in the number of reproductive complexes per proglottid within the genus Ligula, the genus Digramma was erected. However, the validity of Digramma has been previously questioned due to a low variability in the cox1, nad1 and ITS rDNA sequences between the two genera. We undertook a study to greatly increase the amount of sequence data available for resolution of this question by sequencing and characterizing the complete mitogenomes of Digramma interrupta and Ligula intestinalis.

Results

The circular mtDNA molecules of Digramma interrupta and Ligula intestinalis are 13,685 bp and 13,621 bp in size, respectively, both comprising 12 PCGs, 22 tRNA genes, two rRNA genes, and two mNCRs. Both mitogenomes exhibit the same gene order and share 92.7% nucleotide identity, compared with 85.8–86.5% to the most closely related genus Dibothriocephalus. Each gene from D. interrupta and L. intestinalis is almost of the same size, and the sequence identity ranges from 87.5% (trnD) to 100% (trnH, trnQ and trnV). NCR2 sequences of D. interrupta and L. intestinalis are 249 bp and 183 bp in length, respectively, which contributes to the main difference in length between their complete mitogenomes. A sliding window analysis of the 12 PCGs and two rRNAs indicated nucleotide diversity to be higher in nad5, nad6, nad2, nad4 and cox3, whereas the most conserved genes were rrnL and rrnS. Lower sequence identity was also found in nad2, nad4, nad5, nad6 and cox3 genes between the two diphyllobothriids. Within the Diphyllobothriidae, phylogenetic analysis indicated Ligula and Digramma to be most closely related to one another, forming a sister group with Dibothriocephalus.

Conclusions

Owing to higher nucleotide diversity, the genes nad2, nad4, nad5, nad6 and cox3 should be considered optimal candidates to use as molecular markers for population genetics and species identification between the two closely related species. The phylogenetic results in combination with the comparative analysis of the two mitogenomes, consistently support the congeneric status of L. intestinalis and D. interrupta.

Electronic supplementary material

The online version of this article (10.1186/s13071-018-2910-9) contains supplementary material, which is available to authorized users.

An ancient FMRFamide-related peptide-receptor pair induces defence behaviour in a brachiopod larva

Animal behaviour often comprises spatially separated sub-reactions and even ciliated larvae are able to coordinate sub-reactions of complex behaviours (metamorphosis, feeding). How these sub-reactions are coordinated is currently not well understood. Neuropeptides are potential candidates for triggering larval behaviour. However, although their immunoreactivity has been widely analysed, their function in trochozoan larvae has only been studied for a few cases. Here, we investigate the role of neuropeptides in the defence behaviour of brachiopod larvae. When mechanically disturbed, the planktonic larvae of Terebratalia transversa protrude their stiff chaetae and sink down slowly. We identified endogenous FLRFamide-type neuropeptides (AFLRFamide and DFLRFamide) in T. transversa larvae and show that the protrusion of the chaetae as well as the sinking reaction can both be induced by each of these peptides. This also correlates with the presence of FLRFamidergic neurons in the apical lobe and adjacent to the trunk musculature. We deorphanized the AFLRFamide/DFLRFamide receptor and detected its expression in the same tissues. Furthermore, the ability of native and modified FLRFamide-type peptides to activate this receptor was found to correspond with their ability to trigger behavioural responses. Our results show how FLRFamide-type neuropeptides can induce two coherent sub-reactions in a larva with a simple nervous system.

Genetic diversity of bat orthohepadnaviruses in China and a proposed new nomenclature

The orthohepadnaviruses, which include the major human pathogen hepatitis B virus, exist in a wide range of hosts. Since 2013, a large group of orthohepadnaviruses has been identified in bats worldwide and classified as 4 species within the genus Orthohepadnavirus. To further investigate orthohepadnaviruses in the Chinese bat population, 554 archived bat samples from 20 colonies covering 3 southern provinces were screened with results showing that 9 (1.6%) were positive. A systematic phylogenetic analysis has indicated the need for a new nomenclature for bat hepatitis B virus-like viruses: BtHBV, with the addition of 3 new species, one being divided into 6 genotypes. Viruses identified here shared 9.0–19.2% full genome divergence and classified into 3 different genotypes. This study illustrates the genetic diversity of orthohepadnaviruses in the Chinese bat population, and emphasizes need for further investigation of their public health significance.

•

Three new orthohepadnaviral lineages were identified in Chinese bats.

•

A new nomenclature was proposed for bat hepatitis B virus-like viruses.

•

This study indicates genetic diversity of orthohepadnaviruses in Chinese bats.

Virus survey in populations of two subspecies of bent-winged bats (Miniopterus orianae bassanii and oceanensis) in south-eastern Australia reveals a high prevalence of diverse herpesviruses

While bats are often viewed as carriers of infectious disease agents, little research has been conducted on the effects these potential 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 been declining for 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, 351 southern bent-winged bats from four locations were captured, and oral swabs were collected and tested for the presence of potentially pathogenic viruses. Results were compared with those obtained from 116 eastern bent-winged bats (Miniopterus orianae oceanensis) from three different locations. The eastern bent-winged bat is a related but more common and widespread subspecies whose geographical range overlaps partly with southern bent-winged bats. Herpesviruses were detected in bent-winged bats from all seven locations. At least six novel herpesviruses (five betaherpesviruses and one gammaherpesvirus) were identified. The prevalence of herpesvirus infection was higher in eastern bent-winged bats (44%, 51/116), compared to southern bent-winged bats (27%, 95/351), although this varied across the locations and sampling periods. Adenoviruses and a range of different RNA viruses (lyssaviruses, filoviruses, coronaviruses and henipaviruses) were also tested for but not detected. The detected herpesviruses did not appear to be associated with obvious ill health, and may thus not be playing a role in the population decline of the southern bent-winged bat. The detection of multiple novel herpesviruses at a high prevalence of infection is consistent with our understanding of bats as hosts to a rich diversity of viruses.

Population Structure and Local Adaptation of MAC Lung Disease Agent Mycobacterium avium subsp. hominissuis

Mycobacterium avium subsp. hominissuis (MAH) is one of the most common nontuberculous mycobacterial species responsible for chronic lung disease in humans. Despite increasing worldwide incidence, little is known about the genetic mechanisms behind the population evolution of MAH. To elucidate the local adaptation mechanisms of MAH, we assessed genetic population structure, the mutual homologous recombination, and gene content for 36 global MAH isolates, including 12 Japanese isolates sequenced in the present study. We identified five major MAH lineages and found that extensive mutual homologous recombination occurs among them. Two lineages (MahEastAsia1 and MahEastAsia2) were predominant in the Japanese isolates. We identified alleles unique to these two East Asian lineages in the loci responsible for trehalose biosynthesis (treS and mak) and in one mammalian cell entry operon, which presumably originated from as yet undiscovered mycobacterial lineages. Several genes and alleles unique to East Asian strains were located in the fragments introduced via recombination between East Asian lineages, suggesting implication of recombination in local adaptation. These patterns of MAH genomes are consistent with the signature of distribution conjugative transfer, a mode of sexual reproduction reported for other mycobacterial species.

First detection of foot-and-mouth disease virus O/ME-SA/Ind2001 in China

Foot-and-mouth disease (FMD) is endemic in China and is predominantly due to foot-and-mouth disease virus (FMDV) serotype O Mya-98 lineage. In recent years, FMDV O/ME-SA/Ind2001 lineage has spread from the Indian subcontinent to South-East Asia, Middle East and Africa, which may pose potential threats for future trans-regional livestock movements. In this study, we identified the appearance of FMDV O/ME-SA/Ind2001 in China; the first time that this virus lineage has been found there. Sequencing and phylogenetic analysis of VP1 sequences revealed that this newly determined strain belongs to O/ME-SA/Ind2001 sublineage d and is closely related to strains that have caused recent outbreaks of FMD in Nepal, Myanmar, Russia and South Korea. The results suggest extensive movements of the current O/ME-SA/Ind2001 sublineage d viruses and provide essential information for an effective national FMDV control programme in China.

Daphnia diversity on the Tibetan Plateau measured by DNA taxonomy

Daphnia on the Tibetan Plateau has been little studied, and information on species diversity and biogeography is lacking. Here, we conducted a 4-year survey using the barcoding fragment of the mitochondrial COI gene to determine the distribution and diversity of Daphnia species found across the Plateau. Our results show that species richness is higher than previously thought, with total described and provisional species number doubling from 5 to 10. Six of the taxonomic units recovered by DNA taxonomy agreed well with morphology, but DNA barcoding distinguished three clades each for the D. longispina (D. galeata, D. dentifera, and D. longispina) and D. pulex (D. pulex, D. cf. tenebrosa, and D. pulicaria) complexes. The sequence divergence between congeneric species varied within a large range, from 9.25% to 30.71%. The endemic D. tibetana was the most common and widespread species, occurring in 12 hyposaline to mesosaline lakes. The lineage of D. longispina is the first confirmed occurrence in west Tibet.

Diversity of marine planktonic ostracods in South China Sea: a DNA taxonomy approach

Ostracods (Crustacea, Ostracoda) are small bivalved crustaceans, contributing over 200 described species to the marine zooplankton community. They are widely distributed and are relatively abundant components of the mesozooplankton, playing an important role in the transport of organic matter to deep layers. However, identification of ostracods based on micro-morphological characters is extremely difficult and time-consuming. Previous fragmentary taxonomic studies of ostracods in the South China Sea (SCA), were based solely on morphology. Here, by analysing the mitochondrial COI gene, we explore the taxa across the SCA using molecular tools for the first time. Our results show that sequence divergence among species varies within a large range, from 12.93% to 35.82%. Sixteen of the taxonomic units recovered by DNA taxonomy agree well with morphology, but Paraconchoecia oblonga, Conchoecia magna and Halocypris brevirostris split into two clades each, each of which contains cryptic species.

Co-evolutionary interactions between host resistance and pathogen avirulence genes in rice-Magnaporthe oryzae pathosystem

Rice and Magnaporthe oryzae constitutes an ideal pathosystem for studying host-pathogen interaction in cereals crops. There are two alternative hypotheses, viz. Arms race and Trench warfare, which explain the co-evolutionary dynamics of hosts and pathogens which are under continuous confrontation. Arms race proposes that both R- and Avr- genes of host and pathogen, respectively, undergo positive selection. Alternatively, trench warfare suggests that either R- or Avr- gene in the pathosystem is under balanced selection intending to stabilize the genetic advantage gained over the opposition. Here, we made an attempt to test the above-stated hypotheses in rice-M. oryzae pathosystem at loci of three R-Avr gene pairs, Piz-t-AvrPiz-t, Pi54-AvrPi54 and Pita-AvrPita using allele mining approach. Allele mining is an efficient way to capture allelic variants existing in the population and to study the selective forces imposed on the variants during evolution. Results of nucleotide diversity, neutrality statistics and phylogenetic analyses reveal that Piz-t, Pi54 and AvrPita are diversified and under positive selection at their corresponding loci, while their counterparts, AvrPiz-t, AvrPi54 and Pita are conserved and under balancing selection, in nature. These results imply that rice-M. oryzae populations are engaged in a trench warfare at least at the three R/Avr loci studied. It is a maiden attempt to study the co-evolution of three R-Avr gene pairs in this pathosystem. Knowledge gained from this study will help in understanding the evolutionary dynamics of host-pathogen interaction in a better way and will also aid in developing new durable blast resistant rice varieties in future.

Genomic and epigenomic immunity in common bean: the unusual features of NB-LRR gene family

In plants, a key class of genes comprising most of disease resistance (R) genes encodes Nucleotide-binding leucine-rich repeat (NL) proteins. Access to common bean (Phaseolus vulgaris) genome sequence provides unparalleled insight into the organization and evolution of this large gene family (∼400 NL) in this important crop. As observed in other plant species, most common bean NL are organized in cluster of genes. However, a particularity of common bean is that these clusters are often located in subtelomeric regions close to terminal knobs containing the satellite DNA khipu. Phylogenetically related NL are spread between different chromosome ends, suggesting frequent exchanges between non-homologous chromosomes. NL peculiar location, in proximity to heterochromatic regions, led us to study their DNA methylation status using a whole-genome cytosine methylation map. In common bean, NL genes displayed an unusual body methylation pattern since half of them are methylated in the three contexts, reminiscent of the DNA methylation pattern of repeated sequences. Moreover, 90 NL were also abundantly targeted by 24 nt siRNA, with 90% corresponding to methylated NL genes. This suggests the existence of a transcriptional gene silencing mechanism of NL through the RdDM (RNA-directed DNA methylation) pathway in common bean that has not been described in other plant species.

Dropping Hints: Estimating the diets of livestock in rangelands using DNA metabarcoding of faeces

The introduction of domesticated animals into new environments can lead to considerable ecological disruption, and it can be difficult to predict their impact on the new ecosystem. In this study, we use faecal metabarcoding to characterize the diets of three ruminant taxa in the rangelands of south-western New South Wales, Australia. Our study organisms included goats (Capraaegagrushircus) and two breeds of sheep (Ovisaries): Merinos, which have been present in Australia for over two hundred years, and Dorpers, which were introduced in the 1990s. We used High-Throughput Sequencing methods to sequence therbcLand ITS2 genes of plants in the faecal samples, and identified the samples using the GenBank and BOLD online databases, as well as a reference collection of sequences from plants collected in the study area. We found that the diets of all three taxa were dominated by the family Malvaceae, and that the Dorper diet was more diverse than the Merino diet at both the family and the species level. We conclude that Dorpers, like Merinos, are potentially a threat to some vulnerable species in the rangelands of New South Wales.

Genome-wide identification and characterization of the RIO atypical kinase family in plants

Members of the right open reading frame (RIO) atypical kinase family are present in all three domains of life. In eukaryotes, three subfamilies have been identified: RIO1, RIO2, and RIO3. Studies have shown that the yeast and human RIO1 and RIO2 kinases are essential for the biogenesis of small ribosomal subunits. Thus far, RIO3 has been found only in multicellular eukaryotes. In this study, we systematically identified members of the RIO gene family in 37 species representing the major evolutionary lineages in Viridiplantae. A total of 84 RIO genes were identified; among them, 41 were classified as RIO1 and 43 as RIO2. However, no RIO3 gene was found in any of the species examined. Phylogenetic trees constructed for plant RIO1 and RIO2 proteins were generally congruent with the species phylogeny. Subcellular localization analyses showed that the plant RIO proteins were localized mainly in the nucleus and/or cytoplasm. Expression profile analysis of rice, maize, and Arabidopsis RIO genes in different tissues revealed similar expression patterns between RIO1 and RIO2 genes, and their expression levels were high in certain tissues. In addition, the expressions of plant RIO genes were regulated by two drugs: mycophenolic acid and actinomycin D. Function prediction using genome-wide coexpression analysis revealed that most plant RIO genes may be involved in ribosome biogenesis. Our results will be useful for the evolutionary analysis of the ancient RIO kinase family and provide a basis for further functional characterization of RIO genes in plants.

Ancient genomes revisit the ancestry of domestic and Przewalski’s horses

The Eneolithic Botai culture of the Central Asian steppes provides the earliest archaeological evidence for horse husbandry, ~5500 years ago, but the exact nature of early horse domestication remains controversial. We generated 42 ancient-horse genomes, including 20 from Botai. Compared to 46 published ancient- and modern-horse genomes, our data indicate that Przewalski’s horses are the feral descendants of horses herded at Botai and not truly wild horses. All domestic horses dated from ~4000 years ago to present only show ~2.7% of Botai-related ancestry. This indicates that a massive genomic turnover underpins the expansion of the horse stock that gave rise to modern domesticates, which coincides with large-scale human population expansions during the Early Bronze Age.

The Impact of the Tree Prior on Molecular Dating of Data Sets Containing a Mixture of Inter- and Intraspecies Sampling

In Bayesian phylogenetic analyses of genetic data, prior probability distributions need to be specified for the model parameters, including the tree. When Bayesian methods are used for molecular dating, available tree priors include those designed for species-level data, such as the pure-birth and birth-death priors, and coalescent-based priors designed for population-level data. However, molecular dating methods are frequently applied to data sets that include multiple individuals across multiple species. Such data sets violate the assumptions of both the speciation and coalescent-based tree priors, making it unclear which should be chosen and whether this choice can affect the estimation of node times. To investigate this problem, we used a simulation approach to produce data sets with different proportions of within- and between-species sampling under the multispecies coalescent model. These data sets were then analyzed under pure-birth, birth-death, constant-size coalescent, and skyline coalescent tree priors. We also explored the ability of Bayesian model testing to select the best-performing priors. We confirmed the applicability of our results to empirical data sets from cetaceans, phocids, and coregonid whitefish. Estimates of node times were generally robust to the choice of tree prior, but some combinations of tree priors and sampling schemes led to large differences in the age estimates. In particular, the pure-birth tree prior frequently led to inaccurate estimates for data sets containing a mixture of inter- and intraspecific sampling, whereas the birth-death and skyline coalescent priors produced stable results across all scenarios. Model testing provided an adequate means of rejecting inappropriate tree priors. Our results suggest that tree priors do not strongly affect Bayesian molecular dating results in most cases, even when severely misspecified. However, the choice of tree prior can be significant for the accuracy of dating results in the case of data sets with mixed inter- and intraspecies sampling. [Bayesian phylogenetic methods; model testing; molecular dating; node time; tree prior.].

The complete mitochondrial genome of Huananpotamon lichuanense (Decapoda: Brachyura) with phylogenetic implications for freshwater crabs

In the present study, we determined the complete mitochondrial genome of Huananpotamon lichuanense (Decapoda: Brachyura) for the first time. The genome is 15,380bp in length and typically consists of 37 genes. When the gene order was compared to the ancestral crustacean type, two tRNA genes (tRNA^His and tRNA^Gln) were rearranged in H. lichuanense, and the translocation of tRNA^Gln appeared only in Potamoidea crabs, such as Geothelphusa dehaani and Sinopotamon xiushuiense, supporting the monophyly of the Potamoidea superfamily. Thirteen protein-coding genes and 2 rRNA genes were divided into five complexes to perform the phylogenetic analysis, and the results showed that the trees constructed by complex I (ND1-ND6 and ND4L), complex IV (COX1-COX3) and rRNA genes better accord with the morphological classification system, suggesting that molecular markers of higher-level phylogeny can be developed in these three complexes in the future. The estimated divergence time for freshwater crabs is approximately 133.58Ma, and G. dehaani from Japan diverged from the freshwater crabs of mainland China approximately 60.66Ma. A selective pressure analysis based on current data revealed obviously increasing dN/dS ratios (except for ATP6 and ND4L) of freshwater crabs, and the accumulation of nonsynonymous mutations suggests that terrestrial habitats provide a relatively relaxed selective pressure environment for this group.

Cameroonian fruit bats harbor divergent viruses, including rotavirus H, bastroviruses, and picobirnaviruses using an alternative genetic code

Most human emerging infectious diseases originate from wildlife and bats are a major reservoir of viruses, a few of which have been highly pathogenic to humans. In some regions of Cameroon, bats are hunted and eaten as a delicacy. This close proximity between human and bats provides ample opportunity for zoonotic events. To elucidate the viral diversity of Cameroonian fruit bats, we collected and metagenomically screened eighty-seven fecal samples of Eidolon helvum and Epomophorus gambianus fruit bats. The results showed a plethora of known and novel viruses. Phylogenetic analyses of the eleven gene segments of the first complete bat rotavirus H genome, showed clearly separated clusters of human, porcine, and bat rotavirus H strains, not indicating any recent interspecies transmission events. Additionally, we identified and analyzed a bat bastrovirus genome (a novel group of recently described viruses, related to astroviruses and hepatitis E viruses), confirming their recombinant nature, and provide further evidence of additional recombination events among bat bastroviruses. Interestingly, picobirnavirus-like RNA-dependent RNA polymerase gene segments were identified using an alternative mitochondrial genetic code, and further principal component analyses suggested that they may have a similar lifestyle to mitoviruses, a group of virus-like elements known to infect the mitochondria of fungi. Although identified bat coronavirus, parvovirus, and cyclovirus strains belong to established genera, most of the identified partitiviruses and densoviruses constitute putative novel genera in their respective families. Finally, the results of the phage community analyses of these bats indicate a very diverse geographically distinct bat phage population, probably reflecting different diets and gut bacterial ecosystems.

Genome-Wide Identification, Evolutionary and Expression Analyses of the GALACTINOL SYNTHASE Gene Family in Rapeseed and Tobacco

Galactinol synthase (GolS) is a key enzyme in raffinose family oligosaccharide (RFO) biosynthesis. The finding that GolS accumulates in plants exposed to abiotic stresses indicates RFOs function in environmental adaptation. However, the evolutionary relationships and biological functions of GolS family in rapeseed (Brassica napus) and tobacco (Nicotiana tabacum) remain unclear. In this study, we identified 20 BnGolS and 9 NtGolS genes. Subcellular localization predictions showed that most of the proteins are localized to the cytoplasm. Phylogenetic analysis identified a lost event of an ancient GolS copy in the Solanaceae and an ancient duplication event leading to evolution of GolS4/7 in the Brassicaceae. The three-dimensional structures of two GolS proteins were conserved, with an important DxD motif for binding to UDP-galactose (uridine diphosphate-galactose) and inositol. Expression profile analysis indicated that BnGolS and NtGolS genes were expressed in most tissues and highly expressed in one or two specific tissues. Hormone treatments strongly induced the expression of most BnGolS genes and homologous genes in the same subfamilies exhibited divergent-induced expression. Our study provides a comprehensive evolutionary analysis of GolS genes among the Brassicaceae and Solanaceae as well as an insight into the biological function of GolS genes in hormone response in plants.

Characterization of HIV-1 Near Full-Length Proviral Genome Quasispecies from Patients with Undetectable Viral Load Undergoing First-Line HAART Therapy

Increased access to highly active antiretroviral therapy (HAART) by human immunodeficiency virus postive (HIV⁺) individuals has become a reality worldwide. In Brazil, HAART currently reaches over half of HIV-infected subjects. In the context of a remarkable HIV-1 genetic variability, highly related variants, called quasispecies, are generated. HIV quasispecies generated during infection can influence virus persistence and pathogenicity, representing a challenge to treatment. However, the clinical relevance of minority quasispecies is still uncertain. In this study, we have determined the archived proviral sequences, viral subtype and drug resistance mutations from a cohort of HIV⁺ patients with undetectable viral load undergoing HAART as first-line therapy using next-generation sequencing for near full-length virus genome (NFLG) assembly. HIV-1 consensus sequences representing NFLG were obtained for eleven patients, while for another twelve varying genome coverage rates were obtained. Phylogenetic analysis showed the predominance of subtype B (83%; 19/23). Considering the minority variants, 18 patients carried archived virus harboring at least one mutation conferring antiretroviral resistance; for six patients, the mutations correlated with the current ARVs used. These data highlight the importance of monitoring HIV minority drug resistant variants and their clinical impact, to guide future regimen switches and improve HIV treatment success.

Central European Woolly Mammoth Population Dynamics: Insights from Late Pleistocene Mitochondrial Genomes

The population dynamics of the Pleistocene woolly mammoth (Mammuthus primigenius) has been the subject of intensive palaeogenetic research. Although a large number of mitochondrial genomes across Eurasia have been reconstructed, the available data remains geographically sparse and mostly focused on eastern Eurasia. Thus, population dynamics in other regions have not been extensively investigated. Here, we use a multi-method approach utilising proteomic, stable isotope and genetic techniques to identify and generate twenty woolly mammoth mitochondrial genomes, and associated dietary stable isotopic data, from highly fragmentary Late Pleistocene material from central Europe. We begin to address region-specific questions regarding central European woolly mammoth populations, highlighting parallels with a previous replacement event in eastern Eurasia ten thousand years earlier. A high number of shared derived mutations between woolly mammoth mitochondrial clades are identified, questioning previous phylogenetic analysis and thus emphasizing the need for nuclear DNA studies to explicate the increasingly complex genetic history of the woolly mammoth.

Inferring the evolution of the major histocompatibility complex of wild pigs and peccaries using hybridisation DNA capture-based sequencing

The major histocompatibility complex (MHC) is a key genomic model region for understanding the evolution of gene families and the co-evolution between host and pathogen. To date, MHC studies have mostly focused on species from major vertebrate lineages. The evolution of MHC classical (Ia) and non-classical (Ib) genes in pigs has attracted interest because of their antigen presentation roles as part of the adaptive immune system. The pig family Suidae comprises over 18 extant species (mostly wild), but only the domestic pig has been extensively sequenced and annotated. To address this, we used a DNA-capture approach, with probes designed from the domestic pig genome, to generate MHC data for 11 wild species of pigs and their closest living family, Tayassuidae. The approach showed good efficiency for wild pigs (~80% reads mapped, ~87× coverage), compared to tayassuids (~12% reads mapped, ~4× coverage). We retrieved 145 MHC loci across both families. Phylogenetic analyses show that the class Ia and Ib genes underwent multiple duplications and diversifications before suids and tayassuids diverged from their common ancestor. The histocompatibility genes mostly form orthologous groups and there is genetic differentiation for most of these genes between Eurasian and sub-Saharan African wild pigs. Tests of selection showed that the peptide-binding region of class Ib genes was under positive selection. These findings contribute to better understanding of the evolutionary history of the MHC, specifically, the class I genes, and provide useful data for investigating the immune response of wild populations against pathogens.

A novel type I cystatin of parasite origin with atypical legumain-binding domain

Parasite inhibitors of cysteine peptidases are known to influence a vast range of processes linked to a degradation of either the parasites' own proteins or proteins native to their hosts. We characterise a novel type I cystatin (stefin) found in a sanguinivorous fish parasite Eudiplozoon nipponicum (Platyhelminthes: Monogenea). We have identified a transcript of its coding gene in the transcriptome of adult worms. Its amino acid sequence is similar to other stefins except for containing a legumain-binding domain, which is in this type of cystatins rather unusual. As expected, the recombinant form of E. nipponicum stefin (rEnStef) produced in Escherichia coli inhibits clan CA peptidases - cathepsins L and B of the worm - via the standard papain-binding domain. It also blocks haemoglobinolysis by cysteine peptidases in the worm's excretory-secretory products and soluble extracts. Furthermore, we had confirmed its ability to inhibit clan CD asparaginyl endopeptidase (legumain). The presence of a native EnStef in the excretory-secretory products of adult worms, detected by mass spectrometry, suggests that this protein has an important biological function at the host-parasite interface. We discuss the inhibitor's possible role in the regulation of blood digestion, modulation of antigen presentation, and in the regeneration of host tissues.

Producing Designer Oils in Industrial Microalgae by Rational Modulation of Co-evolving Type-2 Diacylglycerol Acyltransferases

Microalgal oils, depending on their degree of unsaturation, can be utilized as either nutritional supplements or fuels; thus, a feedstock with genetically designed and tunable degree of unsaturation is desirable to maximize process efficiency and product versatility. Systematic profiling of ex vivo (in yeast), in vitro, and in vivo activities of type-2 diacylglycerol acyltransferases in Nannochloropsis oceanica (NoDGAT2s or NoDGTTs), via reverse genetics, revealed that NoDGAT2A prefers saturated fatty acids (SFAs), NoDGAT2D prefers monounsaturated fatty acids (MUFAs), and NoDGAT2C exhibits the strongest activity toward polyunsaturated fatty acids (PUFAs). As NoDGAT2A, 2C, and 2D originated from the green alga, red alga, and eukaryotic host ancestral participants of secondary endosymbiosis, respectively, a mechanistic model of oleaginousness was unveiled, in which the indigenous and adopted NoDGAT2s formulated functional complementarity and specific transcript abundance ratio that underlie a rigid SFA:MUFA:PUFA hierarchy in triacylglycerol (TAG). By rationally modulating the ratio of NoDGAT2A:2C:2D transcripts, a bank of N. oceanica strains optimized for nutritional supplement or fuel production with a wide range of degree of unsaturation were created, in which proportion of SFAs, MUFAs, and PUFAs in TAG varied by 1.3-, 3.7-, and 11.2-fold, respectively. This established a novel strategy to simultaneously improve productivity and quality of oils from industrial microalgae.

Mitogenome evolution in the last surviving woolly mammoth population reveals neutral and functional consequences of small population size

The onset of the Holocene was associated with a global temperature increase, which led to a rise in sea levels and isolation of the last surviving population of woolly mammoths on Wrangel Island. Understanding what happened with the population's genetic diversity at the time of the isolation and during the ensuing 6000 years can help clarify the effects of bottlenecks and subsequent limited population sizes in species approaching extinction. Previous genetic studies have highlighted questions about how the Holocene Wrangel population was established and how the isolation event affected genetic diversity. Here, we generated high‐quality mitogenomes from 21 radiocarbon‐dated woolly mammoths to compare the ancestral large and genetically diverse Late Pleistocene Siberian population and the small Holocene Wrangel population. Our results indicate that mitogenome diversity was reduced to one single haplotype at the time of the isolation, and thus that the Holocene Wrangel Island population was established by a single maternal lineage. Moreover, we show that the ensuing small effective population size coincided with fixation of a nonsynonymous mutation, and a comparative analysis of mutation rates suggests that the evolutionary rate was accelerated in the Holocene population. These results suggest that isolation on Wrangel Island led to an increase in the frequency of deleterious genetic variation, and thus are consistent with the hypothesis that strong genetic drift in small populations leads to purifying selection being less effective in removing deleterious mutations.

Genomic, Network, and Phylogenetic Analysis of the Oomycete Effector Arsenal

The oomycetes are a class of microscopic, filamentous eukaryotes and include ecologically significant animal and plant pathogens. Oomycetes secrete large arsenals of effector proteins that degrade host cell components, manipulate host immune responses, and induce necrosis, enabling parasitic colonization. In this study, we catalogued the number and evolution of effectors in 37 oomycete species whose genomes have been completely sequenced. Large expansions of effector protein families in Phytophthora species, including glycoside hydrolases, pectinases, and necrosis-inducing proteins, were observed. Species-specific expansions were detected, including chitinases in Aphanomyces astaci and Pythium oligandrum. Novel effectors which may be involved in suppressing animal immune responses were identified in Ap. astaci and Py. oligandrum. Type 2 necrosis-inducing proteins with an unusual phylogenetic history were also located. This work represents an up-to-date in silico catalogue of the effector arsenal of the oomycetes based on the 37 genomes currently available.

The oomycetes are a class of microscopic, filamentous eukaryotes within the stramenopiles-alveolate-Rhizaria (SAR) supergroup and include ecologically significant animal and plant pathogens. Oomycetes secrete large arsenals of effector proteins that degrade host cell components, manipulate host immune responses, and induce necrosis, enabling parasitic colonization. This study investigated the expansion and evolution of effectors in 37 oomycete species in 4 oomycete orders, including Albuginales, Peronosporales, Pythiales, and Saprolegniales species. Our results highlight the large expansions of effector protein families, including glycoside hydrolases, pectinases, and necrosis-inducing proteins, in Phytophthora species. Species-specific expansions, including expansions of chitinases in Aphanomyces astaci and Pythium oligandrum, were detected. Novel effectors which may be involved in suppressing animal immune responses in Ap. astaci and Py. insidiosum were also identified. Type 2 necrosis-inducing proteins with an unusual phylogenetic history were also located in a number of oomycete species. We also investigated the "RxLR" effector complement of all 37 species and, as expected, observed large expansions in Phytophthora species numbers. Our results provide in-depth sequence information on all putative RxLR effectors from all 37 species. This work represents an up-to-date in silico catalogue of the effector arsenal of the oomycetes based on the 37 genomes currently available.

IMPORTANCE The oomycetes are a class of microscopic, filamentous eukaryotes and include ecologically significant animal and plant pathogens. Oomycetes secrete large arsenals of effector proteins that degrade host cell components, manipulate host immune responses, and induce necrosis, enabling parasitic colonization. In this study, we catalogued the number and evolution of effectors in 37 oomycete species whose genomes have been completely sequenced. Large expansions of effector protein families in Phytophthora species, including glycoside hydrolases, pectinases, and necrosis-inducing proteins, were observed. Species-specific expansions were detected, including chitinases in Aphanomyces astaci and Pythium oligandrum. Novel effectors which may be involved in suppressing animal immune responses were identified in Ap. astaci and Py. oligandrum. Type 2 necrosis-inducing proteins with an unusual phylogenetic history were also located. This work represents an up-to-date in silico catalogue of the effector arsenal of the oomycetes based on the 37 genomes currently available.

The Stone Age Plague and Its Persistence in Eurasia

Yersinia pestis, the etiologic agent of plague, is a bacterium associated with wild rodents and their fleas. Historically it was responsible for three pandemics: the Plague of Justinian in the 6^th century AD, which persisted until the 8^th century [1]; the renowned Black Death of the 14^th century [2, 3], with recurrent outbreaks until the 18^th century [4]; and the most recent 19^th century pandemic, in which Y. pestis spread worldwide [5] and became endemic in several regions [6]. The discovery of molecular signatures of Y. pestis in prehistoric Eurasian individuals and two genomes from Southern Siberia suggest that Y. pestis caused some form of disease in humans prior to the first historically documented pandemic [7]. Here, we present six new European Y. pestis genomes spanning the Late Neolithic to the Bronze Age (LNBA; 4,800 to 3,700 calibrated years before present). This time period is characterized by major transformative cultural and social changes that led to cross-European networks of contact and exchange [8, 9]. We show that all known LNBA strains form a single putatively extinct clade in the Y. pestis phylogeny. Interpreting our data within the context of recent ancient human genomic evidence that suggests an increase in human mobility during the LNBA, we propose a possible scenario for the early spread of Y. pestis: the pathogen may have entered Europe from Central Eurasia following an expansion of people from the steppe, persisted within Europe until the mid-Bronze Age, and moved back toward Central Eurasia in parallel with human populations.

The complete mitochondrial genome of parasitic nematode Camallanus cotti: extreme discontinuity in the rate of mitogenomic architecture evolution within the Chromadorea class

BACKGROUND

Complete mitochondrial genomes are much better suited for the taxonomic identification and phylogenetic studies of nematodes than morphology or traditionally-used molecular markers, but they remain unavailable for the entire Camallanidae family (Chromadorea). As the only published mitogenome in the Camallanina suborder (Dracunculoidea superfamily) exhibited a unique gene order, the other objective of this research was to study the evolution of mitochondrial architecture in the Spirurida order. Thus, we sequenced the complete mitogenome of the Camallanus cotti fish parasite and conducted structural and phylogenomic comparative analyses with all available Spirurida mitogenomes.

RESULTS

The mitogenome is exceptionally large (17,901 bp) among the Chromadorea and, with 46 (pseudo-) genes, exhibits a unique architecture among nematodes. Six protein-coding genes (PCGs) and six tRNAs are duplicated. An additional (seventh) tRNA (Trp) was probably duplicated by the remolding of tRNA-Ser2 (missing). Two pairs of these duplicated PCGs might be functional; three were incomplete and one contained stop codons. Apart from Ala and Asp, all other duplicated tRNAs are conserved and probably functional. Only 19 unique tRNAs were found. Phylogenomic analysis included Gnathostomatidae (Spirurina) in the Camallanina suborder.

CONCLUSIONS

Within the Nematoda, comparable PCG duplications were observed only in the enoplean Mermithidae family, but those result from mitochondrial recombination, whereas characteristics of the studied mitogenome suggest that likely rearrangement mechanisms are either a series of duplications, transpositions and random loss events, or duplication, fragmentation and subsequent reassembly of the mitogenome. We put forward a hypothesis that the evolution of mitogenomic architecture is extremely discontinuous, and that once a long period of stasis in gene order and content has been punctuated by a rearrangement event, such a destabilised mitogenome is much more likely to undergo subsequent rearrangement events, resulting in an exponentially accelerated evolutionary rate of mitogenomic rearrangements. Implications of this model are particularly important for the application of gene order similarity as an additive source of phylogenetic information. Chromadorean nematodes, and particularly Camallanina clade (with C. cotti as an example of extremely accelerated rate of rearrangements), might be a good model to further study this discontinuity in the dynamics of mitogenomic evolution.

Gene Acquisition by a Distinct Phyletic Group within Streptococcus pneumoniae Promotes Adhesion to the Ocular Epithelium

Streptococcus pneumoniae (pneumococcus) displays broad tissue tropism and infects multiple body sites in the human host. However, infections of the conjunctiva are limited to strains within a distinct phyletic group with multilocus sequence types ST448, ST344, ST1186, ST1270, and ST2315. In this study, we sequenced the genomes of six pneumococcal strains isolated from eye infections. The conjunctivitis isolates are grouped in a distinct phyletic group together with a subset of nasopharyngeal isolates. The keratitis (infection of the cornea) and endophthalmitis (infection of the vitreous body) isolates are grouped with the remainder of pneumococcal strains. Phenotypic characterization is consistent with morphological differences associated with the distinct phyletic group. Specifically, isolates from the distinct phyletic group form aggregates in planktonic cultures and chain-like structures in biofilms grown on abiotic surfaces. To begin to investigate the association between genotype and epidemiology, we focused on a predicted surface-exposed adhesin (SspB) encoded exclusively by this distinct phyletic group. Phylogenetic analysis of the gene encoding SspB in the context of a streptococcal species tree suggests that sspB was acquired by lateral gene transfer from Streptococcus suis. Furthermore, an sspB deletion mutant displays decreased adherence to cultured cells from the ocular epithelium compared to the isogenic wild-type and complemented strains. Together these findings suggest that acquisition of genes from outside the species has contributed to pneumococcal tissue tropism by enhancing the ability of a subset of strains to infect the ocular epithelium causing conjunctivitis. IMPORTANCE Changes in the gene content of pathogens can modify their ability to colonize and/or survive in different body sites in the human host. In this study, we investigate a gene acquisition event and its role in the pathogenesis of Streptococccus pneumoniae (pneumococcus). Our findings suggest that the gene encoding the predicted surface protein SspB has been transferred from Streptococcus suis (a distantly related streptococcal species) into a distinct set of pneumococcal strains. This group of strains distinguishes itself from the remainder of pneumococcal strains by extensive differences in genomic composition and by the ability to cause conjunctivitis. We find that the presence of sspB increases adherence of pneumococcus to the ocular epithelium. Thus, our data support the hypothesis that a subset of pneumococcal strains has gained genes from neighboring species that enhance their ability to colonize the epithelium of the eye, thus expanding into a new niche.

A virus carries a gene encoding juvenile hormone acid methyltransferase, a key regulatory enzyme in insect metamorphosis

Microbial parasitism, infection, and symbiosis in animals often modulate host endocrine systems, resulting in alterations of phenotypic traits of the host that can have profound effects on the ecology and evolution of both the microorganisms and their hosts. Information about the mechanisms and genetic bases of such modulations by animal parasites is available from studies of steroid hormones. However, reports involving other hormones are scarce. We found that an insect virus, a betaentomopoxvirus, encodes a juvenile hormone acid methyltransferase that can synthesize an important insect hormone, the sesquiterpenoid juvenile hormone. Phylogenetic analysis suggested that this gene is of bacterial origin. Our study challenges the conventional view that functional enzymes in the late phase of the juvenile hormone biosynthesis pathway are almost exclusive to insects or arthropods, and shed light on juvenoid hormone synthesis beyond Eukaryota. This striking example demonstrates that even animal parasites having no metabolic pathways for molecules resembling host hormones can nevertheless influence the synthesis of such hormones, and provides a new context for studying animal parasite strategies in diverse systems such as host-parasite, host-symbiont or host-vector-parasite.

Bite it forward … bite it better? Incisor procumbency and mechanical advantage in the chisel-tooth and scratch-digger genus Ctenomys (Caviomorpha, Rodentia)

The subterranean genus Ctenomys (∼60 species, ∼100-1000g) constructs its burrows by using both forefeet and teeth throughout a wide range of habitats in South America. They show a high variation in the incisors' angle of attack (procumbency) and a mostly conserved skull morphology, not only amongst their congeners, but within the caviomorph rodents. Traditionally, procumbency has been largely related to tooth-digging. Looking for the possible influence of incisor procumbency on the mechanical advantage (MA) of each of the seven jaw adductor muscles in the genus, we examined 165 skulls representing 24 species. We also evaluated the role of two other relevant traits - i.e. mandibular width and diastema length - in jaw biomechanics and the existence of a relationship between procumbency angle and soil hardness. The in- and out-lever arms (L_i and L_o) of the involve muscles were determined based on their insertion's 3D-coordinates and integrated to calculate their MA. Interspecific scaling relationships for skull and muscle measurements were analyzed through reduced major axis regression performed with phylogenetically independent standardized contrasts. Although the procumbency angle ranged between ∼92.5° (C. mendocinus) and ∼107.2 (C. occultus), we found that it was not significantly correlated with the MA of any jaw adductor muscle. This study also showed that the incisor procumbency variation was not associated with the relative rostral length or soil hardness. This result contradicts previous generalizations about a correlation between habitat conditions and the procumbency of the incisors in subterranean rodents. In sum, our results suggest that, within Ctenomys, possessing more procumbent incisors may not represent a biomechanical advantage, but might be beneficial in other aspects related to chisel-tooth digging or food processing behaviors.

Membrane Trafficking Modulation during Entamoeba Encystation

Entamoeba histolytica is an intestinal parasite that infects 50-100 million people and causes up to 55,000 deaths annually. The transmissive form of E. histolytica is the cyst, with a single infected individual passing up to 45 million cysts per day, making cyst production an attractive target for infection control. Lectins and chitin are secreted to form the cyst wall, although little is known about the underlying membrane trafficking processes supporting encystation. As E. histolytica does not readily form cysts in vitro, we assessed membrane trafficking gene expression during encystation in the closely related model Entamoeba invadens. Genes involved in secretion are up-regulated during cyst formation, as are some trans-Golgi network-to-endosome trafficking genes. Furthermore, endocytic and general trafficking genes are up-regulated in the mature cyst, potentially preserved as mRNA in preparation for excystation. Two divergent dynamin-related proteins found in Entamoeba are predominantly expressed during cyst formation. Phylogenetic analyses indicate that they are paralogous to, but quite distinct from, classical dynamins found in human, suggesting that they may be potential drug targets to block encystation. The membrane-trafficking machinery is clearly regulated during encystation, providing an additional facet to understanding this crucial parasitic process.

Genetic Mapping and Biochemical Basis of Yellow Feather Pigmentation in Budgerigars

Parrot feathers contain red, orange, and yellow polyene pigments called psittacofulvins. Budgerigars are parrots that have been extensively bred for plumage traits during the last century, but the underlying genes are unknown. Here we use genome-wide association mapping and gene-expression analysis to map the Mendelian blue locus, which abolishes yellow pigmentation in the budgerigar. We find that the blue trait maps to a single amino acid substitution (R644W) in an uncharacterized polyketide synthase (MuPKS). When we expressed MuPKS heterologously in yeast, yellow pigments accumulated. Mass spectrometry confirmed that these yellow pigments match those found in feathers. The R644W substitution abolished MuPKS activity. Furthermore, gene-expression data from feathers of different bird species suggest that parrots acquired their colors through regulatory changes that drive high expression of MuPKS in feather epithelia. Our data also help formulate biochemical models that may explain natural color variation in parrots. VIDEO ABSTRACT.

A biologist's guide to Bayesian phylogenetic analysis

Bayesian methods have become very popular in molecular phylogenetics due to the availability of user-friendly software implementing sophisticated models of evolution. However, Bayesian phylogenetic models are complex, and analyses are often carried out using default settings, which may not be appropriate. Here, we summarize the major features of Bayesian phylogenetic inference and discuss Bayesian computation using Markov chain Monte Carlo (MCMC), the diagnosis of an MCMC run, and ways of summarising the MCMC sample. We discuss the specification of the prior, the choice of the substitution model, and partitioning of the data. Finally, we provide a list of common Bayesian phylogenetic software and provide recommendations as to their use.

Identification and preliminary characterization of a putative C-type lectin receptor-like protein in the T. cruzi tomato lectin endocytic-enriched proteome

Trypanosoma cruzi, the etiological agent of the Chagas' disease in Latin America undergoes a complex life cycle involving two hosts, a mammalian host and a reduviid insect vector (triatomine). In the insect midgut the parasite multiplies as epimastigote forms, which rely on endocytosis for their energy requirement. We recently showed that posttranslational modification of endocytic N-glycoproteins by tomato lectin (TL) binding-N-glycans is crucial for receptor-mediated endocytosis (RME) in epimastigote forms. In an attempt to characterize the endocytic proteome we used a TL affinity chromatography, which significantly enriched glycoproteins of the trypanosomal endocytic pathway. In addition to various lysosomal hydrolases, we found an endosomal C-type lectin-like protein, which displays some structural and topological characteristics of the mammalian lectin receptor superfamily. This lectin encoding a large transmembrane protein of around 375kDa contained three putative extracellular N-terminal C-type lectin domains (CTLD) and located inside the flagellar pocket (FP)/cytostome and endosomal compartments of the insect stage of the parasite and on the surface of the plasma membrane of intracellular amastigote parasites. Noteworthy, this endogenous lectin displayed similar sugar-binding specificity to that of TL and therefore could be important in either the N-glycan mediated endocytosis or parasite adhesion to host cells. We postulated that during the evolution of trypanosomatids, genes encoding lectin harboring 3 CTDLs represent an old acquisition present in free-living, monoxenic and heteroxenic trypanosomatids, which would have been secondarily lost in extracellular parasites from the T. brucei clade.

Nucleotide sequence and results of test of adaptive evolution in the α-globin gene of octodontoid rodents

The data presented in this article are related to the research article entitled “Molecular adaptive convergence in the α-globin gene in subterranean octodontid rodents” (Tomasco et al., 2017) [1]. This article shows the nucleotide sequences of α-globin subunit gene of hemoglobin of several South American caviomorph rodents, including subterranean and fossorial species. These sequences are deposited in Genbank, with accession numbers ranging from MF169881 to MF169898. Of a total of 429 nucleotides analyzed (143 codons), 100 variable sites and 43 amino acid replacements were reported. In this article we also show the results of TreeSaap (Woolley et al., 2003) [2] and MEME (Murrell et al., 2012) [3], that identified some replacement changes as interesting for future studies of adaptive evolution in this large rodent radiation.

Genome-wide identification of the MADS-box transcription factor family in pear (Pyrus bretschneideri) reveals evolution and functional divergence

MADS-box transcription factors play significant roles in plant developmental processes such as floral organ conformation, flowering time, and fruit development. Pear (Pyrus), as the third-most crucial temperate fruit crop, has been fully sequenced. However, there is limited information about the MADS family and its functional divergence in pear. In this study, a total of 95 MADS-box genes were identified in the pear genome, and classified into two types by phylogenetic analysis. Type I MADS-box genes were divided into three subfamilies and type II genes into 14 subfamilies. Synteny analysis suggested that whole-genome duplications have played key roles in the expansion of the MADS family, followed by rearrangement events. Purifying selection was the primary force driving MADS-box gene evolution in pear, and one gene pairs presented three codon sites under positive selection. Full-scale expression information for PbrMADS genes in vegetative and reproductive organs was provided and proved by transcriptional and reverse transcription PCR analysis. Furthermore, the PbrMADS11(12) gene, together with partners PbMYB10 and PbbHLH3 was confirmed to activate the promoters of the structural genes in anthocyanin pathway of red pear through dual luciferase assay. In addition, the PbrMADS11 and PbrMADS12 were deduced involving in the regulation of anthocyanin synthesis response to light and temperature changes. These results provide a solid foundation for future functional analysis of PbrMADS genes in different biological processes, especially of pigmentation in pear.

Conservation and diversification of small RNA pathways within flatworms

Background

Small non-coding RNAs, including miRNAs, and gene silencing mediated by RNA interference have been described in free-living and parasitic lineages of flatworms, but only few key factors of the small RNA pathways have been exhaustively investigated in a limited number of species. The availability of flatworm draft genomes and predicted proteomes allowed us to perform an extended survey of the genes involved in small non-coding RNA pathways in this phylum.

Results

Overall, findings show that the small non-coding RNA pathways are conserved in all the analyzed flatworm linages; however notable peculiarities were identified. While Piwi genes are amplified in free-living worms they are completely absent in all parasitic species. Remarkably all flatworms share a specific Argonaute family (FL-Ago) that has been independently amplified in different lineages. Other key factors such as Dicer are also duplicated, with Dicer-2 showing structural differences between trematodes, cestodes and free-living flatworms. Similarly, a very divergent GW182 Argonaute interacting protein was identified in all flatworm linages. Contrasting to this, genes involved in the amplification of the RNAi interfering signal were detected only in the ancestral free living species Macrostomum lignano. We here described all the putative small RNA pathways present in both free living and parasitic flatworm lineages.

Conclusion

These findings highlight innovations specifically evolved in platyhelminths presumably associated with novel mechanisms of gene expression regulation mediated by small RNA pathways that differ to what has been classically described in model organisms. Understanding these phylum-specific innovations and the differences between free living and parasitic species might provide clues to adaptations to parasitism, and would be relevant for gene-silencing technology development for parasitic flatworms that infect hundreds of million people worldwide.

Electronic supplementary material

The online version of this article (10.1186/s12862-017-1061-5) contains supplementary material, which is available to authorized users.

Complete Chloroplast Genome of Pinus massoniana (Pinaceae): Gene Rearrangements, Loss of ndh Genes, and Short Inverted Repeats Contraction, Expansion

The chloroplast genome (CPG) of Pinus massoniana belonging to the genus Pinus (Pinaceae), which is a primary source of turpentine, was sequenced and analyzed in terms of gene rearrangements, ndh genes loss, and the contraction and expansion of short inverted repeats (IRs). P. massoniana CPG has a typical quadripartite structure that includes large single copy (LSC) (65,563 bp), small single copy (SSC) (53,230 bp) and two IRs (IRa and IRb, 485 bp). The 108 unique genes were identified, including 73 protein-coding genes, 31 tRNAs, and 4 rRNAs. Most of the 81 simple sequence repeats (SSRs) identified in CPG were mononucleotides motifs of A/T types and located in non-coding regions. Comparisons with related species revealed an inversion (21,556 bp) in the LSC region; P. massoniana CPG lacks all 11 intact ndh genes (four ndh genes lost completely; the five remained truncated as pseudogenes; and the other two ndh genes remain as pseudogenes because of short insertions or deletions). A pair of short IRs was found instead of large IRs, and size variations among pine species were observed, which resulted from short insertions or deletions and non-synchronized variations between "IRa" and "IRb". The results of phylogenetic analyses based on whole CPG sequences of 16 conifers indicated that the whole CPG sequences could be used as a powerful tool in phylogenetic analyses.

Re-analyses of "Algal" Genes Suggest a Complex Evolutionary History of Oomycetes

The spread of photosynthesis is one of the most important but constantly debated topics in eukaryotic evolution. Various hypotheses have been proposed to explain the plastid distribution in extant eukaryotes. Notably, the chromalveolate hypothesis suggested that multiple eukaryotic lineages were derived from a photosynthetic ancestor that had a red algal endosymbiont. As such, genes of plastid/algal origin in aplastidic chromalveolates, such as oomycetes, were considered to be important supporting evidence. Although the chromalveolate hypothesis has been seriously challenged, some of its supporting evidence has not been carefully investigated. In this study, we re-evaluate the "algal" genes from oomycetes with a larger sampling and careful phylogenetic analyses. Our data provide no conclusive support for a common photosynthetic ancestry of stramenopiles, but show that the initial estimate of "algal" genes in oomycetes was drastically inflated due to limited genome data available then for certain eukaryotic lineages. These findings also suggest that the evolutionary histories of these "algal" genes might be attributed to complex scenarios such as differential gene loss, serial endosymbioses, or horizontal gene transfer.