Lessons learned from the dog genome

Extensive genetic resources and a high-quality genome sequence position the dog as an important model species for understanding genome evolution, population genetics and genes underlying complex phenotypic traits. Newly developed genomic resources have expanded our understanding of canine evolutionary history and dog origins. Domestication involved genetic contributions from multiple populations of gray wolves probably through backcrossing. More recently, the advent of controlled breeding practices has segregated genetic variability into distinct dog breeds that possess specific phenotypic traits. Consequently, genome-wide association and selective sweep scans now allow the discovery of genes underlying breed-specific characteristics. The dog is finally emerging as a novel resource for studying the genetic basis of complex traits, including behavior.     

Mapping four genes from human chromosome 4 to porcine chromosome 8 further develops the comparative map for an economically important chromosome of the swine genome

Because porcine chromosome (SSC) 8 has become the focal point of many efforts aimed at identifying quantitative trait loci affecting ovulation rate, genes distributed across human chromosome (HSA) 4 were physically mapped in the pig. A more refined comparative map of this region for these two species was produced. In this study, four genes were selected based on their location in the human genome, the availability of nucleotide sequence and their genomic organization. The genes selected were fibroblast growth factor basic (FGF2; HSA 4q25-27), gonadotropin releasing hormone receptor (GNRHR; HSA 4q13), phosphodiesterase 6 B (PDE6B; HSA 4p16.3) and aminopeptidase S (PEPS; HSA 4p11-q12). Genomic libraries were screened via PCR and clones were physically assigned using fluorescence in situ hybridization (FISH). These four genes from HSA 4 were physically mapped to SSC 8p2.3 (PDE6B), 8p1.1 (PEPS), 8q1.1-1.2 (GNRHR) and 8q2.2-2.4 (FGF2). These assignments provide additional benchmarks for the comparative map and help define the level of gene order conserved between HSA 4 and SSC 8.     

Whole genome sequencing and methylome analysis of the wild guinea pig

Background

DNA methylation is a heritable mechanism that acts in response to environmental changes, lifestyle and diseases by influencing gene expression in eukaryotes. Epigenetic studies of wild organisms are mandatory to understand their role in e.g. adaptational processes in the great variety of ecological niches. However, strategies to address those questions on a methylome scale are widely missing. In this study we present such a strategy and describe a whole genome sequence and methylome analysis of the wild guinea pig.

Results

We generated a full Wild guinea pig (Cavia aperea) genome sequence with enhanced coverage of methylated regions, benefiting from the available sequence of the domesticated relative Cavia porcellus. This new genome sequence was then used as reference to map the sequence reads of bisulfite treated Wild guinea pig sequencing libraries to investigate DNA-methylation patterns at nucleotide-specific level, by using our here described method, named ‘DNA-enrichment-bisulfite-sequencing’ (MEBS). The results achieved using MEBS matched those of standard methods in other mammalian model species. The technique is cost efficient, and incorporates both methylation enrichment results and a nucleotide-specific resolution even without a whole genome sequence available. Thus MEBS can be easily applied to extend methylation enrichment studies to a nucleotide-specific level.

Conclusions

The approach is suited to study methylomes of not yet sequenced mammals at single nucleotide resolution. The strategy is transferable to other mammalian species by applying the nuclear genome sequence of a close relative. It is therefore of interest for studies on a variety of wild species trying to answer evolutionary, adaptational, ecological or medical questions by epigenetic mechanisms.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1036) contains supplementary material, which is available to authorized users.     

Yeast genome sequencing: the power of comparative genomics

