Plastid phylogenomics resolves infrafamilial relationships of the Styracaceae and sheds light on the backbone relationships of the Ericales

Relationships among the genera of the small, woody family Styracaceae and among families of the large, diverse order Ericales have resisted complete resolution with sequences from one or a few genes. We used plastome sequencing to attempt to resolve the backbone relationships of Styracaceae and Ericales and to explore plastome structural evolution. Complete plastomes for 23 species are newly reported here, including 18 taxa of Styracaceae and five of Ericales (including species of Sapotaceae, Clethraceae, Symplocaceae, and Diapensiaceae). Combined with publicly available complete plastome data, this resulted in a data set of 60 plastomes, including 11 of the 12 genera of Styracaceae and 12 of 22 families of Ericales. Styracaceae plastomes were found to possess the quadripartite structure typical of angiosperms, with sizes ranging from 155 to 159 kb. Most of the plastomes were found to possess the full complement of typical angiosperm plastome genes. Unusual structural features were detected in plastomes of Alniphyllum and Bruinsmia, including the presence of a large 20-kb inversion (14 genes) in the Large Single-Copy region, the loss or pseudogenization of the clpP and accD genes in Bruinsmia, and the loss of the first exon of rps16 in B. styracoides. Likewise, the second intron from clpP was found to be lost in Alniphyllum and Huodendron. Phylogenomic analyses including all 79 plastid protein-coding genes provided improved resolution for relationships among the genera of Styracaceae and families of Ericales. Styracaceae was strongly supported as monophyletic, with Styrax, Huodendron, and a clade of Alniphyllum + Bruinsmia successively sister to the remainder of the family, all with strong support. All genera of Styracaceae were recovered as monophyletic, except for Halesia and Pterostyrax, which were each recovered as polyphyletic with strong support. Within Ericales, all families were recovered as monophyletic with strong support, with Balsaminaceae sister to remaining Ericales. Most relationships recovered in plastome analyses are congruent with previous analyses based on smaller data sets. Our results demonstrate the power of plastid phylogenomics to improve phylogenetic hypotheses among genera and families, and provide new insight into plastome evolution across Ericales.

Mouse Obox and Crxos modulate preimplantation transcriptional profiles revealing similarity between paralogous mouse and human homeobox genes

BACKGROUND

ETCHbox genes are eutherian-specific homeobox genes expressed during preimplantation development at a time when the first cell lineage decisions are being made. The mouse has an unusual repertoire of ETCHbox genes with several gene families lost in evolution and the remaining two, Crxos and Obox, greatly divergent in sequence and number. Each has undergone duplication to give a double homeodomain Crxos locus and a large cluster of over 60 Obox loci. The gene content differences between species raise important questions about how evolution can tolerate loss of genes implicated in key developmental events.

RESULTS

We find that Crxos internal duplication occurred in the mouse lineage, while Obox duplication was stepwise, generating subgroups with distinct sequence and expression. Ectopic expression of three Obox genes and a Crxos transcript in primary mouse embryonic cells followed by transcriptome sequencing allowed investigation into their functional roles. We find distinct transcriptomic influences for different Obox subgroups and Crxos, including modulation of genes related to zygotic genome activation and preparation for blastocyst formation. Comparison with similar experiments performed using human homeobox genes reveals striking overlap between genes downstream of mouse Crxos and genes downstream of human ARGFX.

CONCLUSIONS

Mouse Crxos and human ARGFX homeobox genes are paralogous rather than orthologous, yet they have evolved to regulate a common set of genes. This suggests there was compensation of function alongside gene loss through co-option of a different locus. Functional compensation by non-orthologous genes with dissimilar sequences is unusual but may indicate underlying distributed robustness. Compensation may be driven by the strong evolutionary pressure for successful early embryo development.

Avian genomics lends insights into endocrine function in birds

The genomics era has brought along the completed sequencing of a large number of bird genomes that cover a broad range of the avian phylogenetic tree (>30 orders), leading to major novel insights into avian biology and evolution. Among recent findings, the discovery that birds lack a large number of protein coding genes that are organized in highly conserved syntenic clusters in other vertebrates is very intriguing, given the physiological importance of many of these genes. A considerable number of them play prominent endocrine roles, suggesting that birds evolved compensatory genetic or physiological mechanisms that allowed them to survive and thrive in spite of these losses. While further studies are needed to establish the exact extent of avian gene losses, these findings point to birds as potentially highly relevant model organisms for exploring the genetic basis and possible therapeutic approaches for a wide range of endocrine functions and disorders.

Comparative analyses of plastid genomes from fourteen Cornales species: inferences for phylogenetic relationships and genome evolution

Background

