Survey sequencing and comparative analysis of the elephant shark (Callorhinchus milii) genome

Owing to their phylogenetic position, cartilaginous fishes (sharks, rays, skates, and chimaeras) provide a critical reference for our understanding of vertebrate genome evolution. The relatively small genome of the elephant shark, Callorhinchus milii, a chimaera, makes it an attractive model cartilaginous fish genome for whole-genome sequencing and comparative analysis. Here, the authors describe survey sequencing (1.4× coverage) and comparative analysis of the elephant shark genome, one of the first cartilaginous fish genomes to be sequenced to this depth. Repetitive sequences, represented mainly by a novel family of short interspersed element–like and long interspersed element–like sequences, account for about 28% of the elephant shark genome. Fragments of approximately 15,000 elephant shark genes reveal specific examples of genes that have been lost differentially during the evolution of tetrapod and teleost fish lineages. Interestingly, the degree of conserved synteny and conserved sequences between the human and elephant shark genomes are higher than that between human and teleost fish genomes. Elephant shark contains putative four Hox clusters indicating that, unlike teleost fish genomes, the elephant shark genome has not experienced an additional whole-genome duplication. These findings underscore the importance of the elephant shark as a critical reference vertebrate genome for comparative analysis of the human and other vertebrate genomes. This study also demonstrates that a survey-sequencing approach can be applied productively for comparative analysis of distantly related vertebrate genomes.     

The Role of Gene Duplication and Unconstrained Selective Pressures in the Melanopsin Gene Family Evolution and Vertebrate Circadian Rhythm Regulation

Melanopsin is a photosensitive cell protein involved in regulating circadian rhythms and other non-visual responses to light. The melanopsin gene family is represented by two paralogs, OPN4x and OPN4m, which originated through gene duplication early in the emergence of vertebrates. Here we studied the melanopsin gene family using an integrated gene/protein evolutionary approach, which revealed that the rhabdomeric urbilaterian ancestor had the same amino acid patterns (DRY motif and the Y and E conterions) as extant vertebrate species, suggesting that the mechanism for light detection and regulation is similar to rhabdomeric rhodopsins. Both OPN4m and OPN4x paralogs are found in vertebrate genomic paralogons, suggesting that they diverged following this duplication event about 600 million years ago, when the complex eye emerged in the vertebrate ancestor. Melanopsins generally evolved under negative selection (ω = 0.171) with some minor episodes of positive selection (proportion of sites = 25%) and functional divergence (θ_I = 0.349 and θ_II = 0.126). The OPN4m and OPN4x melanopsin paralogs show evidence of spectral divergence at sites likely involved in melanopsin light absorbance (200F, 273S and 276A). Also, following the teleost lineage-specific whole genome duplication (3R) that prompted the teleost fish radiation, type I divergence (θ_I = 0.181) and positive selection (affecting 11% of sites) contributed to amino acid variability that we related with the photo-activation stability of melanopsin. The melanopsin intracellular regions had unexpectedly high variability in their coupling specificity of G-proteins and we propose that Gq/11 and Gi/o are the two G-proteins most-likely to mediate the melanopsin phototransduction pathway. The selection signatures were mainly observed on retinal-related sites and the third and second intracellular loops, demonstrating the physiological plasticity of the melanopsin protein group. Our results provide new insights on the phototransduction process and additional tools for disentangling and understanding the links between melanopsin gene evolution and the specializations observed in vertebrates, especially in teleost fish.     

Multiple Chromosomal Rearrangements Structured the Ancestral Vertebrate Hox-Bearing Protochromosomes

