Structure of the catfish IGH locus: analysis of the region including the single functional IGHM gene

The catfish IGH locus is large (∼1 Mb) and complex, having undergone multiple internal duplications and transpositions. To define the structure of the locus that contains the single expressed IGHM gene, two overlapping bacterial-artificial-chromosome (BAC) clones spanning the most 3′ end of the channel catfish immunoglobulin heavy (IGH) chain locus have been completely sequenced. The analyses created a contig of 257,153 bp containing 55 VH, 6 D, 12 JH genes and the IGH constant region genes encoding the functional secreted and membrane forms of IgM and the membrane form of IgD. This analysis revealed three major features. First, no C-region genes were found aside from the previously described IGHM1 and IGHD1, with the latter gene being the most 3′ C-region gene of the catfish IGH locus. There was no evidence in the region sequenced for genes that could encode an Ig class similar to the IgZ/IgT described in zebrafish, trout and pufferfish. Second, there are a high number of VH pseudogenes, 28 out of 55 (51%). In contrast, the entire zebrafish IGH locus has 40 functional VH genes and eight pseudogenes (17%). Third, an internal duplication of a 52.4-kb block of VH genes has occurred. These observations suggest that the IGH locus of teleost fish varies significantly from species to species in the diversity of C-region genes as well as the numbers of genes encoding V regions.     

The western painted turtle genome,a model for the evolution of extreme physiological adaptations in a slowly evolving lineage

Background

We describe the genome of the western painted turtle, Chrysemys picta bellii, one of the most widespread, abundant, and well-studied turtles. We place the genome into a comparative evolutionary context, and focus on genomic features associated with tooth loss, immune function, longevity, sex differentiation and determination, and the species'' physiological capacities to withstand extreme anoxia and tissue freezing.

Results

Our phylogenetic analyses confirm that turtles are the sister group to living archosaurs, and demonstrate an extraordinarily slow rate of sequence evolution in the painted turtle. The ability of the painted turtle to withstand complete anoxia and partial freezing appears to be associated with common vertebrate gene networks, and we identify candidate genes for future functional analyses. Tooth loss shares a common pattern of pseudogenization and degradation of tooth-specific genes with birds, although the rate of accumulation of mutations is much slower in the painted turtle. Genes associated with sex differentiation generally reflect phylogeny rather than convergence in sex determination functionality. Among gene families that demonstrate exceptional expansions or show signatures of strong natural selection, immune function and musculoskeletal patterning genes are consistently over-represented.

Conclusions

Our comparative genomic analyses indicate that common vertebrate regulatory networks, some of which have analogs in human diseases, are often involved in the western painted turtle''s extraordinary physiological capacities. As these regulatory pathways are analyzed at the functional level, the painted turtle may offer important insights into the management of a number of human health disorders.     

Effects of acclimation and egg-incubation temperature on selected temperature by hatchling western painted turtles (Chrysemys picta bellii)

The ability of hatchling turtles to detect environmental temperature differences and to effectively select preferred temperature is a function that critically impacts survival. In some turtle species, temperature preference may be influenced by embryonic and post-hatching conditions, such as egg-incubation and acclimation temperature. We tested for effects of embryonic incubation temperature (27.5 °C, 30 °C) and acclimation temperature (20 °C, 25 °C) on the selected temperature and movement patterns of 32 Chrysemys picta bellii (Reptilia: Emydidae) hatchlings in an aquatic thermal gradient of 14-34 °C and in single-temperature (20 °C, 25 °C) control tests. Among 10-11 month old hatchlings, acclimation temperature and egg-incubation temperature influenced temperature selection and movement patterns. Acclimation temperature affected activity and movement: in thermal gradient and single-temperature control tests, 25 °C-acclimated turtles relocated between chambers significantly more frequently than individuals acclimated to 20 °C. Acclimation temperature also affected temperature selection: 20 °C-acclimated turtles selected a specific temperature during gradient tests, but 25 °C-acclimated turtles did not. Among 20 °C-acclimated turtles, egg-incubation temperature was inversely related to selected temperature: hatchling turtles incubated at 27.5 °C selected the warmest temperature available (34 °C); individuals incubated at 30 °C selected the coldest temperature (14 °C). These results suggest that interactions of environmental conditions may influence post-hatching thermoregulatory behavior in C. picta bellii, a factor that ultimately affects fitness.     