The Cornales is the basal lineage of the asterids, the largest angiosperm clade. Phylogenetic relationships within the order were previously not fully resolved. Fifteen plastid genomes representing 14 species, ten genera and seven families of Cornales were newly sequenced for comparative analyses of genome features, evolution, and phylogenomics based on different partitioning schemes and filtering strategies.

Results

All plastomes of the 14 Cornales species had the typical quadripartite structure with a genome size ranging from 156,567 bp to 158,715 bp, which included two inverted repeats (25,859–26,451 bp) separated by a large single-copy region (86,089–87,835 bp) and a small single-copy region (18,250–18,856 bp) region. These plastomes encoded the same set of 114 unique genes including 31 transfer RNA, 4 ribosomal RNA and 79 coding genes, with an identical gene order across all examined Cornales species. Two genes (rpl22 and ycf15) contained premature stop codons in seven and five species respectively. The phylogenetic relationships among all sampled species were fully resolved with maximum support. Different filtering strategies (none, light and strict) of sequence alignment did not have an effect on these relationships. The topology recovered from coding and noncoding data sets was the same as for the whole plastome, regardless of filtering strategy. Moreover, mutational hotspots and highly informative regions were identified.

Conclusions

Phylogenetic relationships among families and intergeneric relationships within family of Cornales were well resolved. Different filtering strategies and partitioning schemes do not influence the relationships. Plastid genomes have great potential to resolve deep phylogenetic relationships of plants.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4319-9) contains supplementary material, which is available to authorized users.

Rhodopsin gene copies in Japanese eel originated in a teleost-specific genome duplication

BACKGROUND

Gene duplication is considered important to increasing the genetic diversity in animals. In fish, visual pigment genes are often independently duplicated, and the evolutionary significance of such duplications has long been of interest. Eels have two rhodopsin genes (rho), one of which (freshwater type, fw-rho) functions in freshwater and the other (deep-sea type, ds-rho) in marine environments. Hence, switching of rho expression in retinal cells is tightly linked with eels' unique life cycle, in which they migrate from rivers or lakes to the sea. These rho genes are apparently paralogous, but the timing of their duplication is unclear due to the deep-branching phylogeny. The aim of the present study is to elucidate the evolutionary origin of the two rho copies in eels using comparative genomics methods.

RESULTS

In the present study, we sequenced the genome of Japanese eel Anguilla japonica and reconstructed two regions containing rho by de novo assembly. We found a single corresponding region in a non-teleostean primitive ray-finned fish (spotted gar) and two regions in a primitive teleost (Asian arowana). The order of ds-rho and the neighboring genes was highly conserved among the three species. With respect to fw-rho, which was lost in Asian arowana, the neighboring genes were also syntenic between Japanese eel and Asian arowana. In particular, the pattern of gene losses in ds-rho and fw-rho regions was the same as that in Asian arowana, and no discrepancy was found in any of the teleost genomes examined. Phylogenetic analysis supports mutual monophyly of these two teleostean synteny groups, which correspond to the ds-rho and fw-rho regions.

CONCLUSIONS

Syntenic and phylogenetic analyses suggest that the duplication of rhodopsin gene in Japanese eel predated the divergence of eel (Elopomorpha) and arowana (Osteoglossomorpha). Thus, based on the principle of parsimony, it is most likely that the rhodopsin paralogs were generated through a whole genome duplication in the ancestor of teleosts, and have remained till the present in eels with distinct functional roles. Our result indicates, for the first time, that teleost-specific genome duplication may have contributed to a gene innovation involved in eel-specific migratory life cycle.

Progressive erosion of the Relaxin1 gene in bovids

The relaxin/insulin-like (RLN/INSL) gene family is a group of genes that encode peptide hormones involved in a variety of physiological functions related to reproduction. Previous studies have shown that relaxin plays a key role in widening of the pubic bone during labor and in gamete maturation. Because of these functions, studying the evolution of RLN1, the gene encoding for relaxin, is relevant in livestock species, most of which belong in the group Laurasiatheria, which includes cow, pig, horse, goat, and sheep in addition to bats, cetaceans and carnivores. Experimental evidence suggests that cows do not synthesize relaxin, but respond to it, and sheep apparently have a truncated RLN1 gene. Thus, we made use of genome sequence data to characterize the genomic locus of the RLN1 gene in Laurasiatherian mammals to better understand how cows lost the ability to synthesize this peptide. We found that all ruminants in our study (cow, giraffe, goat, sheep and Tibetan antelope) lack a functional RLN1 gene, and document the progressive loss of RLN1 in the lineage leading to cows. Our analyses indicate that 1 - all ruminants have lost all key regulatory elements upstream of the first exon, 2 - giraffe, goat, sheep and Tibetan antelope have multiple inactivating mutations in the RLN1 pseudogene, and 3 - the cow genome has lost all traces of RLN1. The 5' regulatory sequence plays a key role in activating expression, and the loss of this sequence would impair synthesis of mRNA. Our results suggest that changes in regulatory sequence preceded mutations in coding sequence and highlight the importance of these regions in maintaining proper gene function. In addition, we found that all bovids examined posses copies of the relaxin receptors, which explains why they are able to respond to relaxin despite their inability to produce it.

