Non-invasive age estimation based on faecal DNA using methylation-sensitive high-resolution melting for Indo-Pacific bottlenose dolphins.

Age is necessary information for the study of life history of wild animals. A general method to estimate the age of odontocetes is counting dental growth layer groups (GLGs). However, this method is highly invasive as it requires the capture and handling of individuals to collect their teeth. Recently, the development of DNA-based age estimation methods has been actively studied as an alternative to such invasive methods, of which many have relied on used biopsy samples. However, if DNA-based age estimation can be developed from faecal samples, age estimation can be performed entirely non-invasively. We developed an age estimation model using the methylation rate of two gene regions, GRIA2 and CDKN2A, measured through methylation-sensitive high-resolution melting (MS-HRM) from faecal samples of wild Indo-Pacific bottlenose dolphins (Tursiops aduncus). The age of individuals was known through conducting longitudinal individual identification surveys underwater. Methylation rates were quantified from 36 samples collected from 30 individuals. Both gene regions showed a significant correlation between age and methylation rate. The age estimation model was constructed based on the methylation rates of both genes which achieved sufficient accuracy (after LOOCV: MAE = 5.08, R2 = 0.33) for the ecological studies of the Indo-Pacific bottlenose dolphins, with a lifespan of 40-50 years. This is the first study to report the use of non-invasive faecal samples to estimate the age of marine mammals.

A wild Indo-Pacific bottlenose dolphin adopts a socially and genetically distant neonate.

Alloparental behaviour and adoption have been reported in many mammals and birds. Such behaviours are energetically costly, and their causes and functions remain unclear. We observed the adoption behaviour of a wild Indo-Pacific bottlenose dolphin (Tursiops aduncus) near Mikura Island, Japan. A calf was seen with its mother on six observation days. Following the mother's death, the calf was observed with a sub-adult female on all 18 observation days from May to September 2012. On three days, the calf was observed swimming with this female in the suckling position and milk was seen leaking from the female's mammary slit. A five-year dataset revealed no significant social or kin relationships between the biological mother and allomother, indicating that kinship and social relationships did not play an important role in the observed adoption.

Porcine MHC classical class I genes are coordinately expressed in superantigen-activated mononuclear cells.

The expression of the major histocompatibility complex (MHC) classical class I genes is important for the adaptive immune response to target virus-infected cells and cancer cells. The up-regulation of the MHC is achieved by hormonal/cytokine signals including IFN-γ-inducible elements. The swine leukocyte antigen (SLA), the MHC class I region of pigs, consists of the duplicated classical class I genes, SLA-1, SLA-2 and SLA-3, but the molecular mechanisms involved in their up-regulation after T cell stimulation have not been fully elucidated. In order to better understand some of the putative regulatory mechanisms of SLA class I gene expression in activated T cells, we examined the coordinated expression of the SLA classical class I, IFN-γ and interferon regulatory factor-1 (IRF-1) genes in the peripheral blood mononuclear cells (PBMCs) of SLA homozygous Clawn miniature swine stimulated for 72 h with either IFN-γ or an enterotoxin produced by Staphylococcus aureus. This enterotoxin, toxic shock syndrome-1 (TSST-1), is known to act as a superantigen (sAG) to activate the T cells in various vertebrate species. We showed by using mAbs and flow cytometry that the CD4(+)CD25(+) cell number of swine PBMCs was also increased by TSST-1 and to a lesser degree by IFN-γ. Time course analyses of the expression of the IFN-γ, IRF-1 and the three classical class I genes, SLA-1, SLA-2, and SLA-3, in PBMCs by quantitative real-time PCR revealed a transitory response to TSST-1 or IFN-γ stimulation. The IFN-γ mRNA levels in the PBMCs were continuously up-regulated over the first 48 h by TSST-1 or IFN-γ. In contrast, SLA class I expression moderately increased at 24h and then decreased to a baseline level or less at 72 h of IFN-γ or TSST-1 stimulation. The three classical SLA class I genes showed similar expression kinetics, although SLA-3 mRNA level was consistently lower than those of SLA-1 and -2. The expression of IRF-1, a modulator of SLA expression, showed similar kinetics to those of the three classical SLA class I genes. The expression profiles detected by flow cytometry of the SLA molecules on the cell surface of PBMCs were maintained at a consistently high level during cell stimulation with either TSST-1 or IFN-γ, which was distinct from the kinetics of mRNA expression. These results showed that miniature swine SLA class I mRNA expression was effectively and equally up-regulated among the three loci and coordinately with IRF-1 gene expression after stimulation of T cell activation by sAG or IFN-γ.