For decades, unicellular yeasts have been general models to help understand the eukaryotic cell and also our own biology. Recently, over a dozen yeast genomes have been sequenced, providing the basis to resolve several complex biological questions. Analysis of the novel sequence data has shown that the minimum number of genes from each species that need to be compared to produce a reliable phylogeny is about 20. Yeast has also become an attractive model to study speciation in eukaryotes, especially to understand molecular mechanisms behind the establishment of reproductive isolation. Comparison of closely related species helps in gene annotation and to answer how many genes there really are within the genomes. Analysis of non-coding regions among closely related species has provided an example of how to determine novel gene regulatory sequences, which were previously difficult to analyse because they are short and degenerate and occupy different positions. Comparative genomics helps to understand the origin of yeasts and points out crucial molecular events in yeast evolutionary history, such as whole-genome duplication and horizontal gene transfer(s). In addition, the accumulating sequence data provide the background to use more yeast species in model studies, to combat pathogens and for efficient manipulation of industrial strains.     

Criteria for gene identification and features of genome organization: analysis of 6.5 Mb of DNA sequence from human chromosome 21

To establish criteria for and the limitations of novel gene identification, to identify novel genes of potential relevance to Down Syndrome and to investigate features of genome organization, 6. 550kb. In total, 41 novel gene models were predicted, and for a subset of these, RT-PCR experiments helped to verify and refine the models, and were used to assess expression in early development and in adult brain regions of potential relevance to Down syndrome. Results suggest generally low and/or restricted patterns of expression, and also reveal examples of complex alternative processing, especially in brain, that may have important implications for regulation of protein function. Analysis of complete gene structures of the known genes identified a number of very large introns, a number of very short intergenic distances, and at least one potentially bi-directional promoter. At least 3/4 of known genes and 1/2 of predicted genes are associated with CpG islands. For novel genes, three cases of overlapping genes are predicted. Results of these analyses illustrate some of the complexities inherent in mammalian genome organization and some of the limitations of current sequence analysis technologies. They also doubled the number of potential genes within the region.     

A high-resolution comparative map between pig chromosome 17 and human chromosomes 4, 8, and 20: identification of synteny breakpoints

We report on the construction of a high-resolution comparative map of porcine chromosome 17 (SSC17) focusing on evolutionary breakpoints with human chromosomes. The comparative map shows high homology with human chromosome 20 but suggests more limited homologies with other human chromosomes. SSC17 is of particular interest in studies of chromosomal organization due to the presence of QTLs that affect meat quality and carcass composition. A total of 158 pig ESTs available in databases or developed by the Sino-Danish Pig Genome Sequencing Consortium were mapped using the INRA-University of Minnesota porcine radiation hybrid panel. The high-resolution map was further anchored by fluorescence in situ hybridization. This study confirmed the extensive conservation between SSC17 and HSA20 and enabled the gene order to be determined. The homology of the SSC17 pericentromeric region was extended to other human chromosomes (HSA4, HSA8) and the chromosomal breakpoint boundaries were accurately defined. In total 15 breakpoints were identified.     

The structure of a conserved region of porcine genome, represented in human genome by chromosome 17

Radiation mapping of nine genes (H3F3B, HLR1, MYL4, STAT5B, THRA1, TOP2A, MCP1, NF1, and MPO) to porcine chromosome 12 was carried out. Also, subchromosomal location of the NF1 gene along with the two loci containing the DNA sequences homologous to the DNA of the two human BAC clones was determined. The NF1 position was ascertained via microdissection of chromosome 12 with subsequent PCR amplification of the gene fragment with specific primers. BAC clones were mapped using FISH. Comparative analysis of the gene order in porcine chromosome 12 and in the homologous human chromosome 17 was performed. It was demonstrated that the gene orders in these chromosomes differed relative to the position of the MPO gene.     

Genomic sequencing and comparative analysis of Epstein-Barr virus genome isolated from primary nasopharyngeal carcinoma biopsy