Complete plastome sequencing of both living species of Circaeasteraceae (Ranunculales) reveals unusual rearrangements and the loss of the ndh gene family

BACKGROUND

Among the 13 families of early-diverging eudicots, only Circaeasteraceae (Ranunculales), which consists of the two monotypic genera Circaeaster and Kingdonia, lacks a published complete plastome sequence. In addition, the phylogenetic position of Circaeasteraceae as sister to Lardizabalaceae has only been weakly or moderately supported in previous studies using smaller data sets. Moreover, previous plastome studies have documented a number of novel structural rearrangements among early-divergent eudicots. Hence it is important to sequence plastomes from Circaeasteraceae to better understand plastome evolution in early-diverging eudicots and to further investigate the phylogenetic position of Circaeasteraceae.

RESULTS

Using an Illumina HiSeq 2000, complete plastomes were sequenced from both living members of Circaeasteraceae: Circaeaster agrestis and Kingdonia uniflora . Plastome structure and gene content were compared between these two plastomes, and with those of other early-diverging eudicot plastomes. Phylogenetic analysis of a 79-gene, 99-taxon data set including exemplars of all families of early-diverging eudicots was conducted to resolve the phylogenetic position of Circaeasteraceae. Both plastomes possess the typical quadripartite structure of land plant plastomes. However, a large ~49 kb inversion and a small ~3.5 kb inversion were found in the large single-copy regions of both plastomes, while Circaeaster possesses a number of other rearrangements, particularly in the Inverted Repeat. In addition, infA was found to be a pseudogene and accD was found to be absent within Circaeaster, whereas all ndh genes, except for ndhE and ndhJ, were found to be either pseudogenized (ΨndhA, ΨndhB, ΨndhD, ΨndhH and ΨndhK) or absent (ndhC, ndhF, ndhI and ndhG) in Kingdonia. Circaeasteraceae was strongly supported as sister to Lardizabalaceae in phylogenetic analyses.

CONCLUSION

The first plastome sequencing of Circaeasteraceae resulted in the discovery of several unusual rearrangements and the loss of ndh genes, and confirms the sister relationship between Circaeasteraceae and Lardizabalaceae. This research provides new insight to characterize plastome structural evolution in early-diverging eudicots and to better understand relationships within Ranunculales .

Comparative genomics analyses of alpha-keratins reveal insights into evolutionary adaptation of marine mammals

Background

Diversity of hair in marine mammals was suggested as an evolutionary innovation to adapt aquatic environment, yet its genetic basis remained poorly explored. We scanned α-keratin genes, one major structural components of hair, in 16 genomes of mammalian species, including seven cetaceans, two pinnipeds, polar bear, manatee and five terrestrial species.

Results

Extensive gene loss and high pseudogenization rate of α-keratin genes were identified in cetaceans when compared to terrestrial artiodactylans (average number of α-keratins 37.29 vs. 58.33; pseudogenization rate 29.89% vs. 8.00%), especially of hair follicle-specific keratin genes (average pseudogenization rate in cetaceans of 43.88% relative to 3.80% artiodactylian average). Compared to toothed whale, the much more number of intact functional α-keratin genes was examined in the baleen whale that had specific keratinized baleen. In contrast, the number of keratin genes in pinnipeds, polar bear and manatee were comparable to those of their respective terrestrial relatives. Additionally, four keratin genes (K39, K9, K42, and K74) were found to be pseudogenes or lost uniquely in cetaceans and manatees.

Conclusions

Species-specific evolution of α-keratin gene family identified in the marine mammals might be responsible for their different hair characteristics. Increased gene loss and pseudogenization rate identified in cetacean lineages was likely to contribute to hair-less phenotype to adaptation for complete aquatic environment. However, the fully aquatic manatee still remained the comparable number of intact genes to its terrestrial relative, probably due to its perioral bristles and bristle-like hairs on the oral disk. By contrast, similar evolution pattern of α-keratin gene repertoire in the pinnipeds, polar bear and their terrestrial relatives was likely due to abundant hair to keep warm when they went ashore. Interestingly, some keratin genes were exclusively lost in cetaceans and manatees, likely as a result of convergent hair-loss phenotype to inhabit completely aquatic environment in both groups.

Electronic supplementary material

The online version of this article (doi:10.1186/s12983-017-0225-x) contains supplementary material, which is available to authorized users.

