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

Analysis of shark NCR3 family genes reveals primordial features of vertebrate NKp30

Natural killer (NK) cells play major roles in innate immunity against viruses and cancer. Natural killer receptors (NKR) expressed by NK cells recognize foreign- or self-ligands on infected and transformed cells as well as healthy cells. NKR genes are the most rapidly evolving loci in vertebrates, and it is generally difficult to detect orthologues in different taxa. The unique exception is NKp30, an activating NKR in mammals that binds to the self-ligand B7H6. The NKp30-encoding gene, NCR3, has been found in most vertebrates including sharks, the oldest vertebrates with human-type adaptive immunity. NCR3 has a special, non-rearranging VJ-type immunoglobulin superfamily (IgSF) domain that predates the emergence of the rearranging antigen receptors. Herein we show that NCR3 loci are linked to the shark major histocompatibility complex (MHC), proving NCR3's primordial association with the MHC. We identified eight subtypes of differentially expressed highly divergent shark NCR3 family genes. Using in situ hybridization, we detected one subtype, NS344823, to be expressed by predominantly single cells outside of splenic B cell zones. The expression by non-B cells was also confirmed by PCR in peripheral blood lymphocytes. Surprisingly, high expression of NS344823 was detected in the thymic cortex, demonstrating NS344823 expression in developing T cells. Finally, we show for the first time that shark T cells are found as single cells or in small clusters in the splenic red pulp, also unassociated with the large B cell follicles we previously identified.

Ancient Use of Ig Variable Domains Contributes Significantly to the TCRδ Repertoire

The loci encoding B and T cell Ag receptors are generally distinct in commonly studied mammals, with each receptor's gene segments limited to intralocus, cis chromosomal rearrangements. The nurse shark (Ginglymostoma cirratum) represents the oldest vertebrate class, the cartilaginous fish, with adaptive immunity provided via Ig and TCR lineages, and is one species among a growing number of taxa employing Ig-TCRδ rearrangements that blend these distinct lineages. Analysis of the nurse shark Ig-TCRδ repertoire found that these rearrangements possess CDR3 characteristics highly similar to canonical TCRδ rearrangements. Furthermore, the Ig-TCRδ rearrangements are expressed with TCRγ, canonically found in the TCRδ heterodimer. We also quantified BCR and TCR transcripts in the thymus for BCR (IgHV-IgHC), chimeric (IgHV-TCRδC), and canonical (TCRδV-TCRδC) transcripts, finding equivalent expression levels in both thymus and spleen. We also characterized the nurse shark TCRαδ locus with a targeted bacterial artifical chromosome sequencing approach and found that the TCRδ locus houses a complex of V segments from multiple lineages. An IgH-like V segment, nestled within the nurse shark TCRδ translocus, grouped with IgHV-like rearrangements we found expressed with TCRδ (but not IgH) rearrangements in our phylogenetic analysis. This distinct lineage of TCRδ-associated IgH-like V segments was termed "TAILVs." Our data illustrate a dynamic TCRδ repertoire employing TCRδVs, NARTCRVs, bona fide trans-rearrangements from shark IgH clusters, and a novel lineage in the TCRδ-associated Ig-like V segments.

Molecular Origin of the Stability Difference in Four Shark IgNAR Constant Domains

Improving the stability of antibodies for manufacture and shelf life is one of the main focuses of antibody engineering. One stabilization strategy is to perform specific mutations in human antibodies based on highly stable antibodies in other species. To identify the key residues for mutagenesis, it is necessary to understand the roles of these residues in stabilizing the antibody. Here, we use molecular dynamics simulations to study the molecular origin of the four shark immunoglobulin new antigen receptors constant domains (C1-C4). According to the unfolding pathways and the conformational free energy surfaces in 8 M urea at 380 K, the C2 domain is the most stable, followed by C4, C1, and C3, which agrees with the experimental findings. The C1 and C3 domains follow a common unfolding pathway in which the unfolding starts from the edge strands, particularly strand g, and then gradually progresses to the inner strands. Detailed structural analysis of the C2 domain reveals a "sandwich-like" R339-E322-R341 salt-bridge cluster on strand g, which grants ultrahigh stability to the C2 domain. We further design two sets of mutations by mutating E322 to alanine or setting all atomic charges in E322 to zero to break the salt-bridge cluster in the C2 domain, which confirms the importance of the salt bridges in stability. In the C4 domain, the D80-K104 salt bridge on strand g also strengthens the stability. On the other hand, in the C1 and C3 domains, there is no salt bridge on strand g. In addition to the salt bridges, the overall hydrophobicity score of the hydrophobic core is also positively correlated with the domain stability. Our findings provide a detailed microscopic picture of the molecular origin of the four shark immunoglobulin new antigen receptors constant domains that not only explains the differences in their structural stability but also provides important insights into future antibody design.

