Abstract To verify the “nearly neutral theory (NNT),” the ratio of nonsynonymous to synonymous substitutions (dN/dS) was compared among populations of different species. To determine the validity of NNT, however, populations that are genetically isolated from each other but share the same selection agents and differ in size should be compared. Genetically different lineages of obligate asexual Daphnia pulex invading Japan from North America are an ideal example as they satisfy these prerequisites. Therefore, we analyzed the whole-genome sequences of 18 genotypes, including those of the two independently invaded D. pulex lineages (JPN1 and JPN2) and compared the dN/dS ratio between the lineages. The base substitution rate of each genotype demonstrated that the JPN1 lineage having a larger distribution range diverged earlier and thus was older than the JPN2 lineage. Comparisons of the genotypes within lineages revealed that changes in dN/dS occurred after the divergence and were larger in the younger lineage, JPN2. These results imply that the JPN1 lineage has been more effectively subjected to purification selections, while slightly deteriorating mutations are less purged in JPN2 with smaller population size. Altogether, the lineage-specific difference in the dN/dS ratio for the obligate asexual D. pulex was well explained by the NNT.

Abstract The enhanced green fluorescent protein (EGFP) is considered to be a harmless protein because the critical expression level that causes growth defects is higher than that of other proteins. Here, we found that overexpression of EGFP, but not a glycolytic protein Gpm1, triggered the cell elongation phenotype in the budding yeast Saccharomyces cerevisiae. By the morphological analysis of the cell overexpressing fluorescent protein and glycolytic enzyme variants, we revealed that cysteine content was associated with the cell elongation phenotype. The abnormal cell morphology triggered by overexpression of EGFP was also observed in the fission yeast Schizosaccharomyces pombe. Overexpression of cysteine-containing protein was toxic, especially at high-temperature, while the toxicity could be modulated by additional protein characteristics. Investigation of protein aggregate formation, morphological abnormalities in mutants, and transcriptomic changes that occur upon overexpression of EGFP variants suggested that perturbation of the proteasome by the exposed cysteine of the overexpressed protein causes cell elongation. Overexpression of proteins with relatively low folding properties, such as EGFP, was also found to promote the formation of SHOTA (Seventy kDa Heat shock protein-containing, Overexpression-Triggered Aggregates), an intracellular aggregate that incorporates Hsp70/Ssa1, which induces a heat shock response, while it was unrelated to cell elongation. Evolutionary analysis of duplicated genes showed that cysteine toxicity may be an evolutionary bias to exclude cysteine from highly expressed proteins. The overexpression of cysteine-less moxGFP, the least toxic protein revealed in this study, would be a good model system to understand the physiological state of protein burden triggered by ultimate overexpression of harmless proteins.

Recent studies have identified numerous RNAs with both coding and noncoding functions. However, the sequence characteristics that determine this bifunctionality remain largely unknown. In the present study, we develop and test the open reading frame (ORF) dominance score, which we define as the fraction of the longest ORF in the sum of all putative ORF lengths. This score correlates with translation efficiency in coding transcripts and with translation of noncoding RNAs. In bacteria and archaea, coding and noncoding transcripts have narrow distributions of high and low ORF dominance, respectively, whereas those of eukaryotes show relatively broader ORF dominance distributions, with considerable overlap between coding and noncoding transcripts. The extent of overlap positively and negatively correlates with the mutation rate of genomes and the effective population size of species, respectively. Tissue-specific transcripts show higher ORF dominance than ubiquitously expressed transcripts, and the majority of tissue-specific transcripts are expressed in mature testes. These data suggest that the decrease in population size and the emergence of testes in eukaryotic organisms allowed for the evolution of potentially bifunctional RNAs.

Abstract Recent progress in high throughput single cell RNA-seq (scRNA-seq) has activated the development of data-driven inferring methods of gene regulatory networks. Most network estimations assume that perturbations produce downstream effects. However, the effects of gene perturbations are sometimes compensated by a gene with redundant functionality (functional compensation). In order to avoid functional compensation, previous studies constructed double gene deletions, but its vast nature of gene combinations was not suitable for comprehensive network estimation. We hypothesized that functional compensation may emerge as a noise change without mean change (noise-only change) due to varying physical properties and strong compensation effects. Here, we show compensated interactions, which are not detected by mean change, are captured by noise-only change quantified from scRNA-seq. We investigated whether noise-only change genes caused by a single deletion of STP1 and STP2, which have strong functional compensation, are enriched in redundantly regulated genes. As a result, noise-only change genes are enriched in their redundantly regulated genes. Furthermore, novel downstream genes detected from noise change are enriched in “transport”, which is related to known downstream genes. Herein, we suggest the noise difference comparison has the potential to be applied as a new strategy for network estimation that capture even compensated interaction.

Wild rice species have long awns at their seed tips, but this trait has been lost through rice domestication. Awn loss mitigates harvest and seed storage; further, awnlessness increases the grain number and, subsequently, improves grain yield in Asian cultivated rice, highlighting the contribution of the loss of awn to modern rice agriculture. Therefore, identifying the genes regulating awn development would facilitate the elucidation of a part of the domestication process in rice and increase our understanding of the complex mechanism in awn morphogenesis. To identify the novel loci regulating awn development and understand the conservation of genes in other wild rice relatives belonging to the AA genome group, we analyzed the chromosome segment substitution lines (CSSL). In this study, we compared a number of CSSL sets derived by crossing wild rice species in the AA genome group with the cultivated species Oryza sativa ssp. japonica. Two loci on chromosomes 7 and 11 were newly discovered to be responsible for awn development. We also found wild relatives that were used as donor parents of the CSSLs carrying the functional alleles responsible for awn elongation, REGULATOR OF AWN ELONGATION 1 (RAE1) and RAE2. To understand the conserveness of RAE1 and RAE2 in wild rice relatives, we analyzed RAE1 and RAE2 sequences of 175 accessions among diverse AA genome species retrieved from the sequence read archive (SRA) database. Comparative sequence analysis demonstrated that most wild rice AA genome species maintained functional RAE1 and RAE2, whereas most Asian rice cultivars have lost either or both functions. In addition, some different loss-of-function alleles of RAE1 and RAE2 were found in Asian cultivated species. These findings suggest that different combinations of dysfunctional alleles of RAE1 and RAE2 were selected after the speciation of O. sativa, and that two-step loss of function in RAE1 and RAE2 contributed to awnlessness in Asian cultivated rice.