The phytochelatin synthase gene in date palm (Phoenix dactylifera L.): Phylogeny, evolution and expression

We studied date palm phytochelatin synthase type I (PdPCS1), which catalyzes the cytosolic synthesis of phytochelatins (PCs), a heavy metal binding protein, in plant cells. The gene encoding PdPCS1 (Pdpcs) consists of 8 exons and 7 introns and encodes a protein of 528 amino acids. PCs gene history was studied using Notung phylogeny. During evolution, gene loss from several lineages was predicted including Proteobacteria, Bilateria and Brassicaceae. In addition, eleven gene duplication events appeared toward interior nodes of the reconciled tree and four gene duplication events appeared toward the external nodes. These latter sequences belong to species with a second copy of PCs suggesting that this gene evolved through subfunctionalization. Pdpcs1 gene expression was measured in seedling hypocotyls exposed to Cd, Cu and Cr using quantitative real-time polymerase chain reaction (qPCR). A Pdpcs1 overexpression was evidenced in P. dactylifera seedlings exposed to metals suggesting that 1-the Pdpcs1 gene is functional, 2-there is an implication of the enzyme in metal detoxification mechanisms. Additionally, the structure of PdPCS1 was predicted using its homologue from Nostoc (cyanobacterium, NsPCS) as a template in Discovery studio and PyMol software. These analyses allowed us to identify the phytochelatin synthase type I enzyme in date palm (PdPCS1) via recognition of key consensus amino acids involved in the catalytic mechanism, and to propose a hypothetical binding and catalytic site for an additional substrate binding cavity.

Statistical analysis of fractionation resistance by functional category and expression

Background

The current literature establishes the importance of gene functional category and expression in promoting or suppressing duplicate gene loss after whole genome doubling in plants, a process known as fractionation. Inspired by studies that have reported gene expression to be the dominating factor in preventing duplicate gene loss, we analyzed the relative effect of functional category and expression.

Methods

We use multivariate methods to study data sets on gene retention, function and expression in rosids and asterids to estimate effects and assess their interaction.

Results

Our results suggest that the effect on duplicate gene retention fractionation by functional category and expression are independent and have no statistical interaction.

Conclusion

In plants, functional category is the more dominant factor in explaining duplicate gene loss.

Evolutionary origin and function of NOX4-art, an arthropod specific NADPH oxidase

BACKGROUND

NADPH oxidases (NOX) are ROS producing enzymes that perform essential roles in cell physiology, including cell signaling and antimicrobial defense. This gene family is present in most eukaryotes, suggesting a common ancestor. To date, only a limited number of phylogenetic studies of metazoan NOXes have been performed, with few arthropod genes. In arthropods, only NOX5 and DUOX genes have been found and a gene called NOXm was found in mosquitoes but its origin and function has not been examined. In this study, we analyzed the evolution of this gene family in arthropods. A thorough search of genomes and transcriptomes was performed enabling us to browse most branches of arthropod phylogeny.

RESULTS

We have found that the subfamilies NOX5 and DUOX are present in all arthropod groups. We also show that a NOX gene, closely related to NOX4 and previously found only in mosquitoes (NOXm), can also be found in other taxonomic groups, leading us to rename it as NOX4-art. Although the accessory protein p22-phox, essential for NOX1-4 activation, was not found in any of the arthropods studied, NOX4-art of Aedes aegypti encodes an active protein that produces H₂O₂. Although NOX4-art has been lost in a number of arthropod lineages, it has all the domains and many signature residues and motifs necessary for ROS production and, when silenced, H₂O₂ production is considerably diminished in A. aegypti cells.

CONCLUSIONS

Combining all bioinformatic analyses and laboratory work we have reached interesting conclusions regarding arthropod NOX gene family evolution. NOX5 and DUOX are present in all arthropod lineages but it seems that a NOX2-like gene was lost in the ancestral lineage leading to Ecdysozoa. The NOX4-art gene originated from a NOX4-like ancestor and is functional. Although no p22-phox was observed in arthropods, there was no evidence of neo-functionalization and this gene probably produces H₂O₂ as in other metazoan NOX4 genes. Although functional and present in the genomes of many species, NOX4-art was lost in a number of arthropod lineages.

Plastome phylogeny and early diversification of Brassicaceae

Background

The family Brassicaceae encompasses diverse species, many of which have high scientific and economic importance. Early diversifications and phylogenetic relationships between major lineages or clades remain unclear. Here we re-investigate Brassicaceae phylogeny with complete plastomes from 51 species representing all four lineages or 5 of 6 major clades (A, B, C, E and F) as identified in earlier studies.

Results