While the proposal that large-scale genome expansions occurred early in vertebrate evolution is widely accepted, the exact mechanisms of the expansion—such as a single or multiple rounds of whole genome duplication, bloc chromosome duplications, large-scale individual gene duplications, or some combination of these—is unclear. Gene families with a single invertebrate member but four vertebrate members, such as the Hox clusters, provided early support for Ohno's hypothesis that two rounds of genome duplication (the 2R-model) occurred in the stem lineage of extant vertebrates. However, despite extensive study, the duplication history of the Hox clusters has remained unclear, calling into question its usefulness in resolving the role of large-scale gene or genome duplications in early vertebrates. Here, we present a phylogenetic analysis of the vertebrate Hox clusters and several linked genes (the Hox “paralogon”) and show that different phylogenies are obtained for Dlx and Col genes than for Hox and ErbB genes. We show that these results are robust to errors in phylogenetic inference and suggest that these competing phylogenies can be resolved if two chromosomal crossover events occurred in the ancestral vertebrate. These results resolve conflicting data on the order of Hox gene duplications and the role of genome duplication in vertebrate evolution and suggest that a period of genome reorganization occurred after genome duplications in early vertebrates.     

Problem-Elephant Translocation: Translocating the Problem and the Elephant?

Human-elephant conflict (HEC) threatens the survival of endangered Asian elephants (Elephas maximus). Translocating “problem-elephants” is an important HEC mitigation and elephant conservation strategy across elephant range, with hundreds translocated annually. In the first comprehensive assessment of elephant translocation, we monitored 16 translocations in Sri Lanka with GPS collars. All translocated elephants were released into national parks. Two were killed within the parks where they were released, while all the others left those parks. Translocated elephants showed variable responses: “homers” returned to the capture site, “wanderers” ranged widely, and “settlers” established home ranges in new areas soon after release. Translocation caused wider propagation and intensification of HEC, and increased elephant mortality. We conclude that translocation defeats both HEC mitigation and elephant conservation goals.     

Comparative analysis of vertebrate EIF2AK2 (PKR) genes and assignment of the equine gene to ECA15q24–q25 and the bovine gene to BTA11q12–q15

The structures of the canine, rabbit, bovine and equine EIF2AK2 genes were determined. Each of these genes has a 5'' non-coding exon as well as 15 coding exons. All of the canine, bovine and equine EIF2AK2 introns have consensus donor and acceptor splice sites. In the equine EIF2AK2 gene, a unique single nucleotide polymorphism that encoded a Tyr329Cys substitution was detected. Regulatory elements predicted in the promoter region were conserved in ungulates, primates, rodents, Afrotheria (elephant) and Insectifora (shrew). Western clawed frog and fugu EIF2AK2 gene sequences were detected in the USCS Genome Browser and compared to those of other vertebrate EIF2AK2 genes. A comparison of EIF2AK2 protein domains in vertebrates indicates that the kinase catalytic domains were evolutionarily more conserved than the nucleic acid-binding motifs. Nucleotide substitution rates were uniform among the vertebrate sequences with the exception of the zebrafish and goldfish EIF2AK2 genes, which showed substitution rates about 20% higher than those of other vertebrates. FISH was used to physically assign the horse and cattle genes to chromosome locations, ECA15q24–q25 and BTA11q12–15, respectively. Comparative mapping data confirmed conservation of synteny between ungulates, humans and rodents.     

Survey Sequencing and Comparative Analysis of the Elephant Shark (Callorhinchus milii) Genome

Owing to their phylogenetic position, cartilaginous fishes (sharks, rays, skates, and chimaeras) provide a critical reference for our understanding of vertebrate genome evolution. The relatively small genome of the elephant shark, Callorhinchus milii, a chimaera, makes it an attractive model cartilaginous fish genome for whole-genome sequencing and comparative analysis. Here, the authors describe survey sequencing (1.4× coverage) and comparative analysis of the elephant shark genome, one of the first cartilaginous fish genomes to be sequenced to this depth. Repetitive sequences, represented mainly by a novel family of short interspersed element–like and long interspersed element–like sequences, account for about 28% of the elephant shark genome. Fragments of approximately 15,000 elephant shark genes reveal specific examples of genes that have been lost differentially during the evolution of tetrapod and teleost fish lineages. Interestingly, the degree of conserved synteny and conserved sequences between the human and elephant shark genomes are higher than that between human and teleost fish genomes. Elephant shark contains putative four Hox clusters indicating that, unlike teleost fish genomes, the elephant shark genome has not experienced an additional whole-genome duplication. These findings underscore the importance of the elephant shark as a critical reference vertebrate genome for comparative analysis of the human and other vertebrate genomes. This study also demonstrates that a survey-sequencing approach can be applied productively for comparative analysis of distantly related vertebrate genomes.     

