Search

The application of current genome-wide sequencing techniques on human populations helps elucidate the considerable gene flow among genusHomo, which includes modern and archaic humans. Gene flow among current human populations has been studied using frequencies of single nucleotide polymorphisms. Unlike single nucleotide polymorphism frequency data, haplotype data are suitable for identifying and tracing rare evolutionary events. Haplotype data can also conveniently detect genomic location and estimate molecular function that may be a target of selection. We analyzed eight loci of the human genome using the same procedure for each locus to infer human haplotype diversity and reevaluate past explanations of the evolutionary mechanisms that affected these loci. These loci have been recognized by separate studies because of their unusual gene genealogy and geographic distributions that are inconsistent with the recent out-of-Africa model. For each locus, we constructed genealogies for haplotypes using sequence data of the 1000 Genomes Project. Then, we performed S* analysis to estimate distinct gene flow events other than out-of-Africa events. Furthermore, we also estimated unevenness of selective pressure between haplotypes by Extended Haplotype Homozygosity analysis. Based on the patterns of results obtained by this combination of analyses, we classified the examined loci without using a specific population model. This simple method helped clarify evolutionary events for each locus, including rare evolutionary events such as introgression, incomplete lineage sorting, selection, and haplotype recombination that may be hard to discriminate from each other.

Felsenstein's PHYLIP package of molecular phylogeny tools has been used globally since 1980. The programs are receiving renewed attention because of their character-based user interface, which has the advantage of being scriptable for use with large-scale data studies based on super-computers or massively parallel computing clusters. However, occasionally we found, the PHYLIP Consense program output text file displays two or more divided bootstrap values for the same cluster in its result table, and when this happens the output Newick tree file incorrectly assigns only the last value to that cluster that disturbs correct estimation of a consensus tree. We ascertained the cause of this aberrant behavior in the bootstrapping calculation. Our rewrite of the Consense program source code outputs bootstrap values, without redundancy, in its result table, and a Newick tree file with appropriate, corresponding bootstrap values. Furthermore, we developed an add-on program and shell script, add_bootstrap.pl and fasta2tre_bs.bsh, to generate a Newick tree containing the topology and branch lengths inferred from the original data along with valid bootstrap values, and to actualize the automated inference of a phylogenetic tree containing the originally inferred topology and branch lengths with bootstrap values, from multiple unaligned sequences, respectively. These programs can be downloaded at: https://github.com/ShimadaMK/PHYLIP_enhance/.

According to the length distribution of human introns, there is a large population of short introns with a threshold of 65 nucleotides (nt) and a peak at 85 nt. Using human genome and transcriptome databases, we investigated the introns shorter than 66 nt, termed ultra-short introns, the identities of which are scarcely known. Here, we provide for the first time a list of bona fide human ultra-short introns, which have never been characterized elsewhere. By conducting BLAST searches of the databases, we screened 22 introns (37–65 nt) with conserved lengths and sequences among closely related species. We then provide experimental and bioinformatic evidence for the splicing of 15 introns, of which 12 introns were remarkably G-rich and 9 introns contained completely inefficient splice sites and/or branch sites. These unorthodox characteristics of ultra-short introns suggest that there are unknown splicing mechanisms that differ from the well-established mechanism.

Maternal kinship is important in primate societies because it affects individual behaviour as well as the sustainability of populations. All members of the Bossou chimpanzee community are descended from 8 individuals (herein referred to as original adults) who were already adults or subadults when field observations were initiated in 1976 and whose genetic relationships were unknown. Sequencing of the control region on the maternally inherited mtDNA revealed that 4 (1 male and 3 females) of the 8 original adults shared an identical haplotype. We investigated the effects of the skewed distribution of mtDNA haplotypes on the following two outcomes. First, we demonstrated that the probability of mtDNA haplotype extinction would be increased under such a skewed composition in a small community. Second, the ratio of potential mating candidates to competitors is likely to decrease if chimpanzees become aware of maternal kinship and avoid incest. We estimated that the magnitude of the decrease in the ratio is 10 times greater in males than in females. Here we demonstrate a scenario in which this matrilineal skewness in a small community accelerates extinction of mtDNA haplotype, which will make it more difficult to find a suitable mate within the community.