Bayesian and maximum likelihood phylogenetic analyses using a partitioned supermatrix of 77 protein coding genes resulted in nearly identical tree topologies exemplified by highly supported relationships between clades. All four lineages were well identified and interrelationships between them were resolved. The previously defined Clade C was found to be paraphyletic (the genus Megadenia formed a separate lineage), while the remaining clades were monophyletic. Clade E (lineage III) was sister to clades B + C rather than to all core Brassicaceae (clades A + B + C or lineages I + II), as suggested by a previous transcriptome study. Molecular dating based on plastome phylogeny supported the origin of major lineages or clades between late Oligocene and early Miocene, and the following radiative diversification across the family took place within a short timescale. In addition, gene losses in the plastomes occurred multiple times during the evolutionary diversification of the family.

Conclusions

Plastome phylogeny illustrates the early diversification of cruciferous species. This phylogeny will facilitate our further understanding of evolution and adaptation of numerous species in the model family Brassicaceae.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3555-3) contains supplementary material, which is available to authorized users.

Plastome phylogeny and early diversification of Brassicaceae

BACKGROUND

The family Brassicaceae encompasses diverse species, many of which have high scientific and economic importance. Early diversifications and phylogenetic relationships between major lineages or clades remain unclear. Here we re-investigate Brassicaceae phylogeny with complete plastomes from 51 species representing all four lineages or 5 of 6 major clades (A, B, C, E and F) as identified in earlier studies.

RESULTS

Bayesian and maximum likelihood phylogenetic analyses using a partitioned supermatrix of 77 protein coding genes resulted in nearly identical tree topologies exemplified by highly supported relationships between clades. All four lineages were well identified and interrelationships between them were resolved. The previously defined Clade C was found to be paraphyletic (the genus Megadenia formed a separate lineage), while the remaining clades were monophyletic. Clade E (lineage III) was sister to clades B + C rather than to all core Brassicaceae (clades A + B + C or lineages I + II), as suggested by a previous transcriptome study. Molecular dating based on plastome phylogeny supported the origin of major lineages or clades between late Oligocene and early Miocene, and the following radiative diversification across the family took place within a short timescale. In addition, gene losses in the plastomes occurred multiple times during the evolutionary diversification of the family.

CONCLUSIONS

Plastome phylogeny illustrates the early diversification of cruciferous species. This phylogeny will facilitate our further understanding of evolution and adaptation of numerous species in the model family Brassicaceae.

Evolution of gremlin 2 in cetartiodactyl mammals: gene loss coincides with lack of upper jaw incisors in ruminants

Understanding the processes that give rise to genomic variability in extant species is an active area of research within evolutionary biology. With the availability of whole genome sequences, it is possible to quantify different forms of variability such as variation in gene copy number, which has been described as an important source of genetic variability and in consequence of phenotypic variability. Most of the research on this topic has been focused on understanding the biological significance of gene duplication, and less attention has been given to the evolutionary role of gene loss. Gremlin 2 is a member of the DAN gene family and plays a significant role in tooth development by blocking the ligand-signaling pathway of BMP2 and BMP4. The goal of this study was to investigate the evolutionary history of gremlin 2 in cetartiodactyl mammals, a group that possesses highly divergent teeth morphology. Results from our analyses indicate that gremlin 2 has experienced a mixture of gene loss, gene duplication, and rate acceleration. Although the last common ancestor of cetartiodactyls possessed a single gene copy, pigs and camels are the only cetartiodactyl groups that have retained gremlin 2. According to the phyletic distribution of this gene and synteny analyses, we propose that gremlin 2 was lost in the common ancestor of ruminants and cetaceans between 56.3 and 63.5 million years ago as a product of a chromosomal rearrangement. Our analyses also indicate that the rate of evolution of gremlin 2 has been accelerated in the two groups that have retained this gene. Additionally, the lack of this gene could explain the high diversity of teeth among cetartiodactyl mammals; specifically, the presence of this gene could act as a biological constraint. Thus, our results support the notions that gene loss is a way to increase phenotypic diversity and that gremlin 2 is a dispensable gene, at least in cetartiodactyl mammals.

The genome of the Gulf pipefish enables understanding of evolutionary innovations

Background

Evolutionary origins of derived morphologies ultimately stem from changes in protein structure, gene regulation, and gene content. A well-assembled, annotated reference genome is a central resource for pursuing these molecular phenomena underlying phenotypic evolution. We explored the genome of the Gulf pipefish (Syngnathus scovelli), which belongs to family Syngnathidae (pipefishes, seahorses, and seadragons). These fishes have dramatically derived bodies and a remarkable novelty among vertebrates, the male brood pouch.

Results