Functional conservation of a forebrain enhancer from the elephant shark (<Emphasis Type="Italic">Callorhinchus milii</Emphasis>) in zebrafish and mice

Background  

The phylogenetic position of the elephant shark (Callorhinchus milii ) is particularly relevant to study the evolution of genes and gene regulation in vertebrates. Here we examine the evolution of Dlx homeobox gene regulation during vertebrate embryonic development with a particular focus on the forebrain. We first identified the elephant shark sequence orthologous to the URE2 cis -regulatory element of the mouse Dlx1/Dlx2 locus (herein named CmURE2). We then conducted a comparative study of the sequence and enhancer activity of CmURE2 with that of orthologous regulatory sequences from zebrafish and mouse.     

Defining synphenotype groups in Xenopus tropicalis by use of antisense morpholino oligonucleotides

To identify novel genes involved in early development, and as proof-of-principle of a large-scale reverse genetics approach in a vertebrate embryo, we have carried out an antisense morpholino oligonucleotide (MO) screen in Xenopus tropicalis, in the course of which we have targeted 202 genes expressed during gastrula stages. MOs were designed to complement sequence between −80 and +25 bases of the initiating AUG codons of the target mRNAs, and the specificities of many were tested by (i) designing different non-overlapping MOs directed against the same mRNA, (ii) injecting MOs differing in five bases, and (iii) performing “rescue” experiments. About 65% of the MOs caused X. tropicalis embryos to develop abnormally (59% of those targeted against novel genes), and we have divided the genes into “synphenotype groups,” members of which cause similar loss-of-function phenotypes and that may function in the same developmental pathways. Analysis of the expression patterns of the 202 genes indicates that members of a synphenotype group are not necessarily members of the same synexpression group. This screen provides new insights into early vertebrate development and paves the way for a more comprehensive MO-based analysis of gene function in X. tropicalis.     

Identification and Comparative Analysis of the Protocadherin Cluster in a Reptile,the Green Anole Lizard

Background

The vertebrate protocadherins are a subfamily of cell adhesion molecules that are predominantly expressed in the nervous system and are believed to play an important role in establishing the complex neural network during animal development. Genes encoding these molecules are organized into a cluster in the genome. Comparative analysis of the protocadherin subcluster organization and gene arrangements in different vertebrates has provided interesting insights into the history of vertebrate genome evolution. Among tetrapods, protocadherin clusters have been fully characterized only in mammals. In this study, we report the identification and comparative analysis of the protocadherin cluster in a reptile, the green anole lizard (Anolis carolinensis).

Methodology/Principal Findings

We show that the anole protocadherin cluster spans over a megabase and encodes a total of 71 genes. The number of genes in the anole protocadherin cluster is significantly higher than that in the coelacanth (49 genes) and mammalian (54–59 genes) clusters. The anole protocadherin genes are organized into four subclusters: the δ, α, β and γ. This subcluster organization is identical to that of the coelacanth protocadherin cluster, but differs from the mammalian clusters which lack the δ subcluster. The gene number expansion in the anole protocadherin cluster is largely due to the extensive gene duplication in the γb subgroup. Similar to coelacanth and elephant shark protocadherin genes, the anole protocadherin genes have experienced a low frequency of gene conversion.

Conclusions/Significance

