Genetic diversity and signatures of selection in BoHuai goat revealed by whole-genome sequencing.

BACKGROUND: Cross breeding is an important way to improve livestock performance. As an important livestock and poultry resource in Henan Province of China, Bohuai goat was formed by crossing Boer goat and Huai goat. After more than 20 years of breeding, BoHuai goats showed many advantages, such as fast growth, good reproductive performance, and high meat yield. In order to better develop and protect Bohuai goats, we sequenced the whole genomes of 30 BoHuai goats and 5 Huai goats to analyze the genetic diversity, population structure and genomic regions under selection of BoHuai goat. Furthermore, we used 126 published genomes of world-wide goat to characterize the genomic variation of BoHuai goat. RESULTS: The results showed that the nucleotide diversity of BoHuai goats was lower and the degree of linkage imbalance was higher than that of other breeds. The analysis of population structure showed that BoHuai goats have obvious differences from other goat breeds. In addition, the BoHuai goat is more closely related to the Boer goat than the Huai goat and is highly similar to the Boer goat. Group by selection signal in the BoHuai goat study, we found that one region on chromosome 7 shows a very strong selection signal, which suggests that it could well be the segment region under the intense artificial selection results. Through selective sweeps, we detected some genes related to important traits such as lipid metabolism (LDLR, STAR, ANGPTL8), fertility (STAR), and disease resistance (CD274, DHPS, PDCD1LG2). CONCLUSION: In this paper, we elucidated the genomic variation, ancestry composition, and selective signals related to important economic traits in BoHuai goats. Our studies on the genome of BoHuai goats will not only help to understand the characteristics of the crossbred but also provide a basis for the improvement of cross-breeding programs.

De Novo Assembly of 20 Chicken Genomes Reveals the Undetectable Phenomenon for Thousands of Core Genes on Microchromosomes and Subtelomeric Regions.

The gene numbers and evolutionary rates of birds were assumed to be much lower than those of mammals, which is in sharp contrast to the huge species number and morphological diversity of birds. It is, therefore, necessary to construct a complete avian genome and analyze its evolution. We constructed a chicken pan-genome from 20 de novo assembled genomes with high sequencing depth, and identified 1,335 protein-coding genes and 3,011 long noncoding RNAs not found in GRCg6a. The majority of these novel genes were detected across most individuals of the examined transcriptomes but were seldomly measured in each of the DNA sequencing data regardless of Illumina or PacBio technology. Furthermore, different from previous pan-genome models, most of these novel genes were overrepresented on chromosomal subtelomeric regions and microchromosomes, surrounded by extremely high proportions of tandem repeats, which strongly blocks DNA sequencing. These hidden genes were proved to be shared by all chicken genomes, included many housekeeping genes, and enriched in immune pathways. Comparative genomics revealed the novel genes had 3-fold elevated substitution rates than known ones, updating the knowledge about evolutionary rates in birds. Our study provides a framework for constructing a better chicken genome, which will contribute toward the understanding of avian evolution and the improvement of poultry breeding.

Assessing genomic diversity and signatures of selection in Pinan cattle using whole-genome sequencing data.

BACKGROUND: Crossbreeding is an important way to improve production beef cattle performance. Pinan cattle is a new hybrid cattle obtained from crossing Piedmontese bulls with Nanyang cows. After more than 30 years of cross-breeding, Pinan cattle show a variety of excellent characteristics, including fast growth, early onset of puberty, and good meat quality. In this study, we analyzed the genetic diversity, population structure, and genomic region under the selection of Pinan cattle based on whole-genome sequencing data of 30 Pinan cattle and 169 published cattle genomic data worldwide.  RESULTS: Estimating ancestry composition analysis showed that the composition proportions for our Pinan cattle were mainly Piedmontese and a small amount of Nanyang cattle. The analyses of nucleotide diversity and linkage disequilibrium decay indicated that the genomic diversity of Pinan cattle was higher than that of European cattle and lower than that of Chinese indigenous cattle. De-correlated composite of multiple selection signals, which combines four different statistics including θπ, CLR, FST, and XP-EHH, was computed to detect the signatures of selection in the Pinan cattle genome. A total of 83 genes were identified, affecting many economically important traits. Functional annotation revealed that these selected genes were related to immune (BOLA-DQA2, BOLA-DQB, LSM14A, SEC13, and NAALADL2), growth traits (CYP4A11, RPL26, and MYH10), embryo development (REV3L, NT5E, CDX2, KDM6B, and ADAMTS9), hornless traits (C1H21orf62), and climate adaptation (ANTXR2). CONCLUSION: In this paper, we elucidated the genomic characteristics, ancestry composition, and selective signals related to important economic traits in Pinan cattle. These results will provide the basis for further genetic improvement of Pinan cattle and reference for other hybrid cattle related studies.

