A network approach to investigating the key microbes and stability of gut microbial communities in a mouse neuropathic pain model.

BACKGROUND: Neuropathic pain is an abnormally increased sensitivity to pain, especially from mechanical or thermal stimuli. To date, the current pharmacological treatments for neuropathic pain are still unsatisfactory. The gut microbiota reportedly plays important roles in inducing neuropathic pain, so probiotics have also been used to treat it. However, the underlying questions around the interactions in and stability of the gut microbiota in a spared nerve injury-induced neuropathic pain model and the key microbes (i.e., the microbes that play critical roles) involved have not been answered. We collected 66 fecal samples over 2 weeks (three mice and 11 time points in spared nerve injury-induced neuropathic pain and Sham groups). The 16S rRNA gene was polymerase chain reaction amplified, sequenced on a MiSeq platform, and analyzed using a MOTHUR- UPARSE pipeline. RESULTS: Here we show that spared nerve injury-induced neuropathic pain alters gut microbial diversity in mice. We successfully constructed reliable microbial interaction networks using the Metagenomic Microbial Interaction Simulator (MetaMIS) and analyzed these networks based on 177,147 simulations. Interestingly, at a higher resolution, our results showed that spared nerve injury-induced neuropathic pain altered both the stability of the microbial community and the key microbes in a gut micro-ecosystem. Oscillospira, which was classified as a low-abundance and core microbe, was identified as the key microbe in the Sham group, whereas Staphylococcus, classified as a rare and non-core microbe, was identified as the key microbe in the spared nerve injury-induced neuropathic pain group. CONCLUSIONS: In summary, our results provide novel experimental evidence that spared nerve injury-induced neuropathic pain reshapes gut microbial diversity, and alters the stability and key microbes in the gut.

Characterization of benzimidazole resistance in Haemonchus contortus: integration of phenotypic, genotypic and proteomic approaches.

The field strain of Haemonchus contortus has a long history of anthelmintic resistance. To understand this phenomenon, the benzimidazole resistance profile was characterized from the Malaysian field-resistant strain by integrating phenotypic, genotypic and proteomic approaches. The faecal egg count reduction test (FECRT) demonstrated that benzimidazole resistance was at a critical level in the studied strain. The primary resistance mechanism was attributed to F200Y mutation in the isotype 1 ß-tubulin gene as revealed by AS-PCR and direct sequencing. Furthermore, the protein response of the resistant strain towards benzimidazole (i.e., albendazole) treatment was investigated via two-dimensional difference gel electrophoresis (2D-DIGE) and tandem liquid chromatography-mass spectrometry (LC-MS/MS). These investigations illustrated an up-regulation of antioxidant (i.e., ATP-binding region and heat-shock protein 90, superoxide dismutase) and metabolic (i.e., glutamate dehydrogenase) enzymes and down-regulation of glutathione S-transferase, malate dehydrogenase, and other structural and cytoskeletal proteins (i.e., actin, troponin T). Findings from this study are pivotal in updating the current knowledge on anthelmintic resistance and providing new insights into the defence mechanisms of resistant nematodes towards drug treatment.

A novel molecular and chromosomal lineage of the anthropophilic Simulium (Simulium) rufibasis subgroup (Diptera: Simuliidae) in Taiwan.

The Simulium rufibasis subgroup is one of three subgroups of the Simulium (Simulium) tuberosum species-group; it is characterized by a pair of clustered stout hairs on the ventral surface of female abdominal segment 7. A member of the S. rufibasis subgroup in Taiwan was investigated morphologically and genetically using the universal cytochrome c oxidase subunit I (COI) barcoding gene and polytene chromosomal banding pattern. The Taiwanese material is morphologically similar to S. rosliramlii Takaoka & Chen from Vietnam and represents the second species of the S. rufibasis subgroup known from Taiwan. It also represents a novel molecular lineage that is distinct from three other primary lineages identified as S. doipuiense, S. doipuiense/S. rufibasis, and S. weji previously reported from Thailand. The mitochondrial evidence for a distinct lineage in Taiwan is supported by chromosomal analysis, which revealed unique sex chromosomes. For nomenclatural stability, we associate the name S. arisanum Shiraki with the Taiwanese entity. Originally described from females from Taiwan, S. arisanum until now has remained an enigmatic species.

MetaMIS: a metagenomic microbial interaction simulator based on microbial community profiles.

