Retrotransposons evolution and impact on lncRNA and protein coding genes in pigs

Population genetic analysis based on the polymorphisms mediated by transposons in the genomes of pig

Transposable elements (TEs) mobility is capable of generating a large number of structural variants (SVs), which can have considerable potential as molecular markers for genetic analysis and molecular breeding in livestock. Our results showed that the pig genome contains mainly TE-SVs generated by short interspersed nuclear elements (51,873/76.49%), followed by long interspersed nuclear elements (11,131/16.41%), and more than 84% of the common TE-SVs (Minor allele frequency, MAF > 0.10) were validated to be polymorphic. Subsequently, we utilized the identified TE-SVs to gain insights into the population structure, resulting in clear differentiation among the three pig groups and facilitating the identification of relationships within Chinese local pig breeds. In addition, we investigated the frequencies of TEs in the gene coding regions of different pig groups and annotated the respective TE types, related genes, and functional pathways. Through genome-wide comparisons of Large White pigs and Chinese local pigs utilizing the Beijing Black pigs, we identified TE-mediated SVs associated with quantitative trait loci and observed that they were mainly involved in carcass traits and meat quality traits. Lastly, we present the first documented evidence of TE transduction in the pig genome.

Evolution of Endogenous Retroviruses in the Subfamily of Caprinae

The interest in endogenous retroviruses (ERVs) has been fueled by their impact on the evolution of the host genome. In this study, we used multiple pipelines to conduct a de novo exploration and annotation of ERVs in 13 species of the Caprinae subfamily. Through analyses of sequence identity, structural organization, and phylogeny, we defined 28 ERV groups within Caprinae, including 19 gamma retrovirus groups and 9 beta retrovirus groups. Notably, we identified four recent and potentially active groups prevalent in the Caprinae genomes. Additionally, our investigation revealed that most long noncoding genes (lncRNA) and protein-coding genes (PC) contain ERV-derived sequences. Specifically, we observed that ERV-derived sequences were present in approximately 75% of protein-coding genes and 81% of lncRNA genes in sheep. Similarly, in goats, ERV-derived sequences were found in approximately 74% of protein-coding genes and 75% of lncRNA genes. Our findings lead to the conclusion that the majority of ERVs in the Caprinae genomes can be categorized as fossils, representing remnants of past retroviral infections that have become permanently integrated into the genomes. Nevertheless, the identification of the Cap_ERV_20, Cap_ERV_21, Cap_ERV_24, and Cap_ERV_25 groups indicates the presence of relatively recent and potentially active ERVs in these genomes. These particular groups may contribute to the ongoing evolution of the Caprinae genome. The identification of putatively active ERVs in the Caprinae genomes raises the possibility of harnessing them for future genetic marker development.

Retroviral Insertion Polymorphism (RIP) of Porcine Endogenous Retroviruses (PERVs) in Pig Genomes

Endogenous retroviruses (ERVs) are one of the superfamilies of long terminal repeat retrotransposons (LTRs) in mice and humans. Approximately 8% of the pig genome is composed of sequences derived from LTRs. While the majority of ERVs in pigs have decayed, a small number of full-length copies can still mobilize within the genome. This study investigated the unexplored retroviral insertion polymorphisms (RIPs) generated by the mobilization of full-length ERVs (Fl-ERVs), and evaluated their impact on phenotypic variation to gain insights into the biological role of Fl-ERVs in pigs. Overall, 39 RIPs (insertions or deletions relative to the pig reference genome) generated by Fl-ERVs were predicted by comparative genomic analysis, and 18 of them were confirmed by PCR detection. Four RIP sites (D5, D14, D15, and D18) were further evaluated by population analysis, and all of them displayed polymorphisms in multiple breeds. The RIP site of ERV-D14, which is a Fl-ERV inserted in the STAB2-like gene, was further confirmed by sequencing. Population analysis of the polymorphic site of ERV-D14 reveals that it presents moderate polymorphism information in the Large White pig breed, and the association analysis reveals that the RIP of ERV-D14 is associated with age variations at 30 kg body weight (p < 0.05) and 100 kg body weight (p < 0.01) in the population of Large White pigs (N = 480). Furthermore, the ERV-D14 RIP is associated with changes in the expression of the target gene STAB2-like in the liver, backfat, and leaf fat in Sushan pigs. These data suggest that some Fl-ERVs are still mobilizing in the pig's genome, and contribute to genomic and phenotypic variations.

