Pigs have been one of the most important sources of meat for humans, and their productivity has been substantially improved by recent strong selection. Here, we present whole-genome resequencing analyses of 55 pigs of five breeds representing Korean native pigs, wild boar and three European origin breeds. 1,673.1 Gb of sequence reads were mapped to the Swine reference assembly, covering ∼99.2% of the reference genome, at an average of ∼11.7-fold coverage. We detected 20,123,573 single-nucleotide polymorphisms (SNPs), of which 25.5% were novel. We extracted 35,458 of non-synonymous SNPs in 9,904 genes, which may contribute to traits of interest. The whole SNP sets were further used to access the population structures of the breeds, using multiple methodologies, including phylogenetic, similarity matrix, and population structure analysis. They showed clear population clusters with respect to each breed. Furthermore, we scanned the whole genomes to identify signatures of selection throughout the genome. The result revealed several promising loci that might underlie economically important traits in pigs, such as the CLDN1 and TWIST1 genes. These discoveries provide useful genomic information for further study of the discrete genetic mechanisms associated with economically important traits in pigs.

Female fertility is a highly regulated process involving the synchronized activities of multiple tissues. The underlying genomic regulation of the tissue synchronization is poorly understood. To understand this better we investigated the transcriptomes of the porcine ovary, endometrium, and oviduct at days 0, 3, 6, 9, 12, 15, or 18 of the oestrous cycle. We analysed the transcriptome profiles of the individual tissues and focus on the bridging genes shared by two or more tissues. The three tissue-networks were connected forming a triangular shape. We identified 65 bridging genes with a high level of connectivity to all other genes in the network. The expression levels showed negative correlations between the ovary and the other two tissues, and low correlations between endometrium and oviduct. The main functional annotations involved biosynthesis of steroid hormones, cell-to-cell adhesion, and cell apoptosis, suggesting that regulation of steroid hormone synthesis and tissue viability are major regulatory mechanisms.

This study provides high-quality variation data of diverse radish genotypes. Genome-wide SNP comparison along with RNA-seq analysis identified candidate genes related to domestication that have potential as trait-related markers for genetics and breeding of radish. Radish (Raphanus sativus L.) is an annual root vegetable crop that also encompasses diverse wild species. Radish has a long history of domestication, but the origins and selective sweep of cultivated radishes remain controversial. Here, we present comprehensive whole-genome resequencing analysis of radish to explore genomic variation between the radish genotypes and to identify genetic bottlenecks due to domestication in Asian cultivars. High-depth resequencing and multi-sample genotyping analysis of ten cultivated and seven wild accessions obtained 4.0 million high-quality homozygous single-nucleotide polymorphisms (SNPs)/insertions or deletions. Variation analysis revealed that Asian cultivated radish types are closely related to wild Asian accessions, but are distinct from European/American cultivated radishes, supporting the notion that Asian cultivars were domesticated from wild Asian genotypes. SNP comparison between Asian genotypes identified 153 candidate domestication regions (CDRs) containing 512 genes. Network analysis of the genes in CDRs functioning in plant signaling pathways and biochemical processes identified group of genes related to root architecture, cell wall, sugar metabolism, and glucosinolate biosynthesis. Expression profiling of the genes during root development suggested that domestication-related selective advantages included a main taproot with few branched lateral roots, reduced cell wall rigidity and favorable taste. Overall, this study provides evolutionary insights into domestication-related genetic selection in radish as well as identification of gene candidates with the potential to act as trait-related markers for background selection of elite lines in molecular breeding.

DNA methylation is an epigenetic mark that plays an essential role in regulating gene expression. CpG islands are DNA methylations regions in promoters known to regulate gene expression through transcriptional silencing of the corresponding gene. DNA methylation at CpG islands is crucial for gene expression and tissue-specific processes. At the current time, a limited number of studies have reported on gene expression associated with DNA methylation in diverse adult tissues at the genome-wide level. Expression levels are rarely affected by DNA methylation in normal adult tissues; however, statistical differences in gene expression level correlated with DNA methylation have recently been revealed. In this study, we examined 20 pairs of DNA methylomes and transcriptomes from RNA-seq and reduced representation bisulfite sequencing (RRBS) data using adult Ogye chicken tissues. A total of 3,133 CpG islands were identified from 20 tissue data in a single chicken sample which could affect downstream genes. Analyzing these CpG island and gene pairs, 121 significant units were statistically correlated. Among them, six genes ( CLDN3, DECR2, EVA1B, NME4, NTSR1 , and XPNPEP2 ) were highly significantly changed by altered DNA methylation. Finally, our data demonstrated how DNA methylation correlated to gene expression in normal adult tissues. Our source codes can be found at https://github.com/wjlim/correlation-between-rna-seq-and-RRBS .