BACKGROUND: The complexity and dynamics of microbial communities are major factors in the ecology of a system. With the NGS technique, metagenomics data provides a new way to explore microbial interactions. Lotka-Volterra models, which have been widely used to infer animal interactions in dynamic systems, have recently been applied to the analysis of metagenomic data. RESULTS: In this paper, we present the Lotka-Volterra model based tool, the Metagenomic Microbial Interacticon Simulator (MetaMIS), which is designed to analyze the time series data of microbial community profiles. MetaMIS first infers underlying microbial interactions from abundance tables for operational taxonomic units (OTUs) and then interprets interaction networks using the Lotka-Volterra model. We also embed a Bray-Curtis dissimilarity method in MetaMIS in order to evaluate the similarity to biological reality. MetaMIS is designed to tolerate a high level of missing data, and can estimate interaction information without the influence of rare microbes. For each interaction network, MetaMIS systematically examines interaction patterns (such as mutualism or competition) and refines the biotic role within microbes. As a case study, we collect a human male fecal microbiome and show that Micrococcaceae, a relatively low abundance OTU, is highly connected with 13 dominant OTUs and seems to play a critical role. MetaMIS is able to organize multiple interaction networks into a consensus network for comparative studies; thus we as a case study have also identified a consensus interaction network between female and male fecal microbiomes. CONCLUSIONS: MetaMIS provides an efficient and user-friendly platform that may reveal new insights into metagenomics data. MetaMIS is freely available at: https://sourceforge.net/projects/metamis/ .

Functional analysis for gut microbes of the brown tree frog (Polypedates megacephalus) in artificial hibernation.

BACKGROUND: Annual hibernation is an adaptation that helps many animals conserve energy during food shortage in winter. This natural cycle is also accompanied by a remodeling of the intestinal immune system, which is an aspect of host biology that is both influenced by, and can itself influence, the microbiota. In amphibians, the bacteria in the intestinal tract show a drop in bacterial counts. The proportion of pathogenic bacteria is greater in hibernating frogs than that found in nonhibernating frogs. This suggests that some intestinal gut microbes in amphibians can be maintained and may contribute to the functions in this closed ecosystem during hibernation. However, these results were derived from culture-based approaches that only covered a small portion of bacteria in the intestinal tract. METHODS: In this study, we use a more comprehensive analysis, including bacterial appearance and functional prediction, to reveal the global changes in gut microbiota during artificial hibernation via high-throughput sequencing technology. RESULTS: Our results suggest that artificial hibernation in the brown tree frog (Polypedates megacephalus) could reduce microbial diversity, and artificially hibernating frogs tend to harbor core operational taxonomic units that are rarely distributed among nonhibernating frogs. In addition, artificial hibernation increased significantly the relative abundance of the red-leg syndrome-related pathogenic genus Citrobacter. Furthermore, functional predictions via PICRUSt and Tax4Fun suggested that artificial hibernation has effects on metabolism, disease, signal transduction, bacterial infection, and primary immunodeficiency. CONCLUSIONS: We infer that artificial hibernation may impose potential effects on primary immunodeficiency and increase the risk of bacterial infections in the brown tree frog.

Examining the condition-specific antisense transcription in S. cerevisiae and S. paradoxus.

BACKGROUND: Recent studies have demonstrated that antisense transcription is pervasive in budding yeasts and is conserved between Saccharomyces cerevisiae and S. paradoxus. While studies have examined antisense transcripts of S. cerevisiae for inverse expression in stationary phase and stress conditions, there is a lack of comprehensive analysis of the conditional specific evolutionary characteristics of antisense transcription between yeasts. Here we attempt to decipher the evolutionary relationship of antisense transcription of S. cerevisiae and S. paradoxus cultured in mid log, early stationary phase, and heat shock conditions. RESULTS: Massively parallel sequencing of sequence strand-specific cDNA library was performed from RNA isolated from S. cerevisiae and S. paradoxus cells at mid log, stationary phase and heat shock conditions. We performed this analysis using a stringent set of sense ORF transcripts and non-coding antisense transcripts that were expressed in all the three conditions, as well as in both species. We found the divergence of the condition-specific anti-sense transcription levels is higher than that in condition-specific sense transcription levels, suggesting that antisense transcription played a potential role in adapting to different conditions. Furthermore, 43% of sense-antisense pairs demonstrated inverse expression in either stationary phase or heat shock conditions relative to the mid log conditions. In addition, a large part of sense-antisense pairs (67%), which demonstrated inverse expression, were highly conserved between the two species. Our results were also concordant with known functional analyses from previous studies and with the evidence from mechanistic experiments of role of individual genes. CONCLUSIONS: By performing a genome-scale computational analysis, we have tried to evaluate the role of antisense transcription in mediating sense transcription under different environmental conditions across and in two related yeast species. Our findings suggest that antisense regulation could control expression of the corresponding sense transcript via inverse expression under a range of different circumstances.