Taming transposable elements in livestock and poultry: a review of their roles and applications

Livestock and poultry play a significant role in human nutrition by converting agricultural by-products into high-quality proteins. To meet the growing demand for safe animal protein, genetic improvement of livestock must be done sustainably while minimizing negative environmental impacts. Transposable elements (TE) are important components of livestock and poultry genomes, contributing to their genetic diversity, chromatin states, gene regulatory networks, and complex traits of economic value. However, compared to other species, research on TE in livestock and poultry is still in its early stages. In this review, we analyze 72 studies published in the past 20 years, summarize the TE composition in livestock and poultry genomes, and focus on their potential roles in functional genomics. We also discuss bioinformatic tools and strategies for integrating multi-omics data with TE, and explore future directions, feasibility, and challenges of TE research in livestock and poultry. In addition, we suggest strategies to apply TE in basic biological research and animal breeding. Our goal is to provide a new perspective on the importance of TE in livestock and poultry genomes.

SINE Insertion in the Pig Carbonic Anhydrase 5B (CA5B) Gene Is Associated with Changes in Gene Expression and Phenotypic Variation

Transposons are genetic elements that are present in mammalian genomes and occupy a large proportion of the pig genome, with retrotransposons being the most abundant. In a previous study, it was found that a SINE retrotransposon was inserted in the 1^st intron of the CA5B gene in pigs, and the present study aimed to investigate the SINE insertion polymorphism in this gene in different pig breeds. Polymerase chain reaction (PCR) was used to confirm the polymorphism in 11 pig breeds and wild boars), and it was found that there was moderate polymorphism information content in 9 of the breeds. Further investigation in cell experiments revealed that the 330 bp SINE insertion in the RIP-CA5B site promoted expression activity in the weak promoter region of this site. Additionally, an enhancer verification vector experiment showed that the 330 bp SINE sequence acted as an enhancer on the core promoter region upstream of the CA5B gene region. The expression of CA5B in adipose tissue (back fat and leaf fat) in individuals with the (SINE^+/+) genotype was significantly higher than those with (SINE^+/-) and (SINE^-/-) genotypes. The association analysis revealed that the (SINE^+/+) genotype was significantly associated with a higher back fat thickness than the (SINE^-/-) genotype. Moreover, it was observed that the insertion of SINE at the RIP-CA5B site carried ATTT repeats, and three types of (ATTT) repeats were identified among different individuals/breeds (i.e., (ATTT)₄, (ATTT)₆ and (ATTT)₉). Overall, the study provides insights into the genetic basis of adipose tissue development in pigs and highlights the role of a SINE insertion in the CA5B gene in this process.

Whole-genome identification of transposable elements reveals the equine repetitive element insertion polymorphism in Chinese horses

Transposable elements (TEs) are diverse, abundant, and complicated in genomes. They not only can drive the genome evolution process but can also act as special resources for adaptation. However, little is known about the evolutionary processes that shaped horses. In this work, 126 horse assemblages involved in most horse breeds in China were used to investigate the patterns of TE variation for the first time. By using RepeatMasker and melt software, we found that the horse-specific short interspersed repetitive elements family, equine repetitive elements (ERE1), exhibited polymorphisms in horse genomes. Phylogenetic analysis based on these ERE1 loci (minor allele frequency ≥0.05) revealed three major horse groups, namely, those in northern China, southern China, and Qinghai-Tibetan, which mirrors the result determined by SNPs to some extent. The present ERE1 family emerged ~0.26 to 1.77 Mya ago, with an activity peak at ~0.49 Mya, which matches the early stage of the horse lineage and decreases after the divergence of Equus caballus and Equus ferus przewalskii. To detect the functional ERE1(s) associated with adaptation, locus-specific branch length, genome-wide association study, and absolute allele frequency difference analyses were conducted and resulted in two common protein-coding genes annotated by candidate ERE1s. They were clustered into the vascular smooth muscle contraction (p = 0.01, EDNRA) and apelin signalling pathways (p = 0.02, NRF1). Notably, ERE1 insertion into the EDNRA gene showed a higher association with adaptation among southern China horses and other horses in 15 populations and 451 individuals (p = 4.55 e-8). Our results provide a comprehensive understanding of TE variations to analyse the phylogenetic relationships and traits relevant to adaptive evolution in horses.