We produce a reference genome, condensed into chromosomes, for the Gulf pipefish. Gene losses and other changes have occurred in pipefish hox and dlx clusters and in the tbx and pitx gene families, candidate mechanisms for the evolution of syngnathid traits, including an elongated axis and the loss of ribs, pelvic fins, and teeth. We measure gene expression changes in pregnant versus non-pregnant brood pouch tissue and characterize the genomic organization of duplicated metalloprotease genes (patristacins) recruited into the function of this novel structure. Phylogenetic inference using ultraconserved sequences provides an alternative hypothesis for the relationship between orders Syngnathiformes and Scombriformes. Comparisons of chromosome structure among percomorphs show that chromosome number in a pipefish ancestor became reduced via chromosomal fusions.

Conclusions

The collected findings from this first syngnathid reference genome open a window into the genomic underpinnings of highly derived morphologies, demonstrating that de novo production of high quality and useful reference genomes is within reach of even small research groups.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-016-1126-6) contains supplementary material, which is available to authorized users.

The loss of photosynthetic pathways in the plastid and nuclear genomes of the non-photosynthetic mycoheterotrophic eudicot Monotropa hypopitys

BACKGROUND

Chloroplasts of most plants are responsible for photosynthesis and contain a conserved set of about 110 genes that encode components of housekeeping gene expression machinery and photosynthesis-related functions. Heterotrophic plants obtaining nutrients from other organisms and their plastid genomes represent model systems in which to study the effects of relaxed selective pressure on photosynthetic function. The most evident is a reduction in the size and gene content of the plastome, which correlates with the loss of genes encoding photosynthetic machinery which become unnecessary. Transition to a non-photosynthetic lifestyle is expected also to relax the selective pressure on photosynthetic machinery in the nuclear genome, however, the corresponding changes are less known.

RESULTS

Here we report the complete sequence of the plastid genome of Monotropa hypopitys, an achlorophyllous obligately mycoheterotrophic plant belonging to the family Ericaceae. The plastome of M. hypopitys is greatly reduced in size (35,336 bp) and lacks the typical quadripartite structure with two single-copy regions and an inverted repeat. Only 45 genes remained presumably intact- those encoding ribosomal proteins, ribosomal and transfer RNA and housekeeping genes infA, matK, accD and clpP. The clpP and accD genes probably remain functional, although their sequences are highly diverged. The sets of genes for ribosomal protein and transfer RNA are incomplete relative to chloroplasts of a photosynthetic plant. Comparison of the plastid genomes of two subspecies-level isolates of M. hypopitys revealed major structural rearrangements associated with repeat-driven recombination and the presence of isolate-specific tRNA genes. Analysis of the M. hypopitys transcriptome by RNA-Seq showed the absence of expression of nuclear-encoded components of photosystem I and II reaction center proteins, components of cytochrome b6f complex, ATP synthase, ribulose bisphosphate carboxylase components, as well as chlorophyll from protoporphyrin IX biosynthesis pathway.

CONCLUSIONS

With the complete loss of genes related to photosynthesis, NADH dehydrogenase, plastid-encoded RNA polymerase and ATP synthase, the M. hypopitys plastid genome is among the most functionally reduced ones characteristic of obligate non-photosynthetic parasitic species. Analysis of the M. hypopitys transcriptome revealed coordinated evolution of the nuclear and plastome genomes and the loss of photosynthesis-related functions in both genomes.

Pinpointing genes underlying annual/perennial transitions with comparative genomics

BACKGROUND

Transitions between perennial and an annual life history occur often in plant lineages, but the genes that control whether a plant is an annual or perennial are largely unknown. To identify genes that confer differences between annuals and perennials we compared the gene content of four pairs of sister lineages (Arabidopsis thaliana/Arabidopsis lyrata, Arabis montbretiana/Arabis alpina, Arabis verna/Aubrieta parviflora and Draba nemorosa/Draba hispanica) in the Brassicaceae in which each pair contains one annual and one perennial, plus one extra annual species (Capsella rubella).

RESULTS

After sorting all genes in all nine species into gene families, we identified five families in which well-annotated genes are present in the perennials A. lyrata and A. alpina, but are not present in any of the annual species. For the eleven genes in perennials in these families, an orthologous pseudogene or otherwise highly diverged gene was found in the syntenic region of the annual species in six cases. The five candidate families identified encode: a kinase, an oxidoreductase, a lactoylglutathione lyase, a F-box protein and a zinc finger protein. By comparing the active gene in the perennial to the pseudogene or heavily altered gene in the annual, dN and dS were calculated. The low dN/dS values in one kinase suggest that it became pseudogenized more recently, while the other kinase, F-box, oxidoreductase and zinc-finger became pseudogenized closer to the divergence between the annual-perennial pair.

CONCLUSIONS

We identified five gene families that may be involved in the life history switch from perennial to annual. Considering the dN and dS data and whether syntenic pseudogenes were found and the potential functions of the genes, the F-box family is considered the most promising candidate for future functional studies to determine if it affects life history.