Identification of activated cryptic 5' splice sites using structure profiles and odds measure.

The activation of cryptic 5' splice sites (5' SSs) is often related to human hereditary diseases. The DNA-based mutation screening strategies are commonly used to recognize the cryptic 5' SSs, because features of the local DNA sequence can influence the choice of cryptic 5' SSs. To improve the identification of the cryptic 5' SSs, we developed a structure-based method, named SPO (structure profiles and odds measure), which combines two parameters, the structural feature derived from hydroxyl radical cleavage pattern and odds measure, to assess the likelihood of a cryptic 5' SS activation in competing with its paired authentic 5' SS. Compared to the current tools for identifying activated cryptic 5' SSs, the SPO algorithm achieves higher prediction accuracy than the other methods, including MaxEnt, MDD, Markov model, weight matrix model, Shapiro and Senapathy matrix, R(i) and ΔG. In addition, the predicted ΔSPO scores from the SPO algorithm exhibited a greater degree of correlation with the strength of cryptic 5' SS activation than that measured from the other seven methods. In conclusion, the SPO algorithm provides an optimal identification of cryptic 5' SSs, can be applied in designing mutagenesis experiments for various splicing events and may be helpful to investigate the relationship between structural variants and human hereditary diseases.

Role of antisense RNAs in evolution of yeast regulatory complexity.

Antisense RNAs (asRNAs) are known to regulate gene expression. However, a genome-wide mechanism of asRNA regulation is unclear, and there is no good explanation why partial asRNAs are not functional. To explore its regulatory role, we investigated asRNAs using an evolutionary approach, as genome-wide experimental data are limited. We found that the percentage of genes coupling with asRNAs in Saccharomyces cerevisiae is negatively associated with regulatory complexity and evolutionary age. Nevertheless, asRNAs evolve more slowly when their sense genes are under more complex regulation. Older genes coupling with asRNAs are more likely to demonstrate inverse expression, reflecting the role of these asRNAs as repressors. Our analyses provide novel evidence, suggesting a minor contribution of asRNAs in developing regulatory complexity. Although our results support the leaky hypothesis for asRNA transcription, our evidence also suggests that partial asRNAs may have evolved as repressors. Our study deepens the understanding of asRNA regulatory evolution.

Evolution of cis-regulatory elements in yeast de novo and duplicated new genes.

BACKGROUND: New genes that originate from non-coding DNA rather than being duplicated from parent genes are called de novo genes. Their short evolution time and lack of parent genes provide a chance to study the evolution of cis-regulatory elements in the initial stage of gene emergence. Although a few reports have discussed cis-regulatory elements in new genes, knowledge of the characteristics of these elements in de novo genes is lacking. Here, we conducted a comprehensive investigation to depict the emergence and establishment of cis-regulatory elements in de novo yeast genes. RESULTS: In a genome-wide investigation, we found that the number of transcription factor binding sites (TFBSs) in de novo genes of S. cerevisiae increased rapidly and quickly became comparable to the number of TFBSs in established genes. This phenomenon might have resulted from certain characteristics of de novo genes; namely, a relatively frequent gain of TFBSs, an unexpectedly high number of preexisting TFBSs, or lower selection pressure in the promoter regions of the de novo genes. Furthermore, we identified differences in the promoter architecture between de novo genes and duplicated new genes, suggesting that distinct regulatory strategies might be employed by genes of different origin. Finally, our functional analyses of the yeast de novo genes revealed that they might be related to reproduction. CONCLUSIONS: Our observations showed that de novo genes and duplicated new genes possess mutually distinct regulatory characteristics, implying that these two types of genes might have different roles in evolution.

MetaABC--an integrated metagenomics platform for data adjustment, binning and clustering.

SUMMARY: MetaABC is a metagenomic platform that integrates several binning tools coupled with methods for removing artifacts, analyzing unassigned reads and controlling sampling biases. It allows users to arrive at a better interpretation via series of distinct combinations of analysis tools. After execution, MetaABC provides outputs in various visual formats such as tables, pie and bar charts as well as clustering result diagrams. AVAILABILITY: MetaABC source code and documentation are available at http://bits2.iis.sinica.edu.tw/MetaABC/ CONTACT: dywang@gate.sinica.edu.tw; hktsai@iis.sinica.edu.tw SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Evidence of association between nucleosome occupancy and the evolution of transcription factor binding sites in yeast.