Loss of genetic diversity increases the risk of extinction of a population through inbreeding depression. In addition, the number of deleterious genetic variations, which might accumulate in a small population through genetic drift, can also make the population vulnerable. Rhododendron boninense is a critically endangered azalea species, which grows on the Bonin Islands, a UNESCO World Heritage Site; its population size in the wild dropped to only a single individual. Although a governmental program for species conservation has been conducted, no successful natural regeneration or population growth was achieved. We hypothesized that the reduced genetic diversity and accumulated deleterious genetic variation might cause the vulnerability of the species, and conducted comparative genetic and genomic analyses between R. boninense and its congeners. Genetic analysis using microsatellite markers indicated that the genetic diversity of R. boninense was ultimately low; microsatellite loci of this species were all fixed, whereas the congeners maintain high allelic diversity. Based on comprehensive transcriptome analysis, the amount of deleterious variation of R. boninense was significantly greater than that of R. amanoi. Repeated generation transitions of the small population of R. boninense on the oceanic islands is likely to have resulted in low genetic diversity as well as more deleterious variations in the genome, and we speculated that the accumulated deleterious variation might result in vulnerability of R. boninense. This study indicates the possibility of evaluating the feasibility of effective conservation programs based on the genetic diversity and the amount of accumulated deleterious variation in the genome of critically endangered species.

BACKGROUND: Copy number variants (CNVs) have been reported to be associated with diseases, traits, and evolution. However, it is hard to determine which gene should have priority as a target for further functional experiments if a CNV is rare or a singleton. In this study, we attempted to overcome this issue by using two approaches: by assessing the influences of gene dosage sensitivity and gene expression sensitivity. Dosage sensitive genes derived from two-round whole-genome duplication in previous studies. In addition, we proposed a cross-sectional omics approach that utilizes open data from GTEx to assess the effect of whole-genome CNVs on gene expression. METHODS: Affymetrix Genome-Wide SNP Array 6.0 was used to detect CNVs by PennCNV and CNV Workshop. After quality controls for population stratification, family relationship and CNV detection, 287 patients with narcolepsy, 133 patients with essential hypersomnia, 380 patients with panic disorders, 164 patients with autism, 784 patients with Alzheimer disease and 1280 healthy individuals remained for the enrichment analysis. RESULTS: Overall, significant enrichment of dosage sensitive genes was found across patients with narcolepsy, panic disorders and autism. Particularly, significant enrichment of dosage-sensitive genes in duplications was observed across all diseases except for Alzheimer disease. For deletions, less or no enrichment of dosage-sensitive genes with deletions was seen in the patients when compared to the healthy individuals. Interestingly, significant enrichments of genes with expression sensitivity in brain were observed in patients with panic disorder and autism. While duplications presented a higher burden, deletions did not cause significant differences when compared to the healthy individuals. When we assess the effect of sensitivity to genome dosage and gene expression at the same time, the highest ratio of enrichment was observed in the group including dosage-sensitive genes and genes with expression sensitivity only in brain. In addition, shared CNV regions among the five neuropsychiatric diseases were also investigated. CONCLUSIONS: This study contributed the evidence that dosage-sensitive genes are associated with CNVs among neuropsychiatric diseases. In addition, we utilized open data from GTEx to assess the effect of whole-genome CNVs on gene expression. We also investigated shared CNV region among neuropsychiatric diseases.

Loss of genetic diversity is known to decrease the fitness of species and is a critical factor that increases extinction risk. However, there is little evidence for higher vulnerability and extinction risk in endangered species based on genomic differences between endangered and non-endangered species. This is true even in the case of functional loci, which are more likely to relate to the fitness of species than neutral loci. Here, we compared the genome-wide genetic diversity, proportion of duplicated genes (PD), and accumulation of deleterious variations of endangered island endemic (EIE) plants from four genera with those of their non-endangered (NE) widespread congeners. We focused on exhaustive sequences of expressed genes obtained by RNA sequencing. Most EIE species exhibited significantly lower genetic diversity and PD than NE species. Additionally, all endangered species accumulated deleterious variations. Our findings provide new insights into the genomic traits of EIE species.

Thermal tolerances of organisms play a role in defining geographic ranges and occurrence of species. In Cuba, three sympatric species of Anolis lizards (Anolis allogus, Anolis homolechis and Anolis sagrei) inhabit different thermal microhabitats. A previous study found that these species showed distinct gene expression patterns in response to temperature stimuli, suggesting the genetically distinct thermal physiology among species. To investigate whether the Anolis species inhabiting locally distinct thermal habitats diverge their thermal tolerances, we first conducted behavioural experiments to analyse the temperatures at which the three Anolis species escape from heat source. Then, for each of the three species, we isolated cDNA encoding a putative molecular heat sensor, transient receptor potential ion channel ankyrin 1 (TRPA1), which has been suggested to play a role on eliciting behavioural responses to heat stimuli. We performed electrophysiological analysis to quantify activation temperature of Anolis TRPA1 to see whether the pattern of divergence in TRPA1 responses is congruent with that of divergence in behavioural responses. We found that temperatures triggering behavioural and TRPA1 responses were significantly lower for shade-dwelling species (A. allogus) than for sun-dwelling species (A. homolechis and A. sagrei). The ambient temperature of shade habitats where A. allogus occurs stays relatively cool compared to that of open habitats where A. homolechis and A. sagrei occur and bask. The high temperature thresholds of A. homolechis and A. sagrei may reflect their heat tolerances that would benefit these species to inhabit the open habitats.

A fraction of genetic variants segregating in any population are deleterious, which negatively impacts individual fitness. The domestication of animals and plants is associated with population bottlenecks and artificial selection, which are predicted to increase the proportion of deleterious variants. However, the extent to which this is a general feature of domestic species is unclear. Here, we examine the effects of domestication on the prevalence of deleterious variation using pooled whole-genome resequencing data from five domestic animal species (dog, pig, rabbit, chicken, and silkworm) and two domestic plant species (rice and soybean) compared with their wild ancestors. We find significantly reduced genetic variation and increased proportion of nonsynonymous amino acid changes in all but one of the domestic species. These differences are observable across a range of allele frequencies, both common and rare. We find proportionally more single nucleotide polymorphisms in highly conserved elements in domestic species and a tendency for domestic species to harbor a higher proportion of changes classified as damaging. Our findings most likely reflect an increased incidence of deleterious variants in domestic species, which is most likely attributable to population bottlenecks that lead to a reduction in the efficacy of selection. An exception to this pattern is displayed by European domestic pigs, which do not show traces of a strong population bottleneck and probably continued to exchange genes with wild boar populations after domestication. The results presented here indicate that an elevated proportion of deleterious variants is a common, but not ubiquitous, feature of domestic species.

Invaded species often can rapidly expand and establish in novel environments through adaptive evolution, resulting in devastating effects on native communities. However, it is unclear if genetic variation at whole-genomic levels is actually reduced in the introduced populations and which genetic changes have occurred responding to adaptation to new environments. In the 1960s, Anolis carolinensis was introduced onto one of the Ogasawara Islands, Japan, and subsequently expanded its range rapidly throughout two of the islands. Morphological comparison showed that lower hindlimb length in the introduced populations tended to be longer than those in its native Florida populations. Using re-sequenced whole genomic data, we estimated that the effective population size at the time of introduction was actually small (less than 50). We also inferred putative genomic regions subject to natural selection after this introduction event using SweeD and a method based on Tajima's D, pi and F-ST. Five candidate genes that were potentially subject to selection were estimated by both methods. The results suggest that there were standing variations that could potentially contribute to adaptation to nonnative environments despite the founder population being small.