Characterization of the horse (<Emphasis Type="Italic">Equus caballus</Emphasis>)<Emphasis Type="Italic"> IGHA</Emphasis> gene

Nucleotide sequences of the immunoglobulin constant heavy chain genes of the horse have been described for IGHM, IGHG and IGHE genes, but not for IGHA. Here, we provide the nucleotide sequence of the genomic IGHA gene of the horse (Equus caballus), including its secretion region and the transmembrane exon. The equine IGHA gene shows the typical structure of a mammalian IGHA gene, with only three exons, separated by two introns of similar size. The hinge exon is located at the 5 end of the CH2 exon and encodes a hinge region of 11 amino acids, which contains five proline residues. The coding nucleotide sequence of the secreted form of the equine IGHA gene shares around 72% identity with the human IGHA1 and IGHA2 genes, as well as the bovine, ovine, porcine and canine IGHA genes, without distinct preference for any of these species. The same species also cluster together in a phylogenetic tree of the IGHA coding regions of various mammals, whereas rodent, rabbit, marsupial and monotreme IGHA genes each build a separate cluster.The nucleotide sequences reported in this paper have been assigned the EMBL/GenBank accession numbers AY247966 and AY351982     

Heritable variation in heat shock gene expression: a potential mechanism for adaptation to thermal stress in embryos of sea turtles

The capacity of species to respond adaptively to warming temperatures will be key to their survival in the Anthropocene. The embryos of egg-laying species such as sea turtles have limited behavioural means for avoiding high nest temperatures, and responses at the physiological level may be critical to coping with predicted global temperature increases. Using the loggerhead sea turtle (Caretta caretta) as a model, we used quantitative PCR to characterise variation in the expression response of heat-shock genes (hsp60, hsp70 and hsp90; molecular chaperones involved in cellular stress response) to an acute non-lethal heat shock. We show significant variation in gene expression at the clutch and population levels for some, but not all hsp genes. Using pedigree information, we estimated heritabilities of the expression response of hsp genes to heat shock and demonstrated both maternal and additive genetic effects. This is the first evidence that the heat-shock response is heritable in sea turtles and operates at the embryonic stage in any reptile. The presence of heritable variation in the expression of key thermotolerance genes is necessary for sea turtles to adapt at a molecular level to warming incubation environments.     

Organization of the variable region of the immunoglobulin heavy-chain gene locus of the rat

We have mapped and annotated the variable region of the immunoglobulin heavy (IGH) gene locus of the Brown Norway (BN) rat (assembly V3.4; Rat Genomic Sequence Consortium). In addition to known variable region genes, we found 12 novel previously unidentified functional IGHV genes and 1 novel functional IGHD gene. In total, the variable region of the rat IGH locus is composed of at least 353 unique IGHV genes, 21 IGHD genes, and 5 IGHJ genes, of which 131, 14, and 4 are potentially functional genes, respectively. Of all species studied so far, the rat seems to have the highest number of functional IGHV genes in the genome. Rat IGHV genes can be classified into 13 IGHV families based on nucleotide sequence identity. The variable region of the BN rat spans a total length of approximately 4.9 Mb and is organized in a typical translocon organization. Like the mouse, members of the various IGHV gene families are more or less grouped together on the genome, albeit some members of IGHV gene families are found intermingled with each other. In the rat, the largest IGHV gene families are IGHV1, IGHV2, and IGHV5. The overall conclusion is that the genomic organization of the variable region of the rat IGH locus is strikingly similar to that of the mouse, illustrating the close evolutionary relationship between these two species.     