Acquisition and loss of virulence-associated factors during genome evolution and speciation in three clades of Bordetella species

Background

The genus Bordetella consists of nine species that include important respiratory pathogens such as the ‘classical’ species B. bronchiseptica, B. pertussis and B. parapertussis and six more distantly related and less extensively studied species. Here we analyze sequence diversity and gene content of 128 genome sequences from all nine species with focus on the evolution of virulence-associated factors.

Results

Both genome-wide sequence-based and gene content-based phylogenetic trees divide the genus into three species clades. The phylogenies are congruent between species suggesting genus-wide co-evolution of sequence diversity and gene content, but less correlated within species, mainly because of strain-specific presence of many different prophages. We compared the genomes with focus on virulence-associated genes and identified multiple clade-specific, species-specific and strain-specific events of gene acquisition and gene loss, including genes encoding O-antigens, protein secretion systems and bacterial toxins. Gene loss was more frequent than gene gain throughout the evolution, and loss of hundreds of genes was associated with the origin of several species, including the recently evolved human-restricted B. pertussis and B. holmesii, B. parapertussis and the avian pathogen B. avium.

Conclusions

Acquisition and loss of multiple genes drive the evolution and speciation in the genus Bordetella, including large scale gene loss associated with the origin of several species. Recent loss and functional inactivation of genes, including those encoding pertussis vaccine components and bacterial toxins, in individual strains emphasize ongoing evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3112-5) contains supplementary material, which is available to authorized users.

One pedigree we all may have come from - did Adam and Eve have the chromosome 2 fusion?

Background

In contrast to Great Apes, who have 48 chromosomes, modern humans and likely Neandertals and Denisovans have and had, respectively, 46 chromosomes. The reduction in chromosome number was caused by the head-to-head fusion of two ancestral chromosomes to form human chromosome 2 (HSA2) and may have contributed to the reproductive barrier with Great Apes.

Results

Next generation sequencing and molecular clock analyses estimated that this fusion arose prior to our last common ancestor with Neandertal and Denisovan hominins ~ 0.74 - 4.5 million years ago.

Hypotheses

I propose that, unlike recurrent Robertsonian translocations in humans, the HSA2 fusion was a single nonrecurrent event that spread through a small polygamous clan population bottleneck. Its heterozygous to homozygous conversion, fixation, and accumulation in the succeeding populations was likely facilitated by an evolutionary advantage through the genomic loss rather than deregulation of expression of the gene(s) flanking the HSA2 fusion site at 2q13.

Conclusions

The origin of HSA2 might have been a critical evolutionary event influencing higher cognitive functions in various early subspecies of hominins. Next generation sequencing of Homo heidelbergensis and Homo erectus genomes and complete reconstruction of DNA sequence of the orthologous subtelomeric chromosomes in Great Apes should enable more precise timing of HSA2 formation and better understanding of its evolutionary consequences.

Catalase in Leishmaniinae: With me or against me?

The catalase gene is a virtually ubiquitous component of the eukaryotic genomes. It is also present in the monoxenous (i.e. parasitizing solely insects) trypanosomatids of the subfamily Leishmaniinae, which have acquired the enzyme by horizontal gene transfer from a bacterium. However, as shown here, the catalase gene was secondarily lost from the genomes of all Leishmania sequenced so far. Due to the potentially key regulatory role of hydrogen peroxide in the inter-stagial transformation of Leishmania spp., this loss seems to be a necessary prerequisite for the emergence of a complex life cycle of these important human pathogens. Hence, in this group of protists, the advantages of keeping catalase were uniquely outweighed by its disadvantages.

Evolution of RLSB, a nuclear-encoded S1 domain RNA binding protein associated with post-transcriptional regulation of plastid-encoded rbcL mRNA in vascular plants

Background

RLSB, an S-1 domain RNA binding protein of Arabidopsis, selectively binds rbcL mRNA and co-localizes with Ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) within chloroplasts of C₃ and C₄ plants. Previous studies using both Arabidopsis (C₃) and maize (C₄) suggest RLSB homologs are post-transcriptional regulators of plastid-encoded rbcL mRNA. While RLSB accumulates in all Arabidopsis leaf chlorenchyma cells, in C₄ leaves RLSB-like proteins accumulate only within Rubisco-containing bundle sheath chloroplasts of Kranz-type species, and only within central compartment chloroplasts in the single cell C₄ plant Bienertia. Our recent evidence implicates this mRNA binding protein as a primary determinant of rbcL expression, cellular localization/compartmentalization, and photosynthetic function in all multicellular green plants. This study addresses the hypothesis that RLSB is a highly conserved Rubisco regulatory factor that occurs in the chloroplasts all higher plants.

Results