Application of high-resolution, massively parallel pyrosequencing for estimation of haplotypes and gene expression levels of swine leukocyte antigen (SLA) class I genes.

The swine is an important animal model for allo- and xeno-transplantation donor studies, which necessitates an extensive characterization of the expression and sequence variations within the highly polygenic and polymorphic swine leukocyte antigen (SLA) region. Massively parallel pyrosequencing is potentially an effective new 2ndGen method for simultaneous high-throughput genotyping and detection of SLA class I gene expression levels. In this study, we compared the 2ndGen method using the Roche Genome Sequencer 454 FLX with the conventional method using sub-cloning and Sanger sequencing to genotype SLA class I genes in five pigs of the Clawn breed and four pigs of the Landrace breed. We obtained an average of 10.4 SLA class I sequences per pig by the 2ndGen method, consistent with the inheritance data, and an average of only 6.0 sequences by the conventional method. We also performed a correlation analysis between the sequence read numbers obtained by the 2ndGen method and the relative expression values obtained by quantitative real-time PCR analysis at the allele level. A significant correlation coefficient (r = 0.899, P < 0.01) was observed between the sequence read numbers and the relative quantitative values for the expressed classical SLA class I genes SLA-1, SLA-2, and SLA-3, suggesting that the sequence read numbers closely reflect the gene expression levels in white blood cells. Overall, five novel class I sequences, different haplotype-specific expression patterns and a splice variant for one of the SLA class I genes were identified by the 2ndGen method at greater efficiency and sensitivity than the conventional method.

Comparative genome analysis of the major histocompatibility complex (MHC) class I B/C segments in primates elucidated by genomic sequencing in common marmoset (Callithrix jacchus).

Common marmoset monkeys (Callithrix jacchus) have emerged as important animal models for biomedical research, necessitating a more extensive characterization of their major histocompatibility complex (MHC) region. However, the genomic information of the marmoset MHC (Caja) is still lacking. The MHC-B/C segment represents the most diverse MHC region among primates. Therefore, in this paper, to elucidate the detailed gene organization and evolutionary processes of the Caja class I B (Caja-B) segment, we determined two parts of the Caja-B sequences with 1,079 kb in total, ranging from H6orf15 to BAT1 and compared the structure and phylogeny with that of other primates. This segment contains 54 genes in total, nine Caja-B genes (Caja-B1 to Caja-B9), two MIC genes (MIC1 and MIC2), eight non-MHC genes, two non-coding genes, and 33 non-MHC pseudogenes that have not been observed in other primate MHC-B/C segments. Caja-B3, Caja-B4, and Caja-B7 encode proper MHC class I proteins according to amino acid structural characteristics. Phylogenetic relationships based on 48 MHC-I nucleotide sequences in primates suggested (1) species-specific divergence for Caja, Mamu, and HLA/Patr/Gogo lineages, (2) independent generation of the "seven coding exon" type MHC-B genes in Mamu and HLA/Patr/Gogo lineages from an ancestral "eight coding exon" type MHC-I gene, (3) parallel correlation with the long and short segmental duplication unit length in Caja and Mamu lineages. These findings indicate that the MHC-B/C segment has been under permanent selective pressure in the evolution of primates.

Study of cynomolgus monkey (Macaca fascicularis) DRA polymorphism in four populations.

To describe the polymorphism of the DRA gene in Macaca fascicularis, we have studied 141 animals either at cDNA level (78 animals from Mauritius, the Philippines, and Vietnam) or genomic level (63 animals from the Philippines, Indonesia, and Vietnam). In total, we characterized 22 cDNA DRA alleles, 13 of which had not been described until now. In the Mauritius population, we confirmed the presence of three DRA alleles. In the Philippine and Vietnam populations, we observed 11 and 14 DRA alleles, respectively. Only two alleles were present in all three populations. All DRA alleles but one differ from the consensus sequence by one to three mutations, most being synonymous; so, only seven DR alpha proteins were deduced from the 22 cDNA alleles. One DRA cDNA allele, Mafa-DRA*02010101, differs from all other alleles by 11 to 14 mutations of which only four are non-synonymous. The two amino acid changes inside the peptide groove of Mafa-DRA*02010101 are highly conservative. The very low proportion of non-synonymous/synonymous mutations is compatible with a purifying selection which is comparable to all previous observations concerning the evolution of the DRA gene in mammals. Homologues of the allele Mafa-DRA*02010101 are also found in two other Asian macaques (Macaca mulatta and Macaca nemestrina). The forces able to maintain this highly divergent allele in three different macaque species remain hypothetical.

