Construction and analysis of a human-chimpanzee comparative clone map

Generation of chimpanzee induced pluripotent stem cell lines for cross-species comparisons

As humans' closest living relatives, chimpanzees offer valuable insights into human evolution. However, technical and ethical limitations hinder investigations into the molecular and cellular foundations that distinguish chimpanzee and human traits. Recently, induced pluripotent stem cells (iPSCs) have emerged as a novel model for functional comparative studies and provided a non-invasive alternative for studying embryonic phenomena. In this study, we generated five new chimpanzee iPSC lines from peripheral blood cells and skin fibroblasts with SeV vectors carrying four reprogramming factors (human OCT3/4, SOX2, KLF4, and L-MYC) and characterized their pluripotency and differentiation potential. We also examined the expression of a human-specific non-coding RNA, HSTR1, which is predicted to be involved in human brain development. Our results show that the chimpanzee iPSCs possess pluripotent characteristics and can differentiate into various cell lineages. Moreover, we found that HSTR1 is expressed in human iPSCs and their neural derivatives but not in chimpanzee counterparts, supporting its possible role in human-specific brain development. As iPSCs are inherently variable due to genetic and epigenetic differences in donor cells or reprogramming procedures, it is essential to expand the number of chimpanzee iPSC lines to comprehensively capture the molecular and cellular properties representative of chimpanzees. Hence, our cells provide a valuable resource for investigating the function and regulation of human-specific transcripts such as HSTR1 and for understanding human evolution more generally.

The KIR repertoire of a West African chimpanzee population is characterized by limited gene, allele, and haplotype variation

Introduction

The killer cell immunoglobulin-like receptors (KIR) play a pivotal role in modulating the NK cell responses, for instance, through interaction with major histocompatibility complex (MHC) class I molecules. Both gene systems map to different chromosomes but co-evolved during evolution. The human KIR gene family is characterized by abundant allelic polymorphism and copy number variation. In contrast, our knowledge of the KIR repertoire in chimpanzees is limited to 39 reported alleles, with no available population data. Only three genomic KIR region configurations have been mapped, and seventeen additional ones were deduced by genotyping.

Methods

Previously, we documented that the chimpanzee MHC class I repertoire has been skewed due to an ancient selective sweep. To understand the depth of the sweep, we set out to determine the full-length KIR transcriptome - in our MHC characterized pedigreed West African chimpanzee cohort - using SMRT sequencing (PacBio). In addition, the genomic organization of 14 KIR haplotypes was characterized by applying a Cas9-mediated enrichment approach in concert with long-read sequencing by Oxford Nanopore Technologies.

Results

In the cohort, we discovered 35 undescribed and 15 already recorded Patr-KIR alleles, and a novel hybrid KIR gene. Some KIR transcripts are subject to evolutionary conserved alternative splicing events. A detailed insight on the KIR region dynamics (location and order of genes) was obtained, however, only five new KIR region configurations were detected. The population data allowed to investigate the distribution of the MHC-C1 and C2-epitope specificity of the inhibitory lineage III KIR repertoire, and appears to be skewed towards C2.

Discussion

Although the KIR region is known to evolve fast, as observed in other primate species, our overall conclusion is that the genomic architecture and repertoire in West African chimpanzees exhibit only limited to moderate levels of variation. Hence, the ancient selective sweep that affected the chimpanzee MHC class I region may also have impacted the KIR system.

Swine as the Animal Model for Testing New Formulations of Anti-Inflammatory Drugs: Carprofen Pharmacokinetics and Bioavailability of the Intramuscular Route

Carprofen (CP) is a non-steroidal anti-inflammatory drug (NSAID) frequently used to treat respiratory diseases in numerous small animals, but also in large species. CP is a formidable candidate for further therapeutic research of human inflammatory diseases using the pig as an animal model. However, CP administration in swine is very uncommon and respective pharmacokinetics/bioavailability studies are scarce. A simultaneous population pharmacokinetic analysis after CP intravenous and intramuscular administrations in pigs has shown high extent and rate of absorption and a similar distribution profile with respect to man and other mammals. However, clearance and half-life values found in swine suggest a slower elimination process than that observed in man and some other animal species. Although not reported in other species, liver and kidney concentrations achieved at 48 h post-intramuscular administration in pigs were ten times lower than those found in plasma. Simulations pointed to 4 mg/kg every 24 h as the best dosage regimen to achieve similar therapeutic levels to those observed in other animal species. All these findings support the use of pig as an animal model to study the anti-inflammatory effects of CP in humans.

The Genomic Organization of the LILR Region Remained Largely Conserved Throughout Primate Evolution: Implications for Health And Disease

The genes of the leukocyte immunoglobulin-like receptor (LILR) family map to the leukocyte receptor complex (LRC) on chromosome 19, and consist of both activating and inhibiting entities. These receptors are often involved in regulating immune responses, and are considered to play a role in health and disease. The human LILR region and evolutionary equivalents in some rodent and bird species have been thoroughly characterized. In non-human primates, the LILR region is annotated, but a thorough comparison between humans and non-human primates has not yet been documented. Therefore, it was decided to undertake a comprehensive comparison of the human and non-human primate LILR region at the genomic level. During primate evolution the organization of the LILR region remained largely conserved. One major exception, however, is provided by the common marmoset, a New World monkey species, which seems to feature a substantial contraction of the number of LILR genes in both the centromeric and the telomeric region. Furthermore, genomic analysis revealed that the killer-cell immunoglobulin-like receptor gene KIR3DX1, which maps in the LILR region, features one copy in humans and great ape species. A second copy, which might have been introduced by a duplication event, was observed in the lesser apes, and in Old and New World monkey species. The highly conserved gene organization allowed us to standardize the LILR gene nomenclature for non-human primate species, and implies that most of the receptors encoded by these genes likely fulfill highly preserved functions.

Overview of Nonhuman Primate Models of SARS-CoV-2 Infection

COVID-19, the disease caused by the SARS-CoV-2 betacoronavirus, was declared a pandemic by the World Health Organization on March 11, 2020. Since then, SARS-CoV-2 has triggered a devastating global health and economic emergency. In response, a broad range of preclinical animal models have been used to identify effective therapies and vaccines. Current animal models do not express the full spectrum of human COVID-19 disease and pathology, with most exhibiting mild to moderate disease without mortality. NHPs are physiologically, genetically, and immunologically more closely related to humans than other animal species; thus, they provide a relevant model for SARS-CoV-2 investigations. This overview summarizes NHP models of SARS-CoV-2 and their role in vaccine and therapeutic development.

Elucidation of the speciation history of three sister species of crown-of-thorns starfish (Acanthaster spp.) based on genomic analysis

The crown-of-thorns starfish (COTS) is a coral predator that is widely distributed in Indo-Pacific Oceans. A previous phylogenetic study using partial mitochondrial sequences suggested that COTS had diverged into four distinct species, but a nuclear genome-based analysis to confirm this was not conducted. To address this, COTS species nuclear genome sequences were analysed here, sequencing Northern Indian Ocean (NIO) and Red Sea (RS) species genomes for the first time, followed by a comparative analysis with the Pacific Ocean (PO) species. Phylogenetic analysis and ADMIXTURE analysis revealed clear divergences between the three COTS species. Furthermore, within the PO species, the phylogenetic position of the Hawaiian sample was further away from the other Pacific-derived samples than expected based on the mitochondrial data, suggesting that it may be a PO subspecies. The pairwise sequentially Markovian coalescent model showed that the trajectories of the population size diverged by region during the Mid-Pleistocene transition when the sea-level was dramatically decreased, strongly suggesting that the three COTS species experienced allopatric speciation. Analysis of the orthologues indicated that there were remarkable genes with species-specific positive selection in the genomes of the PO and RS species, which suggested that there may be local adaptations in the COTS species.

A comprehensive study of a 29-capsid AAV library in a non-human primate central nervous system

