BAC libraries and comparative genomics of aquatic chordate species

Reconstruction of the repetitive antifreeze glycoprotein genomic loci in the cold-water gadids Boreogadus saida and Microgadus tomcod

Antifreeze glycoproteins (AFGPs) are a novel evolutionary innovation in members of the northern cod fish family (Gadidae), crucial in preventing death from inoculative freezing by environmental ice in their frigid Arctic and sub-Arctic habitats. However, the genomic origin and molecular mechanism of evolution of this novel life-saving adaptive genetic trait remained to be definitively determined. To this end, we constructed large insert genomic DNA BAC (bacterial artificial chromosome) libraries for two AFGP-bearing gadids, the high-Arctic polar cod Boreogadus saida and the cold-temperate Atlantic tomcod Microgadus tomcod, to isolate and sequence their AFGP genomic regions for fine resolution evolutionary analyses. The BAC library construction encountered poor cloning efficiency initially, which we resolved by pretreating the agarose-embedded erythrocyte DNA with a cationic detergent, a method that may be of general use to BAC cloning for teleost species and/or where erythrocytes are the source of input DNA. The polar cod BAC library encompassed 92,160 clones with an average insert size of 94.7 kbp, and the Atlantic tomcod library contained 73,728 clones with an average insert size of 89.6 kbp. The genome sizes of B. saida and M. tomcod were estimated by cell flow cytometry to be 836 Mbp and 645 Mbp respectively, thus their BAC libraries have approximately 10- and 9.7-fold genome coverage respectively. The inclusiveness and depth of coverage were empirically confirmed by screening the libraries with three housekeeping genes. The BAC clones that mapped to the AFGP genomic loci of the two gadids were then isolated by screening the BAC libraries with gadid AFGP gene probes. Eight minimal tiling path (MTP) clones were identified for B. saida, sequenced, and assembled. The B. saida AFGP locus reconstruction produced both haplotypes, and the locus comprises three distinct AFGP gene clusters, containing a total of 16 AFGP genes and spanning a combined distance of 512 kbp. The M. tomcod AFGP locus is much smaller at approximately 80 kbp, and contains only three AFGP genes. Fluorescent in situ hybridization with an AFGP gene probe showed the AFGP locus in both species occupies a single chromosomal location. The large AFGP locus with its high gene dosage in B. saida is consistent with its chronically freezing high Arctic habitats, while the small gene family in M. tomcod correlates with its milder habitats in lower latitudes. The results from this study provided the data for fine resolution sequence analyses that would yield insight into the molecular mechanisms and history of gadid AFGP gene evolution driven by northern hemisphere glaciation.

A functional Bucky ball-GFP transgene visualizes germ plasm in living zebrafish

In many animals, the germline is specified by maternal RNA-granules termed germ plasm. The correct localization of germ plasm during embryogenesis is therefore crucial for the specification of germ cells. In zebrafish, we previously identified Bucky ball (Buc) as a key regulator of germ plasm formation. Here, we used a Buc antibody to describe its continuous germ plasm localization. Moreover, we generated a transgenic Buc-GFP line for live imaging, which visualizes germ plasm from its assembly during oogenesis up to the larval stages. Live imaging of Buc-GFP generated stunning movies, as they highlighted the dynamic details of germ plasm movements. Moreover, we discovered that Buc was still detected in primordial germ cells 2 days after fertilization. Interestingly, the transgene rescued buc mutants demonstrating genetically that the Buc-GFP fusion protein is functional. These results show that Buc-GFP exerts all biochemical interactions essential for germline development and highlight the potential of this line to analyze the molecular regulation of germ plasm formation.

Protein genes in repetitive sequence-antifreeze glycoproteins in Atlantic cod genome

Background

Highly repetitive sequences are the bane of genome sequence assembly, and the short read lengths produced by current next generation sequencing technologies further exacerbates this obstacle. An adopted practice is to exclude repetitive sequences in genome data assembly, as the majority of repeats lack protein-coding genes. However, this could result in the exclusion of important genotypes in newly sequenced non-model species. The absence of the antifreeze glycoproteins (AFGP) gene family in the recently sequenced Atlantic cod genome serves as an example.

Results