Short Tandem Repeats (STRs) comprise repeats of one to several base pairs. Because of the high mutability due to strand slippage during DNA synthesis, rapid evolutionary change in the number of repeating units directly shapes the range of repeat-number variation according to selection pressure. However, the remaining questions include: Why are STRs causing repeat expansion diseases maintained in the human population; and why are these limited to neurodegenerative diseases? By evaluating the genome-wide selection pressure on STRs using the database we constructed, we identified two different patterns of relationship in repeat-number polymorphisms between DNA and amino-acid sequences, although both patterns are evolutionary consequences of avoiding the formation of harmful long STRs. First, a mixture of degenerate codons is represented in poly-proline (poly-P) repeats. Second, long poly-glutamine (poly-Q) repeats are favored at the protein level; however, at the DNA level, STRs encoding long poly-Qs are frequently divided by synonymous SNPs. Furthermore, significant enrichments of apoptosis and neurodevelopment were biological processes found specifically in genes encoding poly-Qs with repeat polymorphism. This suggests the existence of a specific molecular function for polymorphic and/or long poly-Q stretches. Given that the poly-Qs causing expansion diseases were longer than other poly-Qs, even in healthy subjects, our results indicate that the evolutionary benefits of long and/or polymorphic poly-Q stretches outweigh the risks of long CAG repeats predisposing to pathological hyper-expansions. Molecular pathways in neurodevelopment requiring long and polymorphic poly-Q stretches may provide a clue to understanding why poly-Q expansion diseases are limited to neurodegenerative diseases.

The population genetic history of a 10.1-kbp noncoding region of the human X chromosome was studied using the males of the HGDP-CEPH Human Genome Diversity Panel (672 individuals from 52 populations). The geographic distribution of patterns of variation was roughly consistent with previous studies, with the major exception that 1 highly divergent haplotype (haplotype X, hX) was observed at low frequency in widely scattered non-African populations and not at all observed in sub-Saharan African populations. Microsatellite (short tandem repeat) variation within the sequenced region was low among copies of hX, even though the estimated time of ancestry of hX and other sequences was 1.44 Myr. The estimated age of the common ancestor of all hX copies was 5,230 years (95% consistency index: 2,000-75,480 years). To further address the presence of hX in Africa, additional samples from Chad and Tanzania were screened. Five additional copies of hX were observed, consistent with a history in which hX was present in Africa prior to the migration of modern humans out of Africa and with eastern Africa being the source of non-African modern human populations. Taken together, these features of hX-that it is much older than other haplotypes and uncommon and patchily distributed throughout Africa, Europe, and Asia-present a cautionary tale for interpretations of human history.

Chromosome rearrangement has been considered to be important in the evolutionary process. Here, we demonstrate the evolutionary relationship of the rearranged human chromosome 12 and the corresponding chromosome XII in apes (chimpanzee, bonobo, gorilla, orangutan, and gibbon) by examining PCR products derived from the breakpoints of inversions and by conducting shotgun sequencing of a gorilla fosmid clone containing the breakpoint and a “duplicated segment” (duplicon). We confirmed that a pair of 23-kb duplicons flank the breakpoints of inversions on the long and short arms of chimpanzee chromosome XII. Although only the 23-kb duplicon on the long arm of chimpanzee chromosome XII and its telomeric flanking sequence are found to be conserved among the hominoids (human, great apes, and gibbons), the duplicon on the short arm of chimpanzee chromosome XII is suggested to be the result of a duplication from that on the long arm. Furthermore, the shotgun sequencing of a gorilla fosmid indicated that the breakpoint on the long arm of the gorilla is located at a different position 1.9 kb from that of chimpanzee. The region is flanked by a sequence homologous to that of human chromosome 6q22. Our findings and sequence analysis suggest a close relationship between segmental duplication and chromosome rearrangement (or breakpoint of inversion) in Hominoidea. The role of the chromosome rearrangement in speciation is also discussed based on our new results.