Non-human primates (NHPs) are a preferred animal model for optimizing adeno-associated virus (AAV)-mediated CNS gene delivery protocols before clinical trials. In spite of its inherent appeal, it is challenging to compare different serotypes, delivery routes, and disease indications in a well-powered, comprehensive, multigroup NHP experiment. Here, a multiplex barcode recombinant AAV (rAAV) vector-tracing strategy has been applied to a systemic analysis of 29 distinct, wild-type (WT), AAV natural isolates and engineered capsids in the CNS of eight macaques. The report describes distribution of each capsid in 15 areas of the macaques' CNS after intraparenchymal (putamen) injection, or cerebrospinal fluid (CSF)-mediated administration routes (intracisternal, intrathecal, or intracerebroventricular). To trace the vector biodistribution (viral DNA) and targeted tissues transduction (viral mRNA) of each capsid in each of the analyzed CNS areas, quantitative next-generation sequencing analysis, assisted by the digital-droplet PCR technology, was used. The report describes the most efficient AAV capsid variants targeting specific CNS areas after each route of administration using the direct side-by-side comparison of WT AAV isolates and a new generation of rationally designed capsids. The newly developed bioinformatics and visualization algorithms, applicable to the comparative analysis of several mammalian brain models, have been developed and made available in the public domain.

Multiple Genomic Events Altering Hominin SIGLEC Biology and Innate Immunity Predated the Common Ancestor of Humans and Archaic Hominins

Human-specific pseudogenization of the CMAH gene eliminated the mammalian sialic acid (Sia) Neu5Gc (generating an excess of its precursor Neu5Ac), thus changing ubiquitous cell surface “self-associated molecular patterns” that modulate innate immunity via engagement of CD33-related-Siglec receptors. The Alu-fusion-mediated loss-of-function of CMAH fixed ∼2–3 Ma, possibly contributing to the origins of the genus Homo. The mutation likely altered human self-associated molecular patterns, triggering multiple events, including emergence of human-adapted pathogens with strong preference for Neu5Ac recognition and/or presenting Neu5Ac-containing molecular mimics of human glycans, which can suppress immune responses via CD33-related-Siglec engagement. Human-specific alterations reported in some gene-encoding Sia-sensing proteins suggested a “hotspot” in hominin evolution. The availability of more hominid genomes including those of two extinct hominins now allows full reanalysis and evolutionary timing. Functional changes occur in 8/13 members of the human genomic cluster encoding CD33-related Siglecs, all predating the human common ancestor. Comparisons with great ape genomes indicate that these changes are unique to hominins. We found no evidence for strong selection after the Human–Neanderthal/Denisovan common ancestor, and these extinct hominin genomes include almost all major changes found in humans, indicating that these changes in hominin sialobiology predate the Neanderthal–human divergence ∼0.6 Ma. Multiple changes in this genomic cluster may also explain human-specific expression of CD33rSiglecs in unexpected locations such as amnion, placental trophoblast, pancreatic islets, ovarian fibroblasts, microglia, Natural Killer(NK) cells, and epithelia. Taken together, our data suggest that innate immune interactions with pathogens markedly altered hominin Siglec biology between 0.6 and 2 Ma, potentially affecting human evolution.

Animal models for SARS-CoV-2 research: A comprehensive literature review

Emerging and re-emerging viral diseases can create devastating effects on human lives and may also lead to economic crises. The ongoing COVID-19 pandemic due to the novel coronavirus (nCoV), severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which originated in Wuhan, China, has caused a global public health emergency. To date, the molecular mechanism of transmission of SARS-CoV-2, its clinical manifestations and pathogenesis is not completely understood. The global scientific community has intensified its efforts in understanding the biology of SARS-CoV-2 for development of vaccines and therapeutic interventions to prevent the rapid spread of the virus and to control mortality and morbidity associated with COVID-19. To understand the pathophysiology of SARS-CoV-2, appropriate animal models that mimic the biology of human SARS-CoV-2 infection are urgently needed. In this review, we outline animal models that have been used to study previous human coronaviruses (HCoVs), including severe acute respiratory syndrome coronavirus (SARS-CoV) and middle east respiratory syndrome coronavirus (MERS-CoV). Importantly, we discuss models that are appropriate for SARS-CoV-2 as well as the advantages and disadvantages of various available methods.

The Faster-Onset Antidepressant Effects of Hypidone Hydrochloride (YL-0919) in Monkeys Subjected to Chronic Unpredictable Stress

Given the limited monkey models of depression available to date, as well as the procedural complexity and time investments that they involve, the ability to test the efficacy and time course of antidepressants in monkey models is greatly restricted. The present study attempted to build a simple and feasible monkey model of depression with chronic unpredictable stress (CUS) and evaluate the antidepressant effect and onset time of fluoxetine hydrochloride (FLX) and the new drug hypidone hydrochloride (YL-0919), a potent and selective 5-HT reuptake inhibitor, 5-HT_1A receptor partial agonist and 5-HT₆ receptor full agonist. Female cynomolgus monkeys with low social status in their colonies were selected and subjected to CUS for 8 weeks by means of food and water deprivation, space restriction, loud noise, strobe light, and intimidation with fake snakes. Huddling, self-clasping, locomotion and environmental exploration were monitored to evaluate behavioral changes. In addition, the window-opening test was used to evaluate the exploratory interest of the monkeys. The present results revealed that CUS-exposed monkeys displayed significant depression-like behaviors, including significant decreases in exploratory interest, locomotion, and exploration as well as significant increases in huddling and self-clasping behavior and the level of fecal cortisol after 8 weeks of CUS. Treatment with FLX (2.4 mg/kg, i. g.) or YL-0919 (1.2 mg/kg, i. g.) markedly reversed the depression-like behaviors caused by CUS, producing significant antidepressant effects. YL-0919 (once daily for 9 days) had a faster-onset antidepressant effect, compared with FLX (once daily for 17 days). In summary, the present study first established a CUS model using female cynomolgus monkeys with low social status and then successfully evaluated the onset time of 5-HTergic antidepressants. The results suggested that monkeys exposed to CUS displayed significant depression-like behaviors, and both FLX and YL-0919 produced antidepressant effects in this model. Moreover, YL-0919 appeared to act faster than FLX. The present study provides a promising prospect for the evaluation of fast-onset antidepressant drugs based on a CUS monkey model.

Comparative genetics of the major histocompatibility complex in humans and nonhuman primates

The major histocompatibility complex (MHC) is one of the most gene-dense regions of the mammalian genome. Multiple genes within the human MHC (HLA) show extensive polymorphism, and currently, more than 26,000 alleles divided over 39 different genes are known. Nonhuman primate (NHP) species are grouped into great and lesser apes and Old and New World monkeys, and their MHC is studied mostly because of their important role as animal models in preclinical research or in connection with conservation biology purposes. The evolutionary equivalents of many of the HLA genes are present in NHP species, and these genes may also show abundant levels of polymorphism. This review is intended to provide a comprehensive comparison relating to the organization and polymorphism of human and NHP MHC regions.

The HLA A03 Supertype and Several Pan Species Major Histocompatibility Complex Class I A Allotypes Share a Preference for Binding Positively Charged Residues in the F Pocket: Implications for Controlling Retroviral Infections

The major histocompatibility complex (MHC) class I region of humans, chimpanzees (Pan troglodytes), and bonobos (Pan paniscus) is highly similar, and orthologues of HLA-A, -B, and -C are present in both Pan species. Based on functional characteristics, the different HLA-A allotypes are classified into different supertypes. One of them, the HLA A03 supertype, is widely distributed among different human populations. All contemporary known chimpanzee and bonobo MHC class I A allotypes cluster genetically into one of the six HLA-A families, the HLA-A1/A3/A11/A30 family. We report here that the peptide-binding motif of the Patr-A*05:01 allotype, which is commonly present in a cohort of western African chimpanzees, has a strong preference for binding peptides with basic amino acids at the carboxyl terminus. This phenomenon is shared with the family members of the HLA A03 supertype. Based on the chemical similarities in the peptide-binding pocket, we inferred that the preference for binding peptides with basic amino acids at the carboxyl terminus is widely present among the human, chimpanzee, and bonobo MHC-A allotypes. Subsequent in silico peptide-binding predictions illustrated that these allotypes have the capacity to target conserved parts of the proteome of human immunodeficiency virus type 1 (HIV-1) and the simian immunodeficiency virus SIVcpz.IMPORTANCE Most experimentally infected chimpanzees seem to control an HIV-1 infection and are therefore considered to be relatively resistant to developing AIDS. Contemporary free-ranging chimpanzees may carry SIVcpz, and there is evidence for AIDS-like symptoms in these free-ranging animals, whereas SIV infections in bonobos appear to be absent. In humans, the natural control of an HIV-1 infection is strongly associated with the presence of particular HLA class I allotypes. The ancestor of the contemporary living chimpanzees and bonobos survived a selective sweep targeting the MHC class I repertoire. We have put forward a hypothesis that this may have been caused by an ancestral retroviral infection similar to SIVcpz. Characterization of the relevant MHC allotypes may contribute to understanding the shaping of their immune repertoire. The abundant presence of MHC-A allotypes that prefer peptides with basic amino acids at the C termini suggests that these molecules may contribute to the control of retroviral infections in humans, chimpanzees, and bonobos.