MHC class I A loci polymorphism and diversity in three Southeast Asian populations of cynomolgus macaque.

Cynomolgus macaques (Macaca fascicularis, Mafa) have emerged as important animal models for biomedical research, necessitating a more extensive characterization of their major histocompatibility complex polymorphic regions. The current information on the polymorphism or diversity of the polygenetic Mafa class I A loci is limited in comparison to the more commonly studied rhesus macaque Mafa class I A loci. Therefore, in this paper, to better elucidate the degree and types of polymorphisms and genetic differences of Mafa-A1 among three native Southeast Asian populations (Indonesian, Vietnamese, and Filipino) and to investigate how the allele differences between macaques and humans might have evolved to affect their respective immune responses, we identified 83 Mafa-A loci-derived alleles by DNA sequencing of which 66 are newly described. Most alleles are unique to each population, but seven of the most frequent alleles were identical in sequence to some alleles in other macaque species. We also revealed (1) the large and dynamic genetic and structural differences and similarities in allelic variation by analyzing the population allele frequencies, Hardy-Weinberg's equilibrium, heterozygosity, nucleotide diversity profiles, and phylogeny, (2) the difference in genetic structure of populations by Wright's FST statistic and hierarchical analysis of molecular variance, and (3) the different demographic and selection pressures on the three populations by performing Tajima's D test of neutrality. The large level of diversity and polymorphism at the Mafa-A1 was less evident in the Filipino than in the Vietnam or the Indonesian populations, which may have important implications in animal capture, selection, and breeding for medical research.

HLA alleles are associated with postherpetic neuralgia but not with herpes zoster.

OBJECTIVE: In some herpes zoster (HZ) patients, one symptom is problematic; postherpetic neuralgia (PHN), a persistence of pain such as spontaneous pain and stimulus-evoked pain, allodynia or hyperpathia, for more than 6 months after healing of the vesicular eruptions. In our previous study, we reported the association between HLA alleles, HLA-A*33 and B*44, and PHN patients. In this study, we aim to refine the association of these alleles with PHN or HZ using higher-resolution HLA typing technique with an increased sample size. METHODS: HLA allele frequencies were compared in 70 PHN patients, 80 HZ patients, 52 HZ(-) patients, and 140 Japanese controls using HLA typing kits of SSOP protocols. RESULT AND CONCLUSION: The allele frequencies of HLA-A*3303, B*4403, and DRB1*1302 in PHN(+) patients were significantly higher than those in Japanese controls (P=0.00007, P=0.000002, and P=0.0003, respectively). The frequencies of above alleles in PHN(+) patients were also significantly higher than those in PHN(-) patients (P=0.03, P=0.006 and P=0.03, respectively). However, no association was found in comparison of HZ(+) patients and HZ(-). And the frequency of HLA-B*5101 in HZ(-) patients was significantly higher than those in HZ(+) and PHN(+) patients (P=0.003 and P=0.009, respectively) indicating that HLA-B*5101 functions as a protective allele to HZ.

A BAC-based contig map of the cynomolgus macaque (Macaca fascicularis) major histocompatibility complex genomic region.

The construction of a cynomolgus macaque (Macaca fascicularis, Mafa) BAC library for genomic comparison between rhesus and cynomolgus macaques is necessary to promote the cynomolgus macaque as one of the important experimental animals for future medical and biological research. In this paper, we constructed a cynomolgus macaque BAC library and a map of the MHC (Mafa) genomic region for comparison of the genomic organization and nucleotide similarities between the human, the chimpanzee, and the rhesus macaque. The BAC library consists of 221,184 clones with an average insert size of 83 kb, providing a sixfold coverage of the haploid genome. A total of 114 BAC clones and 54 PCR primer sets were used to construct a 4.3-Mb contig of the MHC region. Diversity analysis of genomic sequence from selected subregions of the MHC revealed that the cynomolgus sequence varied compared to rhesus macaque, human, and chimpanzee sequences by 0.48, 4.15, and 4.10%, respectively. From these findings, we conclude that the BAC library and Mafa genomic map are useful tools for genome analysis and will have important applications for comparative genomics and identifying regions of consequence in medical research.