The Atlantic cod (Gadus morhua) genome was assembled entirely from Roche 454 short reads, demonstrating the feasibility of this approach. However, a well-known major adaptive trait, the AFGP, essential for survival in frigid Arctic marine habitats was absent in the annotated genome. To assess whether this resulted from population difference, we performed Southern blot analysis of genomic DNA from multiple individuals from the North East Arctic cod population that the sequenced cod belonged, and verified that the AFGP genotype is indeed present. We searched the raw assemblies of the Atlantic cod using our G. morhua AFGP gene, and located partial AFGP coding sequences in two sequence scaffolds. We found these two scaffolds constitute a partial genomic AFGP locus through comparative sequence analyses with our newly assembled genomic AFGP locus of the related polar cod, Boreogadus saida. By examining the sequence assembly and annotation methodologies used for the Atlantic cod genome, we deduced the primary cause of the absence of the AFGP gene family from the annotated genome was the removal of all repetitive Roche 454 short reads before sequence assembly, which would exclude most of the highly repetitive AFGP coding sequences. Secondarily, the model teleost genomes used in projection annotation of the Atlantic cod genome have no antifreeze trait, perpetuating the unawareness that the AFGP gene family is missing.

Conclusions

We recovered some of the missing AFGP coding sequences and reconstructed a partial AFGP locus in the Atlantic cod genome, bringing to light that not all repetitive sequences lack protein coding information. Also, reliance on genomes of model organisms as reference for annotating protein-coding gene content of a newly sequenced non-model species could lead to omission of novel genetic traits.

Isolation of a 97-kb minimal essential MHC B locus from a new reverse-4D BAC library of the golden pheasant

The bacterial artificial chromosome (BAC) system is widely used in isolation of large genomic fragments of interest. Construction of a routine BAC library requires several months for picking clones and arraying BACs into superpools in order to employ 4D-PCR to screen positive BACs, which might be time-consuming and laborious. The major histocompatibility complex (MHC) is a cluster of genes involved in the vertebrate immune system, and the classical avian MHC-B locus is a minimal essential one, occupying a 100-kb genomic region. In this study, we constructed a more effective reverse-4D BAC library for the golden pheasant, which first creates sub-libraries and then only picks clones of positive sub-libraries, and identified several MHC clones within thirty days. The full sequencing of a 97-kb reverse-4D BAC demonstrated that the golden pheasant MHC-B locus contained 20 genes and showed good synteny with that of the chicken. The notable differences between these two species were the numbers of class II B loci and NK genes and the inversions of the TAPBP gene and the TAP1-TAP2 region. Furthermore, the inverse TAP2-TAP1 was unique in the golden pheasant in comparison with that of chicken, turkey, and quail. The newly defined genomic structure of the golden pheasant MHC will give an insight into the evolutionary history of the avian MHC.

Sequencing three crocodilian genomes to illuminate the evolution of archosaurs and amniotes

The International Crocodilian Genomes Working Group (ICGWG) will sequence and assemble the American alligator (Alligator mississippiensis), saltwater crocodile (Crocodylus porosus) and Indian gharial (Gavialis gangeticus) genomes. The status of these projects and our planned analyses are described.

Antarctic notothenioid fishes: genomic resources and strategies for analyzing an adaptive radiation

The perciform suborder Notothenoidei provides a compelling opportunity to study the adaptive radiation of a marine species-flock in the cold Southern Ocean that surrounds Antarctica. To facilitate genome-level studies of the diversification of these fishes, we present estimates of the genome sizes of 11 Antarctic species and describe the production of high-quality bacterial artificial chromosome (BAC) libraries for two, the red-blooded notothen Notothenia coriiceps and the white-blooded icefish Chaenocephalus aceratus. Our results indicate that evolution of phylogenetically derived notothenioid families (e.g., the crown group Channichthyidae [icefishes]), was accompanied by genome expansion. Six species from the basal family Nototheniidae had C-values between 0.98 and 1.20 pg, a range that is consistent with the genome sizes of proposed outgroups (e.g., percids) of the notothenioid suborder. In contrast, four icefishes had C-values in the range 1.66-1.83 pg. The BAC libraries VMRC-19 (N. coriiceps) and VMRC-21 (C. aceratus) comprise 12× and 10× coverage of the respective genomes and have average insert sizes of 138 and 168 kb. Paired BAC-end reads representing ∼0.1% of each genome showed that the repetitive element landscapes of the two genomes (13.4% of the N. coriiceps genome and 14.5% for C. aceratus) were similar. The availability of these high-quality and well-characterized BAC libraries sets the stage for targeted genomic analyses of the unusual anatomical and physiological adaptations of the notothenioids, some of which mimic human diseases. Here we consider the evolution of secondary pelagicism by various taxa of the group and illustrate the utility of Antarctic icefishes as an evolutionary-mutant model of human osteopenia (low-mineral density of bones).