Cells cope with temperature elevations, which cause protein misfolding, by expressing heat shock proteins (HSPs). This adaptive response is called the heat shock response (HSR), and it is regulated mainly by heat shock transcription factor (HSF). Among the four HSF family members in vertebrates, HSF1 is a master regulator of HSP expression during proteotoxic stress including heat shock in mammals, whereas HSF3 is required for the HSR in birds. To examine whether only one of the HSF family members possesses the potential to induce the HSR in vertebrate animals, we isolated cDNA clones encoding lizard and frog HSF genes. The reconstructed phylogenetic tree of vertebrate HSFs demonstrated that HSF3 in one species is unrelated with that in other species. We found that the DNA-binding activity of both HSF1 and HSF3 in lizard and frog cells was induced in response to heat shock. Unexpectedly, overexpression of lizard and frog HSF3 as well as HSF1 induced HSP70 expression in mouse cells during heat shock, indicating that the two factors have the potential to induce the HSR. Furthermore, knockdown of either HSF3 or HSF1 markedly reduced HSP70 induction in lizard cells and resistance to heat shock. These results demonstrated that HSF1 and HSF3 cooperatively regulate the HSR at least in lizards, and suggest complex mechanisms of the HSR in lizards as well as frogs.

The green anole Anolis carolinensis invaded the Ogasawara Islands in Japan, drove various native species to extinction, and its distribution expanded 14years after initial establishment. A. carolinensis invaded Okinawa Island, but it has not expanded its distribution in more than 25years, although its density is extremely high in the southern region. To determine whether A. carolinensis has the potential to expand its distribution on Okinawa Island, we performed phylogenetic analysis of mitochondrial ND2 DNA sequences to study the origin of A. carolinensis that invaded Okinawa Island. We further used a species distribution model (MaxEnt) based on the distribution of native populations in North America to identify ecologically suitable areas on Okinawa Island. Nucleotide sequence analysis shows that the invader A. carolinensis originated in the western part of the Gulf Coast and inland areas of the United States and that a portion of the anoles on Okinawa was not introduced via the Ogasawara Islands. The MaxEnt predictions indicate that most areas in Okinawa Island are suitable for A. carolinensis. Therefore, A.carolinensis may have the potential to expand its distribution in Okinawa Island. The predictions indicate that habitat suitability is high in areas of high annual mean temperature and urbanized areas. The values of precipitation in summer in the northern region of Okinawa Island were higher compared with those of North America, which reduced the habitat suitability in Okinawa Island. Adaptation to low temperatures, an increase in the mean temperature through global warming, and an increase in open environments through land development will likely expand the distribution of A. carolinensis in Okinawa Island. Therefore, we must continue to monitor the introduced populations and be alert to the possibility that city planning that increases open environments may cause their range to expand.

Phaeomelanin is a common pigment that confers a reddish colour to animals. Since phaeomelanogenesis requires the sulfhydryl group from cysteine or glutathione (GSH), which is an important antioxidant, this pigmentation and the associated coloration may be an honest signal, whereby only high-quality individuals (e.g. with lower oxidative stress) are able to develop showy plumage. The present study tested the mechanisms underlying the honest signal hypothesis using nestling barn swallows, Hirundo rustica gutturalis, which exhibit phaeomelanic throat plumage patches. We examined the relationship between phaeomelanin pigmentation levels and physiological condition during trait development, and the expression of the phaeomelanin-related gene agouti-signalling protein (ASIP) and the GSH-related gene glutathione S-transferase (GST) in throat feather follicles. We found that during phaeomelanogenesis, heavier nestlings produced more pigmented feathers, indicating that nestlings with high phaeomelanin concentrations are in better condition. We also found that phaeomelanin concentration was negatively correlated with total GSH level, but not significantly related with measures of oxidative stress. Among the GST genes, GSTM3 exhibited the highest expression in the developing feathers during phaeomelanogenesis. The expression levels of ASIP were positively associated with the amount of phaeomelanin deposition and negatively associated with the expression of GSTM3, reducing the amount of GSH that was available as an antioxidant. These findings suggest that high-quality individuals produce high concentrations of phaeomelanin in their plumage without experiencing increased oxidative stress, despite phaeomelanin production, which is triggered by ASIP, potentially actively consuming the sulfhydryl group from GSH.

As most genes are shared between females and males, DNA methylation is assumed to play a crucial role in sex-biased gene expression. DNA methylation exclusively occurs at CpG dinucleotides, and therefore, we would expect that CpG density around transcription start sites (TSSs) relate to sex-biased gene expression. Here we investigated the relationship between CpG densities around TSSs and the ratio of gene expression levels between sexes in the guppy (Poecilia reticulata), which displays remarkable sexual dimorphisms. We found that genes with sex-biased gene expression had different CpG densities downstream of TSSs compared with genes lacking sex-biased gene expression. Intriguingly, male-biased expression genes with intermediate CpG density downstream of TSSs exhibited greater differences in gene expression between sexes in the gonad and tail. Our findings suggested the possibility that CpGs around TSSs, especially in the downstream regions, play a crucial role in sex-biased gene expression through DNA methylation.

Differential gene expression can play an important role in phenotypic evolution and divergent adaptation. Although differential gene expression can be caused by both local- and distant-regulatory changes, we know little about their relative contribution to transcriptome evolution in natural populations. Here, we conducted expression quantitative trait loci (eQTL) analysis to investigate the genetic architecture underlying transcriptome divergence between marine and stream ecotypes of threespine sticklebacks (Gasterosteus aculeatus). We identified both local and distant eQTLs, some of which constitute hotspots, regions with a disproportionate number of significant eQTLs relative to the genomic background. The majority of local eQTLs including those in the hotspots caused expression changes consistent with the direction of transcriptomic divergence between ecotypes. Genome scan analysis showed that many local eQTLs overlapped with genomic regions of high differentiation. In contrast, nearly half of the distant eQTLs including those in the hotspots caused opposite expression changes, and few overlapped with regions of high differentiation, indicating that distant eQTLs may act as a constraint of transcriptome evolution. Finally, a comparison between two salinity conditions revealed that nearly half of eQTL hotspots were environment specific, suggesting that analysis of genetic architecture in multiple conditions is essential for predicting response to selection.

Degeneration of Y chromosomes is a common evolutionary path of XY sex chromosome systems. Recent genomic studies in flies and plants have revealed that even young neo-sex chromosomes with the age of a few million years show signs of Y degeneration, such as the accumulation of nonsense and frameshift mutations. However, it remains unclear whether neo-Y chromosomes also show rapid degeneration in fishes, which often have homomorphic sex chromosomes. Here, we investigated whether a neo-Y chromosome of Japan Sea stickleback (Gasterosteus nipponicus), which was formed by a Y-autosome fusion within the last 2 million years, accumulates deleterious mutations. Our previous genomic analyses did not detect excess nonsense and frameshift mutations on the Japan Sea stickleback neo-Y. In the present study, we found that the nonrecombining region of the neo-Y near the fusion end has accumulated nonsynonymous mutations altering amino acids of evolutionarily highly conserved residues. Enrichment of gene ontology terms related to protein phosphorylation and cellular protein modification process was found in the genes with potentially deleterious mutations on the neo-Y. These results suggest that the neo-Y of the Japan Sea stickleback has already accumulated mutations that may impair protein functions.