Various canine breeds are remarkably different from each other not only in their sizes and shapes but also in behavioral traits, suggesting that some of them are under genetic control. Although dopaminergic neurotransmission system is considered to affect animal behavior, little is known about related genes in canine. Relations between specific alleles in polymorphic regions of the dopamine receptor D4 gene (DRD4) and personality or psychiatric disorders have been reported in humans, and we first found polymorphism in exon III region of the gene in 4 canine breeds. In this study we surveyed allele frequency distribution in 23 breeds including a total of 1,535 unrelated individuals. In exon III, 8 alleles including a novel allele were identified. A group of breeds in which the alleles 447b, 498 and 549 were frequent tended toward high scores in aggression-related behavioral traits than that with frequent alleles 435 and 447a. Moreover, a polymorphism based on 24 bp insertion/deletion was found in exon I region for the first time in dogs. This information may be of use for candidate gene studies of behavioral variation in dogs.

Simian T-cell leukemia viruses (STLVs) are the simian counterparts of human T-cell leukemia viruses (HTLVs). A novel, divergent type of STLV (STLV-L) from captive baboons was reported in 1994, but its natural prevalence remained unclear. We investigated the prevalence of STLV-L in 519 blood samples from wild-living nonhuman primates in Ethiopia. Seropositive monkeys having cross-reactive antibodies against HTLV were found among 22 out of 40 hamadryas baboons, 8 of 96 anubis baboons, 24 of 50 baboons that are hybrids between hamadryas and anubis baboons, and 41 of 177 grivet monkeys, but not in 156 gelada baboons. A Western blotting assay showed that sera obtained from seropositive hamadryas and hybrid baboons exhibited STLV-L-like reactivity. A PCR assay successfully amplified STLV sequences, which were subsequently sequenced and confirmed as being closely related to STLV-L. Surprisingly, further PCR showed that nearly half of the hamadryas (20 out of 40) and hybrid (19 out of 50) baboons had STLV-L DNA sequences. In contrast, most of the seropositive anubis baboons and grivet monkeys carried typical STLV-1 but not STLV-L. These observations demonstrate that STLV-L naturally prevails among hamadryas and hybrid baboons at significantly high rates. STLV-1 and -2, the close relative of STLV-L, are believed to have jumped across simian-human barriers, which resulted in widespread infection of HTLV-1 and -2. Further studies are required to know if STLV-L is spreading into human populations.

Cercopithecus aethiops can be classified into four subspecies by morphology and by geographic distribution. However, the phylogenetic relationship between these subspecies is unclear. We previously found five distinct haplogroups of mitochondrial DNA (mtDNA) in the subspecies C. aethiops aethiops at the restriction fragment length polymorphism (RFLP) level, and found that those haplogroups are parapatrically distributed in their habitat. To determine the relationship between subspeciation and haplogroup formation in a subspecies, we compared mtDNA control region and 12S rRNA gene sequences (approximately 700 bp) in C. a. aethiops, two other subspecies of C. aethiops, and two species of Cercopithecus. The diversity between haplogroups in C. a. aethiops was almost the same as that between subspecies. This similar level of diversification between and within haplogroups may explain why a previously obtained mtDNA tree did not show monophyletic branching according to subspecies.

Eight Argonne premium coals and an Australian brown coal (Morwell) were extracted in a flowing stream of either tetralin or 1-methylnaphthalene under 10 MPa at 200–400°C to examine the dissolution/depolymerization behavior of coal in liquid phase. The flowing solvent was used to minimize the secondary interaction between the extract and the macromolecular network of coal. Four bituminous coals were extracted with negligible decomposition below 350°C. The extract yield reached 65–80% for the bituminous coals at 350°C, and the extract was separated into about 25–40% of soluble fraction at room temperature (soluble) and about 40% of solid precipitated at room temperature (deposit). The detailed analyses of the soluble, the deposit and the residue through ultimate analysis, 13 C solid NMR, GPC, and matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOFMS) enabled to examine the structure of Pittsburgh #8 (PITT) in detail. The raw PITT coal was found to contain small molecules ( M W =130–600) and large molecules ( M W =1300–3500) in its network. The soluble consisted of only small molecules, but both small and large molecules were involved in the deposit and the residue. The extraction method employed in this work was found to be effective to separate coal into different molecule size fractions without decomposition for bituminous coals. The soluble, being obtained in 20–40% yield and consisting of molecules smaller than 600 of M W , and the deposit, being obtained in 40% yield at maximum and containing little ash, would be well utilized as either clean fuel or raw materials for the subsequent conversion.