BACKGROUND: Divergence of transcription factor binding sites is considered to be an important source of regulatory evolution. The associations between transcription factor binding sites and phenotypic diversity have been investigated in many model organisms. However, the understanding of other factors that contribute to it is still limited. Recent studies have elucidated the effect of chromatin structure on molecular evolution of genomic DNA. Though the profound impact of nucleosome positions on gene regulation has been reported, their influence on transcriptional evolution is still less explored. With the availability of genome-wide nucleosome map in yeast species, it is thus desirable to investigate their impact on transcription factor binding site evolution. Here, we present a comprehensive analysis of the role of nucleosome positioning in the evolution of transcription factor binding sites. RESULTS: We compared the transcription factor binding site frequency in nucleosome occupied regions and nucleosome depleted regions in promoters of old (orthologs among Saccharomycetaceae) and young (Saccharomyces specific) genes; and in duplicate gene pairs. We demonstrated that nucleosome occupied regions accommodate greater binding site variations than nucleosome depleted regions in young genes and in duplicate genes. This finding was confirmed by measuring the difference in evolutionary rates of binding sites in sensu stricto yeasts at nucleosome occupied regions and nucleosome depleted regions. The binding sites at nucleosome occupied regions exhibited a consistently higher evolution rate than those at nucleosome depleted regions, corroborating the difference in the selection constraints at the two regions. Finally, through site-directed mutagenesis experiment, we found that binding site gain or loss events at nucleosome depleted regions may cause more expression differences than those in nucleosome occupied regions. CONCLUSIONS: Our study indicates the existence of different selection constraint on binding sites at nucleosome occupied regions than at the nucleosome depleted regions. We found that the binding sites have a different rate of evolution at nucleosome occupied and depleted regions. Finally, using transcription factor binding site-directed mutagenesis experiment, we confirmed the difference in the impact of binding site changes on expression at these regions. Thus, our work demonstrates the importance of composite analysis of chromatin and transcriptional evolution.

Reanalyze unassigned reads in Sanger based metagenomic data using conserved gene adjacency.

BACKGROUND: Investigation of metagenomes provides greater insight into uncultured microbial communities. The improvement in sequencing technology, which yields a large amount of sequence data, has led to major breakthroughs in the field. However, at present, taxonomic binning tools for metagenomes discard 30-40% of Sanger sequencing data due to the stringency of BLAST cut-offs. In an attempt to provide a comprehensive overview of metagenomic data, we re-analyzed the discarded metagenomes by using less stringent cut-offs. Additionally, we introduced a new criterion, namely, the evolutionary conservation of adjacency between neighboring genes. To evaluate the feasibility of our approach, we re-analyzed discarded contigs and singletons from several environments with different levels of complexity. We also compared the consistency between our taxonomic binning and those reported in the original studies. RESULTS: Among the discarded data, we found that 23.7 ± 3.9% of singletons and 14.1 ± 1.0% of contigs were assigned to taxa. The recovery rates for singletons were higher than those for contigs. The Pearson correlation coefficient revealed a high degree of similarity (0.94 ± 0.03 at the phylum rank and 0.80 ± 0.11 at the family rank) between the proposed taxonomic binning approach and those reported in original studies. In addition, an evaluation using simulated data demonstrated the reliability of the proposed approach. CONCLUSIONS: Our findings suggest that taking account of conserved neighboring gene adjacency improves taxonomic assignment when analyzing metagenomes using Sanger sequencing. In other words, utilizing the conserved gene order as a criterion will reduce the amount of data discarded when analyzing metagenomes.

Roles of trans and cis variation in yeast intraspecies evolution of gene expression.