Our results suggest that similar to the protocadherin clusters in other vertebrates, the evolution of anole protocadherin cluster is driven mainly by lineage-specific gene duplications and degeneration. Our analysis also shows that loss of the protocadherin δ subcluster in the mammalian lineage occurred after the divergence of mammals and reptiles. We present a model for the evolutionary history of the protocadherin cluster in tetrapods.     

“Candidatus Accumulibacter” Population Structure in Enhanced Biological Phosphorus Removal Sludges as Revealed by Polyphosphate Kinase Genes

We investigated the fine-scale population structure of the “Candidatus Accumulibacter” lineage in enhanced biological phosphorus removal (EBPR) systems using the polyphosphate kinase 1 gene (ppk1) as a genetic marker. We retrieved fragments of “Candidatus Accumulibacter” 16S rRNA and ppk1 genes from one laboratory-scale and several full-scale EBPR systems. Phylogenies reconstructed using 16S rRNA genes and ppk1 were largely congruent, with ppk1 granting higher phylogenetic resolution and clearer tree topology and thus serving as a better genetic marker than 16S rRNA for revealing population structure within the “Candidatus Accumulibacter” lineage. Sequences from at least five clades of “Candidatus Accumulibacter” were recovered by ppk1-targeted PCR, and subsequently, specific primer sets were designed to target the ppk1 gene for each clade. Quantitative real-time PCR (qPCR) assays using “Candidatus Accumulibacter”-specific 16S rRNA and “Candidatus Accumulibacter” clade-specific ppk1 primers were developed and conducted on three laboratory-scale and nine full-scale EBPR samples and two full-scale non-EBPR samples to determine the abundance of the total “Candidatus Accumulibacter” lineage and the relative distributions and abundances of the five “Candidatus Accumulibacter” clades. The qPCR-based estimation of the total “Candidatus Accumulibacter” fraction as a proportion of the bacterial community as measured using 16S rRNA genes was not significantly different from the estimation measured using ppk1, demonstrating the power of ppk1 as a genetic marker for detection of all currently defined “Candidatus Accumulibacter” clades. The relative distributions of “Candidatus Accumulibacter” clades varied among different EBPR systems and also temporally within a system. Our results suggest that the “Candidatus Accumulibacter” lineage is more diverse than previously realized and that different clades within the lineage are ecologically distinct.     

Gene Splicing of an Invertebrate Beta Subunit (LCavβ) in the N-Terminal and HOOK Domains and Its Regulation of LCav1 and LCav2 Calcium Channels

The accessory beta subunit (Ca_vβ) of calcium channels first appear in the same genome as Ca_v1 L-type calcium channels in single-celled coanoflagellates. The complexity of this relationship expanded in vertebrates to include four different possible Ca_vβ subunits (β₁, β₂, β₃, β₄) which associate with four Ca_v1 channel isoforms (Ca_v1.1 to Ca_v1.4) and three Ca_v2 channel isoforms (Ca_v2.1 to Ca_v2.3). Here we assess the fundamentally-shared features of the Ca_vβ subunit in an invertebrate model (pond snail Lymnaea stagnalis) that bears only three homologous genes: (LCa_v1, LCa_v2, and LCa_vβ). Invertebrate Ca_vβ subunits (in flatworms, snails, squid and honeybees) slow the inactivation kinetics of Ca_v2 channels, and they do so with variable N-termini and lacking the canonical palmitoylation residues of the vertebrate β2a subunit. Alternative splicing of exon 7 of the HOOK domain is a primary determinant of a slow inactivation kinetics imparted by the invertebrate LCa_vβ subunit. LCa_vβ will also slow the inactivation kinetics of LCa_v3 T-type channels, but this is likely not physiologically relevant in vivo. Variable N-termini have little influence on the voltage-dependent inactivation kinetics of differing invertebrate Ca_vβ subunits, but the expression pattern of N-terminal splice isoforms appears to be highly tissue specific. Molluscan LCa_vβ subunits have an N-terminal “A” isoform (coded by exons: 1a and 1b) that structurally resembles the muscle specific variant of vertebrate β1a subunit, and has a broad mRNA expression profile in brain, heart, muscle and glands. A more variable “B” N-terminus (exon 2) in the exon position of mammalian β3 and has a more brain-centric mRNA expression pattern. Lastly, we suggest that the facilitation of closed-state inactivation (e.g. observed in Ca_v2.2 and Ca_vβ₃ subunit combinations) is a specialization in vertebrates, because neither snail subunit (LCa_v2 nor LCa_vβ) appears to be compatible with this observed property.     