Construction of a Genomic Bacterial Artificial Chromosome (BAC) Library for the Prawn Macrobrachium rosenbergii and Initial Analysis of ZW Chromosome-Derived BAC Inserts

Knowledge on sex determination has proven valuable for commercial production of the prawn Macrobrachium rosenbergii due to sex dimorphism of the male and female individuals. Previous studies indicated that prawn sex is determined by a ZW-ZZ chromosomal system, but no genomic information is available for the sex chromosome. Herein, we constructed a genomic bacterial artificial chromosome (BAC) library and identified the ZW-derived BAC clones for initial analysis of the sex chromosomal DNA sequence. The arrayed BAC library contains 200,448 clones with average insert size of 115.4 kb, corresponding to ∼ 4× coverage of the estimated 5.38 Gb genome. Based on a short female-specific marker, a Z- and a W-fragment were retrieved with the genomic walking method. Screening the BAC library using a ZW-specific marker as probe resulted in 12 positive clones. From these, a Z-derived (P331M17) and a W-derived (P122G2) BAC clones were randomly selected and sequenced by PacBio method. We report the construction of a large insert, deep-coverage, and high-quality BAC library for M. rosenbergii that provides a useful resource for positional cloning of target genes, genomic organization, and comparative genomics analysis. Our study not only confirmed the ZW/ZZ system but also discovered sex-linked genes on ZW chromosomes for the first time, contributing to a comprehensive understanding of the genomic structure of sex chromosomes in M. rosenbergii.

Phylotranscriptomics suggests the jawed vertebrate ancestor could generate diverse helper and regulatory T cell subsets

BACKGROUND

The cartilaginous fishes diverged from other jawed vertebrates ~ 450 million years ago (mya). Despite this key evolutionary position, the only high-quality cartilaginous fish genome available is for the elephant shark (Callorhinchus milii), a chimaera whose ancestors split from the elasmobranch lineage ~ 420 mya. Initial analysis of this resource led to proposals that key components of the cartilaginous fish adaptive immune system, most notably their array of T cell subsets, was primitive compared to mammals. This proposal is at odds with the robust, antigen-specific antibody responses reported in elasmobranchs following immunization. To explore this discrepancy, we generated a multi-tissue transcriptome for small-spotted catshark (Scyliorhinus canicula), a tractable elasmobranch model for functional studies. We searched this, and other newly available sequence datasets, for CD4+ T cell subset-defining genes, aiming to confirm the presence or absence of each subset in cartilaginous fishes.

RESULTS

We generated a new transcriptome based on a normalised, multi-tissue RNA pool, aiming to maximise representation of tissue-specific and lowly expressed genes. We utilized multiple transcriptomic datasets and assembly variants in phylogenetic reconstructions to unambiguously identify several T cell subset-specific molecules in cartilaginous fishes for the first time, including interleukins, interleukin receptors, and key transcription factors. Our results reveal the inability of standard phylogenetic reconstruction approaches to capture the site-specific evolutionary processes of fast-evolving immune genes but show that site-heterogeneous mixture models can adequately do so.

CONCLUSIONS

Our analyses reveal that cartilaginous fishes are capable of producing a range of CD4+ T cell subsets comparable to that of mammals. Further, that the key molecules required for the differentiation and functioning of these subsets existed in the jawed vertebrate ancestor. Additionally, we highlight the importance of considering phylogenetic diversity and, where possible, utilizing multiple datasets for individual species, to accurately infer gene presence or absence at higher taxonomic levels.

Biased Immunoglobulin Light Chain Gene Usage in the Shark