Whether certain Epstein-Barr virus (EBV) strains are associated with pathogenesis of nasopharyngeal carcinoma (NPC) is still an unresolved question. In the present study, EBV genome contained in a primary NPC tumor biopsy was amplified by Polymerase Chain Reaction (PCR), and sequenced using next-generation (Illumina) and conventional dideoxy-DNA sequencing. The EBV genome, designated HKNPC1 (Genbank accession number JQ009376) is a type 1 EBV of approximately 171.5 kb. The virus appears to be a uniform strain in line with accepted monoclonal nature of EBV in NPC but is heterogeneous at 172 nucleotide positions. Phylogenetic analysis with the four published EBV strains, B95-8, AG876, GD1, and GD2, indicated HKNPC1 was more closely related to the Chinese NPC patient-derived strains, GD1 and GD2. HKNPC1 contains 1,589 single nucleotide variations (SNVs) and 132 insertions or deletions (indels) in comparison to the reference EBV sequence (accession number NC007605). When compared to AG876, a strain derived from Ghanaian Burkitt's lymphoma, we found 322 SNVs, of which 76 were non-synonymous SNVs and were shared amongst the Chinese GD1, GD2 and HKNPC1 isolates. We observed 88 non-synonymous SNVs shared only by HKNPC1 and GD2, the only other NPC tumor-derived strain reported thus far. Non-synonymous SNVs were mainly found in the latent, tegument and glycoprotein genes. The same point mutations were found in glycoprotein (BLLF1 and BALF4) genes of GD1, GD2 and HKNPC1 strains and might affect cell type specific binding. Variations in LMP1 and EBNA3B epitopes and mutations in Cp (11404 C>T) and Qp (50134 G>C) found in GD1, GD2 and HKNPC1 could potentially affect CD8(+) T cell recognition and latent gene expression pattern in NPC, respectively. In conclusion, we showed that whole genome sequencing of EBV in NPC may facilitate discovery of previously unknown variations of pathogenic significance.     

Mapping of the chicken cleft primary palate mutation on chromosome 11 and sequencing of the 4.9 Mb linked region

An embryonic lethal mutation in chicken named cleft primary palate (cpp) is inherited in an autosomal recessive mode and results in a severely truncated upper beak. In this study, genotyping and sequencing techniques were employed to advance our genetic and genomic knowledge of the mutation’s chromosomal location, candidate region and possible causative element using a congenic inbred line. Herein, the candidate region for the cpp developmental mutation was established as a ca. 5.1 Mb region of chicken chromosome 11 (GGA 11) through the use of a 600K Affymetrix SNP array. The SNPs identified from this array linked to cpp were used to genotype individuals from the congenic inbred line over several generations and thereby fine-map the causative region resulting in an approximately 200 kb size reduction. This candidate region (4.9 Mb) was sequenced via capture array in a cohort of 24 individuals, including carriers, mutants and their wild type (wt) siblings. Interestingly, the GGA 11 region for cpp encompasses the predicted centromere location and is thus unlikely to be highly disrupted by further recombination. Here we report on the variation unique to the cpp mutation, i.e. single-nucleotide variants and insertions or deletions. Although the candidate region contains several genes of interest with regard to the cpp phenotype, only one cpp-linked variant was predicted to have a significant physiological effect by causing a frameshift mutation in ESRP2, which has a role in tissue-specific splicing during development.     

Complete genome sequencing and comparative analysis of the linezolid-resistant Enterococcus faecalis strain DENG1

Genome level analysis of bacterial strains provides information on genetic composition and resistance mechanisms to clinically relevant antibiotics. To date, whole genome characterization of linezolid-resistant Enterococcus faecalis isolated in the clinic is lacking. In this study, we report the entire genome sequence, genomic characteristics and virulence factors of a pathogenic E. faecalis strain, DENG1. Our results showed considerable differences in genomic characteristics and virulence factors compared with other E. faecalis strains (V583 and OG1RF). The genome of this LZD-resistant E. faecalis strain can be used as a reference to study the mechanism of LZD resistance and the phylogenetic relationship of E. faecalis strains worldwide.     

Survey sequencing and comparative analysis of the elephant shark (Callorhinchus milii) genome