A Comprehensive Analysis of the Phylogeny,Genomic Organization and Expression of Immunoglobulin Light Chain Genes in Alligator sinensis,an Endangered Reptile Species

Crocodilians are evolutionarily distinct reptiles that are distantly related to lizards and are thought to be the closest relatives of birds. Compared with birds and mammals, few studies have investigated the Ig light chain of crocodilians. Here, employing an Alligator sinensis genomic bacterial artificial chromosome (BAC) library and available genome data, we characterized the genomic organization of the Alligator sinensis IgL gene loci. The Alligator sinensis has two IgL isotypes, λ and κ, the same as Anolis carolinensis. The Igλ locus contains 6 C_λ genes, each preceded by a J_λ gene, and 86 potentially functional V_λ genes upstream of (J_λ-C_λ)_n. The Igκ locus contains a single C_κ gene, 6 J_κs and 62 functional V_κs. All V_L genes are classified into a total of 31 families: 19 V_λ families and 12 V_κ families. Based on an analysis of the chromosomal location of the light chain genes among mammals, birds, lizards and frogs, the data further confirm that there are two IgL isotypes in the Alligator sinensis: Igλ and Igκ. By analyzing the cloned Igλ/κ cDNA, we identified a biased usage pattern of V families in the expressed V_λ and V_κ. An analysis of the junctions of the recombined VJ revealed the presence of N and P nucleotides in both expressed λ and κ sequences. Phylogenetic analysis of the V genes revealed V families shared by mammals, birds, reptiles and Xenopus, suggesting that these conserved V families are orthologous and have been retained during the evolution of IgL. Our data suggest that the Alligator sinensis IgL gene repertoire is highly diverse and complex and provide insight into immunoglobulin gene evolution in vertebrates.     

CRISPR Diversity in E. coli Isolates from Australian Animals,Humans and Environmental Waters

Seventy four SNP genotypes and 54 E. coli genomes from kangaroo, Tasmanian devil, reptile, cattle, dog, horse, duck, bird, fish, rodent, human and environmental water sources were screened for the presence of the CRISPR 2.1 loci flanked by cas2 and iap genes. CRISPR 2.1 regions were found in 49% of the strains analysed. The majority of human E. coli isolates lacked the CRISPR 2.1 locus. We described 76 CRISPR 2.1 positive isolates originating from Australian animals and humans, which contained a total of 764 spacer sequences. CRISPR arrays demonstrated a long history of phage attacks especially in isolates from birds (up to 40 spacers). The most prevalent spacer (1.6%) was an ancient spacer found mainly in human, horse, duck, rodent, reptile and environmental water sources. The sequence of this spacer matched the intestinal P7 phage and the pO111 plasmid of E. coli.     

Digital Gene Expression Profiling reveals transcriptional responses to acute cold stress in Chinese soft-shelled turtle Pelodiscus sinensis juveniles

Turtles are well known for their stress tolerance, including an ability to deal with temperature extremes or rapid thermal change. To know more about the comprehensive molecular basis of thermal stress responses in turtles, we assessed differentially expressed genes (DEGs) in the brain, liver and kidney of juvenile soft-shelled turtles, Pelodiscus sinensis, after acute cold stress (28 °C–8 °C acute transfer and held for 12 h) and following recovery (back to 28 °C and held for 24 h) by digital gene expression profiling. Selected DEGs were also validated via real-time PCR. We found the fewest DEGs in the brain, only one-tenth of the number seen in liver, indicating a tissue-specific gene expression pattern. The DEGs indicated the potential activation of several important functions in response to cold stress and recovery in P. sinensis. This included response to oxidative stress or regulation of reactive oxygen species metabolism in the brain and liver, cerebral inositol metabolism, hepatic monosaccharide metabolism, hepatic complement system, renal DNA repair mechanisms, and TNF and PI3K-Akt signaling pathways in the kidney. These functions likely responded to cold stress in different tissues of P. sinensis to help minimize or repair cell damage as well as enhance innate immunity. The outcomes of this study provide some fundamental insight into the tissue specific complex mechanisms underlining cold stress responses in the soft-shelled turtle P. sinensis.     

