Makoto K. Shimada

High Prevalence of Simian T-Lymphotropic Virus Type L in Wild Ethiopian Baboons

ABSTRACT 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.

Distribution and effects of nonsense polymorphisms in human genes.

Background A great amount of data has been accumulated on genetic variations in the human genome, but we still do not know much about how the genetic variations affect gene function. In particular, little is known about the distribution of nonsense polymorphisms in human genes despite their drastic effects on gene products. Methodology/Principal Findings To detect polymorphisms affecting gene function, we analyzed all publicly available polymorphisms in a database for single nucleotide polymorphisms (dbSNP build 125) located in the exons of 36,712 known and predicted protein-coding genes that were defined in an annotation project of all human genes and transcripts (H-InvDB ver3.8). We found a total of 252,555 single nucleotide polymorphisms (SNPs) and 8,479 insertion and deletions in the representative transcripts in these genes. The SNPs located in ORFs include 40,484 synonymous and 53,754 nonsynonymous SNPs, and 1,258 SNPs that were predicted to be nonsense SNPs or read-through SNPs. We estimated the density of nonsense SNPs to be 0.85×10−3 per site, which is lower than that of nonsynonymous SNPs (2.1×10−3 per site). On average, nonsense SNPs were located 250 codons upstream of the original termination codon, with the substitution occurring most frequently at the first codon position. Of the nonsense SNPs, 581 were predicted to cause nonsense-mediated decay (NMD) of transcripts that would prevent translation. We found that nonsense SNPs causing NMD were more common in genes involving kinase activity and transport. The remaining 602 nonsense SNPs are predicted to produce truncated polypeptides, with an average truncation of 75 amino acids. In addition, 110 read-through SNPs at termination codons were detected. Conclusion/Significance Our comprehensive exploration of nonsense polymorphisms showed that nonsense SNPs exist at a lower density than nonsynonymous SNPs, suggesting that nonsense mutations have more severe effects than amino acid changes. The correspondence of nonsense SNPs to known pathological variants suggests that phenotypic effects of nonsense SNPs have been reported for only a small fraction of nonsense SNPs, and that nonsense SNPs causing NMD are more likely to be involved in phenotypic variations. These nonsense SNPs may include pathological variants that have not yet been reported. These data are available from Transcript View of H-InvDB and VarySysDB (http://h-invitational.jp/varygene/).

Geographic Distribution of Mitochondrial DNA Variations among Grivet (Cercopithecus aethiops aethiops) Populations in Central Ethiopia

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.

VarySysDB: a human genetic polymorphism database based on all H-InvDB transcripts

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/.

Genetic variation of blood proteins within and between local populations of grivet monkey (Cercopithecus aethiops aethiops) in central Ethiopia

A population genetic study by blood protein electrophoresis revealed that populations of wild grivet monkeys in central Ethiopia show a comparatively low level of variability and less differentiation among local populations over broad geographical areas. This is evaluated by comparison with other wild primate studies using the same electrophoretic technique. A total of 196 blood samples, collected from 10 local populations comprising 11 troops distributed along approximately 600 km of the Awash River, were examined for 33 genetic loci. The low level of variability was indicated by the proportion of polymorphic loci (Ppoly), which was on average 11% with an average heterozygosity (H) of 3%. A tendency for lower genetic differentiation among local populations was shown by theGST value of 0.09, an averageFST of 0.08, andNeis genetic distance; ranging from 0.002 to 0.023. Considering paleoclimatological studies of the area and ecological traits of this species, these findings can be explained as a consequence of a comparatively recent and repeated series of rapid habitat expansions following severe climatic conditions.

A comprehensive survey of human polymorphisms at conserved splice dinucleotides and its evolutionary relationship with alternative splicing

BackgroundAlternative 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.ResultsWe 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.ConclusionssdSNPs 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.

Skewed Matrilineal Genetic Composition in a Small Wild Chimpanzee Community

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.

Selection pressure on human STR loci and its relevance in repeat expansion disease

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.

A modification of the PHYLIP program: A solution for the redundant cluster problem, and an implementation of an automatic bootstrapping on trees inferred from original data

Felsensteins 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/.

Identification and Validation of Evolutionarily Conserved Unusually Short Pre-mRNA Introns in the Human Genome.

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.