This study of a large family of κ L chain clusters in nurse shark completes the characterization of its classical Ig gene content (two H chain isotypes, μ and ω, and four L chain isotypes, κ, λ, σ, and σ-2). The shark κ clusters are minigenes consisting of a simple VL-JL-CL array, where V to J recombination occurs over an ~500-bp interval, and functional clusters are widely separated by at least 100 kb. Six out of ~39 κ clusters are prerearranged in the germline (germline joined). Unlike the complex gene organization and multistep assembly process of Ig in mammals, each shark Ig rearrangement, somatic or in the germline, appears to be an independent event localized to the minigene. This study examined the expression of functional, nonproductive, and sterile transcripts of the κ clusters compared with the other three L chain isotypes. κ cluster usage was investigated in young sharks, and a skewed pattern of split gene expression was observed, one similar in functional and nonproductive rearrangements. These results show that the individual activation of the spatially distant κ clusters is nonrandom. Although both split and germline-joined κ genes are expressed, the latter are prominent in young animals and wane with age. We speculate that, in the shark, the differential activation of the multiple isotypes can be advantageously used in receptor editing.

SkateBase, an elasmobranch genome project and collection of molecular resources for chondrichthyan fishes

Chondrichthyan fishes are a diverse class of gnathostomes that provide a valuable perspective on fundamental characteristics shared by all jawed and limbed vertebrates. Studies of phylogeny, species diversity, population structure, conservation, and physiology are accelerated by genomic, transcriptomic and protein sequence data. These data are widely available for many sarcopterygii (coelacanth, lungfish and tetrapods) and actinoptergii (ray-finned fish including teleosts) taxa, but limited for chondrichthyan fishes.  In this study, we summarize available data for chondrichthyes and describe resources for one of the largest projects to characterize one of these fish, Leucoraja erinacea, the little skate.  SkateBase ( http://skatebase.org) serves as the skate genome project portal linking data, research tools, and teaching resources.

What the shark immune system can and cannot provide for the expanding design landscape of immunotherapy

INTRODUCTION

Sharks have successfully lived in marine ecosystems, often atop food chains as apex predators, for nearly one and a half billion years. Throughout this period they have benefitted from an immune system with the same fundamental components found in terrestrial vertebrates like man. Additionally, sharks have some rather extraordinary immune mechanisms which mammals lack.

AREAS COVERED

In this review the author briefly orients the reader to sharks, their adaptive immunity, and their important phylogenetic position in comparative immunology. The author also differentiates some of the myths from facts concerning these animals, their cartilage, and cancer. From thereon, the author explores some of the more remarkable capabilities and products of shark lymphocytes. Sharks have an isotype of light chain-less antibodies that are useful tools in molecular biology and are moving towards translational use in the clinic. These special antibodies are just one of the several tricks of shark lymphocyte antigen receptor systems.

EXPERT OPINION

While shark cartilage has not helped oncology patients, shark immunoglobulins and T cell receptors do offer exciting novel possibilities for immunotherapeutics. Much of the clinical immunology developmental pipeline has turned from traditional vaccines to passively delivered monoclonal antibody-based drugs for targeted depletion, activation, blocking and immunomodulation. The immunogenetic tools of shark lymphocytes, battle-tested since the dawn of our adaptive immune system, are well poised to expand the design landscape for the next generation of immunotherapy products.

Construction and analysis of Siberian tiger bacterial artificial chromosome library with approximately 6.5-fold genome equivalent coverage

Bacterial artificial chromosome (BAC) libraries are extremely valuable for the genome-wide genetic dissection of complex organisms. The Siberian tiger, one of the most well-known wild primitive carnivores in China, is an endangered animal. In order to promote research on its genome, a high-redundancy BAC library of the Siberian tiger was constructed and characterized. The library is divided into two sub-libraries prepared from blood cells and two sub-libraries prepared from fibroblasts. This BAC library contains 153,600 individually archived clones; for PCR-based screening of the library, BACs were placed into 40 superpools of 10 × 384-deep well microplates. The average insert size of BAC clones was estimated to be 116.5 kb, representing approximately 6.46 genome equivalents of the haploid genome and affording a 98.86% statistical probability of obtaining at least one clone containing a unique DNA sequence. Screening the library with 19 microsatellite markers and a SRY sequence revealed that each of these markers were present in the library; the average number of positive clones per marker was 6.74 (range 2 to 12), consistent with 6.46 coverage of the tiger genome. Additionally, we identified 72 microsatellite markers that could potentially be used as genetic markers. This BAC library will serve as a valuable resource for physical mapping, comparative genomic study and large-scale genome sequencing in the tiger.

Comparative BAC-based physical mapping of Oryza sativa ssp. indica var. 93-11 and evaluation of the two rice reference sequence assemblies

Reference sequences are sequences that are used for public consultation, and therefore must be of high quality. Using the whole-genome shotgun/next-generation sequencing approach, many genome sequences of complex higher plants have been generated in recent years, and are generally considered reference sequences. However, none of these sequences has been experimentally evaluated at the whole-genome sequence assembly level. Rice has a relatively simple plant genome, and the genome sequences for its two sub-species obtained using different sequencing approaches were published approximately 10 years ago. This provides a unique system for a case study to evaluate the qualities and utilities of published plant genome sequences. We constructed a robust BAC physical map embedding a large number of BAC end sequences forrice variety 93-11. Through BAC end sequence alignments and tri-assembly comparisons of the 93-11 physical map and the two reference sequences, we found that the Nipponbare reference sequence generated using the clone-by-clone approach has a high quality but still contains small artifact inversions and missing sequences. In contrast, the 93-11 reference sequence generated using the whole-genome shotgun approach contains many large and varied assembly errors, such as inversions, duplications and translocations, as well as missing sequences. The 93-11 physical map provides an invaluable resource for evaluation and improvements toward completion of both Nipponbare and 93-11 reference sequences.

Global genomic diversity of Oryza sativa varieties revealed by comparative physical mapping

Bacterial artificial chromosome (BAC) physical maps embedding a large number of BAC end sequences (BESs) were generated for Oryza sativa ssp. indica varieties Minghui 63 (MH63) and Zhenshan 97 (ZS97) and were compared with the genome sequences of O. sativa spp. japonica cv. Nipponbare and O. sativa ssp. indica cv. 93-11. The comparisons exhibited substantial diversities in terms of large structural variations and small substitutions and indels. Genome-wide BAC-sized and contig-sized structural variations were detected, and the shared variations were analyzed. In the expansion regions of the Nipponbare reference sequence, in comparison to the MH63 and ZS97 physical maps, as well as to the previously constructed 93-11 physical map, the amounts and types of the repeat contents, and the outputs of gene ontology analysis, were significantly different from those of the whole genome. Using the physical maps of four wild Oryza species from OMAP (http://www.omap.org) as a control, we detected many conserved and divergent regions related to the evolution process of O. sativa. Between the BESs of MH63 and ZS97 and the two reference sequences, a total of 1532 polymorphic simple sequence repeats (SSRs), 71,383 SNPs, 1767 multiple nucleotide polymorphisms, 6340 insertions, and 9137 deletions were identified. This study provides independent whole-genome resources for intra- and intersubspecies comparisons and functional genomics studies in O. sativa. Both the comparative physical maps and the GBrowse, which integrated the QTL and molecular markers from GRAMENE (http://www.gramene.org) with our physical maps and analysis results, are open to the public through our Web site (http://gresource.hzau.edu.cn/resource/resource.html).

A BAC based physical map and genome survey of the rice false smut fungus Villosiclava virens

Background

Rice false smut caused by Villosiclava virens is a devastating fungal disease that spreads in major rice-growing regions throughout the world. However, the genomic information for this fungal pathogen is limited and the pathogenic mechanism of this disease is still not clear. To facilitate genetic, molecular and genomic studies of this fungal pathogen, we constructed the first BAC-based physical map and performed the first genome survey for this species.

Results

High molecular weight genomic DNA was isolated from young mycelia of the Villosiclava virens strain UV-8b and a high-quality, large-insert and deep-coverage Bacterial Artificial Chromosome (BAC) library was constructed with the restriction enzyme HindIII. The BAC library consisted of 5,760 clones, which covers 22.7-fold of the UV-8b genome, with an average insert size of 140 kb and an empty clone rate of lower than 1%. BAC fingerprinting generated successful fingerprints for 2,290 BAC clones. Using the fingerprints, a whole genome-wide BAC physical map was constructed that contained 194 contigs (2,035 clones) spanning 51.2 Mb in physical length. Bidirectional-end sequencing of 4,512 BAC clones generated 6,560 high quality BAC end sequences (BESs), with a total length of 3,030,658 bp, representing 8.54% of the genome sequence. Analysis of the BESs revealed general genome information, including 51.52% GC content, 22.51% repetitive sequences, 376.12/Mb simple sequence repeat (SSR) density and approximately 36.01% coding regions. Sequence comparisons to other available fungal genome sequences through BESs showed high similarities to Metarhizium anisopliae, Trichoderma reesei, Nectria haematococca and Cordyceps militaris, which were generally in agreement with the 18S rRNA gene analysis results.

Conclusion

This study provides the first BAC-based physical map and genome information for the important rice fungal pathogen Villosiclava virens. The BAC clones, physical map and genome information will serve as fundamental resources to accelerate the genetic, molecular and genomic studies of this pathogen, including positional cloning, comparative genomic analysis and whole genome sequencing. The BAC library and physical map have been opened to researchers as public genomic resources (http://gresource.hzau.edu.cn/resource/resource.html).

Genomic resources for gene discovery, functional genome annotation, and evolutionary studies of maize and its close relatives

Maize is one of the most important food crops and a key model for genetics and developmental biology. A genetically anchored and high-quality draft genome sequence of maize inbred B73 has been obtained to serve as a reference sequence. To facilitate evolutionary studies in maize and its close relatives, much like the Oryza Map Alignment Project (OMAP) (www.OMAP.org) bacterial artificial chromosome (BAC) resource did for the rice community, we constructed BAC libraries for maize inbred lines Zheng58, Chang7-2, and Mo17 and maize wild relatives Zea mays ssp. parviglumis and Tripsacum dactyloides. Furthermore, to extend functional genomic studies to maize and sorghum, we also constructed binary BAC (BIBAC) libraries for the maize inbred B73 and the sorghum landrace Nengsi-1. The BAC/BIBAC vectors facilitate transfer of large intact DNA inserts from BAC clones to the BIBAC vector and functional complementation of large DNA fragments. These seven Zea Map Alignment Project (ZMAP) BAC/BIBAC libraries have average insert sizes ranging from 92 to 148 kb, organellar DNA from 0.17 to 2.3%, empty vector rates between 0.35 and 5.56%, and genome equivalents of 4.7- to 8.4-fold. The usefulness of the Parviglumis and Tripsacum BAC libraries was demonstrated by mapping clones to the reference genome. Novel genes and alleles present in these ZMAP libraries can now be used for functional complementation studies and positional or homology-based cloning of genes for translational genomics.

Shark IgW C region diversification through RNA processing and isotype switching

Sharks and skates represent the earliest vertebrates with an adaptive immune system based on lymphocyte Ag receptors generated by V(D)J recombination. Shark B cells express two classical Igs, IgM and IgW, encoded by an early, alternative gene organization consisting of numerous autonomous miniloci, where the individual gene cluster carries a few rearranging gene segments and one C region, μ or ω. We have characterized eight distinct Ig miniloci encoding the nurse shark ω H chain. Each cluster consists of VH, D, and JH segments and six to eight C domain exons. Two interspersed secretory exons, in addition to the 3'-most C exon with tailpiece, provide the gene cluster with the ability to generate at least six secreted isoforms that differ as to polypeptide length and C domain combination. All clusters appear to be functional, as judged by the capability for rearrangement and absence of defects in the deduced amino acid sequence. We previously showed that IgW VDJ can perform isotype switching to μ C regions; in this study, we found that switching also occurs between ω clusters. Thus, C region diversification for any IgW VDJ can take place at the DNA level by switching to other ω or μ C regions, as well as by RNA processing to generate different C isoforms. The wide array of pathogens recognized by Abs requires different disposal pathways, and our findings demonstrate complex and unique pathways for C effector function diversity that evolved independently in cartilaginous fishes.

Construction and preliminary characterization analysis of Wuzhishan miniature pig bacterial artificial chromosome library with approximately 8-fold genome equivalent coverage

Bacterial artificial chromosome (BAC) libraries have been invaluable tools for the genome-wide genetic dissection of complex organisms. Here, we report the construction and characterization of a high-redundancy BAC library from a very valuable pig breed in China, Wuzhishan miniature pig (Sus scrofa), using its blood cells and fibroblasts, respectively. The library contains approximately 153,600 clones ordered in 40 superpools of 10 × 384-deep well microplates. The average insert size of BAC clones was estimated to be 152.3 kb, representing approximately 7.68 genome equivalents of the porcine haploid genome and a 99.93% statistical probability of obtaining at least one clone containing a unique DNA sequence in the library. 19 pairs of microsatellite marker primers covering porcine chromosomes were used for screening the BAC library, which showed that each of these markers was positive in the library; the positive clone number was 2 to 9, and the average number was 7.89, which was consistent with 7.68-fold coverage of the porcine genome. And there were no significant differences of genomic BAC library from blood cells and fibroblast cells. Therefore, we identified 19 microsatellite markers that could potentially be used as genetic markers. As a result, this BAC library will serve as a valuable resource for gene identification, physical mapping, and comparative genomics and large-scale genome sequencing in the porcine.

Dynamic intra-japonica subspecies variation and resource application

We constructed a physical map of O. sativa ssp. japonica cv. ZH11 and compared it and its random sample sequences with the Nipponbare RefSeq derived from the same subspecies. This comparison showed that the two japonica genomes were highly syntenic but revealed substantial differences in terms of structural variations, rates of substitutions and indels, and transposable element content. For example, contractions/expansions as large as 450 kb and repeat sequences that were present in high copy numbers only in ZH11 were detected. In tri-alignment regions using the indica variety 93-11 sequence as an outgroup, we found that: (1) the substitution rates of the two japonica-indica inter-subspecies comparison combinations were close but almost a magnitude higher than the substitution rate between the japonica rice varieties ZH11 and Nipponbare; (2) of the substitutions found between ZH11 and Nipponbare, 47.2% occurred in ZH11 and 52.6% in Nipponbare; (3) of the indels found between ZH11 and Nipponbare, the indels that occurred in ZH11 were 15.8 times of those in Nipponbare. Of the indels that occurred in ZH11, 75.67% were insertions and 24.33% deletions. Of the indels that occurred in Nipponbare, 48.23% were insertions and 51.77% were deletions. The ZH11 comparative map covered four Nipponbare physical gaps, detected assembly errors in the Nipponbare sequence, and was integrated with the FSTs of a large ZH11 T-DNA insertion mutant library. ZH11 BAC clones can be browsed, searched, and obtained at our website, http://GResource.hzau.edu.cn.

A pair of new BAC and BIBAC vectors that facilitate BAC/BIBAC library construction and intact large genomic DNA insert exchange

BACKGROUND

Large-insert BAC and BIBAC libraries are important tools for structural and functional genomics studies of eukaryotic genomes. To facilitate the construction of BAC and BIBAC libraries and the transfer of complete large BAC inserts into BIBAC vectors, which is desired in positional cloning, we developed a pair of new BAC and BIBAC vectors.

RESULTS

The new BAC vector pIndigoBAC536-S and the new BIBAC vector BIBAC-S have the following features: 1) both contain two 18-bp non-palindromic I-SceI sites in an inverted orientation at positions that flank an identical DNA fragment containing the lacZ selection marker and the cloning site. Large DNA inserts can be excised from the vectors as single fragments by cutting with I-SceI, allowing the inserts to be easily sized. More importantly, because the two vectors contain different antibiotic resistance genes for transformant selection and produce the same non-complementary 3' protruding ATAA ends by I-SceI that suppress self- and inter-ligations, the exchange of intact large genomic DNA inserts between the BAC and BIBAC vectors is straightforward; 2) both were constructed as high-copy composite vectors. Reliable linearized and dephosphorylated original low-copy pIndigoBAC536-S and BIBAC-S vectors that are ready for library construction can be prepared from the high-copy composite vectors pHZAUBAC1 and pHZAUBIBAC1, respectively, without the need for additional preparation steps or special reagents, thus simplifying the construction of BAC and BIBAC libraries. BIBAC clones constructed with the new BIBAC-S vector are stable in both E. coli and Agrobacterium. The vectors can be accessed through our website http://GResource.hzau.edu.cn.

CONCLUSIONS

The two new vectors and their respective high-copy composite vectors can largely facilitate the construction and characterization of BAC and BIBAC libraries. The transfer of complete large genomic DNA inserts from one vector to the other is made straightforward.

The multiple shark Ig H chain genes rearrange and hypermutate autonomously

Sharks and skates are representatives of the earliest vertebrates with an immune system based on V(D)J rearrangement. They possess a unique Ig gene organization consisting of 15 to >50 individual IgM loci, each with one VH, two DH, one JH, and one set of constant region exons. The present study attempts to understand how multiple Ig genes are regulated with respect to rearrangement initiation and to targeting during somatic hypermutation. The linkage of three single-copy IgH genes was determined, and single-cell genomic PCR studies in a neonatal animal were used to examine any relationship between relative gene position and likelihood of rearrangement. Our results show that one to three IgH genes are activated independently of linkage or allelic position and the data best fit with a probability model based on the hypothesis that V(D)J rearrangement occurs as a sequence of trials within the B cell. In the neonatal cell set, two closely related IgH, G2A, and G2B, rearranged at similar frequencies, and their membrane forms were expressed at similar levels, like in other young animals. However, older animals displayed a bias in favor of the G2A isotype, which suggests that although rearrangement at G2A and G2B was randomly initiated during primary repertoire generation, the two very similar IgM sequences appear to be differentially expressed with age and exposure to Ag. We performed genomic single-cell PCR on B cells from an immunized individual to study activation-induced cytidine deaminase targeting and found that hypermutation, like V(D)J rearrangement, occurred independently among the many shark IgH.

The 19 genomes of Drosophila: a BAC library resource for genus-wide and genome-scale comparative evolutionary research

The genus Drosophila has been the subject of intense comparative phylogenomics characterization to provide insights into genome evolution under diverse biological and ecological contexts and to functionally annotate the Drosophila melanogaster genome, a model system for animal and insect genetics. Recent sequencing of 11 additional Drosophila species from various divergence points of the genus is a first step in this direction. However, to fully reap the benefits of this resource, the Drosophila community is faced with two critical needs: i.e., the expansion of genomic resources from a much broader range of phylogenetic diversity and the development of additional resources to aid in finishing the existing draft genomes. To address these needs, we report the first synthesis of a comprehensive set of bacterial artificial chromosome (BAC) resources for 19 Drosophila species from all three subgenera. Ten libraries were derived from the exact source used to generate 10 of the 12 draft genomes, while the rest were generated from a strategically selected set of species on the basis of salient ecological and life history features and their phylogenetic positions. The majority of the new species have at least one sequenced reference genome for immediate comparative benefit. This 19-BAC library set was rigorously characterized and shown to have large insert sizes (125-168 kb), low nonrecombinant clone content (0.3-5.3%), and deep coverage (9.1-42.9×). Further, we demonstrated the utility of this BAC resource for generating physical maps of targeted loci, refining draft sequence assemblies and identifying potential genomic rearrangements across the phylogeny.

Primordial linkage of β2-microglobulin to the MHC

β2-Microglobulin (β2M) is believed to have arisen in a basal jawed vertebrate (gnathostome) and is the essential L chain that associates with most MHC class I molecules. It contains a distinctive molecular structure called a constant-1 Ig superfamily domain, which is shared with other adaptive immune molecules including MHC class I and class II. Despite its structural similarity to class I and class II and its conserved function, β2M is encoded outside the MHC in all examined species from bony fish to mammals, but it is assumed to have translocated from its original location within the MHC early in gnathostome evolution. We screened a nurse shark bacterial artificial chromosome library and isolated clones containing β2M genes. A gene present in the MHC of all other vertebrates (ring3) was found in the bacterial artificial chromosome clone, and the close linkage of ring3 and β2M to MHC class I and class II genes was determined by single-strand conformational polymorphism and allele-specific PCR. This study satisfies the long-held conjecture that β2M was linked to the primordial MHC (Ur MHC); furthermore, the apparent stability of the shark genome may yield other genes predicted to have had a primordial association with the MHC specifically and with immunity in general.

Elasmobranch qPCR reference genes: a case study of hypoxia preconditioned epaulette sharks

Background

Elasmobranch fishes are an ancient group of vertebrates which have high potential as model species for research into evolutionary physiology and genomics. However, no comparative studies have established suitable reference genes for quantitative PCR (qPCR) in elasmobranchs for any physiological conditions. Oxygen availability has been a major force shaping the physiological evolution of vertebrates, especially fishes. Here we examined the suitability of 9 reference candidates from various functional categories after a single hypoxic insult or after hypoxia preconditioning in epaulette shark (Hemiscyllium ocellatum).

Results

Epaulette sharks were caught and exposed to hypoxia. Tissues were collected from 10 controls, 10 individuals with single hypoxic insult and 10 individuals with hypoxia preconditioning (8 hypoxic insults, 12 hours apart). We produced sequence information for reference gene candidates and monitored mRNA expression levels in four tissues: cerebellum, heart, gill and eye. The stability of the genes was examined with analysis of variance, geNorm and NormFinder. The best ranking genes in our study were eukaryotic translation elongation factor 1 beta (eef1b), ubiquitin (ubq) and polymerase (RNA) II (DNA directed) polypeptide F (polr2f). The performance of the ribosomal protein L6 (rpl6) was tissue-dependent. Notably, in one tissue the analysis of variance indicated statistically significant differences between treatments for genes that were ranked as the most stable candidates by reference gene software.

Conclusions

Our results indicate that eef1b and ubq are generally the most suitable reference genes for the conditions and tissues in the present epaulette shark studies. These genes could also be potential reference gene candidates for other physiological studies examining stress in elasmobranchs. The results emphasise the importance of inter-group variation in reference gene evaluation.

The evolution of multiple isotypic IgM heavy chain genes in the shark

The IgM H chain gene organization of cartilaginous fishes consists of 15-200 miniloci, each with a few gene segments (V(H)-D1-D2-J(H)) and one C gene. This is a gene arrangement ancestral to the complex IgH locus that exists in all other vertebrate classes. To understand the molecular evolution of this system, we studied the nurse shark, which has relatively fewer loci, and characterized the IgH isotypes for organization, functionality, and the somatic diversification mechanisms that act upon them. Gene numbers differ slightly between individuals ( approximately 15), but five active IgM subclasses are always present. Each gene undergoes rearrangement that is strictly confined within the minilocus; in B cells there is no interaction between adjacent loci located > or =120 kb apart. Without combinatorial events, the shark IgM H chain repertoire is based on junctional diversity and, subsequently, somatic hypermutation. We suggest that the significant contribution by junctional diversification reflects the selected novelty introduced by RAG in the early vertebrate ancestor, whereas combinatorial diversity coevolved with the complex translocon organization. Moreover, unlike other cartilaginous fishes, there are no germline-joined VDJ at any nurse shark mu locus, and we suggest that such genes, when functional, are species-specific and may have specialized roles. With an entire complement of IgM genes available for the first time, phylogenetic analyses were performed to examine how the multiple Ig loci evolved. We found that all domains changed at comparable rates, but V(H) appears to be under strong positive selection for increased amino acid sequence diversity, and surprisingly, so does Cmicro2.

Construction of a BAC library and mapping BAC clones to the linkage map of Barramundi, Lates calcarifer

BACKGROUND

Barramundi (Lates calcarifer) is an important farmed marine food fish species. Its first generation linkage map has been applied to map QTL for growth traits. To identify genes located in QTL responsible for specific traits, genomic large insert libraries are of crucial importance. We reported herein a bacterial artificial chromosome (BAC) library and the mapping of BAC clones to the linkage map.

RESULTS

This BAC library consisted of 49,152 clones with an average insert size of 98 kb, representing 6.9-fold haploid genome coverage. Screening the library with 24 microsatellites and 15 ESTs/genes demonstrated that the library had good genome coverage. In addition, 62 novel microsatellites each isolated from 62 BAC clones were mapped onto the first generation linkage map. A total of 86 BAC clones were anchored on the linkage map with at least one BAC clone on each linkage group.

CONCLUSION

We have constructed the first BAC library for L. calcarifer and mapped 86 BAC clones to the first generation linkage map. This BAC library and the improved linkage map with 302 DNA markers not only supply an indispensable tool to the integration of physical and linkage maps, the fine mapping of QTL and map based cloning genes located in QTL of commercial importance, but also contribute to comparative genomic studies and eventually whole genome sequencing.