Transcriptome-wide study revealed m6A regulation of embryonic muscle development in Dingan goose (Anser cygnoides orientalis).

BACKGROUND: The number of myofiber is determined during the embryonic stage and does not increase during the postnatal period for birds, including goose. Thus, muscle production of adult goose is pre-determined during embryogenesis. Previous studies show N6-methyladenosine (m6A) is an important regulator for skeletal muscle development of birds and miRNAs play as a co-regulator for the skeletal muscle development in birds. Herein, we sequenced m6A and miRNA transcriptomes to investigate the profiles of m6A and their potential mechanism of regulating breast muscle development in Dingan Goose. RESULTS: We selected embryonic 21th day (E21) and embryonic 30th day (E30) to investigate the roles of transcriptome-wide m6A modification combining with mRNAs and miRNAs in goose breast muscle development. In this study, m6A peaks were mainly enriched in coding sequence (CDS) and start codon and397 genes were identified as differentially methylated genes (DMGs). GO and KEGG analysis showed that DMGs were highly related to cellular and metabolic process and that most DMGs were enriched in muscle-related pathways including Wnt signaling pathway, mTOR signaling and FoxO signaling pathway. Interestingly, a negative correlation between m6A methylation level and mRNA abundance was found through the analysis of m6A-RNA and RNA-seq data. Besides, we found 26 muscle-related genes in 397 DMGs. We also detected 228 differentially expressed miRNAs (DEMs), and further found 329 genes shared by the target genes of DEMs and DMGs (m6A-miRNA-genes), suggesting a tightly relationship between DEMs and DMGs. Among the m6A-miRNA-genes, we found 10 genes are related to breast muscle development. We further picked out an m6A-miRNA-gene, PDK3, from the 10 genes to visualize it and the result showed differentially methylated peaks on the mRNA transcript consistent with our m6A-seq results. CONCLUSION: GO and KEGG of DMGs between E21 and E30 showed most DMGs were muscle-related. In total, 228 DEMs were found, and the majority of DMGs were overlapped with the targets of DEGs. The differentially methylated peaks along with an m6A-miRNA-gene, PDK3, showed the similar results with m6A-seq results. Taken together, the results presented here provide a reference for further investigation of embryonic skeletal muscle development mechanism in goose.

Assessing genomic diversity and signatures of selection in Jiaxian Red cattle using whole-genome sequencing data.

BACKGROUND: Native cattle breeds are an important source of genetic variation because they might carry alleles that enable them to adapt to local environment and tough feeding conditions. Jiaxian Red, a Chinese native cattle breed, is reported to have originated from crossbreeding between taurine and indicine cattle; their history as a draft and meat animal dates back at least 30 years. Using whole-genome sequencing (WGS) data of 30 animals from the core breeding farm, we investigated the genetic diversity, population structure and genomic regions under selection of Jiaxian Red cattle. Furthermore, we used 131 published genomes of world-wide cattle to characterize the genomic variation of Jiaxian Red cattle. RESULTS: The population structure analysis revealed that Jiaxian Red cattle harboured the ancestry with East Asian taurine (0.493), Chinese indicine (0.379), European taurine (0.095) and Indian indicine (0.033). Three methods (nucleotide diversity, linkage disequilibrium decay and runs of homozygosity) implied the relatively high genomic diversity in Jiaxian Red cattle. We used θπ, CLR, FST and XP-EHH methods to look for the candidate signatures of positive selection in Jiaxian Red cattle. A total number of 171 (θπ and CLR) and 17 (FST and XP-EHH) shared genes were identified using different detection strategies. Functional annotation analysis revealed that these genes are potentially responsible for growth and feed efficiency (CCSER1), meat quality traits (ROCK2, PPP1R12A, CYB5R4, EYA3, PHACTR1), fertility (RFX4, SRD5A2) and immune system response (SLAMF1, CD84 and SLAMF6). CONCLUSION: We provide a comprehensive overview of sequence variations in Jiaxian Red cattle genomes. Selection signatures were detected in genomic regions that are possibly related to economically important traits in Jiaxian Red cattle. We observed a high level of genomic diversity and low inbreeding in Jiaxian Red cattle. These results provide a basis for further resource protection and breeding improvement of this breed.

m6A and miRNA jointly regulate the development of breast muscles in duck embryonic stages.