CRISPR/Cas9-mediated genome editing in nonhuman primates

Owing to their high similarity to humans, non-human primates (NHPs) provide an exceedingly suitable model for the study of human disease. In this Review, we summarize the history of transgenic NHP models and the progress of CRISPR/Cas9-mediated genome editing in NHPs, from the first proof-of-principle green fluorescent protein-expressing monkeys to sophisticated NHP models of human neurodegenerative disease that accurately phenocopy several complex disease features. We discuss not only the breakthroughs and advantages, but also the potential shortcomings of the application of the CRISPR/Cas9 system to NHPs that have emerged from the expanded understanding of this technology in recent years. Although off-target and mosaic mutations are the main concerns in CRISPR/Cas9-mediated NHP modeling, recent progress in genome editing techniques make it likely that these technical limitations will be overcome soon, bringing excellent prospects to human disease studies.

EB 2017 Article: Changes in advanced glycation end products, beta-defensin-3, and interleukin-17 during diabetic periodontitis development in rhesus monkeys

The bidirectional relationship between diabetes mellitus (DM) and periodontal disease has drawn great attention; however, the mechanisms underlying their association remain unclear. In this study, we aimed to develop a rhesus monkey model of diabetic periodontitis and explore the potential mechanisms by which DM affects the progression of periodontal disease. Three healthy rhesus monkeys were selected as the control group. Five streptozotocin-induced diabetic rhesus monkeys were chosen as the experimental group. Ligature placement was used to induce periodontitis. The changes in the levels of advanced glycation end products (AGEs), beta-defensin-3 (BD-3), and interleukin-17 (IL-17) were measured using enzyme-linked immunosorbent assays (ELISA) and real-time reverse transcription polymerase chain reaction (RT-PCR) at different stages during disease progression. Periodontitis was confirmed by clinical assessment, radiographic images, and histological examination. Significant changes in the levels of AGEs and BD-3 in serum were observed at the periodontitis stage in diabetic rhesus monkeys ( P < 0.05). The expression of BD-3 mRNA in the gingiva of diabetic group at baseline was significantly high ( P < 0.05). Diabetic monkeys exhibited significantly enhanced IL-17 mRNA expression at the periodontitis stage ( P < 0.05). Our findings indicated that the rhesus monkey can serve as an ideal model for exploring the pathogenesis of diabetic periodontitis, and the hyperglycemic environment may accelerate inflammatory response and weaken the defense system in periodontal tissues. Impact statement The mechanism underlying the association between diabetes mellitus (DM) and periodontal disease is not yet fully understood. Hence, there is a need to establish animal models to reveal the effect of DM on the pathogenesis of periodontitis. In this study, we explored the appropriate methods for inducing periodontitis and shortening the modeling time in rhesus monkeys, to investigate the pathogenesis of diabetic periodontitis and develop innovative therapies. Our results suggest that a hyperglycemic environment might lead to the destruction of periodontal tissues by accelerating inflammatory response and weakening the defense system in periodontal tissues. Therefore, this study has significant treatment implications regarding the regulation of the immune response against periodontal diseases in patients with DM.

Successful application of human-based methyl capture sequencing for methylome analysis in non-human primate models

BACKGROUND

The characterization of genomic or epigenomic variation in human and animal models could provide important insight into pathophysiological mechanisms of various diseases, and lead to new developments in disease diagnosis and clinical intervention. The African green monkey (AGM; Chlorocebus aethiops) and cynomolgus monkey (CM; Macaca fascicularis) have long been considered important animal models in biomedical research. However, non-human primate-specific methods applicable to epigenomic analyses in AGM and CM are lacking. The recent development of methyl-capture sequencing (MC-seq) has an unprecedented advantage of cost-effectiveness, and further allows for extending the methylome coverage compared to conventional sequencing approaches.

RESULTS

Here, we used a human probe-designed MC-seq method to assay DNA methylation in DNA obtained from 13 CM and three AGM blood samples. To effectively adapt the human probe-designed target region for methylome analysis in non-human primates, we redefined the target regions, focusing on regulatory regions and intragenic regions with consideration of interspecific sequence homology and promoter region variation. Methyl-capture efficiency was controlled by the sequence identity between the captured probes based on the human reference genome and the AGM and CM genome sequences, respectively. Using reasonable guidelines, 56 and 62% of the human-based capture probes could be effectively mapped for DNA methylome profiling in the AGM and CM genome, respectively, according to numeric global statistics. In particular, our method could cover up to 89 and 87% of the regulatory regions of the AGM and CM genome, respectively.

CONCLUSIONS

Use of human-based MC-seq methods provides an attractive, cost-effective approach for the methylome profiling of non-human primates at the single-base resolution level.

Direct estimation of de novo mutation rates in a chimpanzee parent-offspring trio by ultra-deep whole genome sequencing

Mutations generate genetic variation and are a major driving force of evolution. Therefore, examining mutation rates and modes are essential for understanding the genetic basis of the physiology and evolution of organisms. Here, we aim to identify germline de novo mutations through the whole-genome surveyance of Mendelian inheritance error sites (MIEs), those not inherited through the Mendelian inheritance manner from either of the parents, using ultra-deep whole genome sequences (>150-fold) from a chimpanzee parent-offspring trio. We identified such 889 MIEs and classified them into four categories based on the pattern of inheritance and the sequence read depth: [i] de novo single nucleotide variants (SNVs), [ii] copy number neutral inherited variants, [iii] hemizygous deletion inherited variants, and [iv] de novo copy number variants (CNVs). From de novo SNV candidates, we estimated a germline de novo SNV mutation rate as 1.48 × 10^-8 per site per generation or 0.62 × 10^-9 per site per year. In summary, this study demonstrates the significance of ultra-deep whole genome sequencing not only for the direct estimation of mutation rates but also for discerning various mutation modes including de novo allelic conversion and de novo CNVs by identifying MIEs through the transmission of genomes from parents to offspring.

Hox gene cluster of the ascidian, Halocynthia roretzi, reveals multiple ancient steps of cluster disintegration during ascidian evolution

BACKGROUND

Hox gene clusters with at least 13 paralog group (PG) members are common in vertebrate genomes and in that of amphioxus. Ascidians, which belong to the subphylum Tunicata (Urochordata), are phylogenetically positioned between vertebrates and amphioxus, and traditionally divided into two groups: the Pleurogona and the Enterogona. An enterogonan ascidian, Ciona intestinalis (Ci), possesses nine Hox genes localized on two chromosomes; thus, the Hox gene cluster is disintegrated. We investigated the Hox gene cluster of a pleurogonan ascidian, Halocynthia roretzi (Hr) to investigate whether Hox gene cluster disintegration is common among ascidians, and if so, how such disintegration occurred during ascidian or tunicate evolution.

RESULTS