To clarify the relation between organic sulfur distribution in coal and the emission of sulfur-containing gases, six Japanese standard coals, three Argonne premium coals, and one Chinese coal were pyrolyzed using a thermogravimetric analyzer-mass spectrometer (TG-MS) and a Curie-point pyrolyzer. The changes in sulfur distributions in coal and chars were estimated by a modified controlled-atmosphere-programmed-temperature-oxidation (CAPTO) method. It was clarified that the aliphatic sulfur decomposed below 500 °C and that the aromatic sulfur decomposed at 400−700 °C, irrespective of coal type. The decomposition of thiophenic sulfur was strongly affected by coal type. Finally, the H2S formation rate and the decomposition rate of organic sulfur were analyzed using a new distributed activation energy model (DAEM).

Since mitochondrial DNA (mtDNA) are maternally inherited without recombination, geographic distribution of mtDNA in semiterrestrial cercopithecines is considered to be influenced by female philopatry. I examined the effect of sex difference in migration patterns on geographic distribution in a habitat whose environment has changed frequently. I investigated ten groups (n = 77) of grivets (Cercopithecus aethiops aethiops) along a 600-km stretch of the Awash River, Ethiopia. I examined the mtDNA distribution among natural local populations whose nuclear variation was already shown to have a widely homogeneous distribution. RFLP analysis of whole mtDNA genome using 17 enzymes identified ten haplotypes in five clusters (haplogroups). Sequence divergence within haplogroups ranged from 0.17%–0.38%, while divergence between haplogroups ranged between 1.0%–2.5%. Haplogroups were distributed in blocks which ranged from 120–250 km along the Awash River. The haplotype distribution pattern of males indicated that they migrate between the boundaries of these blocks. Moreover, a clumped distribution pattern suggests the history of matrilineal distribution by group fission and geographic expansion.

It is unknown how very short introns (65 nt; termed 'ultra-short' introns) could be spliced in a massive spliceosome (2.7 MDa) without steric hindrance. By screening an annotated human transcriptome database (H-InvDB), we identified three model ultra-short introns: the 56-nt intron in the HNRNPH1 (hnRNP H1) gene, the 49-nt intron in the NDOR1 (NADPH dependent diflavin oxidoreductase 1) gene, and the 43-nt intron in the ESRP2 (epithelial splicing regulatory protein 2) gene. We verified that these endogenous ultra-short introns are spliced, and also recapitulated this in cultured cells transfected with the corresponding mini-genes. The splicing of these ultra-short introns was repressed by a splicing inhibitor, spliceostatin A, suggesting that SF3b (a U2 snRNP component) is involved in their splicing processes. The 56-nt intron containing a pyrimidine-rich tract was spliced out in a lariat form, and this splicing was inhibited by the disruption of U1, U2, or U4 snRNA. In contrast, the 49- and 43-nt introns were purine-rich overall without any pyrimidine-rich tract, and these lariat RNAs were not detectable. Remarkably, shared G-rich intronic sequences in the 49- and 43-nt introns were required for their splicing, suggesting that these ultra-short introns may recruit a novel auxiliary splicing mechanism linked to G-rich intronic splicing enhancers.

Survival for All: Let’s Share Benefits and Hardships Libraries have played a central role in the inheritance of human knowledge. Thus, publicity and systematic management are keys in accomplishing the library’s mission. In natural science, natural history museums have been functioning as a public domain of knowledge inheritance. Because of this, morphologists and taxonomists have been conducting research by using specimens in museums. Recently, biological databases have provided public access to various types of research sources, such as research papers, sequences of nucleic and amino acids, and medical samples. Although such new public resources for life science have great potential, they have unstable financial bases.

