Kun-Hyang Park

Novel mechanism of conjoined gene formation in the human genome

Recently, conjoined genes (CGs) have emerged as important genetic factors necessary for understanding the human genome. However, their formation mechanism and precise structures have remained mysterious. Based on a detailed structural analysis of 57 human CG transcript variants (CGTVs, discovered in this study) and all (833) known CGs in the human genome, we discovered that the poly(A) signal site from the upstream parent gene region is completely removed via the skipping or truncation of the final exon; consequently, CG transcription is terminated at the poly(A) signal site of the downstream parent gene. This result led us to propose a novel mechanism of CG formation: the complete removal of the poly(A) signal site from the upstream parent gene is a prerequisite for the CG transcriptional machinery to continue transcribing uninterrupted into the intergenic region and downstream parent gene. The removal of the poly(A) signal sequence from the upstream gene region appears to be caused by a deletion or truncation mutation in the human genome rather than post-transcriptional trans-splicing events. With respect to the characteristics of CG sequence structures, we found that intergenic regions are hot spots for novel exon creation during CGTV formation and that exons farther from the intergenic regions are more highly conserved in the CGTVs. Interestingly, many novel exons newly created within the intergenic and intragenic regions originated from transposable element sequences. Additionally, the CGTVs showed tumor tissue-biased expression. In conclusion, our study provides novel insights into the CG formation mechanism and expands the present concepts of the genetic structural landscape, gene regulation, and gene formation mechanisms in the human genome.

Genome analysis of the domestic dog (Korean Jindo) by massively parallel sequencing.

Although pioneering sequencing projects have shed light on the boxer and poodle genomes, a number of challenges need to be met before the sequencing and annotation of the dog genome can be considered complete. Here, we present the DNA sequence of the Jindo dog genome, sequenced to 45-fold average coverage using Illumina massively parallel sequencing technology. A comparison of the sequence to the reference boxer genome led to the identification of 4 675 437 single nucleotide polymorphisms (SNPs, including 3 346 058 novel SNPs), 71 642 indels and 8131 structural variations. Of these, 339 non-synonymous SNPs and 3 indels are located within coding sequences (CDS). In particular, 3 non-synonymous SNPs and a 26-bp deletion occur in the TCOF1 locus, implying that the difference observed in cranial facial morphology between Jindo and boxer dogs might be influenced by those variations. Through the annotation of the Jindo olfactory receptor gene family, we found 2 unique olfactory receptor genes and 236 olfactory receptor genes harbouring non-synonymous homozygous SNPs that are likely to affect smelling capability. In addition, we determined the DNA sequence of the Jindo dog mitochondrial genome and identified Jindo dog-specific mtDNA genotypes. This Jindo genome data upgrade our understanding of dog genomic architecture and will be a very valuable resource for investigating not only dog genetics and genomics but also human and dog disease genetics and comparative genomics.

Analysis of indel variations in the human disease-associated genes CDKN2AIP, WDR66, USP20 and OR7C2 in a Korean population.

Recently, the human genes CDKN2AIP, WDR66, USP20 and OR7C2 have emerged as important genetic factors that could be biologically associated with cancer, haematological diseases and olfactory dysfunction. In this regard, analysis of indel (insertion–deletion) variations in these genes at a population scale is of significant interest. In this study, we performed PCR amplifications and sequencing of the loci of these four genes using genomic DNA from 100 Korean individuals. We analysed the indels in these genes and made predictions about the functional consequences that are likely caused by the indels. We discovered a 3-bp deletion in CDKN2AIP, a 15-bp insertion in WDR66, a 3-bp deletion in USP20 and a 3-bp insertion in OR7C2 with indel allele frequencies of 33%, 79%, 99% and 100%, respectively. The results of 3-D structural predictions and analysis of the proteins encoded by the four genes showed that in-frame amino acid deletions or insertions caused by these indels could result in hindrance of the formation of the optimal functional structures of these proteins, which could affect their functions. In particular, the indels in USP20 and OR7C2 have severely biased frequencies in the investigated Korean population, which may reflect a high susceptibility to certain cancer types and a biased preference to a kind of smell in Korean individuals. Additionally, our results could help to identify therapeutic targets for treating possible genetic diseases in individuals possessing homozygous genotypes for these indels in future studies.

Major chimpanzee-specific structural changes in sperm development-associated genes.

A comprehensive analysis of transcriptional structures of chimpanzee sperm development-associated genes is of significant interest for deeply understanding sperm development and male reproductive process. In this study, we sequenced 7,680 clones from a chimpanzee testis full-length cDNA library and obtained 1,933 nonredundant high-quality full-length cDNA sequences. Comparative analysis between human and chimpanzee showed that 78 sperm development-associated genes, most of which were yet uncharacterized, had undergone severe structural changes (mutations at the start/stop codons, INDELs, alternative splicing variations and fusion forms) on genomic and transcript levels throughout chimpanzee evolution. Specifically, among the 78 sperm development-associated genes, 39 including ODF2, UBC, and CD59 showed markedly chimpanzee-specific structural changes. Through dN/dS analysis, we found that 56 transcripts (including seven sperm development-associated genes) had values of greater than one when comparing human and chimpanzee DNA sequences, whereas the values were less than one when comparing humans and orangutans. Gene ontology annotation and expression profiling showed that the chimpanzee testis transcriptome was enriched with genes that are associated with chimpanzee male germ cell development. Taken together, our study provides the first comprehensive molecular evidence that many chimpanzee sperm development-associated genes had experienced severe structural changes over the course of evolution on genomic and transcript levels.

Gene selection tool (GST): A R-based tool for genetic disorders based on the sliding-window proportion test using whole-exome sequencing data

Whole-exome sequencing (WES) can identify causative mutations in hereditary diseases. However, WES data might have a large candidate variant list, including false positives. Moreover, in families, it is more difficult to select disease-associated variants because many variants are shared among members. To reduce false positives and extract accurate candidates, we used a multilocus variant instead of a single-locus variant (SNV). We set up a specific window to analyze the multilocus variant and devised a sliding-window approach to observe all variants. We developed the gene selection tool (GST) based on proportion tests for linkage analysis using WES data. This tool is R program coded and has high sensitivity. We tested our code to find the gene for hereditary spastic paraplegia using SNVs from a specific family and identified the gene known to cause the disease in a significant gene list. The list identified other genes that might be associated with the disease.

UbC gene allele frequency in Korean population and novel UbC mosaic repeat unit formation

The genomic structural organization of human UbC CDS repeat units could be representative of concerted evolution. The structure of the UbC gene and its repeat unit number frequency at scales of different human ethnic populations remain to be sufficiently determined. In this study, we performed comparative analysis of UbC CDS regions in genomes from 140 Korean individuals. We found that the UbC gene allele types 9, 8 and 7 are present in the Korean population in proportions of 97.1%, 0.4% and 2.5%, respectively. Interestingly, we discovered that the allele types 7 and 8 harbor the novel UbC gene mosaic repeat units 3∧5 (combined between sequence parts derived from standard repeat units 3 and 5) and 8∧9 (combined between sequence parts derived from standard repeat units 8 and 9) within their sequence structures, respectively. Our analysis showed that the novel mosaic repeat unit 3^5 lacks the highly human-specific amino acid S38, implying a functional consequence. These results suggest that the genomic organization of UbC repeat units is still undergoing dynamic structural changes due to concerted evolution through unequal crossing-over. Our results could represent valuable data for future investigations related to treating genetic diseases caused by UbC gene mutations and variations.