Our phylogenetic analysis reveals that the Hr Hox gene complement comprises nine members, including one with a relatively divergent Hox homeodomain sequence. Eight of nine Hr Hox genes were orthologous to Ci-Hox1, 2, 3, 4, 5, 10, 12 and 13. Following the phylogenetic classification into 13 PGs, we designated Hr Hox genes as Hox1, 2, 3, 4, 5, 10, 11/12/13.a, 11/12/13.b and HoxX. To address the chromosomal arrangement of the nine Hox genes, we performed two-color chromosomal fluorescent in situ hybridization, which revealed that the nine Hox genes are localized on a single chromosome in Hr, distinct from their arrangement in Ci. We further examined the order of the nine Hox genes on the chromosome by chromosome/scaffold walking. This analysis suggested a gene order of Hox1, 11/12/13.b, 11/12/13.a, 10, 5, X, followed by either Hox4, 3, 2 or Hox2, 3, 4 on the chromosome. Based on the present results and those previously reported in Ci, we discuss the establishment of the Hox gene complement and disintegration of Hox gene clusters during the course of ascidian or tunicate evolution.

CONCLUSIONS

The Hox gene cluster and the genome must have experienced extensive reorganization during the course of evolution from the ancestral tunicate to Hr and Ci. Nevertheless, some features are shared in Hox gene components and gene arrangement on the chromosomes, suggesting that Hox gene cluster disintegration in ascidians involved early events common to tunicates as well as later ascidian lineage-specific events.

Unique XCI evolution in Tokudaia: initial XCI of the neo-X chromosome in Tokudaia muenninki and function loss of XIST in Tokudaia osimensis

X chromosome inactivation (XCI) is an essential mechanism to compensate gene dosage in mammals. Here, we show that XCI has evolved differently in two species of the genus Tokudaia. The Amami spiny rat, Tokudaia osimensis, has a single X chromosome in males and females (XO/XO). By contrast, the Okinawa spiny rat, Tokudaia muenninki, has XX/XY sex chromosomes like most mammals, although the X chromosome has acquired a neo-X region by fusion with an autosome. BAC clones containing the XIST gene, which produces the long non-coding RNA XIST required for XCI, were obtained by screening of T. osimensis and T. muenninki BAC libraries. Each clone was mapped to the homologous region of the X inactivation center in the X chromosome of the two species by BAC-FISH. XIST RNAs were expressed in T. muenninki females, whereas no expression was observed in T. osimensis. The sequence of the XIST RNA was compared with that of mouse, showing that the XIST gene is highly conserved in T. muenninki. XIST RNAs were localized to the ancestral X region (Xq), to the heterochromatic region (pericentromeric region), and partially to the neo-X region (Xp). The hybridization pattern correlated with LINE-1 accumulation in Xq but not in Xp. Dosage of genes located on the neo-X chromosome was not compensated, suggesting that the neo-X region is in an early state of XCI. By contrast, many mutations were observed in the XIST gene of T. osimensis, indicating its loss of function in the XO/XO species.

Limited MHC class I intron 2 repertoire variation in bonobos

Common chimpanzees (Pan troglodytes) experienced a selective sweep, probably caused by a SIV-like virus, which targeted their MHC class I repertoire. Based on MHC class I intron 2 data analyses, this selective sweep took place about 2-3 million years ago. As a consequence, common chimpanzees have a skewed MHC class I repertoire that is enriched for allotypes that are able to recognise conserved regions of the SIV proteome. The bonobo (Pan paniscus) shared an ancestor with common chimpanzees approximately 1.5 to 2 million years ago. To investigate whether the signature of this selective sweep is also detectable in bonobos, the MHC class I gene repertoire of two bonobo panels comprising in total 29 animals was investigated by Sanger sequencing. We identified 14 Papa-A, 20 Papa-B and 11 Papa-C alleles, of which eight, five and eight alleles, respectively, have not been reported previously. Within this pool of MHC class I variation, we recovered only 2 Papa-A, 3 Papa-B and 6 Papa-C intron 2 sequences. As compared to humans, bonobos appear to have an even more diminished MHC class I intron 2 lineage repertoire than common chimpanzees. This supports the notion that the selective sweep may have predated the speciation of common chimpanzees and bonobos. The further reduction of the MHC class I intron 2 lineage repertoire observed in bonobos as compared to the common chimpanzee may be explained by a founding effect or other subsequent selective processes.

Hematocrit and Hemoglobin Levels of Nonhuman Apes at Moderate Altitudes: A Comparison with Humans

Mortola, Jacopo P. and DeeAnn Wilfong. Hematocrit and hemoglobin levels of nonhuman apes at moderate altitudes: a comparison with humans. High Alt Med Biol. 17:323-335, 2016.-We asked to what extent the hematologic response (increase in hematocrit [Hct] and in blood hemoglobin concentration [Hb]) of humans to altitude hypoxia was shared by our closest relatives, the nonhuman apes. Data were collected from 29 specimens of 7 species of apes at 2073 m altitude (barometric pressure Pb = 598 mm Hg); additional data originated from apes located at a lower altitude (1493 m, Pb = 639 mm Hg). The human altitude profiles of Hct and Hb between sea level and 3000 m were constructed from a compilation of literature sources that (all combined) comprised data sets of 10,000-12,000 subjects for each gender. These human data were binned for 0-250 m altitude (sea level) and for each 500 m of progressively higher altitudes. Values of Hb and Hct of both men and women were significantly higher than at sea level at the 1500 bin (1250-1750 m); hence, the altitude threshold for the human hematological responses must be between 1000 and 1500 m. In the nonhuman apes, no increase in Hct or Hb was apparent at 1500 m; at 2000 m, the increase was significant only for the Hb of females. At either altitude in the group of nonhuman apes, the increase in Hct was much less than in humans, and that of Hb was significantly less at 1500 m. We conclude that lack of, or minimal, hematopoietic response to moderate altitude can occur in mammalian species that are not genetically adapted to high altitudes. Polycythemia is not a common response to altitude hypoxia and, at least at moderate altitudes, the degree of the human response may represent the exception among apes rather than the rule.

No Distinction of Orthology/Paralogy between Human and Chimpanzee Rh Blood Group Genes

On human (Homo sapiens) chromosome 1, there is a tandem duplication encompassing Rh blood group genes (Hosa_RHD and Hosa_RHCE). This duplication occurred in the common ancestor of humans, chimpanzees (Pan troglodytes), and gorillas, after splitting from their common ancestor with orangutans. Although several studies have been conducted on ape Rh blood group genes, the clear genome structures of the gene clusters remain unknown. Here, we determined the genome structure of the gene cluster of chimpanzee Rh genes by sequencing five BAC (Bacterial Artificial Chromosome) clones derived from chimpanzees. We characterized three complete loci (Patr_RHα, Patr_RHβ, and Patr_RHγ). In the Patr_RHβ locus, a short version of the gene, which lacked the middle part containing exons 4-8, was observed. The Patr_RHα and Patr_RHβ genes were located on the locations corresponding to Hosa_RHD and Hosa_RHCE, respectively, and Patr_RHγ was in the immediate vicinity of Patr_RHβ. Sequence comparisons revealed high sequence similarity between Patr_RHβ and Hosa_RHCE, while the chimpanzee Rh gene closest to Hosa_RHD was not Patr_RHα but rather Patr_RHγ. The results suggest that rearrangements and gene conversions frequently occurred between these genes and that the classic orthology/paralogy dichotomy no longer holds between human and chimpanzee Rh blood group genes.

Resequencing of the common marmoset genome improves genome assemblies and gene-coding sequence analysis

The first draft of the common marmoset (Callithrix jacchus) genome was published by the Marmoset Genome Sequencing and Analysis Consortium. The draft was based on whole-genome shotgun sequencing, and the current assembly version is Callithrix_jacches-3.2.1, but there still exist 187,214 undetermined gap regions and supercontigs and relatively short contigs that are unmapped to chromosomes in the draft genome. We performed resequencing and assembly of the genome of common marmoset by deep sequencing with high-throughput sequencing technology. Several different sequence runs using Illumina sequencing platforms were executed, and 181 Gbp of high-quality bases including mate-pairs with long insert lengths of 3, 8, 20, and 40 Kbp were obtained, that is, approximately 60× coverage. The resequencing significantly improved the MGSAC draft genome sequence. The N50 of the contigs, which is a statistical measure used to evaluate assembly quality, doubled. As a result, 51% of the contigs (total length: 299 Mbp) that were unmapped to chromosomes in the MGSAC draft were merged with chromosomal contigs, and the improved genome sequence helped to detect 5,288 new genes that are homologous to human cDNAs and the gaps in 5,187 transcripts of the Ensembl gene annotations were completely filled.