When the shoe doesn’t fit: applying conservation unit concepts to western painted turtles at their northern periphery

As biodiversity continues to be lost at an alarming rate, strategies for prioritizing populations for conservation have become increasingly important. Maintaining intraspecific genetic diversity is of particular importance for preserving evolutionary history and the potential for future adaptation. In order to effectively protect this diversity, units below the species level need to be defined. However, delineation of such units is subject to many challenges, with no one strategy applying universally across taxa. In this study we carried out the first genetic assessment of the western painted turtle (Chrysemys picta bellii) at its northern periphery in British Columbia (BC), Canada, using mitochondrial DNA haplotypic and microsatellite genotypic data to examine population structure and demographic history. We compared the application of evolutionarily significant unit and management unit criteria with Canadian designatable unit guidelines to determine appropriate conservation units. Our results show that BC western painted turtles form a single evolutionarily significant unit, with each occupied site constituting a separate management unit. In contrast, there is evidence for six discrete designatable units. Patterns of genetic variation in BC western painted turtles indicate that the conservation of each region is important to maintaining regional diversity and evolutionary novelty in this widespread species.     

Comparison research and phylogenetic implications of mitochondrial control regions in four soft-shelled turtles of Trionychia (Reptilia, Testudinata)

The mitochondrial control regions (CRs) and flanking sequences of Pelodiscus sinensis, Apalone ferox, Palea steindachneri and Carettochelys insculpta were obtained using Long-PCR with gene-specific primers. The CR lengths of the four species were 1843 bp, 1356 bp, 1725 bp, and 969 bp. The base composition percentages of A+T were 60.5%, 60.7%, 65.7%, 64.7%, respectively. Combined with CR sequences of other three soft-shelled turtles published in GenBank (Pelodiscus sinensis, Korea, AY962573; Dogania subplana, AF366350; Lissemys punctata, EF050073), we compared the CR structures and identified three functional domains (TAS, CD and CSB) in which conserved sequence blocks (TAS, CSB -F, CSB-1, CSB-2 and CSB-3) were also successfully identified according to their sequence similarities to those of other turtles. The variable numbers of tandem repeats (VNTRs 1) with 50–52 bp motif were identified at 5′-end of CR among the five soft-shelled turtles P. sinensis (China), P. sinensis (Korea), A. ferox, P. steindachneri, D. subplana. The copy number of the VNTRs varied from 5 to 15. VNTRs 2 with 2–11 bp motif were identified in the 3′- end of CR among all of the six soft-shelled turtles with variable number of motifs from 4 to 29. Moreover, VNTRs 3 with 6 bp motif were identified between CSB-1 and CSB-2 of CR both in P. sinensis (China) and P. sinensis (Korea), in which the number of motifs varied from 19 to 29. The types and distribution of VNTRs of the six soft-shelled turtles were also discussed. With Alligator mississippiensis as an outgroup, combined with the CR sequences (excluding VNTRs) of other five turtles which were published in GenBank, the molecular phylogenetic trees were constructed using PAUP 4.0b10 and MrBayes ver. 3.0. The results strongly supported the monophyly of Carretochelyidae and Carettochelyidae as sister group to an assemblage of Cryptodira. Our research suggested that the earliest phylogenetic tree splits into three separated basal branches; the Pelomedusidira (Pelomedusa subrufa), the Carettochelyidae (C. insculpta), and an assemblage of Cryptodira and the C. insculpta that might be a representation of distinctive suborder.     

Structural organization,classification and phylogenetic relationship of cytochrome P450 genes in Citrus clementina and Citrus sinensis