Pangenome obtained by long-read sequencing of 11 genomes reveal hidden functional structural variants in pigs

Long-read sequencing (LRS) facilitates both the genome assembly and the discovery of structural variants (SVs). Here, we built a graph-based pig pangenome by incorporating 11 LRS genomes with an average of 94.01% BUSCO completeness score, revealing 206-Mb novel sequences. We discovered 183,352 nonredundant SVs (63% novel), representing 12.12% of the reference genome. By genotyping SVs in an additional 196 short-read sequencing samples, we identified thousands of population stratified SVs. Particularly, we detected 7,568 Tibetan specific SVs, some of which demonstrate significant population differentiation between Tibetan and low-altitude pigs, which might be associated with the high-altitude hypoxia adaptation in Tibetan pigs. Further integrating functional genomic data, the most promising candidate genes within the SVs that might contribute to the high-altitude hypoxia adaptation were discovered. Overall, our study generates a benchmark pangenome resource for illustrating the important roles of SVs in adaptive evolution, domestication, and genetic improvement of agronomic traits in pigs.

Two Retrotransposon Elements in Intron of Porcine BMPR1B Is Associated with Phenotypic Variation

It has been established that through binding to bone morphogenetic proteins (BMPs), bone morphogenetic protein receptor I B (BMPR1B) can mediate transforming growth factor β (TGF-β) signal transduction, and is involved in the regulation of several biological processes, such as bone and muscle formation and homeostasis, as well as folliculogenesis. Also known as FecB, BMPR1B has been reported as the major gene for sheep prolificacy. A number of previous studies have analyzed the relationship between single nucleotide polymorphisms (SNPs) in this gene and its related performance. In recent years, with the illustration of the effect of retrotransposon insertion on the expression of the proximal genes or phenotypic variation, retrotransposon insertion polymorphisms (RIPs) have been used as a novel type of molecular marker in the evaluation of evolution, population structure and breeding of plant and domestic animals. In this study, the RIPs in porcine BMPR1B gene were excavated, and thereafter verified using a comparative genome and polymerase chain reaction (PCR). The potential effects of phenotype, gene expression and functions related to RIPs were also explored. The results showed that 13 distinct RIPs were identified in introns of porcine BMPR1B. Among these, only BMPR1B-SINE-RIP9 and BMPR1B-LINE-RIP13 displayed a close relationship with the growth traits of Large White pigs. Moreover, the total number of BMPR1B-SINE+/+-RIP9 individuals born was found to be significantly higher than that of SINE−/− (p < 0.05). These two RIPs showed an obvious distribution pattern among Chinese indigenous breeds and Western commercial breeds. The expression of BMPR1B in ovaries of adult BMPR1B-SINE+/+-RIP9 Sushan pigs was found to be significantly higher in comparison to those of BMPR1B-SINE−/−-RIP9 (p < 0.05). SINE insertion of BMPR1B-SINE-RIP9 and LINE insertion of BMPR1B-LINE-RIP13 were observed to significantly increase the activity of Octamer binding transcription factor 4 (OCT4) minipromoter in CHO and C2C12 cells (p < 0.01). Therefore, these two RIPs could serve as useful molecular markers for modulating the growth or reproductive traits in assisted selection of pig breeding, while the mechanisms of the insertion function should be studied further.

SINE Insertion May Act as a Repressor to Affect the Expression of Pig LEPROT and Growth Traits

Retrotransposon is an important component of the mammalian genome. Previous studies have shown that the expression of protein-coding genes was affected by the insertion of retrotransposon into the proximal genes, and the phenotype variations would be related to the retrotransposon insertion polymorphisms (RIPs). In this study, leptin (LEP), leptin receptor (LEPR), and leptin receptor overlapping transcript (LEPROT), which play important roles in the regulation of fat synthesis and body weight, were screened to search for the RIPs and their effect on phenotype and gene expression, as well as to further study the function of the insertion. The results showed that three RIPs located in intron 1 of LEPROT and intron 2 and 21 of LEPR were identified, and they were all SINEA1, which was one type of retrotransposon. The SINE insertion at the LEPROT was the dominant allele in native pig breeds. The age of 100 kg body weight of SINE+/+ Large White individuals was significantly higher than those of SINE+/− and SINE−/− individuals (p < 0.05). The LEPROT gene expression in the liver and suet of 30-day-old SINE−/− Sujiang piglets were significantly higher than those of SINE+/+ and SINE+/− piglets (p < 0.01). The dual-luciferase reporter gene assay showed that SINE insertion in PK15 and 3T3-L1 cells significantly reduced the promoter activity of the LEPROT gene (p < 0.01). Therefore, SINE insertion can be a repressor to reduce the expression of LEPROT and could be a useful molecular marker for assisted selection of growth traits in pig breeding.