Development and analysis of a germline BAC resource for the sea lamprey, a vertebrate that undergoes substantial chromatin diminution

Over the last several years, the sea lamprey (Petromyzon marinus) has grown substantially as a model for understanding the evolutionary fundaments and capacity of vertebrate developmental and genome biology. Recent work on the lamprey genome has resulted in a preliminary assembly of the lamprey genome and led to the realization that nearly all somatic cell lineages undergo extensive programmed rearrangements. Here we describe the development of a bacterial artificial chromosome (BAC) resource for lamprey germline DNA and use sequence information from this resource to probe the subchromosomal structure of the lamprey genome. The arrayed germline BAC library represents approximately 10x coverage of the lamprey genome. Analyses of BAC-end sequences reveal that the lamprey genome possesses a high content of repetitive sequences (relative to human), which show strong clustering at the subchromosomal level. This pattern is not unexpected given that the sea lamprey genome is dispersed across a large number of chromosomes (n approximately 99) and suggests a low-copy DNA targeting strategy for efficiently generating informative paired-BAC-end linkages from highly repetitive genomes. This library therefore represents a new and biologically informed resource for understanding the structure of the lamprey genome and the biology of programmed genome rearrangement.

A bacterial artificial chromosome library for the Australian saltwater crocodile (Crocodylus porosus) and its utilization in gene isolation and genome characterization

Background

Crocodilians (Order Crocodylia) are an ancient vertebrate group of tremendous ecological, social, and evolutionary importance. They are the only extant reptilian members of Archosauria, a monophyletic group that also includes birds, dinosaurs, and pterosaurs. Consequently, crocodilian genomes represent a gateway through which the molecular evolution of avian lineages can be explored. To facilitate comparative genomics within Crocodylia and between crocodilians and other archosaurs, we have constructed a bacterial artificial chromosome (BAC) library for the Australian saltwater crocodile, Crocodylus porosus. This is the first BAC library for a crocodile and only the second BAC resource for a crocodilian.

Results

The C. porosus BAC library consists of 101,760 individually archived clones stored in 384-well microtiter plates. NotI digestion of random clones indicates an average insert size of 102 kb. Based on a genome size estimate of 2778 Mb, the library affords 3.7 fold (3.7×) coverage of the C. porosus genome. To investigate the utility of the library in studying sequence distribution, probes derived from CR1a and CR1b, two crocodilian CR1-like retrotransposon subfamilies, were hybridized to C. porosus macroarrays. The results indicate that there are a minimum of 20,000 CR1a/b elements in C. porosus and that their distribution throughout the genome is decidedly non-random. To demonstrate the utility of the library in gene isolation, we probed the C. porosus macroarrays with an overgo designed from a C-mos (oocyte maturation factor) partial cDNA. A BAC containing C-mos was identified and the C-mos locus was sequenced. Nucleotide and amino acid sequence alignment of the C. porosus C-mos coding sequence with avian and reptilian C-mos orthologs reveals greater sequence similarity between C. porosus and birds (specifically chicken and zebra finch) than between C. porosus and squamates (green anole).

Conclusion

We have demonstrated the utility of the Crocodylus porosus BAC library as a tool in genomics research. The BAC library should expedite complete genome sequencing of C. porosus and facilitate detailed analysis of genome evolution within Crocodylia and between crocodilians and diverse amniote lineages including birds, mammals, and other non-avian reptiles.

Evolutionary insights into the unique electromotility motor of mammalian outer hair cells

Prestin (SLC26A5) is the molecular motor responsible for cochlear amplification by mammalian cochlea outer hair cells and has the unique combined properties of energy-independent motility, voltage sensitivity, and speed of cellular shape change. The ion transporter capability, typical of SLC26A members, was exchanged for electromotility function and is a newly derived feature of the therian cochlea. A putative minimal essential motif for the electromotility motor (meEM) was identified through the amalgamation of comparative genomic, evolution, and structural diversification approaches. Comparisons were done among nonmammalian vertebrates, eutherian mammalian species, and the opossum and platypus. The opossum and platypus SLC26A5 proteins were comparable to the eutherian consensus sequence. Suggested from the point-accepted mutation analysis, the meEM motif spans all the transmembrane segments and represented residues 66-503. Within the eutherian clade, the meEM was highly conserved with a substitution frequency of only 39/7497 (0.5%) residues, compared with 5.7% in SLC26A4 and 12.8% in SLC26A6 genes. Clade-specific substitutions were not observed and there was no sequence correlation with low or high hearing frequency specialists. We were able to identify that within the highly conserved meEM motif two regions, which are unique to all therian species, appear to be the most derived features in the SLC26A5 peptide.