Adaptation to different salinities can drive and maintain divergence between populations of aquatic organisms. Anadromous and stream ecotypes of threespine stickleback (Gasterosteus aculeatus) are an excellent model to explore the genetic mechanisms underlying osmoregulation divergence. Using a parapatric pair of anadromous and stream stickleback ecotypes, we employed an integrated genomic approach to identify candidate genes important for adaptation to different salinity environments. Quantitative trait loci (QTL) mapping of plasma sodium concentrations under a seawater challenge experiment identified a significant QTL on chromosome 16. To identify candidate genes within this QTL, we first conducted RNA-seq and microarray analysis on gill tissue to find ecotypic differences in gene expression that were associated with plasma Na+ levels. This resulted in the identification of ten candidate genes. Quantitative PCR analysis on gill tissue of additional Japanese stickleback populations revealed that the majority of the candidate genes showed parallel divergence in expression levels. Second, we conducted whole-genome sequencing and found five genes that are predicted to have functionally important amino acid substitutions. Finally, we conducted genome scan analysis and found that eight of these candidate genes were located in genomic islands of high differentiation, suggesting that they may be under divergent selection. The candidate genes included those involved in ATP synthesis and hormonal signalling, whose expression or amino acid changes may underlie the variation in salinity tolerance. Further functional molecular analysis of these genes will reveal the causative genetic and genomic changes underlying divergent adaptation.

Over-expression of Winged-Eye (WGE) in the Drosophila eye imaginal disc induces an eye-to-wing transformation. Endogenous WGE is required for organ development, and wge-deficient mutants exhibit growth arrest at the larval stage, suggesting that WGE is critical for normal growth. The function of WGE, however, remains unclear. Here, we analyzed the subcellular localization of WGE to gain insight into its endogenous function. Immunostaining showed that WGE localized to specific nuclear foci called the histone locus body (HLB), an evolutionarily conserved nuclear body required for S phase-specific histone mRNA production. Histone mRNA levels and protein levels in cytosolic fractions were aberrantly up-regulated in wge mutant larva, suggesting a role for WGE in regulating histone gene expression. Genetic analyses showed that wge suppresses position-effect variegation, and that WGE and a HLB component Mute appears to be synergistically involved in heterochromatin formation. Further supporting a role in chromatin regulation, wge-deficient mutants showed derepression of retrotransposons and increased H2Av signals, a DNA damage marker. These findings suggest that WGE is a component of HLB in Drosophila with a role in heterochromatin formation and transposon silencing. We propose that WGE at HLB contributes to genomic stability and development by regulating heterochromatin structure via histone gene regulation.

How animals achieve evolutionary adaptation to different thermal environments is an important issue for evolutionary biology as well as for biodiversity conservation in the context of recent global warming. In Cuba, three sympatric species of Anolis lizards (Anolis allogus, A. homolechis and A. sagrei) inhabit different thermal microhabitats, thereby providing an excellent opportunity to examine how they have adapted to different environmental temperatures. Here, we performed RNA-seq on the brain, liver and skin tissues from these three species to analyse their transcriptional responses at two different temperatures. In total, we identified 400, 816 and 781 differentially expressed genes (DEGs) between the two temperatures in A. allogus, A. homolechis and A. sagrei, respectively. Only 62 of these DEGs were shared across the three species, indicating that global transcriptional responses have diverged among these species. Gene ontology (GO) analysis showed that large numbers of ribosomal protein genes were DEGs in the warm-adapted A. homolechis, suggesting that the upregulation of protein synthesis is an important physiological mechanism in the adaptation of this species to hotter environments. GO analysis also showed that GO terms associated with circadian regulation were enriched in all three species. A gene associated with circadian regulation, Nr1d1, was detected as a DEG with opposite expression patterns between the cool-adapted A. allogus and the hot-adapted A. sagrei. Because the environmental temperature fluctuates more widely in open habitats than in forests throughout the day, the circadian thermoregulation could also be important for adaptation to distinct thermal habitats.

Invasive species pose a major threat to biological diversity. Although introduced populations often experience population bottlenecks, some invasive species are thought to be originated from hybridization between multiple populations or species, which can contribute to the maintenance of high genetic diversity. Recent advances in genome sequencing enable us to trace the evolutionary history of invasive species even at whole-genome level and may help to identify the history of past hybridization that may be overlooked by traditional marker-based analysis. Here, we conducted whole-genome sequencing of eight threespine stickleback (Gasterosteus aculeatus) individuals, four from a recently introduced crater lake population and four of the putative source population. We found that both populations have several small genomic regions with high genetic diversity, which resulted from introgression from a closely related species (Gasterosteus nipponicus). The sizes of the regions were too small to be detected with traditional marker-based analysis or even some reduced-representation sequencing methods. Further amplicon sequencing revealed linkage disequilibrium around an introgression site, which suggests the possibility of selective sweep at the introgression site. Thus, interspecies introgression might predate introduction and increase genetic variation in the source population. Whole-genome sequencing of even a small number of individuals can therefore provide higher resolution inference of history of introduced populations.

Background: Understanding the evolutionary forces that influence variation in gene regulatory regions in natural populations is an important challenge for evolutionary biology because natural selection for such variations could promote adaptive phenotypic evolution. Recently, whole-genome sequence analyses have identified regulatory regions subject to natural selection. However, these studies could not identify the relationship between sequence variation in the detected regions and change in gene expression levels. We analyzed sequence variations in core promoter regions, which are critical regions for gene regulation in higher eukaryotes, in a natural population of Drosophila melanogaster, and identified core promoter sequence variations associated with differences in gene expression levels subjected to natural selection. Results: Among the core promoter regions whose sequence variation could change transcription factor binding sites and explain differences in expression levels, three core promoter regions were detected as candidates associated with purifying selection or selective sweep and seven as candidates associated with balancing selection, excluding the possibility of linkage between these regions and core promoter regions. CHKov1, which confers resistance to the sigma virus and related insecticides, was identified as core promoter regions that has been subject to selective sweep, although it could not be denied that selection for variation in core promoter regions was due to linked single nucleotide polymorphisms in the regulatory region outside core promoter regions. Nucleotide changes in core promoter regions of CHKov1 caused the loss of two basal transcription factor binding sites and acquisition of one transcription factor binding site, resulting in decreased gene expression levels. Of nine core promoter regions regions associated with balancing selection, brat, and CG9044 are associated with neuromuscular junction development, and Nmda1 are associated with learning, behavioral plasticity, and memory. Diversity of neural and behavioral traits may have been maintained by balancing selection. Conclusions: Our results revealed the evolutionary process occurring by natural selection for differences in gene expression levels caused by sequence variation in core promoter regions in a natural population. The sequences of core promoter regions were diverse even within the population, possibly providing a source for natural selection.