Retrotransposon Insertion Polymorphisms (RIPs) in Pig Reproductive Candidate Genes

Retrotransposons account for more than one-third of the pig reference genome. On account of the genome variability in different breeds, structural variation (SV) caused by retrotranspos-on-generated deletion or insertion (indel) may have a function in the genome. Litter size is one of the most important reproductive traits and significantly impacts profitability in terms of pig production. We used the method of bioinformatics, genetics, and molecular biology to make an analysis among different pig genomes. Predicted 100 SVs were annotated as retrotransposon indel in 20 genes related to reproductive performance. The PCR detection based on these predicted SVs revealed 20 RIPs in 20 genes, that most RIPs (12) were generated by SINE indel, and eight RIPs were generated by the ERV indel. We selected 12 RIPs to make the second round PCR detection in 24 individuals among nine pig breeds. The PCR detection results revealed that the RIP-A1CF-4 insertion in the breed of Bama, Large White, and Meishan only had the homozygous genotype but low to moderately polymorphisms were present in other breeds. We found that RIP-CWH43-9, RIP-IDO2-9, RIP-PRLR-6, RIP-VMP1-12, and RIP-OPN-1 had a rich polymorphism in the breed of Large White pigs. The statistical analysis revealed that RIP-CWH43-9 had a SINE insertion profitable to the reproductive traits of TNB and NBA but was significantly affected (p < 0.01) and (p < 0.05) in the reproductive traits of litter birthweight (LW) in Large White. On the other hand, the SINE insertion in IDO2-9 may be a disadvantage to the reproductive traits of LW, which was significantly affected (p < 0.05) in Large White. These two RIPs are significant in pig genome research and could be useful molecular markers in the breeding system.

Genetic Evaluation and Population Structure of Jiangsu Native Pigs in China Revealed by SINE Insertion Polymorphisms

Simple Summary

In a previous study, 30 SINE-RIPs were applied for population genetic analysis in 7 Chinese miniature pig populations and approved effectively in the genetic distances and breed-relationships between these populations. There are abundant indigenous pigs famous across the world for their prolificacy in the Jiangsu Province of eastern China, such as Meishan, Erhualian. Since pork production relies on limited commercial breeds such as Landrace, Large White, and Duroc pigs, characterized by maximizing productivity in intensive production systems, these indigenous pigs are nowadays decrease sharply. The genetic characterizations of these resources are essential requirements for the development of conservation, selection, and sustainable utilizations. Therefore, SINE-RIPs were selected to evaluate the genetic variation and population structure of Jiangsu pig populations and the results may assist with the conservation and utilization of these native pig populations.

Abstract

Short interspersed nuclear elements (SINEs), one type of retrotransposon, are considered to be ideal molecular markers due to their wide distribution in the genome, high copy number, and high polymorphism. Preliminary studies have identified more than 35,000 SINE-retrotransposon insertion polymorphisms (RIPs) in the pig genome. In this study, 18 SINE-RIPs were used to evaluate the genetic variation and population structure of seven native pig populations and two crossbreeds in the Jiangsu Province of China. Two commercial pig breeds (Duroc and Large White) and one Italian native breed (Sicilian Black pig) were selected as the control. The results showed that all 18 SINE-RIPs were polymorphic among these pigs. The Jiangsu native pig populations (Erhualian, Fengjing, Middle Meishan, Mi, Shawutou, Small Meishan, and Huai) were shown to be more polymorphic than the crossbreeds (Sushan and Sujiang) and external breeds (Sicilian Black pig, Large White, and Duroc) based on the expected heterozygosity and polymorphic information content values. Some native pigs, including Small Meishan, Mi, Middle Meishan, and Erhualian, had a higher degree of inbreeding according to the F_IS values. Based on the neighbor-joining tree, all of the Jiangsu native pig populations formed one branch, while the three external pig breeds formed the other branches, with the two crossbreeds containing more than 50% external pig ancestry. The Huai pigs were independent of the other Jiangsu native pigs but shared a common ancestor with Sujiang and Mi. The results provide a new perspective on the population structure of these native pig breeds and will assist with the conservation and utilization of Chinese native pigs.