Rhomboids of Mycobacteria: Characterization Using an aarA Mutant of Providencia stuartii and Gene Deletion in Mycobacterium smegmatis

Background

Rhomboids are ubiquitous proteins with unknown roles in mycobacteria. However, bioinformatics suggested putative roles in DNA replication pathways and metabolite transport. Here, mycobacterial rhomboid-encoding genes were characterized; first, using the Providencia stuartii null-rhomboid mutant and then deleted from Mycobacterium smegmatis for additional insight in mycobacteria.

Methodology/Principal Findings

Using in silico analysis we identified in M. tuberculosis genome the genes encoding two putative rhomboid proteins; Rv0110 (referred to as “rhomboid protease 1”) and Rv1337 (“rhomboid protease 2”). Genes encoding orthologs of these proteins are widely represented in all mycobacterial species. When transformed into P. stuartii null-rhomboid mutant (ΔaarA), genes encoding mycobacterial orthologs of “rhomboid protease 2” fully restored AarA activity (AarA is the rhomboid protein of P. stuartii). However, most genes encoding mycobacterial “rhomboid protease 1” orthologs did not. Furthermore, upon gene deletion in M. smegmatis, the ΔMSMEG_4904 single mutant (which lost the gene encoding MSMEG_4904, orthologous to Rv1337, “rhomboid protease 2”) formed the least biofilms and was also more susceptible to ciprofloxacin and novobiocin, antimicrobials that inhibit DNA gyrase. However, the ΔMSMEG_5036 single mutant (which lost the gene encoding MSMEG_5036, orthologous to Rv0110, “rhomboid protease 1”) was not as susceptible. Surprisingly, the double rhomboid mutant ΔMSMEG_4904–ΔMSMEG_5036 (which lost genes encoding both homologs) was also not as susceptible suggesting compensatory effects following deletion of both rhomboid-encoding genes. Indeed, transforming the double mutant with a plasmid encoding MSMEG_5036 produced phenotypes of the ΔMSMEG_4904 single mutant (i.e. susceptibility to ciprofloxacin and novobiocin).

Conclusions/Significance

Mycobacterial rhomboid-encoding genes exhibit differences in complementing aarA whereby it''s only genes encoding “rhomboid protease 2” orthologs that fully restore AarA activity. Additionally, gene deletion data suggests inhibition of DNA gyrase by MSMEG_4904; however, the ameliorated effect in the double mutant suggests occurrence of compensatory mechanisms following deletion of genes encoding both rhomboids.     

Importance of phosphoinositide-dependent signaling pathways in the control of gene expression in resting cells and in response to phytohormones

“Phosphoinositide” refers to phosphorylated forms of phosphatidylinositol, including phosphatidylinositol-4-phosphate and phosphatidylinositol-4,5-bisphosphate. Both of these molecules could be in vivo substrates of plant phospholipase C. These phosphoinositides can also be biologically active “per se,” by directly binding to proteins and thus altering their location and/or activity. The use of pharmacological agents in Arabidopsis suspension cells allowed us to identify genes whose expression was positively or negatively controlled, in the basal state, by products of phosphoinositide-dependent phospholipase C. In this basal state, it seems that no genes exhibit a phosphoinositide-dependent expression “per se.” However, many genes whose expression is altered in the presence of phospholipase C inhibitors appeared to be responsive to salicylic acid. This allowed us to show that salicylic acid acts both by increasing the phosphoinositide pool and by inhibiting the phospholipase C. In response to salicylic acid it is possible to identify genes whose expression is controlled by products of PI-PLC, but also genes whose expression is controlled by phosphoinositides “per se.” Our data highlight the importance of phosphoinositide-dependent pathways in gene expression in resting cells and in response to phytohormones.     

