Splicing Modulators Are Involved in Human Polyglutamine Diversification via Protein Complexes Shuttling between Nucleus and Cytoplasm

Length polymorphisms of polyglutamine (polyQs) in triplet-repeat-disease-causing genes have diversified during primate evolution despite them conferring a risk of human-specific diseases. To explain the evolutionary process of this diversification, there is a need to focus on mechanisms by which rapid evolutionary changes can occur, such as alternative splicing. Proteins that can bind polyQs are known to act as splicing factors and may provide clues about the rapid evolutionary process. PolyQs are also characterized by the formation of intrinsically disordered (ID) regions, so I hypothesized that polyQs are involved in the transportation of various molecules between the nucleus and cytoplasm to regulate mechanisms characteristic of humans such as neural development. To determine target molecules for empirical research to understand the evolutionary change, I explored protein-protein interactions (PPIs) involving the relevant proteins. This study identified pathways related to polyQ binding as hub proteins scattered across various regulatory systems, including regulation via PQBP1, VCP, or CREBBP. Nine ID hub proteins with both nuclear and cytoplasmic localization were found. Functional annotations suggested that ID proteins containing polyQs are involved in regulating transcription and ubiquitination by flexibly changing PPI formation. These findings explain the relationships among splicing complex, polyQ length variations, and modifications in neural development.

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

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

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.

Mechanistic insights into human pre-mRNA splicing of human ultra-short introns: potential unusual mechanism identifies G-rich introns

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.

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

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.

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.

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

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⁻³ per site, which is lower than that of nonsynonymous SNPs (2.1×10⁻³ 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/).

Divergent haplotypes and human history as revealed in a worldwide survey of X-linked DNA sequence variation

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.

Nucleotide sequence comparison of a chromosome rearrangement on human chromosome 12 and the corresponding ape chromosomes

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.

Mitochondrial DNA genealogy of chimpanzees in the Nimba mountains and Bossou, West Africa

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.

Allele frequency distribution of the canine dopamine receptor D4 gene exon III and I in 23 breeds

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.

Polymorphism in the second intron of dopamine receptor D4 gene in humans and apes

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 6bp insertion (allele frequency: 0.002) and 8bp 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.

High prevalence of simian T-lymphotropic virus type L in wild ethiopian baboons

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.

Mitochondrial sequence diversity within a subspecies of savanna monkeys (Cercopithecus aethiops) is similar to that between subspecies

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.