Retrotransposon Insertion Polymorphisms (RIPs) in Pig Coat Color Candidate Genes

The diversity of livestock coat color results from human positive selection and represents an indispensable part of breed identity. As an important biodiversity resource, pigs have many special characteristics, including the most visualized feature, coat color, and excellent adaptation, and the coat color represents an important phenotypic characteristic of the pig breed. Exploring the genetic mechanisms of phenotypic characteristics and the melanocortin system is of considerable interest in domestic animals because their energy metabolism and pigmentation have been under strong selection. In this study, 20 genes related to coat color in mammals were selected, and the structural variations (SVs) in these genic regions were identified by sequence alignment across 17 assembled pig genomes, from representing different types of pigs (miniature, lean, and fat type). A total of 167 large structural variations (>50 bp) of coat-color genes, which overlap with retrotransposon insertions (>50 bp), were obtained and designated as putative RIPs. Finally, 42 RIPs were confirmed by PCR detection. Additionally, eleven RIP sites were further evaluated for their genotypic distributions by PCR in more individuals of eleven domesticated breeds representing different coat color groups. Differential distributions of these RIPs were observed across populations, and some RIPs may be associated with breed differences.

A 192 bp ERV fragment insertion in the first intron of porcine TLR6 may act as an enhancer associated with the increased expressions of TLR6 and TLR1

Background

Toll-like receptors (TLRs) play important roles in building innate immune and inducing adaptive immune responses. Associations of the TLR genes polymorphisms with disease susceptibility, which are the basis of molecular breeding for disease resistant animals, have been reported extensively. Retrotransposon insertion polymorphisms (RIPs), as a new type of molecular markers developed recently, have great potential in population genetics and quantitative trait locus mapping. In this study, bioinformatic prediction combined with PCR-based amplification was employed to screen for RIPs in porcine TLR genes. Their population distribution was examined, and for one RIP the impact on gene activity and phenotype was further evaluated.

Results

Five RIPs, located at the 3' flank of TLR3, 5' flank of TLR5, intron 1 of TLR6, intron 1 of TLR7, and 3' flank of TLR8 respectively, were identified. These RIPs were detected in different breeds with an uneven distribution among them. By using the dual luciferase activity assay a 192 bp endogenous retrovirus (ERV) in the intron 1 of TLR6 was shown to act as an enhancer increasing the activities of TLR6 putative promoter and two mini-promoters. Furthermore, real-time quantitative polymerase chain reaction (qPCR) analysis revealed significant association (p < 0.05) of the ERV insertion with increased mRNA expression of TLR6, the neighboring gene TLR1, and genes downstream in the TLR signaling pathway such as MyD88 (Myeloid differentiation factor 88), Rac1 (Rac family small GTPase 1), TIRAP (TIR domain containing adaptor protein), Tollip (Toll interacting protein) as well as the inflammatory factors IL6 (Interleukin 6), IL8 (Interleukin 8), and TNFα (Tumor necrosis factor alpha) in tissues of 30 day-old piglet. In addition, serum IL6 and TNFα concentrations were also significantly upregulated by the ERV insertion (p < 0.05).

Conclusions

A total of five RIPs were identified in five different TLR loci. The 192 bp ERV insertion in the first intron of TLR6 was associated with higher expression of TLR6, TLR1, and several genes downstream in the signaling cascade. Thus, the ERV insertion may act as an enhancer affecting regulation of the TLR signaling pathways, and can be potentially applied in breeding of disease resistant animals.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13100-021-00248-w.

Distinct Retrotransposon Evolution Profile in the Genome of Rabbit (Oryctolagus cuniculus)