Initiation of recombination suppression and PAR formation during the early stages of neo-sex chromosome differentiation in the Okinawa spiny rat, Tokudaia muenninki

Background

Sex chromosomes of extant eutherian species are too ancient to reveal the process that initiated sex-chromosome differentiation. By contrast, the neo-sex chromosomes generated by sex-autosome fusions of recent origin in Tokudaia muenninki are expected to be evolutionarily ‘young’, and therefore provide a good model in which to elucidate the early phases of eutherian sex chromosome evolution. Here we describe the genomic evolution of T. muenninki in neo-sex chromosome differentiation.

Results

FISH mapping of a T. muenninki male, using 50 BAC clones as probes, revealed no chromosomal rearrangements between the neo-sex chromosomes. Substitution-direction analysis disclosed that sequence evolution toward GC-richness, which positively correlates with recombination activity, occurred in the peritelomeric regions, but not middle regions of the neo-sex chromosomes. In contrast, the sequence evolution toward AT-richness was observed in those pericentromeric regions. Furthermore, we showed genetic differentiation between the pericentromeric regions as well as an accelerated rate of evolution in the neo-Y region through the detection of male-specific substitutions by gene sequencing in multiple males and females, and each neo-sex–derived BAC sequencing.

Conclusions

Our results suggest that recombination has been suppressed in the pericentromeric region of neo-sex chromosomes without chromosome rearrangement, whereas high levels of recombination activity is limited in the peritelomeric region of almost undifferentiated neo-sex chromosomes. We conclude that PAR might have been formed on the peritelomeric region of sex chromosomes as an independent event from spread of recombination suppression during the early stages of sex chromosome differentiation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0514-y) contains supplementary material, which is available to authorized users.

Establishment of cell lines derived from the genus Macaca through controlled expression of cell cycle regulators

Nonhuman primates are useful animal models for the study of human diseases. However, the number of established cell lines from nonhuman primates is quite limited compared with the number established from other experimental animals. The establishment of nonhuman primate cell lines would allow drug testing on those cell lines before moving experiments into primates. In this study, we established nonhuman primate primary cell lines by introducing the genes for CDK4R24C, cyclin D1, and hTERT. These cell lines proliferated more rapidly than primary cells and bypassed cellular senescence. Karyotype analysis showed that the chromosome patterns were intact in the immortalized cell lines. Furthermore, we showed that the expression of introduced genes could be precisely controlled through the Tet-Off system with the addition of doxycycline. The present study shows that introduction of the CDK4R24C, cyclin D1, and hTERT genes are effective methods of establishing nonhuman primate cell lines.

Molecular evolution of the CYP2D subfamily in primates: purifying selection on substrate recognition sites without the frequent or long-tract gene conversion

The human cytochrome P450 (CYP) 2D6 gene is a member of the CYP2D gene subfamily, along with the CYP2D7P and CYP2D8P pseudogenes. Although the CYP2D6 enzyme has been studied extensively because of its clinical importance, the evolution of the CYP2D subfamily has not yet been fully understood. Therefore, the goal of this study was to reveal the evolutionary process of the human drug metabolic system. Here, we investigate molecular evolution of the CYP2D subfamily in primates by comparing 14 CYP2D sequences from humans to New World monkey genomes. Window analysis and statistical tests revealed that entire genomic sequences of paralogous genes were extensively homogenized by gene conversion during molecular evolution of CYP2D genes in primates. A neighbor-joining tree based on genomic sequences at the nonsubstrate recognition sites showed that CYP2D6 and CYP2D8 genes were clustered together due to gene conversion. In contrast, a phylogenetic tree using amino acid sequences at substrate recognition sites did not cluster the CYP2D6 and CYP2D8 genes, suggesting that the functional constraint on substrate specificity is one of the causes for purifying selection at the substrate recognition sites. Our results suggest that the CYP2D gene subfamily in primates has evolved to maintain the regioselectivity for a substrate hydroxylation activity between individual enzymes, even though extensive gene conversion has occurred across CYP2D coding sequences.

Experimental primates and non-human primate (NHP) models of human diseases in China: current status and progress

Non-human primates (NHPs) are phylogenetically close to humans, with many similarities in terms of physiology, anatomy, immunology, as well as neurology, all of which make them excellent experimental models for biomedical research. Compared with developed countries in America and Europe, China has relatively rich primate resources and has continually aimed to develop NHPs resources. Currently, China is a leading producer and a major supplier of NHPs on the international market. However, there are some deficiencies in feeding and management that have hampered China's growth in NHP research and materials. Nonetheless, China has recently established a number of primate animal models for human diseases and achieved marked scientific progress on infectious diseases, cardiovascular diseases, endocrine diseases, reproductive diseases, neurological diseases, and ophthalmic diseases, etc. Advances in these fields via NHP models will undoubtedly further promote the development of China's life sciences and pharmaceutical industry, and enhance China's position as a leader in NHP research. This review covers the current status of NHPs in China and other areas, highlighting the latest developments in disease models using NHPs, as well as outlining basic problems and proposing effective countermeasures to better utilize NHP resources and further foster NHP research in China.

Mutations in the testis-specific enhancer of SOX9 in the SRY independent sex-determining mechanism in the genus Tokudaia

SRY (sex-determining region Y) is widely conserved in eutherian mammals as a sex-determining gene located on the Y chromosome. SRY proteins bind to the testis-specific enhancer of SOX9 (TES) with SF1 to upregulate SOX9 expression in undifferentiated gonads of XY embryos of humans and mice. The core region within TES, named TESCO, is an important enhancer for mammalian sex determination. We show that TESCO of the genus Tokudaia lost enhancer activity caused by mutations in its SRY and SF1 binding sites. Two species of Tokudaia do not have the Y chromosome or SRY, and one species has multiple SRYs located on the neo-Y chromosome consisting of the Y fused with an autosome. The sequence of Tokudaia TESCO exhibited more than 83% identity with mouse TESCO, however, nucleotide substitution(s) were found in two out of three SRY binding sites and in five out of six SF1 binding sites. TESCO of all species showed low enhancer activity in cells co-transfected with SRY and SF1, and SOX9 and SF1 in reporter gene assays. Mutated TESCO, in which nucleotide substitutions found in SRY and SF1 binding sites were replaced with mouse sequence, recovered the activity. Furthermore, SRYs of the SRY-positive species could not activate the mutated TESCO or mouse TESCO, suggesting that SRYs lost function as a sex-determining gene any more. Our results indicate that the SRY dependent sex-determining mechanism was lost in a common ancestor of the genus Tokudaia caused by nucleotide substitutions in SRY and SF1 binding sites after emergence of a new sex-determining gene. We present the first evidence for an intermediate stage of the switchover from SRY to a new sex-determining gene in the evolution of mammalian sex-determining mechanism.

Endocrinology: advances through omics and related technologies

The rapid development of new omics technologies to measure changes at genetic, transcriptomic, proteomic, and metabolomics levels together with the evolution of methods to analyze and integrate the data at a systems level are revolutionizing the study of biological processes. Here we discuss how new approaches using omics technologies have expanded our knowledge especially in nontraditional models. Our increasing knowledge of these interactions and evolutionary pathway conservation facilitates the use of nontraditional species, both invertebrate and vertebrate, as new model species for biological and endocrinology research. The increasing availability of technology to create organisms overexpressing key genes in endocrine function allows manipulation of complex regulatory networks such as growth hormone (GH) in transgenic fish where disregulation of GH production to produce larger fish has also permitted exploration of the role that GH plays in testis development, suggesting that it does so through interactions with insulin-like growth factors. The availability of omics tools to monitor changes at nearly any level in any organism, manipulate gene expression and behavior, and integrate data across biological levels, provides novel opportunities to explore endocrine function across many species and understand the complex roles that key genes play in different aspects of the endocrine function.

Practical and critical instruction for nonhuman primate diabetic models