Modern human populations are known to contain "ancient haplotypes" that originated from archaic humans by hybridization. Some of them had been reported before the development of human genomic diversity databases, such as HapMap. Consequently, some of them have no information about linkage disequilibrium (LD) regions. Because genetic information within LD is tightly linked, to know LD region containing ancient haplotypes will be useful to estimate basic parameters of admixture events, and to infer biological functions that linked with the ancient haplotypes. One of these ancient haplotypes, haplotype X (hX) was found in a 10.1 kb-region located on Xp11.22, which diverged at 1.4 M years ago, with low diversity within the cluster in gene genealogy and worldwide distribution in low frequency. We determined the LD region around the ancient haplotypes using LD information obtained in the HapMap project. The LD determination presents that the LD region surrounding the hX is stable and contains genic regions that may associate with neural and brain functions.

Mitochondrial DNA (mtDNA) is present in much higher copies, 100–1,000 times more, in a cell than nuclear DNA, and sequencing of mtDNA is not subject to the artifacts such as stutter band or allelic dropout that are common in microsatellite genotyping, which is a popular method used in wildlife conservation studies. Although the scope of mtDNA analysis is limited to matrilineal genetic information, it provides haplotypic information (i.e., DNA sequences without interruptions from recombination). This feature of mtDNA allow us to address interesting questions about Bossou and Nimba chimpanzees, including their breeding history, their genetic relationship, social behavior among matrilineal kin, and genetic diversity. Our mtDNA sequencing of Bossou and Nimba chimpanzees found the following: (1) no clear population structure emerged as in other populations of West African chimpanzee (Pan troglodytes verus), suggesting that gene flow was sufficient in the common ancestral population of West African chimpanzees; (2) there is no matrilineal gene exchange between Bossou and Nimba chimpanzees among the current generations; (3) four of the eight original adult members of the Bossou community are likely matrilineal relatives; (4) so far, the matrilineal genetic diversity of the Bossou community is still relatively stable, but there is a high risk of haplotype extinction because of a skewed haplotype distribution and a predicted relaxation of incest avoidance; and (5) finally, hair samples showed an inferior result to other sample types, such as urine, feces, and wadges (chewed-up fruit remains).

Since pre-mRNA splicing is processed by a spliceosome that is a huge complex consists of RNAs and as many as 200 proteins, it is considered that minimum intron size is limited by molecular size of the spliceosome. To understand splicing mechanism, we investigate human intron size distribution using annotated transcriptome database, H-InvDB. Distribution of obtained intron length shows that there is a mode at 83 nt in length with 4049 transcripts and number of introns decreases drastically in shorter than 65 nt, where numbers of transcripts in each locus also decrease. However, some introns less than 65 nt are observed, both shows high Codon Adaptation Index (CAI) and also observed orthologous transcripts in other mammals. These suggests that there is a limitation in minimum length to spliced out in general splicing mechanism of human, but also unknown mechanism may be exist to splice small intron exceeding the limitation.

Background Alternative splicing (AS) is a key molecular process that endows biological functions with diversity and complexity. Generally, functional redundancy leads to the generation of new functions through relaxation of selective pressure in evolution, as exemplified by duplicated genes. It is also known that alternatively spliced exons (ASEs) are subject to relaxed selective pressure. Within consensus sequences at the splice junctions, the most conserved sites are dinucleotides at both ends of introns (splice dinucleotides). However, a small number of single nucleotide polymorphisms (SNPs) occur at splice dinucleotides. An intriguing question relating to the evolution of AS diversity is whether mutations at splice dinucleotides are maintained as polymorphisms and produce diversity in splice patterns within the human population. We therefore surveyed validated SNPs in the database dbSNP located at splice dinucleotides of all human genes that are defined by the H-Invitational Database. Results We found 212 validated SNPs at splice dinucleotides (sdSNPs); these were confirmed to be consistent with the GT-AG rule at either allele. Moreover, 53 of them were observed to neighbor ASEs (AE dinucleotides). No significant differences were observed between sdSNPs at AE dinucleotides and those at constitutive exons (CE dinucleotides) in SNP properties including average heterozygosity, SNP density, ratio of predicted alleles consistent with the GT-AG rule, and scores of splice sites formed with the predicted allele. We also found that the proportion of non-conserved exons was higher for exons with sdSNPs than for other exons. Conclusions sdSNPs are found at CE dinucleotides in addition to those at AE dinucleotides, suggesting two possibilities. First, sdSNPs at CE dinucleotides may be robust against sdSNPs because of unknown mechanisms. Second, similar to sdSNPs at AE dinucleotides, those at CE dinucleotides cause differences in AS patterns because of the arbitrariness in the classification of exons into alternative and constitutive type that varies according to the dataset. Taking into account the absence of differences in sdSNP properties between those at AE and CE dinucleotides, the increased proportion of non-conserved exons found in exons flanked by sdSNPs suggests the hypothesis that sdSNPs are maintained at the splice dinucleotides of newly generated exons at which negative selection pressure is relaxed.