N6-methyladenosine (m6A) is an abundant internal mRNA modification and plays a crucial regulatory role in animal growth and development. In recent years, m6A modification has been found to play a key role in skeletal muscles. However, whether m6A modification contributes to embryonic breast muscle development of Pekin ducks has not been explored. To explore the role of m6A in embryonic breast muscle development of ducks, we performed m6A sequencing and miRNA sequencing for the breast muscle of duck embryos on the 19th (E19) and 27th (E27) days. A total of 12,717 m6A peaks were identified at E19, representing a total of 7,438 gene transcripts. A total of 14,703 m6A peaks were identified, which overlapped with the transcripts of 7,753 genes at E27. Comparing E19 and E27, we identified 2,347 differential m6A peaks, which overlapped with 1,605 m6A-modified genes (MMGs). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses revealed that MMGs were enriched in multiple muscle- or fat-related pathways, which was also revealed from our analysis of differentially expressed genes (DEGs). Conjoint analysis of m6A-seq and RNA-seq data showed that pathways related to ß-oxidation of fatty acids and skeletal muscle development were significantly enriched, suggesting that m6A modification is involved in the regulation of fat deposition and skeletal muscle development. There were 90 upregulated and 102 downregulated miRNAs identified between the E19 and E27 stages. Through overlapping analysis of genes shared by MMGs and DEGs and the targets of differentially expressed miRNAs (DEMs), we identified six m6A-mRNA-regulated miRNAs. Finally, we found that m6A modification can regulate fat deposition and skeletal muscle development. In conclusion, our results suggest that m6A modification is a key regulator for embryonic breast muscle development and fat deposition of ducks by affecting expressions of mRNAs and miRNAs. This is the first study to comprehensively characterize the m6A patterns in the duck transcriptome. These data provide a solid basis for future work aimed at determining the potential functional roles of m6A modification in adipose deposition and muscle growth.

An atlas of CNV maps in cattle, goat and sheep.

Copy number variation (CNV) is the most prevalent type of genetic structural variation that has been recognized as an important source of phenotypic variation in humans, animals and plants. However, the mechanisms underlying the evolution of CNVs and their function in natural or artificial selection remain unknown. Here, we generated CNV region (CNVR) datasets which were diverged or shared among cattle, goat, and sheep, including 886 individuals from 171 diverse populations. Using 9 environmental factors for genome-wide association study (GWAS), we identified a series of candidate CNVRs, including genes relating to immunity, tick resistance, multi-drug resistance, and muscle development. The number of CNVRs shared between species is significantly higher than expected (P<0.00001), and these CNVRs may be more persist than the single nucleotide polymorphisms (SNPs) shared between species. We also identified genomic regions under long-term balancing selection and uncovered the potential diversity of the selected CNVRs close to the important functional genes. This study provides the evidence that balancing selection might be more common in mammals than previously considered, and might play an important role in the daily activities of these ruminant species.

LncRNA-CTS promotes metastasis and epithelial-to-mesenchymal transition through regulating miR-505/ZEB2 axis in cervical cancer.

Cervical carcinoma (CC) is the second most common cancer in females. In order to improve current anti-metastasis strategies for CC, it is important to improve our understanding of the mechanisms involved in epithelial-to-mesenchymal transition (EMT). This study aimed to elucidate the potential role of a novel long non-coding RNA (lncRNA)-CTS and the mechanisms underlying EMT in CC. The expression levels of lncRNA-CTS and miR-505 were detected using quantitative reverse transcriptase polymerase chain reaction in CC specimens and cells (HeLa, SiHa, Ca-Ski, C-33A, and HT-3). Further experiments including wound scratch and transwell invasion assays, Western blotting, immunofluorescence, and luciferase assays were used to investigate the function of lncRNA-CTS/miR-505/ZEB2 in vitro. In addition, a tumor xenograft model was used to assess the effect of lncRNA-CTS in vivo. The expression levels of lncRNA-CTS and miR-505 were correlated with the metastasis-associated clinicopathological features of CC patients. Moreover, lncRNA-CTS was associated with a poor prognosis in CC patients. In vitro and in vivo experiments, along with gain- and loss-of-function studies, showed that lncRNA-CTS enhanced cell migration, invasion, and the transforming growth factor (TGF)-ß1-induced-EMT process. Data also showed that lncRNA-CTS could function as a competing endogenous RNA for miR-505 in CC cells. Further investigations disclosed that ZEB2 was demonstrated as a downstream target of miR-505, and subsequently exerted its metastatic effects via the lncRNA-CTS/miR-505/ZEB2 axis in CC cells. Finally, lncRNA-CTS activated the SMAD/TGF pathway via miR-505 in CC cells. Collectively, our results demonstrate the importance of the lncRNA-CTS/miR-505/ZEB2 axis in CC. LncRNA-CTS can predispose CC patients to metastases and may represent a promising therapeutic target for CC.