The genus Citrus is an important fruit crop and nutritional source for the good health of humans. Cytochrome P450s represent about 1 % of the proteome and mediate diverse biochemical reactions pertaining to both primary and secondary metabolism. Analysis of Citrus genomic resources identified 296 plant cytochrome P450s (CYP) coding genes in Citrus clementina, 272 in double haploid (dh) Citrus sinensis, and 202 in C. sinensis. In C. clementina and dh C. sinensis, CYP genes are distributed into nine clans. In the three genomes, single intron containing CYP genes are predominant in the A-type families. Among non-A-type CYP families, multiple intron containing genes are predominant. More number of genes in CYP A-type families over non-A-type families is attributed to rapid evolution of A-type genes facilitated by their gene organization. Further, complex gene organization of non-A-type genes with the presence of multiple introns might have contributed to the slower evolvement of paralogs. Majority of introns (1,660) from three genomes showed canonical GT-AG splice sites. However, 33 introns showed non-conventional GC… PyAG splice sites and functionality of these splice sites is confirmed by the ESTs lacking this intron. Across the families, gene organization is conserved between the three genomes. In dh C. sinensis, 22 genes were identified to have alternate splicing. Examination of scaffolds in C. clementina revealed that majority of the Citrus CYP genes are solitary and a few of them are in clusters of 3–8 genes. PCR amplification of C. sinensis genomic DNA with gene-specific primers failed to amplify out-grouped genes Ccl-CYP706A16 and Ccl-CYP706B1, confirming that they are specific to C. clementina. Differential number of CYP genes observed between C. clementina and C. sinensis is attributed to the extent of variability between their parents representing ancestral taxa.     

Polymorphism of the IGHA gene in sheep

Genetic variation in immunoglobulin A, the most abundant immunoglobulin in mammalian cells, has not been reported in ruminants. In this study, variation in the immunoglobulin heavy alpha chain constant gene (IGHA) of sheep was investigated by amplification of a fragment that included the hinge coding sequence, followed by single-strand conformational polymorphism (SSCP) analysis and DNA sequencing. Three novel sequences, each characterized by unique SSCP banding patterns, were identified. One or two sequences were detected in individual sheep and all the sequences identified shared high homology to the published ovine and bovine IGHA sequences, suggesting that these sequences represent allelic variants of the IGHA gene in sheep. Sequence alignment showed that these sequences differed mainly in the 3′ end of exon 1 and in the coding sequence of the hinge region. There was either a deletion or an insertion of two codons in the hinge coding region in these allelic variants. Codon usage in the hinge coding region was quite different from that in the non-hinge coding regions of the gene, suggesting different evolution of the IGHA hinge sequence. Three novel amino acid sequences of ovine IGHA were also predicted, and variation in these sequences might not only affect antigen recognition but also susceptibility to cleavage by bacterial or parasitic proteases. Nucleotide sequence data reported in this paper have been submitted to the NCBI GenBank nucleotide sequence database and have been assigned the accession nos. AY956424–AY956426.     

Accuracy and coverage assessment of Oryctolagus cuniculus (rabbit) genes encoding immunoglobulins in the whole genome sequence assembly (OryCun2.0) and localization of the IGH locus to chromosome 20

We report on the analyses of genes encoding immunoglobulin heavy and light chains in the rabbit 6.51× whole genome assembly. This OryCun2.0 assembly confirms previous mapping of the duplicated IGK1 and IGK2 loci to chromosome 2 and the IGL lambda light chain locus to chromosome 21. The most frequently rearranged and expressed IGHV1 that is closest to IG DH and IGHJ genes encodes rabbit VHa allotypes. The partially inbred Thorbecke strain rabbit used for whole-genome sequencing was homozygous at the IGK but heterozygous with the IGHV1a1 allele in one of 79 IGHV-containing unplaced scaffolds and IGHV1a2, IGHM, IGHG, and IGHE sequences in another. Some IGKV, IGLV, and IGHA genes are also in other unplaced scaffolds. By fluorescence in situ hybridization, we assigned the previously unmapped IGH locus to the q-telomeric region of rabbit chromosome 20. An approximately 3-Mb segment of human chromosome 14 including IGH genes predicted to map to this telomeric region based on synteny analysis could not be located on assembled chromosome 20. Unplaced scaffold chrUn0053 contains some of the genes that comparative mapping predicts to be missing. We identified discrepancies between previous targeted studies and the OryCun2.0 assembly and some new BAC clones with IGH sequences that can guide other studies to further sequence and improve the OryCun2.0 assembly. Complete knowledge of gene sequences encoding variable regions of rabbit heavy, kappa, and lambda chains will lead to better understanding of how and why rabbits produce antibodies of high specificity and affinity through gene conversion and somatic hypermutation.     