Although the rabbit genome has already been annotated, it is mobilome remains largely unknown. Here, multiple pipelines were used to de novo mine and annotate the mobilome in rabbit. Four families and 19 subfamilies of LINE1s, two families and nine subfamilies of SINEs, and 12 ERV families were defined in rabbit based on sequence identity, structural organization, and phylogenetic tree. The analysis of insertion age and polymerase chain reaction suggests that a number of families are very young and may remain active, such as L1B, L1D, OcuSINEA, and OcuERV1. RepeatMasker annotation revealed a distinct transposable element landscape within the genome, with approximately two million copies of SINEs, representing the greatest proportion of the genome (19.61%), followed by LINEs (15.44%), and LTRs (4.11%), respectively, considerably different from most other mammal mobilomes except hedgehog and tree shrew, in which LINEs have the highest proportion. Furthermore, a very high rate of insertion polymorphisms (>85%) for the youngest subfamily (OcuSINEA1) was identified by polymerase chain reaction. The majority of retrotransposon insertions overlapped with protein-coding regions (>80%) and lncRNA (90%) genes. Genomic distribution bias was observed for retrotransposons, with those immediately upstream (-1 kb) and downstream (1 kb) of genes significantly depleted. Local GC content in 50-kb widows had significantly negative correlations with LINE (rs=-0.996) and LTR (rs=-0.829) insertions. The current study revealed a distinct mobilome landscape in rabbit, which will assist in the elucidation of the evolution of the genome of lagomorphs, and even other mammals.

SINE jumping contributes to large-scale polymorphisms in the pig genomes

BACKGROUND

Molecular markers based on retrotransposon insertion polymorphisms (RIPs) have been developed and are widely used in plants and animals. Short interspersed nuclear elements (SINEs) exert wide impacts on gene activity and even on phenotypes. However, SINE RIP profiles in livestock remain largely unknown, and not be revealed in pigs.

RESULTS

Our data revealed that SINEA1 displayed the most polymorphic insertions (22.5 % intragenic and 26.5 % intergenic), followed by SINEA2 (10.5 % intragenic and 9 % intergenic) and SINEA3 (12.5 % intragenic and 5.0 % intergenic). We developed a genome-wide SINE RIP mining protocol and obtained a large number of SINE RIPs (36,284), with over 80 % accuracy and an even distribution in chromosomes (14.5/Mb), and 74.34 % of SINE RIPs generated by SINEA1 element. Over 65 % of pig SINE RIPs overlap with genes, most of them (> 95 %) are in introns. Overall, about one forth (23.09 %) of the total genes contain SINE RIPs. Significant biases of SINE RIPs in the transcripts of protein coding genes were observed. Nearly half of the RIPs are common in these pig breeds. Sixteen SINE RIPs were applied for population genetic analysis in 23 pig breeds, the phylogeny tree and cluster analysis were generally consistent with the geographical distributions of native pig breeds in China.

CONCLUSIONS

Our analysis revealed that SINEA1-3 elements, particularly SINEA1, are high polymorphic across different pig breeds, and generate large-scale structural variations in the pig genomes. And over 35,000 SINE RIP markers were obtained. These data indicate that young SINE elements play important roles in creating new genetic variations and shaping the evolution of pig genome, and also provide strong evidences to support the great potential of SINE RIPs as genetic markers, which can be used for population genetic analysis and quantitative trait locus (QTL) mapping in pig.

SINE Insertion in the Intron of Pig GHR May Decrease Its Expression by Acting as a Repressor

Simple Summary

GH/IGF axis genes play a central role in the regulation of skeletal accretion during development and growth, and thus represent candidate genes for growth traits. Retrotransposon insertion polymorphisms are major contributors to structural variations. They tend to generate large effect mutations resulting in variations in target gene activity and phenotype due to the fact that they carry functional elements, such as enhancers, insulators, or promoters. In the present study, RIPs in four GH/IGF axis genes (GH, GHR, IGF1, and IGF1R) were investigated by comparative genomics and PCR. Four RIPs in the GHR gene and one RIP in the IGF1 gene were identified. Further analysis revealed that one RIP in the first intron of GHR might play a role in the regulation of GHR expression by acting as a repressor. These findings contribute to the understanding of the role of RIPs in the genetic variation of GH/IGF axis genes and phenotypic variation in pigs.

Abstract