Creation of a vast variety of proteins is accomplished by genetic variation and a variety of alternative splicing transcripts. Currently, however, the abundant available data on genetic variation and the transcriptome are stored independently and in a dispersed fashion. In order to provide a research resource regarding the effects of human genetic polymorphism on various transcripts, we developed VarySysDB, a genetic polymorphism database based on 187 156 extensively annotated matured mRNA transcripts from 36 073 loci provided by H-InvDB. VarySysDB offers information encompassing published human genetic polymorphisms for each of these transcripts separately. This allows comparisons of effects derived from a polymorphism on different transcripts. The published information we analyzed includes single nucleotide polymorphisms and deletion–insertion polymorphisms from dbSNP, copy number variations from Database of Genomic Variants, short tandem repeats and single amino acid repeats from H-InvDB and linkage disequilibrium regions from D-HaploDB. The information can be searched and retrieved by features, functions and effects of polymorphisms, as well as by keywords. VarySysDB combines two kinds of viewers, GBrowse and Sequence View, to facilitate understanding of the positional relationship among polymorphisms, genome, transcripts, loci and functional domains. We expect that VarySysDB will yield useful information on polymorphisms affecting gene expression and phenotypes. VarySysDB is available at http://h-invitational.jp/varygene/.

The chimpanzee populations of the Bossou and Nimba regions in West Africa were genetically surveyed to 1) reveal the genetic relationship between the Bossou and Nimba populations, and 2) elucidate the evolutionary relationship between the Bossou-Nimba and other West African populations. The chimpanzee group at Bossou is characterized by its small population size, no evidence of contact with neighboring populations, and no female immigration. It is believed that most females and adolescent males emigrate from this population. To reveal the genetic signature of these characteristics, we examined the genetic diversity of Bossou and two neighboring populations (Seringbara and Yealé) in the Nimba Mountains by sequencing approximately 605 bp of the mitochondrial DNA (mtDNA) control region. A total of 20 distinct mtDNA variants were observed from 56 sequences of noninvasively collected, anonymous samples. Nucleotide diversity in the Nimba Mountain populations was 0.03-0.04, and did not differ significantly from that in the Bossou population. Very few mitochondrial variants are shared among the sites sampled, which suggests that there is little gene flow involving mtDNA. Nevertheless, no clear population structures were revealed in either population. A comparison with published sequences from West African chimpanzees (Pan troglodytes verus) indicates that the variants observed in the Bossou and Nimba regions are scattered throughout the subspecies, rather than clustered according to geographic region. This suggests that the Bossou-Nimba populations derived only recently from the common ancestral population of the West African chimpanzees, and did not pass through a bottleneck.

The dopamine receptor D4 ( DRD4 ) has received increasing research attention in behavioral science, psychiatry, and psychopharmacology. However, the number of available genetic markers for primates is still insufficient. We identified a novel variation/polymorphism in the second intron of DRD4 in humans based on the survey of 210 Japanese: a 6 bp insertion (allele frequency: 0.002) and 8 bp deletion (0.024); however, 94 Hungarian Caucasians were found to be monomorphic. Polymorphisms of the homologous region were also found in a survey of 93 specimens from four species of great apes and 51 specimens from seven species of gibbons. The polymorphisms consist of both single nucleotide substitutions and variations in the number of tandem duplications of short GC-rich sequences. Because of usefulness of primates in behavioral science, this polymorphism may be a useful marker for association studies with behavioral traits in both humans and apes.