Phylogenetic analysis has identified RLSB orthologs and paralogs in all major plant groups, from ancient liverworts to recent angiosperms. RLSB homologs were also identified in algae of the division Charophyta, a lineage closely related to land plants. RLSB-like sequences were not identified in any other algae, suggesting that it may be specific to the evolutionary line leading to land plants. The RLSB family occurs in single copy across most angiosperms, although a few species with two copies were identified, seemingly randomly distributed throughout the various taxa, although perhaps correlating in some cases with known ancient whole genome duplications. Monocots of the order Poales (Poaceae and Cyperaceae) were found to contain two copies, designated here as RLSB-a and RLSB-b, with only RLSB-a implicated in the regulation of rbcL across the maize developmental gradient. Analysis of microsynteny in angiosperms revealed high levels of conservation across eudicot species and for both paralogs in grasses, highlighting the possible importance of maintaining this gene and its surrounding genomic regions.

Conclusions

Findings presented here indicate that the RLSB family originated as a unique gene in land plant evolution, perhaps in the common ancestor of charophytes and higher plants. Purifying selection has maintained this as a highly conserved single- or two-copy gene across most extant species, with several conserved gene duplications. Together with previous findings, this study suggests that RLSB has been sustained as an important regulatory protein throughout the course of land plant evolution. While only RLSB-a has been directly implicated in rbcL regulation in maize, RLSB-b could have an overlapping function in the co-regulation of rbcL, or may have diverged as a regulator of one or more other plastid-encoded mRNAs. This analysis confirms that RLSB is an important and unique photosynthetic regulatory protein that has been continuously expressed in land plants as they emerged and diversified from their ancient common ancestor.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0713-1) contains supplementary material, which is available to authorized users.

Comparative analyses of developmental transcription factor repertoires in sponges reveal unexpected complexity of the earliest animals

Developmental transcription factors (DTFs) control development of animals by affecting expression of target genes, some of which are transcription factors themselves. In bilaterians and cnidarians, conserved DTFs are involved in homologous processes such as gastrulation or specification of neurons. The genome of Amphimedon queenslandica, the first sponge to be sequenced, revealed that only a fraction of these conserved DTF families are present in demosponges. This finding was in line with the view that morphological complexity in the animal lineage correlates with developmental toolkit complexity. However, as the phylum Porifera is very diverse, Amphimedon's genome may not be representative of all sponges. The recently sequenced genomes of calcareous sponges Sycon ciliatum and Leucosolenia complicata allowed investigations of DTFs in a sponge lineage evolutionarily distant from demosponges. Surprisingly, the phylogenetic analyses of identified DTFs revealed striking differences between the calcareous sponges and Amphimedon. As these differences appear to be a result of independent gene loss events in the two sponge lineages, the last common ancestor of sponges had to possess a much more diverse repertoire of DTFs than extant sponges. Developmental expression of sponge homologs of genes involved in specification of the Bilaterian endomesoderm and the neurosensory cells suggests that roles of many DTFs date back to the last common ancestor of all animals. Strikingly, even DTFs displaying apparent pan-metazoan conservation of sequence and function are not immune to being lost from individual species genomes. The quest for a comprehensive picture of the developmental toolkit in the last common metazoan ancestor is thus greatly benefitting from the increasing accessibility of sequencing, allowing comparisons of multiple genomes within each phylum.

Lineage-specific loss of FGF17 within the avian orders Galliformes and Passeriformes

The genomic and developmental complexity of vertebrates is commonly attributed to two rounds of whole genome duplications which occurred at the base of the vertebrate radiation. These duplications led to the rise of several, multi-gene families of developmental proteins like the fibroblast growth factors (FGFs); a signaling protein family which functions at various stages of embryonic development. One of the major FGF assemblages arising from these duplications is the FGF8 subfamily, which includes FGF8, FGF17, and FGF18 in tetrapods. While FGF8 and FGF18 are found in all tetrapods and are critical for embryonic survival, genomic analyses suggest putative loss of FGF17 in various lineages ranging from frogs and fish, to the chicken. This study utilizes 27 avian genomes in conjunction with molecular analyses of chicken embryos to confirm the loss of FGF17 in chicken as a true, biological occurrence. FGF17 is also missing in the turkey, black grouse, Japanese quail and northern bobwhite genomes. These species, along with chicken, form a monophyletic clade in the order Galliformes. Four additional species, members of the clade Passeroidea, within the order Passeriformes, are also missing FGF17. Additionally, analysis of intact FGF17 in other avian lineages reveals that it is still under strong purifying selection, despite being seemingly dispensable. Thus, FGF17 likely represents a molecular spandrel arising from a genome duplication event and due to its high connectivity with FGF8/FGF18, and potential for interference with their function, is retained under strong purifying selection, despite itself not having a strong selective advantage.