The genetic diversity of the GH/IGF axis genes and their association with the variation of gene expression and phenotypic traits, principally represented by SNPs, have been extensively reported. Nevertheless, the impact of retrotransposon insertion polymorphisms (RIPs) on the GH/IGF axis gene activity has not been reported. In the present study, bioinformatic prediction and PCR verification were performed to screen RIPs in four GH/IGF axis genes (GH, GHR, IGF1 and IGF1R). In total, five RIPs, including one SINE RIP in intron 3 of IGF1, one L1 RIP in intron 7 of GHR, and three SINE RIPs in intron 1, intron 5 and intron 9 of GHR, were confirmed by PCR, displaying polymorphisms in diverse breeds. Dual luciferase reporter assay revealed that the SINE insertion in intron 1 of GHR significantly repressed the GHR promoter activity in PK15, Hela, C2C12 and 3T3-L1 cells. Furthermore, qPCR results confirmed that this SINE insertion was associated with a decreased expression of GHR in the leg muscle and longissimus dorsi, indicating that it may act as a repressor involved in the regulation of GHR expression. In summary, our data revealed that RIPs contribute to the genetic variation of GH/IGF axis genes, whereby one SINE RIP in the intron 1 of GHR may decrease the expression of GHR by acting as a repressor.

Genetic Diversity and Population Structures in Chinese Miniature Pigs Revealed by SINE Retrotransposon Insertion Polymorphisms, a New Type of Genetic Markers

Simple Summary

Our previous studies suggested that the short interspersed nuclear element (SINE) retrotransposon insertion polymorphisms (RIPs), as a new type of molecular marker developed very recently, are ideal molecular markers and have the potential to be used for population genetic analysis and molecular breeding in pigs and possibly it can be extended to other livestock animals as well. However, no report is available for the application of SINE RIPs in population genetic analysis in livestock, including pigs. Here, we evaluated 30 SINE RIPs in several indigenous Chinese miniature pig breeds, including three subpopulations of Bama pigs (BM-cov, BM-clo, and BM-inb). BM-cov is a subpopulation conserved in the national conservation farm, and BM-clo is a closed population maintained over 30 years with only 2 boars and 14 sows imported from its original area, while BM-inb herd is an 18 generation continuous inbreeding line based on the BM-clo population. To our knowledge, it is the first time to report the genetic diversity, breed differentiation, and population structures for these populations by using SINE RIPs, and which suggests the feasibility of SINE RIPs in pig genetic analysis.

Abstract

RIPs have been developed as effective genetic markers and popularly applied for genetic analysis in plants, but few reports are available for domestic animals. Here, we established 30 new molecular markers based on the SINE RIPs, and applied them for population genetic analysis in seven Chinese miniature pigs. The data revealed that the closed herd (BM-clo), inbreeding herd (BM-inb) of Bama miniature pigs were distinctly different from the BM-cov herds in the conservation farm, and other miniature pigs (Wuzhishan, Congjiang Xiang, Tibetan, and Mingguang small ear). These later five miniature pig breeds can further be classified into two clades based on a phylogenetic tree: one included BM-cov and Wuzhishan, the other included Congjiang Xiang, Tibetan, and Mingguang small ear, which was well-supported by structure analysis. The polymorphic information contents estimated by using SINE RIPs are lower than the predictions based on microsatellites. Overall, the genetic distances and breed-relationships between these populations revealed by 30 SINE RIPs generally agree with their evolutions and geographic distributions. We demonstrated the potential of SINE RIPs as new genetic markers for genetic monitoring and population structure analysis in pigs, which can even be extended to other livestock animals.

Genome-wide DNA methylation dynamics during epigenetic reprogramming in the porcine germline

Background

Prior work in mice has shown that some retrotransposed elements remain substantially methylated during DNA methylation reprogramming of germ cells. In the pig, however, information about this process is scarce. The present study was designed to examine the methylation profiles of porcine germ cells during the time course of epigenetic reprogramming.

Results

Sows were artificially inseminated, and their fetuses were collected 28, 32, 36, 39, and 42 days later. At each time point, genital ridges were dissected from the mesonephros and germ cells were isolated through magnetic-activated cell sorting using an anti-SSEA-1 antibody, and recovered germ cells were subjected to whole-genome bisulphite sequencing. Methylation levels were quantified using SeqMonk software by performing an unbiased analysis, and persistently methylated regions (PMRs) in each sex were determined to extract those regions showing 50% or more methylation. Most genomic elements underwent a dramatic loss of methylation from day 28 to day 36, when the lowest levels were shown. By day 42, there was evidence for the initiation of genomic re-methylation. We identified a total of 1456 and 1122 PMRs in male and female germ cells, respectively, and large numbers of transposable elements (SINEs, LINEs, and LTRs) were found to be located within these PMRs. Twenty-one percent of the introns located in these PMRs were found to be the first introns of a gene, suggesting their regulatory role in the expression of these genes. Interestingly, most of the identified PMRs were demethylated at the blastocyst stage.