The mechanism by which genetic systems affect environmental adaptation is a focus of considerable attention in the fields of ecology, evolution, and conservation. However, the genomic characteristics that constrain adaptive evolution have remained unknown. A recent study showed that the proportion of duplicated genes in whole Drosophila genomes correlated with environmental variability within habitat, but it remains unclear whether the correlation is observed even in vertebrates whose genomes including a large number of duplicated genes generated by whole-genome duplication (WGD). Here, we focus on fully sequenced mammalian genomes that experienced WGD in early vertebrate lineages and show that the proportion of small-scale duplication (SSD) genes in the genome, but not that of WGD genes, is significantly correlated with habitat variability. Moreover, species with low habitat variability have a higher proportion of lost duplicated genes, particularly SSD genes, than those with high habitat variability. These results indicate that species that inhabit variable environments may maintain more SSD genes in their genomes and suggest that SSD genes are important for adapting to novel environments and surviving environmental changes. These insights may be applied to predicting invasive and endangered species.

Sex chromosomes turn over rapidly in some taxonomic groups, where closely related species have different sex chromosomes. Although there are many examples of sex chromosome turnover, we know little about the functional roles of sex chromosome turnover in phenotypic diversification and genomic evolution. The sympatric pair of Japanese threespine stickleback (Gasterosteus aculeatus) provides an excellent system to address these questions: the Japan Sea species has a neo-sex chromosome system resulting from a fusion between an ancestral Y chromosome and an autosome, while the sympatric Pacific Ocean species has a simple XY sex chromosome system. Furthermore, previous quantitative trait locus (QTL) mapping demonstrated that the Japan Sea neo-X chromosome contributes to phenotypic divergence and reproductive isolation between these sympatric species. To investigate the genomic basis for the accumulation of genes important for speciation on the neo-X chromosome, we conducted whole genome sequencing of males and females of both the Japan Sea and the Pacific Ocean species. No substantial degeneration has yet occurred on the neo-Y chromosome, but the nucleotide sequence of the neo-X and the neo-Y has started to diverge, particularly at regions near the fusion. The neo-sex chromosomes also harbor an excess of genes with sex-biased expression. Furthermore, genes on the neo-X chromosome showed higher non-synonymous substitution rates than autosomal genes in the Japan Sea lineage. Genomic regions of higher sequence divergence between species, genes with divergent expression between species, and QTL for inter-species phenotypic differences were found not only at the regions near the fusion site, but also at other regions along the neo-X chromosome. Neo-sex chromosomes can therefore accumulate substitutions causing species differences even in the absence of substantial neo-Y degeneration.

A number of rare copy number variants (CNVs), including both deletions and duplications, have been associated with developmental disorders, including schizophrenia, autism, intellectual disability, and epilepsy. Pathogenicity may derive from dosage sensitivity of one or more genes contained within the CNV locus. To understand pathophysiology, the specific disease-causing gene(s) within each CNV need to be identified. In the present study, we test the hypothesis that ohnologs (genes retained after ancestral whole-genome duplication events, which are frequently dosage sensitive) are overrepresented in pathogenic CNVs. We selected three sets of genes implicated in copy number pathogenicity: (i) genes mapping within rare disease-associated CNVs, (ii) genes within de novo CNVs under negative genetic selection, and (iii) genes identified by clinical array comparative genome hybridization studies as potentially pathogenic. We compared the proportion of ohnologs between these gene sets and control genes, mapping to CNVs not known to be disease associated. We found that ohnologs are significantly overrepresented in genes mapping to pathogenic CNVs, irrespective of how CNVs were identified, with over 90% containing an ohnolog, compared with control CNVs >100 kb, where only about 30% contained an ohnolog. In some CNVs, such as del15p11.2 (CYFIP1) and dup/del16p13.11 (NDE1), the most plausible prior candidate gene was also an ohnolog, as were the genes VIPR2 and NRXN1, each found in short CNVs containing no other genes. Our results support the hypothesis that ohnologs represent critical dosage-sensitive elements of the genome, possibly responsible for some of the deleterious phenotypes observed for pathogenic CNVs and as such are readily identifiable candidate genes for further study.

Most copy number variations are neutral, but some are deleterious and associated with various human diseases. Copy number variations are distributed non-randomly in vertebrate genomes, and it was recently reported that ohnologs, which are duplicated genes derived from whole genome duplication, are refractory to copy number variations. However, it is unclear what genomic factors affect the deleterious effects of copy number variations and the biological significance of the biased genomic distribution of copy number variations remains poorly understood. Here we show that non-ohnologs neighbouring ohnologs are unlikely to have copy number variations, resulting in ohnolog-rich regions in vertebrate genomes being copy number variation deserts. Our results suggest that the genomic location of ohnologs is a determining factor in the retention of copy number variations and that the dosage-balanced ohnologs are likely to cause the deleterious effects of copy number variations in these regions. We propose that investigating copy number variation of genes in regions that are typically copy number variation deserts is an efficient means to find disease-related copy number variations.

Variations in species richness of local assemblages may be explained by local ecological processes or large-scale evolutionary and biogeographical processes. In Anolis lizards, species with different ecomorphs can coexist by occupying different niches. In addition, several species with the same ecomorph (e. g., trunk-ground) can coexist, and the number of trunk-ground anole species varies among local species assemblages. In this study, we assessed the importance of ecological interactions, number of speciation events, and range expansion for local and regional species diversity of these lizards. We examined the species richness and thermal microhabitat partitioning (considered to be a measure of ecological interaction) of 12 trunk-ground anole species in 11 local assemblages in Cuba. The results indicated that the phylogenetic structure of trunk-ground anole lizard assemblages was random. However, there was an overdispersion of preferences for thermal microhabitat use, which indicates that differences in microhabitat use are likely to occur within assemblages. We suggest that the number of speciation events within regions and the number of sympatrically coexisting species increases species richness at the local level. Migration appeared to be limited, leading to the range expansion of only three species with different thermal requirements. The thermal niches of species were conserved within Anolis allogus clade, whereas species within the Anolis homolechis and Anolis sagrei clades tended to change their thermal niches. Our results suggest that the species composition and richness in local assemblages could be explained by evolutionary history (the number of speciation events and limits to range expansion) and ecological processes (habitat partitioning). Of the ecological factors, the number of thermal (microhabitat use) and structural niches (e.g., vegetation) could limit the potential number of coexisting species within a local assemblage.

Gene overexpression beyond a permissible limit causes defects in cellular functions. However, the permissible limits of most genes are unclear. Previously, we developed a genetic method designated genetic tug-of-war (gTOW) to measure the copy number limit of overexpression of a target gene. In the current study, we applied gTOW to the analysis of all protein-coding genes in the budding yeast Saccharomyces cerevisiae. We showed that the yeast cellular system was robust against an increase in the copy number by up to 100 copies in >80% of the genes. After frameshift and segmentation analyses, we isolated 115 dosage-sensitive genes (DSGs) with copy number limits of 10 or less. DSGs contained a significant number of genes involved in cytoskeletal organization and intracellular transport. DSGs tended to be highly expressed and to encode protein complex members. We demonstrated that the protein burden caused the dosage sensitivity of highly expressed genes using a gTOW experiment in which the open reading frame was replaced with GFP. Dosage sensitivities of some DSGs were rescued by the simultaneous increase in the copy numbers of partner genes, indicating that stoichiometric imbalances among complexes cause dosage sensitivity. The results obtained in this study will provide basic knowledge about the physiology of chromosomal abnormalities and the evolution of chromosomal composition.