Comparative genomics in vertebrate evolution and development

The vast quantities of publicly available DNA sequencing data and genome resources are enabling biologists to investigate age-old problems in biology that were not addressable previously. In this review, we discuss how comparative genomics is practiced and how the data can be used to make biological inferences with respect to vertebrate evolution and development. Examples are taken from the well-known HOX clusters, which are always a high-priority target for genomic analyses due to their inferred role in the evolution of metazoans. In addition, we briefly discuss the application of genomic approaches to problems in comparative endocrinology.

Construction of a nurse shark (Ginglymostoma cirratum) bacterial artificial chromosome (BAC) library and a preliminary genome survey

BACKGROUND

Sharks are members of the taxonomic class Chondrichthyes, the oldest living jawed vertebrates. Genomic studies of this group, in comparison to representative species in other vertebrate taxa, will allow us to theorize about the fundamental genetic, developmental, and functional characteristics in the common ancestor of all jawed vertebrates.

AIMS

In order to obtain mapping and sequencing data for comparative genomics, we constructed a bacterial artificial chromosome (BAC) library for the nurse shark, Ginglymostoma cirratum.

RESULTS

The BAC library consists of 313,344 clones with an average insert size of 144 kb, covering ~4.5 x 1010 bp and thus providing an 11-fold coverage of the haploid genome. BAC end sequence analyses revealed, in addition to LINEs and SINEs commonly found in other animal and plant genomes, two new groups of nurse shark-specific repetitive elements, NSRE1 and NSRE2 that seem to be major components of the nurse shark genome. Screening the library with single-copy or multi-copy gene probes showed 6-28 primary positive clones per probe of which 50-90% were true positives, demonstrating that the BAC library is representative of the different regions of the nurse shark genome. Furthermore, some BAC clones contained multiple genes, making physical mapping feasible.

CONCLUSION

We have constructed a deep-coverage, high-quality, large insert, and publicly available BAC library for a cartilaginous fish. It will be very useful to the scientific community interested in shark genomic structure, comparative genomics, and functional studies. We found two new groups of repetitive elements specific to the nurse shark genome, which may contribute to the architecture and evolution of the nurse shark genome.

Ancient divergence of a complex family of immune-type receptor genes

Multigene families of activating/inhibitory receptors belonging to the immunoglobulin superfamily (IgSF) regulate immunological and other cell-cell interactions. A new family of such genes, termed modular domain immune-type receptors (MDIRs), has been identified in the clearnose skate (Raja eglanteria), a phylogenetically ancient vertebrate. At least five different major forms of predicted MDIR proteins are comprised of four different subfamilies of IgSF ectodomains of the intermediate (I)- or C2-set. The predicted number of individual IgSF ectodomains in MDIRs varies from one to six. MDIR1 contains a positively charged transmembrane residue and MDIR2 and MDIR3 each possesses at least one immunoreceptor tyrosine-based inhibitory motif in their cytoplasmic regions. MDIR4 and MDIR5 lack characteristic activating/inhibitory signalling motifs. MDIRs are encoded in a particularly large and complex multigene family. MDIR domains exhibit distant sequence similarity to mammalian CMRF-35-like molecules, polymeric immunoglobulin receptors, triggering receptors expressed on myeloid cells (TREMs), TREM-like transcripts, NKp44 and FcR homologs, as well as to sequences identified in several different vertebrate genomes. Phylogenetic analyses suggest that MDIRs are representative members of an extended family of IgSF genes that diverged before or very early in evolution of the vertebrates and subsequently came to occupy multiple, fully independent distributions in the present day.

Utilization of a zebra finch BAC library to determine the structure of an avian androgen receptor genomic region

The zebra finch (Taeniopygia guttata) is an important model organism for studying behavior, neuroscience, avian biology, and evolution. To support the study of its genome, we constructed a BAC library (TG__Ba) using DNA from livers of females. The BAC library consists of 147,456 clones with 98% containing inserts of an average size of 134 kb and represents 15.5 haploid genome equivalents. By sequencing a whole BAC, a full-length androgen receptor open reading frame was identified, the first in an avian species. Comparison of BAC end sequences and the whole BAC sequence with the chicken genome draft sequence showed a high degree of conserved synteny between the zebra finch and the chicken genome.