Human mitochondrial transcription factor A (TFAM) has been implicated in promoting tumor growth and invasion. TFAM activates mitochondrial DNA (mtDNA) transcription, and affects nuclear gene expression through mitochondrial retrograde signaling. In this study, we investigated the effects of TFAM depletion on the morphology and transcriptome of MKN45 gastric cancer cells. Morphology alteration became visible at 12 h after TFAM knockdown: the proportion of growth-arrested polygonal cells versus oval-shaped cells increased, reaching a half-maximum at 24 h and a near-maximum at 36 h. TFAM knockdown upregulated four genes and downregulated six genes by more than threefold at 24 h and similarly at 48 h. Among them, the knockdown of CFAP65 (cilia and flagella associated protein 65) or PCK1 (cytoplasmic phosphoenolpyruvate carboxykinase) rescued the effects of TFAM depletion on cell morphology and proliferation. PCK1 was found to act downstream of CFAP65 in calcium-mediated retrograde signaling. Furthermore, mtDNA depletion by 2′,3′-dideoxycytidine was sufficient for induction of CFAP65 and PCK1 expression and inhibition of cell proliferation, but oxidative phosphorylation blockade or mitochondrial membrane potential depolarization was not. Thus, the TFAM–mtDNA–calcium–CFAP65–PCK1 axis participates in mitochondrial retrograde signaling, affecting tumor cell differentiation and proliferation.

Indigenous breeds develop their own genomic characteristics by adapting to local environments or cultures over long periods of time. Most of them are not particularly productive in commercial terms, but they have abilities to survive in harsh environments or tolerate to specific diseases. Their adaptive characteristics play an important role as genetic materials for improving commercial breeds. As a step toward this goal, we analyzed the genome of Korean indigenous goats within 10 goat breeds. We collected 136 goat individuals by sequencing 46 new goats and employing 90 publicly available goats. Our whole-genome data was comprised of three indigenous breeds (Korean indigenous goat, Iranian indigenous goat, and Moroccan indigenous goat; n = 29, 18, 20), six commercial breeds (Saanen, Boer, Anglo-Nubian, British Alpine, Alpine, and Korean crossbred; n = 16, 11, 5, 5, 2, 13), and their ancestral species ( Capra aegagrus ; n = 17). We identified that the Iranian indigenous goat and the Moroccan indigenous goat have relatively similar genomic characteristics within a large category of genomic diversity but found that the Korean indigenous goat has unique genomic characteristics distinguished from the other nine breeds. Through population analysis, we confirmed that these characteristics have resulted from a near-isolated environment with strong genetic drift. The Korean indigenous goat experienced a severe genetic bottleneck upon entering the Korean Peninsula about 2,000 years ago, and has subsequently rarely experienced genetic interactions with other goat breeds. From selection analysis and gene-set enrichment analysis, we revealed selection signals for Salmonella infection and cardiomyopathy in the genome of the Korean indigenous goat. These adaptive characteristics were further identified with genomic-based evidence. We uncovered genomic regions of selective sweeps in the LBP and BPI genes ( Salmonella infection) and the TTN and ITGB6 genes (cardiomyopathy), among several candidate genes. Our research presents unique genomic characteristics and distinctive selection signals of the Korean indigenous goat based on the extensive comparison. Although the adaptive traits require further validation through biological experiments, our findings are expected to provide a direction for future biodiversity conservation strategies and to contribute another option to genomic-based breeding programmes for improving the viability of Capra hircus .

Genome-wide association studies (GWAS) have enabled the discovery of candidate markers that play significant roles in various complex traits in plants. Recently, with increased interest in the search for candidate markers, studies on epistatic interactions between single nucleotide polymorphism (SNP) markers have also increased, thus enabling the identification of more candidate markers along with GWAS on single-variant-additive-effect. Here, we focused on the identification of candidate markers associated with flowering time in soybean (Glycine max). A large population of 2,662 cultivated soybean accessions was genotyped using the 180k Axiom ® SoyaSNP array, and the genomic architecture of these accessions was investigated to confirm the population structure. Then, GWAS was conducted to evaluate the association between SNP markers and flowering time. A total of 93 significant SNP markers were detected within 59 significant genes, including E1 and E3, which are the main determinants of flowering time. Based on the GWAS results, multilocus epistatic interactions were examined between the significant and non-significant SNP markers. Two significant and 16 non-significant SNP markers were discovered as candidate markers affecting flowering time via interactions with each other. These 18 candidate SNP markers mapped to 18 candidate genes including E1 and E3, and the 18 candidate genes were involved in six major flowering pathways. Although further biological validation is needed, our results provide additional information on the existing flowering time markers and present another option to marker-assisted breeding programs for regulating flowering time of soybean.

The crab-eating monkey is widely used in biomedical research for pharmacological experiments. Epigenetic regulation in the brain regions of primates involves complex patterns of DNA methylation. Previous studies of methylated CpG-binding domains using microarray technology or peak identification of sequence reads mostly focused on developmental stages or disease, rather than normal brains. To identify correlations between gene expression and DNA methylation levels that may be related to transcriptional regulation, we generated RNA-seq and whole-genome bisulfite sequencing data from seven different brain regions from a single crab-eating monkey. We identified 92 genes whose expression levels were significantly correlated, positively or negatively, with DNA methylation levels. Among them, 11 genes exhibited brain region-specific characteristics, and their expression patterns were strongly correlated with DNA methylation level. Nine genes ( SLC2A5 , MCM5 , DRAM1 , TTC12 , DHX40 , COR01A , LRAT , FLVCR2 , and PTER ) had effects on brain and eye function and development, and two ( LHX6 and MEST ) were previously identified as genes in which DNA methylation levels change significantly in the promoter region and are therefore considered brain epigenetic markers. Furthermore, we characterized DNA methylation of repetitive elements at the whole genome through repeat annotation at single-base resolution. Our results reveal the diverse roles of DNA methylation at single-base resolution throughout the genome and reflect the epigenetic variations in adult brain tissues.