Whole genome duplication (WGD) has made a significant contribution to many eukaryotic genomes including yeast, plants, and vertebrates. Following WGD, some ohnologs (WGD paralogs) remain in the genome arranged in blocks of conserved gene order and content (paralogons). However, the most common outcome is loss of one of the ohnolog pair. It is unclear what factors, if any, govern gene loss from paralogons. Recent studies have reported physical clustering (genetic linkage) of functionally linked (interacting) genes in the human genome and propose a biological significance for the clustering of interacting genes such as coexpression or preservation of epistatic interactions. Here we conduct a novel test of a hypothesis that functionally linked genes in the same paralogon are preferentially retained in cis after WGD. We compare the number of protein-protein interactions (PPIs) between linked singletons within a paralogon (defined as cis-PPIs) with that of PPIs between singletons across paralogon pairs (defined as trans-PPIs). We find that paralogons in which the number of cis-PPIs is greater than that of trans-PPIs are significantly enriched in human and yeast. The trend is similar in plants, but it is difficult to assess statistical significance due to multiple, overlapping WGD events. Interestingly, human singletons participating in cis-PPIs tend to be classified into "response to stimulus.'' We uncover strong evidence of biased gene loss after WGD, which further supports the hypothesis of biologically significant gene clusters in eukaryotic genomes. These observations give us new insight for understanding the evolution of genome structure and of protein interaction networks.

The factors limiting the habitat range of species are crucial in understanding their biodiversity and response to environmental change. Yet the genetic and genomic architectures that produce genetic variation to enable environmental adaptation have remained poorly understood. Here we show that the proportion of duplicated genes (PD) in the whole genomes of fully sequenced Drosophila species is significantly correlated with environmental variability within the habitats measured by the climatic envelope and habitat diversity. Furthermore, species with a low PD tend to lose the duplicated genes owing to their faster evolution. These results indicate that the rapid relaxation of functional constraints on duplicated genes resulted in a low PD for species with lower habitat diversity, and suggest that the maintenance of duplicated genes gives organisms an ecological advantage during evolution. We therefore propose that the PD in a genome is related to adaptation to environmental variation. © 2012 The Author.

On this study, we investigated the evolution of vertebrate tissues by examining the potential association among gene expression, duplication, and base substitution patterns. On particular, we compared whole-genome duplication (WGD) with small-scale duplication (SSD), as well as tissue restricted with ubiquitously expressed genes. All patterns were also analysed in the light of gene evolutionary rates. Among those genes characterized by rapid evolution and expressed in a restricted range of tissues, SSD was represented in a larger proportion than WGD. Conversely, genes with ubiquitous expression were associated with slower evolutionary rates and a larger proportion of WGD. The results also show that evolutionary rates were faster in genes expressed in endodermal tissues and slower in ectodermal genes. Accordingly, the proportion of the SSD and WGD genes was highest in the endoderm and ectoderm, respectively. Therefore, quickly evolving SSD genes might have contributed to the faster evolution of endodermal tissues, whereas the comparatively slowly evolving WGD genes might have functioned to maintain the basic characteristics of ectodermal tissues. Mesenchymal tissues occupied an intermediate position in this regard, whereas the patterns observed for haemocytes were unique. Rapid tissue evolution could be related to a specific gene duplication mode (SSD) and faster molecular evolution in response to exposure to the external environment. These findings reveal general patterns underlying the evolution of tissues and their corresponding genes.

How and why female somatic X-chromosome inactivation (XCI) evolved in mammals remains poorly understood. It has been proposed that XCI is a dosage-compensation mechanism that evolved to equalize expression levels of X-linked genes in females (2X) and males (1X), with a prior twofold increase in expression of X-linked genes in both sexes ("Ohno's hypothesis"). Whereas the parity of X chromosome expression between the sexes has been clearly demonstrated, tests for the doubling of expression levels globally along the X chromosome have returned contradictory results. However, changes in gene dosage during sex-chromosome evolution are not expected to impact on all genes equally, and should have greater consequences for dosage-sensitive genes. We show that, for genes encoding components of large protein complexes (>= 7 members)-a class of genes that is expected to be dosage-sensitive-expression of X-linked genes is similar to that of autosomal genes within the complex. These data support Ohno's hypothesis that XCI acts as a dosage-compensation mechanism, and allow us to refine Ohno's model of XCI evolution. We also explore the contribution of dosage-sensitive genes to X aneuploidy phenotypes in humans, such as Turner (X0) and Klinefelter (XXY) syndromes. X aneuploidy in humans is common and is known to have mild effects because most of the supernumerary X genes are inactivated and not affected by aneuploidy. Only genes escaping XCI experience dosage changes in X-aneuploidy patients. We combined data on dosage sensitivity and XCI to compute a list of candidate genes for X-aneuploidy syndromes.

Nuclear sequence markers are useful tool for the study of the history of populations and adaptation. However, it is not easy to obtain multiple nuclear primers for organisms with poor or no genomic sequence information. Here we used the genomes of organisms that have been fully sequenced to design comprehensive sets of primers to amplify polymorphic genomic fragments of multiple nuclear genes in non-sequenced organisms. First, we identified a large number of candidate polymorphic regions that were flanked on each side by conserved regions in the reference genomes. We then designed primers based on these conserved sequences and examined whether the primers could be used to amplify sequences in target species, montane brown frog (Rana ornativentris), anole lizard (Anolis sagrei), guppy (Poecilia reticulata), and fruit fly (Drosophila melanogaster), for population genetic analysis. We successfully obtained polymorphic markers for all target species studied. In addition, we found that sequence identities of the regions between the primer sites in the reference genomes affected the experimental success of DNA amplification and identification of polymorphic loci in the target genomes, and that exonic primers had a higher success rate than intronic primers in amplifying readable sequences. We conclude that this comparative genomic approach is a time-and cost-effective way to obtain polymorphic markers for non-sequenced organisms, and that it will contribute to the further development of evolutionary ecology and population genetics for non-sequenced organisms, aiding in the understanding of the genetic basis of adaptation.

Sex steroids mediate the expression of sexually dimorphic or sex-specific traits that are important both for mate choice within species and for behavioral isolation between species. We investigated divergence in sex steroid signaling between two sympatric species of threespine stickleback (Gasterosteus aculeatus): the Japan Sea form and the Pacific Ocean form. These sympatric forms diverge in both male display traits and female mate choice behaviors, which together contribute to asymmetric behavioral isolation in sympatry. Here, we found that plasma levels of testosterone and 17 beta-estradiol differed between spawning females of the two sympatric forms. Transcript levels of follicle-stimulating hormone-beta (FSH beta) gene were also higher in the pituitary gland of spawning Japan Sea females than in the pituitary gland of spawning Pacific Ocean females. By contrast, none of the sex steroids examined were significantly different between nesting males of the two forms. However, combining the plasma sex steroid data with testis transcriptome data suggested that the efficiency of the conversion of testosterone into 11-ketotestosterone has likely diverged between forms. Within forms, plasma testosterone levels in males were significantly correlated with male body size, a trait important for female mate choice in the two sympatric species. These results demonstrate that substantial divergence in sex steroid signaling can occur between incipient sympatric species. We suggest that investigation of the genetic and ecological mechanisms underlying divergence in hormonal signaling between incipient sympatric species will provide a better understanding of the mechanisms of speciation in animals.