Both cis and trans mutations contribute to gene expression divergence within and between species. We used Saccharomyces cerevisiae as a model organism to estimate the relative contributions of cis and trans variations to the expression divergence between a laboratory (BY) and a wild (RM) strain of yeast. We examined whether genes regulated by a single transcription factor (TF; single input module, SIM genes) or genes regulated by multiple TFs (multiple input module, MIM genes) are more susceptible to trans variation. Because a SIM gene is regulated by a single immediate upstream TF, the chance for a change to occur in its trans-acting factors would, on average, be smaller than that for a MIM gene. We chose 232 genes that exhibited expression divergence between BY and RM to test this hypothesis. We examined the expression patterns of these genes in a BY-RM coculture system and in a BY-RM diploid hybrid. We found that trans variation is far more important than cis variation for expression divergence between the two strains. However, because in 75% of the genes studied, cis variation has significantly contributed to expression divergence, cis change also plays a significant role in intraspecific expression evolution. Interestingly, we found that the proportion of genes with diverged expression between BY and RM is larger for MIM genes than for SIM genes; in fact, the proportion tends to increase with the number of transcription factors that regulate the gene. Moreover, MIM genes are, on average, subject to stronger trans effects than SIM genes, though the difference between the two types of genes is not conspicuous.

Structure and membership of gut microbial communities in multiple fish cryptic species under potential migratory effects.

The animal gut microbiota evolves quickly towards a complex community and plays crucial roles in its host's health and development. Factors such as host genetics and environmental changes are regarded as important for controlling the dynamics of animal gut microbiota. Migratory animals are an important group for studying how these factors influence gut microbiota because they experience strong environmental perturbations during migration. The commercially important grey mullet, Mugil cephalus, is a cosmopolitan species complex that display reproductive migration behaviour. There are three cryptic species of M. cephalus fish distributed across the Northwest Pacific, and their spawning sites overlap in the Taiwan Strait. This extraordinary natural occurrence makes the grey mullet an ideal model organism for exploring the nature of wild animal-gut microbiota relationships and interactions. This study investigates the diversity and structure of the gut microbial community in three cryptic M. cephalus species using 16S rRNA amplicon sequencing. Gut microbial compositions from adult and juvenile fish samples were analysed. Our results indicate that gut microbial communities within the grey mullet share a core microbiome dominated by Proteobacteria, Firmicutes and Actinobacteria. However, the structures of gut microbial communities were more distinct between adult mullet groups than they were between juvenile ones. Intriguingly, we found that adult fish that migrate to different geographical tracts harbour gut microbiota similar to historical records of seawater microflora, along their respective migration routes. This observation provides new insights into the interaction between aquatic animal gut microbial communities and the environments along their hosts' migratory routes, and thus warrants future study.

Maternal gut microbes shape the early-life assembly of gut microbiota in passerine chicks via nests.

BACKGROUND: Knowledge is growing on how gut microbiota are established, but the effects of maternal symbiotic microbes throughout early microbial successions in birds remain elusive. In this study, we examined the contributions and transmission modes of maternal microbes into the neonatal microbiota of a passerine, the zebra finch (Taeniopygia guttata), based on fostering experiments. RESULTS: Using 16S rRNA amplicon sequencing, we found that zebra finch chicks raised by their biological or foster parents (the society finch Lonchura striata domestica) had gut microbial communities converging with those of the parents that reared them. Moreover, source-tracking models revealed high contribution of zebra finches' oral cavity/crop microbiota to their chicks' early gut microbiota, which were largely replaced by the parental gut microbiota at later stages. The results suggest that oral feeding only affects the early stage of hatchling gut microbial development. CONCLUSIONS: Our study indicates that passerine chicks mainly acquire symbionts through indirect maternal transmission-passive environmental uptake from nests that were smeared with the intestinal and cloacal microbes of parents that raised them. Gut microbial diversity was low in hand-reared chicks, emphasizing the importance of parental care in shaping the gut microbiota. In addition, several probiotics were found in chicks fostered by society finches, which are excellent foster parents for other finches in bird farms and hosts of brood parasitism by zebra finches in aviaries; this finding implies that avian species that can transfer probiotics to chicks may become selectively preferred hosts of brood parasitism in nature. Video Abstract.

A network-based approach to deciphering a dynamic microbiome's response to a subtle perturbation.

Over the past decades, one main issue that has emerged in ecological and environmental research is how losses in biodiversity influence ecosystem dynamics and functioning, and consequently human society. Although biodiversity is a common indicator of ecosystem functioning, it is difficult to measure biodiversity in microbial communities exposed to subtle or chronic environmental perturbations. Consequently, there is a need for alternative bioindicators to detect, measure, and monitor gradual changes in microbial communities against these slight, chronic, and continuous perturbations. In this study, microbial networks before and after subtle perturbations by adding S. acidaminiphila showed diverse topological niches and 4-node motifs in which microbes with co-occurrence patterns played the central roles in regulating and adjusting the intertwined relationships among microorganisms in response to the subtle environmental changes. This study demonstrates that microbial networks are a good bioindicator for chronic perturbation and should be applied in a variety of ecological investigations.