Diabetes mellitus, a disease of metabolic dysregulation, is characterized by inappropriate hyperglycemia resulting from progressive loss of insulin secretion or action. The potential of nonhuman primate (NHP) models in diabetes research has been well understood. NHPs have long been regarded as the "gold standard" for preclinical studies. However, there are persistent, severe obstacles to the development and application of these models. At present, a consensus for standardized strategies of diabetic induction has not been achieved. The different modeling methods of diabetes has led to various characterizations of the pathology of the disease; however, there are deficiencies of systemic evaluation programs for nonhuman primate diabetes models. In this scenario, experimental systemic programs provide the highly required guidelines for NHP diabetic models. Moreover, given the expensive and relatively small population of primates and the fatal diabetic complications, it is imperative to carefully manage the care and use of these animals in biomedical research studies. This article briefly reviews the technical and managerial aspects of NHP diabetes models providing practical and critical instruction on housing and care, routine management, development strategy, modeling diagnosis, evaluation, and disease control, as well as guidelines for model selection for various purposes. The present article sought to provide guidelines for NHP models of diabetes in their development and application. It is not intended to outline mandatory requirements for clinical accreditation.

Whole genome comparative analysis of channel catfish (Ictalurus punctatus) with four model fish species

BACKGROUND

Comparative mapping is a powerful tool to study evolution of genomes. It allows transfer of genome information from the well-studied model species to non-model species. Catfish is an economically important aquaculture species in United States. A large amount of genome resources have been developed from catfish including genetic linkage maps, physical maps, BAC end sequences (BES), integrated linkage and physical maps using BES-derived markers, physical map contig-specific sequences, and draft genome sequences. Application of such genome resources should allow comparative analysis at the genome scale with several other model fish species.

RESULTS

In this study, we conducted whole genome comparative analysis between channel catfish and four model fish species with fully sequenced genomes, zebrafish, medaka, stickleback and Tetraodon. A total of 517 Mb draft genome sequences of catfish were anchored to its genetic linkage map, which accounted for 62% of the total draft genome sequences. Based on the location of homologous genes, homologous chromosomes were determined among catfish and the four model fish species. A large number of conserved syntenic blocks were identified. Analysis of the syntenic relationships between catfish and the four model fishes supported that the catfish genome is most similar to the genome of zebrafish.

CONCLUSION

The organization of the catfish genome is similar to that of the four teleost species, zebrafish, medaka, stickleback, and Tetraodon such that homologous chromosomes can be identified. Within each chromosome, extended syntenic blocks were evident, but the conserved syntenies at the chromosome level involve extensive inter-chromosomal and intra-chromosomal rearrangements. This whole genome comparative map should facilitate the whole genome assembly and annotation in catfish, and will be useful for genomic studies of various other fish species.

Derivation of induced pluripotent stem cells from the baboon: a nonhuman primate model for preclinical testing of stem cell therapies

Development of effective pluripotent stem cell-based therapies will require safety and efficacy testing in a clinically relevant preclinical model such as nonhuman primates (NHPs). Baboons and macaques are equally similar to humans genetically and both have been extensively used for biomedical research. Macaques are preferred for human immunodeficiency virus/acquired immunodeficiency syndrome (HIV/AIDS) research whereas baboons are preferred for transplantation studies because of the greater similarity of their anatomy and immunogenetic system to those of humans. We generated four induced pluripotent stem cell (iPSC) lines from skin cells of the olive baboon (Papio anubis). Each line shows the distinct morphology of primate pluripotent stem cells, including flat colonies with well-defined borders and a high nuclear/cytoplasm ratio. Each is positive for the pluripotency markers OCT4, SOX2, NANOG, and SSEA4. Pluripotency was confirmed in two lines by teratoma formation with representative tissues from each germ layer, whereas a third produced cells from all three germ layers following embryoid body differentiation. Three lines have a normal male karyotype and the fourth is missing the short arm of one copy of chromosome 18. This may serve as an in vitro model for the human developmental disorder 18p-, which impacts 1 in 50,000 births/year. These iPSC lines represent the first step toward establishing the baboon as a NHP model for developing stem cell-based therapies.

Identification of an intact ParaHox cluster with temporal colinearity but altered spatial colinearity in the hemichordate Ptychodera flava

BACKGROUND

ParaHox and Hox genes are thought to have evolved from a common ancestral ProtoHox cluster or from tandem duplication prior to the divergence of cnidarians and bilaterians. Similar to Hox clusters, chordate ParaHox genes including Gsx, Xlox, and Cdx, are clustered and their expression exhibits temporal and spatial colinearity. In non-chordate animals, however, studies on the genomic organization of ParaHox genes are limited to only a few animal taxa. Hemichordates, such as the Enteropneust acorn worms, have been used to gain insights into the origins of chordate characters. In this study, we investigated the genomic organization and expression of ParaHox genes in the indirect developing hemichordate acorn worm Ptychodera flava.

RESULTS

We found that P. flava contains an intact ParaHox cluster with a similar arrangement to that of chordates. The temporal expression order of the P. flava ParaHox genes is the same as that of the chordate ParaHox genes. During embryogenesis, the spatial expression pattern of PfCdx in the posterior endoderm represents a conserved feature similar to the expression of its orthologs in other animals. On the other hand, PfXlox and PfGsx show a novel expression pattern in the blastopore. Nevertheless, during metamorphosis, PfXlox and PfCdx are expressed in the endoderm in a spatially staggered pattern similar to the situation in chordates.

CONCLUSIONS

Our study shows that P. flava ParaHox genes, despite forming an intact cluster, exhibit temporal colinearity but lose spatial colinearity during embryogenesis. During metamorphosis, partial spatial colinearity is retained in the transforming larva. These results strongly suggest that intact ParaHox gene clustering was retained in the deuterostome ancestor and is correlated with temporal colinearity.

Comparative genomic analysis of catfish linkage group 8 reveals two homologous chromosomes in zebrafish and other teleosts with extensive inter-chromosomal rearrangements

Background

Comparative genomics is a powerful tool to transfer genomic information from model species to related non-model species. Channel catfish (Ictalurus punctatus) is the primary aquaculture species in the United States. Its existing genome resources such as genomic sequences generated from next generation sequencing, BAC end sequences (BES), physical maps, linkage maps, and integrated linkage and physical maps using BES-associated markers provide a platform for comparative genomic analysis between catfish and other model teleost fish species. This study aimed to gain understanding of genome organizations and similarities among catfish and several sequenced teleost genomes using linkage group 8 (LG8) as a pilot study.

Results

With existing genome resources, 287 unique genes were identified in LG8. Comparative genome analysis indicated that most of these 287 genes on catfish LG8 are located on two homologous chromosomes of zebrafish, medaka, stickleback, and three chromosomes of green-spotted pufferfish. Large numbers of conserved syntenies were identified. Detailed analysis of the conserved syntenies in relation to chromosome level similarities revealed extensive inter-chromosomal and intra-chromosomal rearrangements during evolution. Of the 287 genes, 35 genes were found to be duplicated in the catfish genome, with the vast majority of the duplications being interchromosomal.

Conclusions

Comparative genome analysis is a powerful tool even in the absence of a well-assembled whole genome sequence. In spite of sequence stacking due to low resolution of the linkage and physical maps, conserved syntenies can be identified although the exact gene order and orientation are unknown at present. Through chromosome-level comparative analysis, homologous chromosomes among teleosts can be identified. Syntenic analysis should facilitate annotation of the catfish genome, which in turn, should facilitate functional inference of genes based on their orthology.

The HIV-1 pandemic: does the selective sweep in chimpanzees mirror humankind's future?

An HIV-1 infection progresses in most human individuals sooner or later into AIDS, a devastating disease that kills more than a million people worldwide on an annual basis. Nonetheless, certain HIV-1-infected persons appear to act as long-term non-progressors, and elite control is associated with the presence of particular MHC class I allotypes such as HLA-B*27 or -B*57. The HIV-1 pandemic in humans arose from the cross-species transmission of SIV_cpz originating from chimpanzees. Chimpanzees, however, appear to be relatively resistant to developing AIDS after HIV-1/SIV_cpz infection. Mounting evidence illustrates that, in the distant past, chimpanzees experienced a selective sweep resulting in a severe reduction of their MHC class I repertoire. This was most likely caused by an HIV-1/SIV-like retrovirus, suggesting that chimpanzees may have experienced long-lasting host-virus relationships with SIV-like viruses. Hence, if natural selection is allowed to follow its course, prospects for the human population may look grim, thus underscoring the desperate need for an effective vaccine.