Infection Dynamics of the Tick-Borne Pathogen “Candidatus Neoehrlichia mikurensis” and Coinfections with Borrelia afzelii in Bank Voles in Southern Sweden

The tick-borne bacterium “Candidatus Neoehrlichia mikurensis” has recently been recognized as a human pathogen. Together with Borrelia afzelii, it is one of the most common pathogens found in the tick Ixodes ricinus. Here, we compared the epidemiologies of “Ca. Neoehrlichia mikurensis” and B. afzelii by longitudinal sampling from May to September in one of their most abundant vertebrate hosts, the bank vole (Myodes glareolus), using real-time PCR for detection and quantification. The prevalences of “Ca. Neoehrlichia mikurensis” and B. afzelii were determined to be 19% (50/261) and 22% (56/261), respectively. The prevalence of “Ca. Neoehrlichia mikurensis” increased significantly during the sampling season. The clearance rate of “Ca. Neoehrlichia mikurensis” was significantly higher than that of B. afzelii. We found a high frequency of double infections; 46% of all samples infected with “Ca. Neoehrlichia mikurensis” also had a coinfection with B. afzelii. The frequency of coinfections was significantly higher than expected from the prevalence of each pathogen. The high level of coinfections can be caused by interactions between the pathogens or might reflect variation in general susceptibility among voles.     

The Peculiarities of Large Intron Splicing in Animals

In mammals a considerable 92% of genes contain introns, with hundreds and hundreds of these introns reaching the incredible size of over 50,000 nucleotides. These “large introns” must be spliced out of the pre-mRNA in a timely fashion, which involves bringing together distant 5′ and 3′ acceptor and donor splice sites. In invertebrates, especially Drosophila, it has been shown that larger introns can be spliced efficiently through a process known as recursive splicing—a consecutive splicing from the 5′-end at a series of combined donor-acceptor splice sites called RP-sites. Using a computational analysis of the genomic sequences, we show that vertebrates lack the proper enrichment of RP-sites in their large introns, and, therefore, require some other method to aid splicing. We analyzed over 15,000 non-redundant, large introns from six mammals, 1,600 from chicken and zebrafish, and 560 non-redundant large introns from five invertebrates. Our bioinformatic investigation demonstrates that, unlike the studied invertebrates, the studied vertebrate genomes contain consistently abundant amounts of direct and complementary strand interspersed repetitive elements (mainly SINEs and LINEs) that may form stems with each other in large introns. This examination showed that predicted stems are indeed abundant and stable in the large introns of mammals. We hypothesize that such stems with long loops within large introns allow intron splice sites to find each other more quickly by folding the intronic RNA upon itself at smaller intervals and, thus, reducing the distance between donor and acceptor sites.     

α-Proteobacterial Symbionts of Marine Bryozoans in the Genus Watersipora

Bacterial symbionts that resembled mollicutes were discovered in the marine bryozoan Watersipora arcuata in the 1980s. In this study, we used PCR and sequencing of 16S rRNA genes, specific fluorescence in situ hybridization, and phylogenetic analysis to determine that the bacterial symbionts of “W. subtorquata” and “W. arcuata” from several locations along the California coast are actually closely related α-Proteobacteria, not mollicutes. We propose the names “Candidatus Endowatersipora palomitas” and “Candidatus Endowatersipora rubus” for the symbionts of “W. subtorquata” and “W. arcuata,” respectively.