Background: There is increasing interest in the evolution of protein-protein interactions because this should ultimately be informative of the patterns of evolution of new protein functions within the cell. One model proposes that the evolution of new protein-protein interactions and protein complexes proceeds through the duplication of self-interacting genes. This model is supported by data from yeast. We examined the relationship between gene duplication and self-interaction in the human genome. Results: We investigated the patterns of self-interaction and duplication among 34808 interactions encoded by 8881 human genes, and show that self-interacting proteins are encoded by genes with higher duplicability than genes whose proteins lack this type of interaction. We show that this result is robust against the system used to define duplicate genes. Finally we compared the presence of self-interactions amongst proteins whose genes have duplicated either through whole-genome duplication (WGD) or small-scale duplication (SSD), and show that the former tend to have more interactions in general. After controlling for age differences between the two sets of duplicates this result can be explained by the time since the gene duplication. Conclusions: Genes encoding self-interacting proteins tend to have higher duplicability than proteins lacking self-interactions. Moreover these duplicate genes have more often arisen through whole-genome rather than small-scale duplication. Finally, self-interacting WGD genes tend to have more interaction partners in general in the PIN, which can be explained by their overall greater age. This work adds to our growing knowledge of the importance of contextual factors in gene duplicability.

About 30% of protein-coding genes in the human genome are related through two whole genome duplication (WGD) events. Although WGD is often credited with great evolutionary importance, the processes governing the retention of these genes and their biological significance remain unclear. One increasingly popular hypothesis is that dosage balance constraints are a major determinant of duplicate gene retention. We test this hypothesis and show that WGD-duplicated genes (ohnologs) have rarely experienced subsequent small-scale duplication (SSD) and are also refractory to copy number variation (CNV) in human populations and are thus likely to be sensitive to relative quantities (i.e., they are dosage-balanced). By contrast, genes that have experienced SSD in the vertebrate lineage are more likely to also display CNV. This supports the hypothesis of biased retention of dosage-balanced genes after WGD. We also show that ohnologs have a strong association with human disease. In particular, Down Syndrome (DS) caused by trisomy 21 is widely assumed to be caused by dosage effects, and 75% of previously reported candidate genes for this syndrome are ohnologs that experienced no other copy number changes. We propose the remaining dosage-balanced ohnologs on chromosome 21 as candidate DS genes. These observations clearly show a persistent resistance to dose changes in genes duplicated by WGD. Dosage balance constraints simultaneously explain duplicate gene retention and essentiality after WGD.

In yeast and worm, duplicate genes overlap in function so that deleting one of a pair from the genome is less likely to be lethal than deleting a singleton gene. By contrast, previous analyses showed that mouse duplicate genes were as essential as singletons. We show that the relationship between gene duplication and essentiality is complex in multicellular organisms, with developmental genes and genes that were duplicated by whole genome duplication being more essential than other duplicated genes.

It is increasingly clear that eukaryotic gene order is nonrandom, being constrained in some instances by expression patterns, expression levels, protein interactions or epistatic effects. The relationship between gene order and function is patchy (not every gene is in a functional cluster) and the factors that influence it are numerous and interconnected. Not surprisingly, then, our knowledge of genome structure is still growing. Here we review reported gene clusters in eukaryote genomes and obstacles such as leaky expression and tandem duplication effects for identification of functional gene clusters. In particular, we show interacting gene clusters, which are identified by protein-protein interactions (PPIs), are robust against the problems. Furthermore, we emphasize that evolutionary analyses of functional gene clusters are very important to assess their biological meaning.

Pterogobius elapoides and Pterogobius zonoleucus are common free-swimming gobies found in rocky and weedy shores along the temperate coast of Japan. We collected individuals of both species from 23 locations around the coast of Japan and compared the mitochondrial nucleotide sequences of two gene regions, CytB and ND2. Phylogenetic trees constructed using the neighbor-joining, maximum parsimony, and maximum likelihood methods consistently indicated that all 125 samples of the two species, which are collected from a variety of locations in Japan, can be clearly divided into the following four clades: "Pacific P. elapoides" (Pa-ela), "Sea of Japan P. elapoides" (SJ-ela), "Pacific P. zonoleucus" (Pa-zon), and "Sea of Japan P. zonoleucus" (SJ-zon). These four monophyletic clades were supported with very high bootstrap values. Although Pa-ela and SJ-ela composed a monophyletic clade, it is noteworthy that the two clades of P. elapoides also formed a monophyletic group together with SJ-zon with a bootstrap value of 95% and 97% by the maximum likelihood and neighbor-joining methods, respectively. We observed several morphological differences between Pa-ela and SJ-ela, including; 1) six dark bands on the body in the former versus seven dark bands in the latter and 2) more pectoral-fin rays numbering 21-24 (mode 22) in the latter compared to the former (19-22, mode 21). Furthermore, the scatter plots of scores on principal components 1 and 2 based on the morphometric characters roughly separated the populations from each other. Moreover, we documented the following morphological differences between Pa-zon and SJ-zon for the first time; 1) six light bands on the body in the former versus five light bands in the latter and 2) the light bands from both eyes forming a complete U-shaped marking on the occipital region occurred in 55% of the specimens in the former versus 16% in the latter. However, no significant differences were found in the morphometric characters between the two populations of P. zonoleucus. The estimated divergence time of the two P. zonoleucus populations was 15.06 +/- 2.72 (mean +/- 1 S.E.) times earlier than that of the two P elapoides populations. However, the morphological differences between the two populations of the former were much smaller than those of the latter. An explanation for this obvious discrepancy between morphological and molecular features is proposed from an evolutionary point of view. (C) 2008 Elsevier B.V. All rights reserved.