Evaluation of IL-28B polymorphisms and serum IP-10 in hepatitis C infected chimpanzees

In humans, clearance of hepatitis C virus (HCV) infection is associated with genetic variation near the IL-28B gene and the induction of interferon-stimulated genes, like IP-10. Also in chimpanzees spontaneous clearance of HCV is observed. To study whether similar correlations exist in these animals, a direct comparison of IP-10 and IL-28B polymorphism between chimpanzees and patients was performed. All chimpanzees studied were monomorphic for the human IL-28B SNPs which are associated with spontaneous and treatment induced HCV clearance in humans. As a result, these particular SNPs cannot be used for clinical association studies in chimpanzees. Although these human SNPs were absent in chimpanzees, gene variation in this region was present however, no correlation was observed between different SNP-genotypes and HCV outcome. Strikingly, IP-10 levels in chimpanzees correlated with HCV-RNA load and γGT, while such correlations were not observed in humans. The correlation between IP-10, γGT and virus load in chimpanzees was not found in patients and may be due to the lack of lifestyle-related confounding factors in chimpanzees. Direct comparison of IP-10 and IL-28B polymorphism between chimpanzees and patients in relation to HCV infection, illustrates that the IFN-pathways are important during HCV infection in both species. The Genbank EMBL accession numbers assigned to chimpanzees specific sequences near the IL-28B gene are HE599784 and HE599785.

Genome physical mapping of polyploids: a BIBAC physical map of cultivated tetraploid cotton, Gossypium hirsutum L

Polyploids account for approximately 70% of flowering plants, including many field, horticulture and forage crops. Cottons are a world-leading fiber and important oilseed crop, and a model species for study of plant polyploidization, cellulose biosynthesis and cell wall biogenesis. This study has addressed the concerns of physical mapping of polyploids with BACs and/or BIBACs by constructing a physical map of the tetraploid cotton, Gossypium hirsutum L. The physical map consists of 3,450 BIBAC contigs with an N50 contig size of 863 kb, collectively spanning 2,244 Mb. We sorted the map contigs according to their origin of subgenome, showing that we assembled physical maps for the A- and D-subgenomes of the tetraploid cotton, separately. We also identified the BIBACs in the map minimal tilling path, which consists of 15,277 clones. Moreover, we have marked the physical map with nearly 10,000 BIBAC ends (BESs), making one BES in approximately 250 kb. This physical map provides a line of evidence and a strategy for physical mapping of polyploids, and a platform for advanced research of the tetraploid cotton genome, particularly fine mapping and cloning the cotton agronomic genes and QTLs, and sequencing and assembling the cotton genome using the modern next-generation sequencing technology.

A comparative physical map reveals the pattern of chromosomal evolution between the turkey (Meleagris gallopavo) and chicken (Gallus gallus) genomes

Background

A robust bacterial artificial chromosome (BAC)-based physical map is essential for many aspects of genomics research, including an understanding of chromosome evolution, high-resolution genome mapping, marker-assisted breeding, positional cloning of genes, and quantitative trait analysis. To facilitate turkey genetics research and better understand avian genome evolution, a BAC-based integrated physical, genetic, and comparative map was developed for this important agricultural species.

Results

The turkey genome physical map was constructed based on 74,013 BAC fingerprints (11.9 × coverage) from two independent libraries, and it was integrated with the turkey genetic map and chicken genome sequence using over 41,400 BAC assignments identified by 3,499 overgo hybridization probes along with > 43,000 BAC end sequences. The physical-comparative map consists of 74 BAC contigs, with an average contig size of 13.6 Mb. All but four of the turkey chromosomes were spanned on this map by three or fewer contigs, with 14 chromosomes spanned by a single contig and nine chromosomes spanned by two contigs. This map predicts 20 to 27 major rearrangements distinguishing turkey and chicken chromosomes, despite up to 40 million years of separate evolution between the two species. These data elucidate the chromosomal evolutionary pattern within the Phasianidae that led to the modern turkey and chicken karyotypes. The predominant rearrangement mode involves intra-chromosomal inversions, and there is a clear bias for these to result in centromere locations at or near telomeres in turkey chromosomes, in comparison to interstitial centromeres in the orthologous chicken chromosomes.

Conclusion

The BAC-based turkey-chicken comparative map provides novel insights into the evolution of avian genomes, a framework for assembly of turkey whole genome shotgun sequencing data, and tools for enhanced genetic improvement of these important agricultural and model species.

Pancreas anatomy and surgical procedure for pancreatectomy in rhesus monkeys

BACKGROUND

The aim of this study was to investigate the pancreas anatomy and surgical procedure for harvesting pancreas for islet isolation while performing pancreatectomy to induce diabetes in rhesus monkeys.

METHODS

The necropsy was performed in three cadaveric monkeys. Two monkeys underwent the total pancreatectomy and four underwent partial pancreatectomy (70-75%).

RESULTS

The greater omentum without ligament to transverse colon, the cystic artery arising from the proper hepatic artery and the branches supplying the paries posterior gastricus from the splenic artery were observed. For pancreatectomy, resected pancreas can be used for islet isolation. Diabetes was not induced in the monkeys undergoing partial pancreatectomy (70-75%).

CONCLUSIONS

Pancreas anatomy in rhesus monkeys is not the same as in human. Diabetes can be induced in rhesus monkeys by total but not partial pancreatectomy (70-75%). Resected pancreas can be used for islet isolation while performing pancreatectomy to induce diabetes.

Genome sequence of an Australian kangaroo, Macropus eugenii, provides insight into the evolution of mammalian reproduction and development

BACKGROUND

We present the genome sequence of the tammar wallaby, Macropus eugenii, which is a member of the kangaroo family and the first representative of the iconic hopping mammals that symbolize Australia to be sequenced. The tammar has many unusual biological characteristics, including the longest period of embryonic diapause of any mammal, extremely synchronized seasonal breeding and prolonged and sophisticated lactation within a well-defined pouch. Like other marsupials, it gives birth to highly altricial young, and has a small number of very large chromosomes, making it a valuable model for genomics, reproduction and development.

RESULTS

The genome has been sequenced to 2 × coverage using Sanger sequencing, enhanced with additional next generation sequencing and the integration of extensive physical and linkage maps to build the genome assembly. We also sequenced the tammar transcriptome across many tissues and developmental time points. Our analyses of these data shed light on mammalian reproduction, development and genome evolution: there is innovation in reproductive and lactational genes, rapid evolution of germ cell genes, and incomplete, locus-specific X inactivation. We also observe novel retrotransposons and a highly rearranged major histocompatibility complex, with many class I genes located outside the complex. Novel microRNAs in the tammar HOX clusters uncover new potential mammalian HOX regulatory elements.

CONCLUSIONS

Analyses of these resources enhance our understanding of marsupial gene evolution, identify marsupial-specific conserved non-coding elements and critical genes across a range of biological systems, including reproduction, development and immunity, and provide new insight into marsupial and mammalian biology and genome evolution.

Using the Acropora digitifera genome to understand coral responses to environmental change

Despite the enormous ecological and economic importance of coral reefs, the keystone organisms in their establishment, the scleractinian corals, increasingly face a range of anthropogenic challenges including ocean acidification and seawater temperature rise. To understand better the molecular mechanisms underlying coral biology, here we decoded the approximately 420-megabase genome of Acropora digitifera using next-generation sequencing technology. This genome contains approximately 23,700 gene models. Molecular phylogenetics indicate that the coral and the sea anemone Nematostella vectensis diverged approximately 500 million years ago, considerably earlier than the time over which modern corals are represented in the fossil record (∼240 million years ago). Despite the long evolutionary history of the endosymbiosis, no evidence was found for horizontal transfer of genes from symbiont to host. However, unlike several other corals, Acropora seems to lack an enzyme essential for cysteine biosynthesis, implying dependency of this coral on its symbionts for this amino acid. Corals inhabit environments where they are frequently exposed to high levels of solar radiation, and analysis of the Acropora genome data indicates that the coral host can independently carry out de novo synthesis of mycosporine-like amino acids, which are potent ultraviolet-protective compounds. In addition, the coral innate immunity repertoire is notably more complex than that of the sea anemone, indicating that some of these genes may have roles in symbiosis or coloniality. A number of genes with putative roles in calcification were identified, and several of these are restricted to corals. The coral genome provides a platform for understanding the molecular basis of symbiosis and responses to environmental changes.