The mitochondrial genome of the gymnosperm Cycas taitungensis contains a novel family of short interspersed elements, Bpu sequences, and abundant RNA editing sites.

The mtDNA of Cycas taitungensis is a circular molecule of 414,903 bp, making it 2- to 6-fold larger than the known mtDNAs of charophytes and bryophytes, but similar to the average of 7 elucidated angiosperm mtDNAs. It is characterized by abundant RNA editing sites (1,084), more than twice the number found in the angiosperm mtDNAs. The A + T content of Cycas mtDNA is 53.1%, the lowest among known land plants. About 5% of the Cycas mtDNA is composed of a novel family of mobile elements, which we designated as "Bpu sequences." They share a consensus sequence of 36 bp with 2 terminal direct repeats (AAGG) and a recognition site for the Bpu 10I restriction endonuclease (CCTGAAGC). Comparison of the Cycas mtDNA with other plant mtDNAs revealed many new insights into the biology and evolution of land plant mtDNAs. For example, the noncoding sequences in mtDNAs have drastically expanded as land plants have evolved, with abrupt increases appearing in the bryophytes, and then in the seed plants. As a result, the genomic organizations of seed plant mtDNAs are much less compact than in other plants. Also, the Cycas mtDNA appears to have been exempted from the frequent gene loss observed in angiosperm mtDNAs. Similar to the angiosperms, the 3 Cycas genes nad1, nad2, and nad5 are disrupted by 5 group II intron squences, which have brought the genes into trans-splicing arrangements. The evolutionary origin and invasion/duplication mechanism of the Bpu sequences in Cycas mtDNA are hypothesized and discussed.

Roles of cis- and trans-changes in the regulatory evolution of genes in the gluconeogenic pathway in yeast.

The yeast Saccharomyces cerevisiae proliferates rapidly in glucose-containing media. As glucose is getting depleted, yeast cells enter the transition from fermentative to nonfermentative metabolism, known as the diauxic shift, which is associated with major changes in gene expression. To understand the expression evolution of genes involved in the diauxic shift and in nonfermentative metabolism within species, a laboratory strain (BY), a wild strain (RM), and a clinical isolate (YJM) were used in this study. Our data showed that the RM strain enters into the diauxic shift approximately 1 h earlier than the BY strain with an earlier, higher induction of many key transcription factors (TFs) involved in the diauxic shift. Our sequence data revealed sequence variations between BY and RM in both coding and promoter regions of the majority of these TFs. The key TF Cat8p, a zinc-finger cluster protein, is required for the expression of many genes in gluconeogenesis under nonfermentative growth, and its derepression is mediated by deactivation of Mig1p. Our kinetic study of CAT8 expression revealed that CAT8 induction corresponded to the timing of glucose depletion in both BY and RM and CAT8 was induced up to 50- to 90-folds in RM, whereas only 20- to 30-folds in BY. In order to decipher the relative importance of cis- and trans-variations in expression divergence in the gluconeogenic pathway during the diauxic shift, we studied the expression levels of MIG1, CAT8, and their downstream target genes in the cocultures and in the hybrid diploids of BY-RM, BY-YJM, and RM-YJM and in strains with swapped promoters. Our data showed that the differences between BY and RM in the expression of MIG1, the upstream regulator of CAT8, were affected mainly by changes in cis-elements, though also by changes in trans-acting factors, whereas those of CAT8 and its downstream target genes were predominantly affected by changes in trans-acting factors.

Co-expression of neighboring genes in the zebrafish (Danio rerio) genome.

Neighboring genes in the eukaryotic genome have a tendency to express concurrently, and the proximity of two adjacent genes is often considered a possible explanation for their co-expression behavior. However, the actual contribution of the physical distance between two genes to their co-expression behavior has yet to be defined. To further investigate this issue, we studied the co-expression of neighboring genes in zebrafish, which has a compact genome and has experienced a whole genome duplication event. Our analysis shows that the proportion of highly co-expressed neighboring pairs (Pearson's correlation coefficient R>0.7) is low (0.24% approximately 0.67%); however, it is still significantly higher than that of random pairs. In particular, the statistical result implies that the co-expression tendency of neighboring pairs is negatively correlated with their physical distance. Our findings therefore suggest that physical distance may play an important role in the co-expression of neighboring genes. Possible mechanisms related to the neighboring genes' co-expression are also discussed.