Here we report the new features and improvements in our latest release of the H-Invitational Database (H-InvDB; http://www.h-invitational.jp/), a comprehensive annotation resource for human genes and transcripts. H-InvDB, originally developed as an integrated database of the human transcriptome based on extensive annotation of large sets of full-length cDNA (FLcDNA) clones, now provides annotation for 120 558 human mRNAs extracted from the International Nucleotide Sequence Databases (INSD), in addition to 54 978 human FLcDNAs, in the latest release H-InvDB_4.6. We mapped those human transcripts onto the human genome sequences (NCBI build 36.1) and determined 34 699 human gene clusters, which could define 34 057 (98.1%) protein-coding and 642 (1.9%) non-protein-coding loci; 858 (2.5%) transcribed loci overlapped with predicted pseudogenes. For all these transcripts and genes, we provide comprehensive annotation including gene structures, gene functions, alternative splicing variants, functional non-protein-coding RNAs, functional domains, predicted sub cellular localizations, metabolic pathways, predictions of protein 3D structure, mapping of SNPs and microsatellite repeat motifs, co-localization with orphan diseases, gene expression profiles, orthologous genes, proteinprotein interactions (PPI) and annotation for gene families. The current H-InvDB annotation resources consist of two main views: Transcript view and Locus view and eight sub-databases: the DiseaseInfo Viewer, H-ANGEL, the Clustering Viewer, G-integra, the TOPO Viewer, Evola, the PPI view and the Gene family/group.

Protein-protein interactions (PPIs) are one of the most important components of biological networks. It is important to understand the evolutionary process of PPIs in order to elucidate how the evolution of biological networks has contributed to diversification of the existent organisms. We focused on the evolutionary rates of proteins involved with PPIs, because it had been shown that for a given protein-coding gene the number of its PPIs in a biological network was one of the important factors in determining the evolutionary rate of the gene. We studied the evolutionary rates of duplicated gene products that were involved with PPIs, reviewing the current situation of this subject. In addition, we focused on how the evolutionary rates of proteins were influenced by the characteristic features of PPIs. We, then, concluded that the evolutionary rates of the proteins in the PPI networks were strongly influenced by their PPI partners. Finally, we emphasized that evolutionary considerations of the PPI proteins were very important for understanding the building up of the current PPI networks. Copyright © 2007 S. Karger AG.

in the network of protein-protein interactions (PPIs), a loss and gain of the partnering proteins can cause drastic changes of network formation during evolution. With the aim of examining the evolutionary effects of the loss and gain of the partnering proteins on PPIs, we examined a relationship between evolutionary rates and losses and/or gains of PPIs for duplicated gene pairs encoding proteins involved in the PPI network. For duplicated pairs, which provided us with a unique opportunity of making fair comparisons of the genes with the same initial condition, we found that the evolutionary rate of the protein with more PPI partners is much slower than that of the other with fewer PPI partners. Moreover, when the ratio of evolutionary rates (faster rate/slower rate) was computed for each of the duplicated pairs, the ratio for the duplicated pair sharing any PPI partners was significantly lower than that for the pair sharing no PPI partners. These results indicate that the duplicated gene pairs differentiate through the losses and/or gains of the PPI partners, resulting in a change in their evolutionary rates. In particular, we point out that the PPI losses for the duplicated gene products that are involved in the functional classes of 'transcription' and 'protein fate' have an impact on their evolutionary rates more than the PPI losses for others. (c) 2006 Elsevier B.V All rights reserved.

Rates of protein evolution are thought to be influenced by features of protein-protein interaction (PPI). However, the most important features of interaction for determining the evolutionary rate are poorly understood. Here, we consider four categories for PPIs in Saccharomyces cerevisiae. Properties we consider are the extent to which proteins interact with proteins of the same function or different function (DF) and the extent to which these interactions involve connections in the dense part or sparse part (SP) of a PPI network. Our findings are that proteins with DF-SP interactions evolve at the slowest rate of all the proteins examined.

The human genome sequence defines our inherent biological potential; the realization of the biology encoded therein requires knowledge of the function of each gene. Currently, our knowledge in this area is still limited. Several lines of investigation have been used to elucidate the structure and function of the genes in the human genome. Even so, gene prediction remains a difficult task, as the varieties of transcripts of a gene may vary to a great extent. We thus performed an exhaustive integrative characterization of 41,118 full-length cDNAs that capture the gene transcripts as complete functional cassettes, providing an unequivocal report of structural and functional diversity at the gene level. Our international collaboration has validated 21,037 human gene candidates by analysis of high-quality full-length cDNA clones through curation using unified criteria. This led to the identification of 5,155 new gene candidates. It also manifested the most reliable way to control the quality of the cDNA clones. We have developed a human gene database, called the H-Invitational Database (H-InvDB; http://www.h-invitational.jp/). It provides the following: integrative annotation of human genes, description of gene structures, details of novel alternative splicing isoforms, non-protein-coding RNAs, functional domains, subcellular localizations, metabolic pathways, predictions of protein three-dimensional structure, mapping of known single nucleotide polymorphisms (SNPs), identification of polymorphic microsatellite repeats within human genes, and comparative results with mouse full-length cDNAs. The H-InvDB analysis has shown that up to 4% of the human genome sequence (National Center for Biotechnology Information build 34 assembly) may contain misassembled or missing regions. We found that 6.5% of the human gene candidates (1,377 loci) did not have a good protein-coding open reading frame, of which 296 loci are strong candidates for nonprotein-coding RNA genes. In addition, among 72,027 uniquely mapped SNPs and insertions/deletions localized within human genes, 13,215 nonsynonymous SNPs, 315 nonsense SNPs, and 452 indels occurred in coding regions. Together with 25 polymorphic microsatellite repeats present in coding regions, they may alter protein structure, causing phenotypic effects or resulting in disease. The H-InvDB platform represents a substantial contribution to resources needed for the exploration of human biology and pathology.

Penicillum sp. 40, which can grow in an extremely acidic medium at pH 2.0 was screened from an acidic soil. This fungus produces xylanases when grown in a medium containing xylan as a sole carbon source. A major xylanase was purified from the culture supernatant of Penicillium sp. 40 and designated XynA. The molecular mass of XynA was estimated to be 25,000 by sodium dodecyl sulfate-polyacrylamide gel electrophoresis. XynA has an optimum pH at 2.0 and is stable in pH 2.0-5.0. Western blot analysis using anit-XynA antibody showed that XynA was induced by xylan and repressed by glucose. Also, its production was increased by an acidic medium. The gene encoding XynA (xynA) was isolated from the genomic library of Penicillium sp. 40. The structural part of xynA was found to be 721 bp. The nucleotide sequence of cDNA amplified by RT-PCR showed that the open reading frame of xynA was interrupted by a single intron which was 58 bp in size and encoded 221 amino acids. Direct N-terminal amino acid sequencing showed that the precursor of XynA had a signal peptide composed of 31 amino acids. The molecular mass caliculated from the deduced amino acid sequence of XynA is 20,713. This is lower than that estimated by gel electrophoresis, suggesting that XynA is a glycoprotein. The predicted amino acid sequence of XynA has strong similarity to other family11 xylanases from fungi.

Promoter activity of the genes encoding Taka-amylase A and phosphoglycerate kinase of Aspergillus oryzae was analyzed using Escherichia coli β-glucuronidase as a reporter in a shoyu-koji mold A. oryzae KBN616. Assay of the β-glucuronidase extracted from the mycelia of the transformants grown in the media containing different carbon sources suggests that expression of the Taka-amylase A gene was not induced by starch in A. oryzae KBN616. Analysis of proteins in the culture supernatant of A. oryzae KBN616 after cultivation in starch medium supports this result. The phosphoglycerate kinase gene of A. oryzae KBN616 was shown to be expressed constitutively in the medium containing glucose and starch.