Positional information resolves structural variations and uncovers an evolutionarily divergent genetic locus in accessions of Arabidopsis thaliana

Genome sequencing of closely related individuals has yielded valuable insights that link genome evolution to phenotypic variations. However, advancement in sequencing technology has also led to an escalation in the number of poor quality–drafted genomes assembled based on reference genomes that can have highly divergent or haplotypic regions. The self-fertilizing nature of Arabidopsis thaliana poses an advantage to sequencing projects because its genome is mostly homozygous. To determine the accuracy of an Arabidopsis drafted genome in less conserved regions, we performed a resequencing experiment on a ∼371-kb genomic interval in the Landsberg erecta (Ler-0) accession. We identified novel structural variations (SVs) between Ler-0 and the reference accession Col-0 using a long-range polymerase chain reaction approach to generate an Illumina data set that has positional information, that is, a data set with reads that map to a known location. Positional information is important for accurate genome assembly and the resolution of SVs particularly in highly duplicated or repetitive regions. Sixty-one regions with misassembly signatures were identified from the Ler-0 draft, suggesting the presence of novel SVs that are not represented in the draft sequence. Sixty of those were resolved by iterative mapping using our data set. Fifteen large indels (>100 bp) identified from this study were found to be located either within protein-coding regions or upstream regulatory regions, suggesting the formation of novel alleles or altered regulation of existing genes in Ler-0. We propose future genome-sequencing experiments to follow a clone-based approach that incorporates positional information to ultimately reveal haplotype-specific differences between accessions.

Treatment and risk factor analysis of hypoglycemia in diabetic rhesus monkeys

In order to anticipate and promptly treat hypoglycemia in diabetic monkeys treated with insulin or other glucose-lowering drugs, the relationships between the incidence and symptoms of hypoglycemia in these animals, and many factors involved in model development and sustainment were analyzed. Different procedures were performed on 22 monkeys for the induction of diabetes. The monkey models were evaluated by blood glucose, insulin, C-peptide levels and intravenous glucose tolerance tests. A glucose treatment program for the diabetic monkeys was administered and laboratory tests were regularly performed. A standard procedure of hypoglycemia treatment was established and the risk factors of hypoglycemia were analyzed by a logistic regression model. Furthermore, the relationships between the four methods of diabetes induction, renal function, glycemic control and hypoglycemia were studied using one-way analysis of variance and t-test. We found that the hypoglycemic conditions of diabetic monkeys were improved rapidly by our treatment. The statistical analysis suggested that the modeling methods, renal function and glycemic control were related to the incidence of hypoglycemia. In detail, the progress of diabetes, effects of glycemic control and, particularly, the severity of the hypoglycemia differed according to the induction strategy used. The models induced by partial pancreatectomy with low-dose streptozotocin were not prone to hypoglycemia and their glycemic controls were stable. However, the models induced by total pancreatectomy were more vulnerable to severe hypoglycemia and their glycemic controls were the most unstable. Moreover, the levels of blood creatinine and triglyceride increased after the development of diabetes, which was related to the occurrence of hypoglycemia. In conclusion, we suggested that total pancreatectomy and renal impairment are two important risk factors for hypoglycemia in diabetic monkeys. More attention should be paid to daily care of diabetic monkeys, particularly monitoring and protecting their renal function.

Divergence, demography and gene loss along the human lineage

Genomic DNA sequences are an irreplaceable source for reconstructing the vanished past of living organisms. Based on updated sequence data, this paper summarizes our studies on species divergence time, ancient population size and functional loss of genes in the primate lineage leading to modern humans (Homo sapiens sapiens). The inter- and intraspecific comparisons of DNA sequences suggest that the human lineage experienced a rather severe bottleneck in the Middle Pleistocene, throughout which period the subdivided African population played a predominant role in shaping the genetic architecture of modern humans. Also, published and newly identified human-specific pseudogenes (HSPs) are enumerated in order to infer their significance for human evolution. Of the 121 candidate genes obtained, authentic HSPs turn out to comprise only 25 olfactory receptor genes, four T cell receptor genes and nine other genes. The fixation of HSPs has been too rare over the past 6–7 Myr to account for species differences between humans and chimpanzees.

Diversification of bitter taste receptor gene family in western chimpanzees

In mammals, bitter taste is mediated by T2R genes, which belong to the large family of seven transmembrane G protein-coupled receptors. Because T2Rs are directly involved in the interaction between mammals and their dietary sources, it is likely that these genes evolved to reflect species' specific diets during mammalian evolution. Here, we investigated the sequences of all 28 putative functional chimpanzee T2R genes (cT2Rs) in 46 western chimpanzees to compare the intraspecies variations in chimpanzees to those already known for all 25 human functional T2R genes (hT2Rs). The numbers of functional genes varied among individuals in western chimpanzees, and most chimpanzees had two or three more functional genes than humans. Similarly to hT2Rs, cT2Rs showed high nucleotide diversity along with a large number of amino acid substitutions. Comparison of the nucleotide substitution patterns in cT2Rs with those in five cT2R pseudogenes and 14 autosomal intergenic noncoding regions among the same individuals revealed that the evolution of cT2R genes was almost identical to that of putative neutral regions with slight but significantly positive Tajima's D values, suggesting that selective constraint on these genes was relaxed with weak balancing selection. These trends have resulted in the occurrence of various divergent alleles of T2Rs within the western chimpanzee populations and in heterozygous individuals who might have the ability to taste a broader range of substances.

A Comparative BAC map for the gilthead sea bream (Sparus aurata L.)

This study presents the first comparative BAC map of the gilthead sea bream (Sparus aurata), a highly valuated marine aquaculture fish species in the Mediterranean. High-throughput end sequencing of a BAC library yielded 92,468 reads (60.6 Mbp). Comparative mapping was achieved by anchoring BAC end sequences to the three-spined stickleback (Gasterosteus aculeatus) genome. BACs that were consistently ordered along the stickleback chromosomes accounted for 14,265 clones. A fraction of 5,249 BACs constituted a minimal tiling path that covers 73.5% of the stickleback chromosomes and 70.2% of the genes that have been annotated. The N50 size of 1,485 “BACtigs” consisting of redundant BACs is 337,253 bp. The largest BACtig covers 2.15 Mbp in the stickleback genome. According to the insert size distribution of mapped BACs the sea bream genome is 1.71-fold larger than the stickleback genome. These results represent a valuable tool to researchers in the field and may support future projects to elucidate the whole sea bream genome.

Comparative genomics in teleost species: Knowledge transfer by linking the genomes of model and non-model fish species

Comparative genomics is a powerful tool to transfer knowledge coming from model fish species to non-model fish species of economic or/and evolutionary interest. Such transfer is of importance as functional studies either are difficult to perform with most non-model species. The first comparative map constructed using the human and the chimpanzee genome allowed the identification of putative orthologues. Although comparative mapping in teleosts is still in its infancy, five model teleost genomes from different orders have been fully sequenced to date and the sequencing of several commercially important species are also underway or near completion. The accessibility of these whole genome sequences and rapid developments in genomics of fish species are paving the way towards new and valuable research in comparative genetics and genomics. With the accumulation of information in model species, the genetic and genomic characterization of non-model, but economically, physiologically or evolutionary important species is now feasible. Furthermore, comparison of low coverage gene maps of non-model fish species against fully sequenced fish species will enhance the efficiency of candidate gene identification projected for quantitative trait loci (QTL) scans for traits of special interest.