Molecular characterization and functional analysis of MyD88 in Chinese soft-shelled turtle Trionyx sinensis

Myeloid differentiation factor 88 (MyD88) is one of the key adaptor proteins to signal transduction that triggers downstream cascades involved in innate immunity. In this study, the MyD88 gene from Chinese soft-shelled turtle (Trionyx sinensis) (tMyD88) was identified, representing the fist example from reptile species. The tMyD88 has a 894-bp ORF and encodes a polypeptide of 297 amino acids including a typical death domain (DD) at the N-terminus and a conservative Toll/IL-1R (TIR) domain at the C-terminus. It was expressed at high levels in spleen, blood, lungs and liver, but marginal in kidneys and intestines of turtles challenged with live cells of Aeromonas hydrophila, as determined by real-time PCR. RAW 264.7 cells transfected with pcDNA-tMyD88 showed higher NF-κB activity than the vector control (673.78 vs 410.72, P < 0.05). Expression of proinflammatory cytokines IL-1β and TNF-α was also significantly higher in RAW 264.7 cells transfected with pcDNA-tMyD88 than those having pcDNA3.1 control vector (P < 0.01). These results indicate that tMyD88 might possess an important role in defense against microbial infection in Chinese soft-shelled turtles similar to that in mammals.     

Isolation of the MAT1-1 mating type idiomorph and evidence for selfing in the Chinese medicinal fungus Ophiocordyceps sinensis

Ophiocordyceps sinensis is one of the most valued medicinal fungi in China. Research on the mating system and sexual development is vitally important to this endangered species. Previous efforts devoted to investigate the mating type (MAT) locus of O. sinensis, however, resulted in an incomplete understanding. In this study, the MAT1-1 locus of O. sinensis was investigated. The conserved α-box and HMG-box regions of the MAT1-1-1 and MAT1-1-3 genes, respectively, and a conserved region of the DNA lyase gene were successfully amplified using degenerate PCR. A combination of TAIL-PCR and long-range PCR were used to connect these genes and obtain the sequence of the MAT1-1 locus. Screening of 22 single spore isolates by PCR demonstrated that both the MAT1-1-1 and MAT1-2-1 genes cooccurred within the same isolate. Additionally, both MAT1-1-1 and MAT1-2-1 are expressed in vegetative mycelia, providing evidence that O. sinensis is likely capable of selfing. DAPI (4,6-diamidino-2-phenylindole) staining of ascospores and hyphae showed that a majority of hyphal compartments are binucleate, suggesting that O. sinensis may be pseudohomothallic. Analyses of sequence diversity showed lower levels of genetic diversity in MAT1-1-1 compared to MAT1-2-1, indicating the possibility that different selective pressures act on the two MAT idiomorphs. The MAT1-1-1 sequences of O. sinensis and Tolypocladium inflatum cluster as a monophyletic group consistent with phylogenetic classification of Ophiocordycipitaceae. Comparison of the structure of the MAT1-1 locus across hypocrealean taxa showed that O. sinensis contains all three mating type genes (MAT1-1-1, MAT1-1-2, and MAT1-1-3) and supported previous observations that of the four families in Hypocreales, MAT1-1-3 has undergone a lineage specific loss only in some members of the Cordycipitaceae.