Owing to their phylogenetic position, cartilaginous fishes (sharks, rays, skates, and chimaeras) provide a critical reference for our understanding of vertebrate genome evolution. The relatively small genome of the elephant shark, Callorhinchus milii, a chimaera, makes it an attractive model cartilaginous fish genome for whole-genome sequencing and comparative analysis. Here, the authors describe survey sequencing (1.4× coverage) and comparative analysis of the elephant shark genome, one of the first cartilaginous fish genomes to be sequenced to this depth. Repetitive sequences, represented mainly by a novel family of short interspersed element–like and long interspersed element–like sequences, account for about 28% of the elephant shark genome. Fragments of approximately 15,000 elephant shark genes reveal specific examples of genes that have been lost differentially during the evolution of tetrapod and teleost fish lineages. Interestingly, the degree of conserved synteny and conserved sequences between the human and elephant shark genomes are higher than that between human and teleost fish genomes. Elephant shark contains putative four Hox clusters indicating that, unlike teleost fish genomes, the elephant shark genome has not experienced an additional whole-genome duplication. These findings underscore the importance of the elephant shark as a critical reference vertebrate genome for comparative analysis of the human and other vertebrate genomes. This study also demonstrates that a survey-sequencing approach can be applied productively for comparative analysis of distantly related vertebrate genomes.     

Lessons from comparative analysis of species-specific imprinted genes

Genomic imprinting is generally believed to be conserved in all mammals except for egg-laying monotremes, suggesting that it is closely related to placental and fetal growth. As expected, the imprinting status of most imprinted genes is conserved between mouse and human, and some are imprinted even in marsupials. On the other hand, a small number of genes were reported to exhibit species-specific imprinting that is not necessarily accounted for by either the placenta or conflict hypotheses. Since mouse and human represent a single, phylogenetically restricted clade in the mammalian class, a much broader comparison including mammals diverged earlier than rodents is necessary to fully understand the species-specificity and variation in evolution of genomic imprinting. Indeed, comparative analysis of a species-specific imprinted gene Impact using a broader range of mammals led us to propose an alternative dosage control hypothesis for the evolution of genomic imprinting.     

Systematic sequencing of the Escherichia coli genome: analysis of the 0-2.4 min region.

A contiguous 111,402-nucleotide sequence corresponding to the 0 to 2.4 min region of the E. coli chromosome was determined as a first step to complete structural analysis of the genome. The resulting sequence was used to predict open reading frames and to search for sequence similarity against the PIR protein database. A number of novel genes were found whose predicted protein sequences showed significant homology with known proteins from various organisms, including several clusters of genes similar to those involved in fatty acid metabolism in bacteria (e.g., betT, baiF) and higher organisms, iron transport (sfuA, B, C) in Serratia marcescens, and symbiotic nitrogen fixation or electron transport (fixA, B, C, X) in Azorhizobium caulinodans. In addition, several genes and IS elements that had been mapped but not sequenced (e.g., leuA, B, C, D) were identified. We estimate that about 90 genes are represented in this region of the chromosome with little spacer.     

Distinct evolutionary patterns of Oryza glaberrima deciphered by genome sequencing and comparative analysis

Here we present the genomic sequence of the African cultivated rice, Oryza glaberrima, and compare these data with the genome sequence of Asian cultivated rice, Oryza sativa. We obtained gene‐enriched sequences of O. glaberrima that correspond to about 25% of the gene regions of the O. sativa (japonica) genome by methylation filtration and subtractive hybridization of repetitive sequences. While patterns of amino acid changes did not differ between the two species in terms of the biochemical properties, genes of O. glaberrima generally showed a larger synonymous–nonsynonymous substitution ratio, suggesting that O. glaberrima has undergone a genome‐wide relaxation of purifying selection. We further investigated nucleotide substitutions around splice sites and found that eight genes of O. sativa experienced changes at splice sites after the divergence from O. glaberrima. These changes produced novel introns that partially truncated functional domains, suggesting that these newly emerged introns affect gene function. We also identified 2451 simple sequence repeats (SSRs) from the genomes of O. glaberrima and O. sativa. Although tri‐nucleotide repeats were most common among the SSRs and were overrepresented in the protein‐coding sequences, we found that selection against indels of tri‐nucleotide repeats was relatively weak in both African and Asian rice. Our genome‐wide sequencing of O. glaberrima and in‐depth analyses provide rice researchers not only with useful genomic resources for future breeding but also with new insights into the genomic evolution of the African and Asian rice species.     