Conclusions

Our findings indicate that methylation reprogramming in pig germ cells follows the general dynamics shown in mice and human, unveiling genomic elements that behave differently between male and female germ cells.

Long Interspersed Nuclear Elements 1 (LINE1): The chimeric transcript L1-MET and its involvement in cancer

Long interspersed nuclear elements 1 (LINE1) are non-LTR retrotransposons that represent the greatest remodeling force of the human genome during evolution. Genomically, LINE1 are constituted by a 5´ untranslated region (UTR), where the promoter regions are located, three open reading frames (ORF0, ORF1, and ORF2) and one 3´UTR, which has a poly(A) tail that harbors the short interspersed nuclear elements (SINEs) Alu and SVA. Although the intrinsic nature of LINE1 is to be copied and inserted into the genome, an increase in their mobility produces genomic instability. In response to this, the cell has "designed" many mechanisms controlling the retrotransposition levels of LINE1; however, alterations in these regulation systems can increase LINE1 mobility and the formation of chimeric genes. Evidence indicates that 988 human genes have LINE1 inserted in their sequence, resulting in the transcriptional control of genes by their own promoters, as well as by the LINE1 antisense promoter (ASP). To date, very little is known about the biologic impact of this and the L1-MET chimera is a more or less studied case. ASP hypomethylation has been observed in all studied cancer types, leading to increased L1-MET expression. In specific types of cancer, this L1-MET increase controls both low and high MET protein levels. It remains to be clarified if this protein product is a chimeric protein.

Development and Deployment of High-Throughput Retrotransposon-Based Markers Reveal Genetic Diversity and Population Structure of Asian Bamboo

Bamboo, a non-timber grass species, known for exceptionally fast growth is a commercially viable crop. Long terminal repeat (LTR) retrotransposons, the main class I mobile genetic elements in plant genomes, are highly abundant (46%) in bamboo, contributing to genome diversity. They play significant roles in the regulation of gene expression, chromosome size and structure as well as in genome integrity. Due to their random insertion behavior, interspaces of retrotransposons can vary significantly among bamboo genotypes. Capitalizing this feature, inter-retrotransposon amplified polymorphism (IRAP) is a high-throughput marker system to study the genetic diversity of plant species. To date, there are no transposon based markers reported from the bamboo genome and particularly using IRAP markers on genetic diversity. Phyllostachys genus of Asian bamboo is the largest of the Bambusoideae subfamily, with great economic importance. We report structure-based analysis of bamboo genome for the LTR-retrotransposon superfamilies, Ty3-gypsy and Ty1-copia, which revealed a total of 98,850 retrotransposons with intact LTR sequences at both the ends. Grouped into 64,281 clusters/scaffold using CD-HIT-EST software, only 13 clusters of retroelements were found with more than 30 LTR sequences and with at least one copy having all intact protein domains such as gag and polyprotein. A total of 16 IRAP primers were synthesized, based on the high copy numbers of conserved LTR sequences. A study using these IRAP markers on genetic diversity and population structure of 58 Asian bamboo accessions belonging to the genus Phyllostachys revealed 3340 amplicons with an average of 98% polymorphism. The bamboo accessions were collected from nine different provinces of China, as well as from Italy and America. A three phased approach using hierarchical clustering, principal components and a model based population structure divided the bamboo accessions into four sub-populations, PhSP1, PhSP2, PhSP3 and PhSP4. All the three analyses produced significant sub-population wise consensus. Further, all the sub-populations revealed admixture of alleles. The analysis of molecular variance (AMOVA) among the sub-populations revealed high intra-population genetic variation (75%) than inter-population. The results suggest that Phyllostachys bamboos are not well evolutionarily diversified, although geographic speciation could have occurred at a limited level. This study highlights the usability of IRAP markers in determining the inter-species variability of Asian bamboos.