一株新的产吡咯喹啉醌甲基营养菌全基因组测序和比较基因组学分析

【背景】由于甲基营养菌被发现的时间较短,而且可以生产吡咯喹啉醌(pyrroloquinoline quinone,PQQ)的甲基杆菌属细菌只有少数菌株的全基因组序列被公布,增加了该类细菌基因组学和生物代谢途径研究的难度。【目的】将本实验室筛选的PQQ生产菌经多种诱变方式处理,用于提高PQQ的发酵产量。对高产突变菌株进行全基因组解析,以探究甲基杆菌PQQ合成的分子机制,为后续分子育种提供序列背景信息。【方法】将野生型PQQ生产菌株进行紫外诱变、亚硝基胍诱变、甲基磺酸乙酯诱变、硫酸二乙酯诱变和紫外-氯化锂复合诱变。将突变菌株利用PromethION三代测序平台和MGISEQ-2000二代测序平台测序,然后进行组装和功能注释。组装得到的全基因组序列与模式菌株扭脱甲基杆菌AM1 (Methylobacterium extorquens AM1)进行比较基因组学分析。【结果】经11轮诱变获得一株突变菌株NI91,其PQQ产量为19.49mg/L,相较原始菌株提高44.91%。突变菌株NI91的基因组由一个5 409 262 bp的染色体组成,共编码4 957个蛋白,与模式菌株M. extorqu...     

Lessons learned from the syndrome of oculopalatal tremor

Journal of Computational Neuroscience - The syndrome of oculopalatal tremor (OPT) featuring the olivo-cerebellar hypersychrony leads to disabling pendular nystagmus and palatal myoclonus. This rare...     

DNA sequence comparative analysis of the 3pter-p26 region of human genome

Most proterminal regions of human chromosomes are GC-rich and gene-rich. Chromosome 3p is an exception. Its proterminal region is GC-poor, and likely to lose heterozy-gosity, thus causing a number of fatal diseases. Except one gap left in the telomeric position, the proterminal region of human chromosome 3p has been completely sequenced. The detailed sequence analysis showed: (i) the GC content of this region was 38.5%, being the lowest among all the human proterminal regions; (ii) this region contained 20 known genes and 22 predicted genes, with an average gene size of 97.5 kb. The previously mapped gene Cntn3 was not found in this region, but instead located in the 74 Mb position of human chromosome 3p; (iii) the interspersed repeats of this region were more active than the average level of the whole human genome, especially (TA)n, the content of which was twice the genome average; (iv) this region had a conserved synteny extending from 104.1 Mb to 112.4 Mb on the mouse chromosome 6, which was 8% larger in size, not in accordance with the whole genome comparison, probably because the 3pter-p26 region was more likely to lose neocleitides and its mouse synteny had more active interspersed repeats.     

DNA sequence comparative analysis of the 3pter-p26 region of human genome

Most proterminal regions of human chromosomes are GC-rich and gene-rich. Chromosome 3p is an exception. Its proterminal region is GC-poor, and likely to lose heterozygosity, thus causing a number of fatal diseases. Except one gap left in the telomeric position, the proterminal region of human chromosome 3p has been completely sequenced. The detailed sequence analysis showed: (i) the GC content of this region was 38.5%, being the lowest among all the human proterminal regions; (ii) this region contained 20 known genes and 22 predicted genes, with an average gene size of 97.5 kb. The previously mapped gene Cntn3 was not found in this region, but instead located in the 74 Mb position of human chromosome 3p; (iii) the interspersed repeats of this region were more active than the average level of the whole human genome, especially (TA)n, the content of which was twice the genome average; (iv) this region had a conserved synteny extending from 104.1 Mb to 112.4 Mb on the mouse chromosome 6, which was 8% larger in size, not in accordance with the whole genome comparison, probably because the 3pter-p26 region was more likely to lose neocleitides and its mouse synteny had more active interspersed repeats.