Whole-genome resequencing reveals new mutations in candidate genes for Beichuan-white goat prolificacya

In this study, we evaluated the copy number variation in the genomes of two groups of Beichuan-white goat populations with large differences in litter size by FST method, and identified 1739 genes and 485 missense mutations in the genes subject to positive selection. Through functional enrichment, ITGAV, LRP4, CDH23, TPRN, RYR2 and CELSR1 genes, involved in embryonic morphogenesis, were essential for litter size trait, which received intensive attention. In addition, some mutation sites of these genes have been proposed (ITGAV: c.38C > T; TPRN: c.133A > T, c.1192A > G, c.1250A > C; CELSR1: c.7640T > C), whose allele frequencies were significantly changed in the high fecundity goat group. Besides, we found that new mutations at these sites altered the hydrophilicity and 3D structure of the protein. Candidate genes related to litter size in this study and their missense mutation sites were identified. These candidate genes are helpful to understand the genetic mechanism of fecundity in Beichuan white goat, and have important significance for future goat breeding.

Deciphering the population structure and genetic basis of growth traits from whole-genome resequencing of the leopard coral grouper ( Plectropomus leopardus)

The leopard coral grouper ( Plectropomus leopardus) is a species of significant economic importance. Although artificial cultivation of P. leopardus has thrived in recent decades, the advancement of selective breeding has been hindered by the lack of comprehensive population genomic data. In this study, we identified over 8.73 million single nucleotide polymorphisms (SNPs) through whole-genome resequencing of 326 individuals spanning six distinct groups. Furthermore, we categorized 226 individuals with high-coverage sequencing depth (≥14×) into eight clusters based on their genetic profiles and phylogenetic relationships. Notably, four of these clusters exhibited pronounced genetic differentiation compared with the other populations. To identify potentially advantageous loci for P. leopardus, we examined genomic regions exhibiting selective sweeps by analyzing the nucleotide diversity ( θπ) and fixation index ( F ST) in these four clusters. Using these high-coverage resequencing data, we successfully constructed the first haplotype reference panel specific to P. leopardus. This achievement holds promise for enabling high-quality, cost-effective imputation methods. Additionally, we combined low-coverage sequencing data with imputation techniques for a genome-wide association study, aiming to identify candidate SNP loci and genes associated with growth traits. A significant concentration of these genes was observed on chromosome 17, which is primarily involved in skeletal muscle and embryonic development and cell proliferation. Notably, our detailed investigation of growth-related SNPs across the eight clusters revealed that cluster 5 harbored the most promising candidate SNPs, showing potential for genetic selective breeding efforts. These findings provide a robust toolkit and valuable insights into the management of germplasm resources and genome-driven breeding initiatives targeting P. leopardus.

Integration of transcriptome and whole-genome re-sequencing analyses reveal growth-related candidate genes in Procambarus clarkii

Growth is a crucial economic trait of all aquaculture species. It is important to explore the molecular regulation on growth, which could help improve the growth rate of species. Mining the growth-related genes is the foundation for revealing its molecular regulation on growth. Presently, the molecular regulation of growth in Procambarus clarkii is not clear, and the study on exploring growth-related genes is limited. In this study, RNA-Seq was used to compare gene expression profiles of the individuals with different growth rates involved in four groups including Big Male (BM), Big Female (BF), Small male (SM), and Small Female (SF) from one P. clarkii family, and the analyses were performed in combination with sex. Meanwhile, whole-genome resequencing data was used to get growth-specific SNP (Single Nucleotide Polymorphism)/InDel (Insertion/Deletion) sites information. Totally, we identified 16,127 genes, of which 9065 were successfully annotated in the GO database. Among these, 1328 DEGs were identified in BM vs. SM, with 357 up-regulated and 971 down-regulated. Additionally, 3507 DEGs were identified in BF vs. SF, with 241 up-regulated and 3266 down-regulated. 96 DEGs were up-regulated and 820 DEGs were down-regulated in Growth-related Group. The expression levels of nine DEGs were validated by RT-qPCR to verify the analysis results of sequencing. 684,040 growth-related SNPs and 182,050 growth-related InDels were obtained after screened. These findings provide candidate growth-related genes and growth-specific SNP/InDel sites for regulation of growth traits in P. clarkii, and new insight into the molecular regulation of P. clarkii growth.

Whole-genome resequencing confirms the genetic effects of dams on an endangered fish Hemibagrus guttatus (Siluriformes: Bagridae): A case study in a tributary of the Pearl River

Dam construction in riverine ecosystems has fragmented natural aquatic habitats and has altered environmental conditions. As a result, damming has been demonstrated to threaten aquatic biodiversity by reducing species distribution ranges and hindering gene exchange, leading to the inability to adapt to environmental changes. Knowledge of the contemporary genetic diversity and genetic structure of fish populations that are separated by dams is vital to developing effective conservation strategies, particularly for endangered fish species. We chose the Lianjiang River, a tributary of the Pearl River, as a case study to assess the effects of dams on the genetic diversity and genetic structure of an endangered fish species, Hemibagrus guttatus, using whole-genome resequencing data from 63 fish samples. The results indicated low levels of genetic diversity, high levels of inbreeding and decreasing trend of effective population size in fragmented H. guttatus populations. In addition, there were significant genetic structure and genetic differentiation among populations, suggesting that the dams might have affected H. guttatus populations. Our findings may benefit management and conservation practices for this endangered species that is currently suffering from the effects of dam construction.

Identification of the genetic characteristics of copy number variations in experimental specific pathogen-free ducks using whole-genome resequencing

BACKGROUND

Specific pathogen-free ducks are a valuable laboratory resource for waterfowl disease research and poultry vaccine development. High throughput sequencing allows the systematic identification of structural variants in genomes. Copy number variation (CNV) can explain the variation of important duck genetic traits. Herein, the genome-wide CNVs of the three experimental duck species in China (Jinding ducks (JD), Shaoxing ducks (SX), and Fujian Shanma ducks (SM)) were characterized using resequencing to determine their genetic characteristics and selection signatures.

RESULTS

We obtained 4,810 CNV regions (CNVRs) by merging 73,012 CNVs, covering 4.2% of the duck genome. Functional analysis revealed that the shared CNVR-harbored genes were significantly enriched for 31 gene ontology terms and 16 Kyoto Encyclopedia of Genes and Genomes pathways (e.g., olfactory transduction and immune system). Based on the genome-wide fixation index for each CNVR, growth (SPAG17 and PTH1R), disease resistance (CATHL3 and DMBT1), and thermoregulation (TRPC4 and SLIT3) candidate genes were identified in strongly selected signatures specific to JD, SM, and SX, respectively.

CONCLUSIONS

In conclusion, we investigated the genome-wide distribution of experimental duck CNVs, providing a reference to establish the genetic basis of different phenotypic traits, thus contributing to the management of experimental animal genetic resources.

Karyotyping of aneuploid and polyploid plants from low coverage whole-genome resequencing

BACKGROUND

Karyotype, as a basic characteristic of species, provides valuable information for fundamental theoretical research and germplasm resource innovation. However, traditional karyotyping techniques, including fluorescence in situ hybridization (FISH), are challenging and low in efficiency, especially when karyotyping aneuploid and polyploid plants. The use of low coverage whole-genome resequencing (lcWGR) data for karyotyping was explored, but existing methods are complicated and require control samples.

RESULTS

In this study, a new protocol for molecular karyotype analysis was provided, which proved to be a simpler, faster, and more accurate method, requiring no control. Notably, our method not only provided the copy number of each chromosome of an individual but also an accurate evaluation of the genomic contribution from its parents. Moreover, we verified the method through FISH and published resequencing data.

CONCLUSIONS

This method is of great significance for species evolution analysis, chromosome engineering, crop improvement, and breeding.

Whole-genome resequencing to explore genome‑wide single nucleotide polymorphisms and genes associated with avian leukosis virus subgroup J infection in chicken

Avian leukosis virus subgroup J (ALV-J) is an oncogenic virus that causes serious economic loss in the poultry industry. Currently, no effective vaccine or drug is available against this virus. Therefore, it is imperative to explore and understand the molecular regulatory mechanisms underlying ALV-J infection. In this study, blood samples from 21 ALV-J-infected and 22 ALV-J-uninfected (DZ) chickens (JZ) were analyzed by whole-genome resequencing (WGR). By combining the fixation index (FST) with the nucleotide diversity (π) ratio based on WGR data, 425 candidate genes were identified. Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis revealed the top 20 enriched pathways, among which 9 pathways were significantly associated with diseases, including endometrial cancer, Chagas disease, PD-L1 expression and PD-1 checkpoint pathway in cancer, colorectal cancer, endocrine resistance, fluid shear stress, atherosclerosis, basal cell carcinoma, non-small cell lung cancer, and melanoma. Fourteen single nucleotide polymorphisms related to twelve genes showed a notable difference between DZ and JZ group chickens. The genes included COMMD3, PPP1CB, VEGFA, GTF2H1, NOTCH2, ITPR1, FGFR4, GNAS, NECTIN1, WNT2B, PPP1CC, and MRC2. These findings may provide a valuable foundation for further exploration of the pathogenesis of ALV-J in chickens.

Genome and whole-genome resequencing of Cinnamomum camphora elucidate its dominance in subtropical urban landscapes

BACKGROUND

Lauraceae is well known for its significant phylogenetic position as well as important economic and ornamental value; however, most evergreen species in Lauraceae are restricted to tropical regions. In contrast, camphor tree (Cinnamomum camphora) is the most dominant evergreen broadleaved tree in subtropical urban landscapes.

RESULTS

Here, we present a high-quality reference genome of C. camphora and conduct comparative genomics between C. camphora and C. kanehirae. Our findings demonstrated the significance of key genes in circadian rhythms and phenylpropanoid metabolism in enhancing cold response, and terpene synthases (TPSs) improved defence response with tandem duplication and gene cluster formation in C. camphora. Additionally, the first comprehensive catalogue of C. camphora based on whole-genome resequencing of 75 accessions was constructed, which confirmed the crucial roles of the above pathways and revealed candidate genes under selection in more popular C. camphora, and indicated that enhancing environmental adaptation is the primary force driving C. camphora breeding and dominance.

CONCLUSIONS

These results decipher the dominance of C. camphora in subtropical urban landscapes and provide abundant genomic resources for enlarging the application scopes of evergreen broadleaved trees.

Integrating genomics and transcriptomics to identify candidate genes for high egg production in Wulong geese (Anser cygnoides orientalis)

BACKGROUND

Wulong geese (Anser cygnoides orientalis) are known for their excellent egg-laying performance. However, they show considerable population differences in egg-laying behavior. This study combined genome-wide selection signal analysis with transcriptome analysis (RNA-seq) to identify the genes related to high egg production in Wulong geese.

RESULTS

A total of 132 selected genomic regions were screened using genome-wide selection signal analysis, and 130 genes related to high egg production were annotated in these regions. These selected genes were enriched in pathways related to egg production, including oocyte meiosis, the estrogen signaling pathway, the oxytocin signaling pathway, and progesterone-mediated oocyte maturation. Furthermore, a total of 890 differentially expressed genes (DEGs), including 340 up-regulated and 550 down-regulated genes, were identified by RNA-seq. Two genes - GCG and FAP - were common to the list of selected genes and DEGs. A non-synonymous single nucleotide polymorphism was identified in an exon of FAP.

CONCLUSIONS

Based on genome-wide selection signal analysis and transcriptome data, GCG and FAP were identified as candidate genes associated with high egg production in Wulong geese. These findings could promote the breeding of Wulong geese with high egg production abilities and provide a theoretical basis for exploring the mechanisms of reproductive regulation in poultry.

Whole-genome resequencing of the native sheep provides insights into the microevolution and identifies genes associated with reproduction traits

BACKGROUND

Sheep genomes undergo numerous genes losses, gains and mutation that generates genome variability among breeds of the same species after long time natural and artificial selection. However, the microevolution of native sheep in northwest China remains elusive. Our aim was to compare the genomes and relevant reproductive traits of four sheep breeds from different climatic environments, to unveil the selection challenges that this species cope with, and the microevolutionary differences in sheep genomes. Here, we resequenced the genomes of 4 representative sheep breeds in northwest China, including Kazakh sheep and Duolang sheep of native breeds, and Hu sheep and Suffolk sheep of exotic breeds with different reproductive characteristics.

RESULTS

We found that these four breeds had a similar expansion experience from ~ 10,000 to 1,000,000 years ago. In the past 10,000 years, the selection intensity of the four breeds was inconsistent, resulting in differences in reproductive traits. We explored the sheep variome and selection signatures by F_ST and θπ. The genomic regions containing genes associated with different reproductive traits that may be potential targets for breeding and selection were detected. Furthermore, non-synonymous mutations in a set of plausible candidate genes and significant differences in their allele frequency distributions across breeds with different reproductive characteristics were found. We identified PAK1, CYP19A1 and PER1 as a likely causal gene for seasonal reproduction in native sheep through qPCR, Western blot and ELISA analyses. Also, the haplotype frequencies of 3 tested gene regions related to reproduction were significantly different among four sheep breeds.

CONCLUSIONS

Our results provide insights into the microevolution of native sheep and valuable genomic information for identifying genes associated with important reproductive traits in sheep.

Whole-genome resequencing reveals recent signatures of selection in five populations of largemouth bass ( Micropterus salmoides)

Largemouth bass ( Micropterus salmoides) is an economically important fish species in North America, Europe, and China. Various genetic improvement programs and domestication processes have modified its genome sequence through selective pressure, leaving nucleotide signals that can be detected at the genomic level. In this study, we sequenced 149 largemouth bass fish, including protospecies (imported from the US) and improved breeds (four domestic breeding populations from China). We detected genomic regions harboring certain genes associated with improved traits, which may be useful molecular markers for practical domestication, breeding, and selection. Subsequent analyses of genetic diversity and population structure revealed that the improved breeds have undergone more rigorous genetic changes. Through selective signal analysis, we identified hundreds of putative selective sweep regions in each largemouth bass line. Interestingly, we predicted 103 putative candidate genes potentially subjected to selection, including several associated with growth (p sst1 and grb10), early development ( klf9, sp4, and sp8), and immune traits ( pkn2, sept2, bcl6, and ripk2). These candidate genes represent potential genomic landmarks that could be used to improve important traits of biological and commercial interest. In summary, this study provides a genome-wide map of genetic variations and selection footprints in largemouth bass, which may benefit genetic studies and accelerate genetic improvement of this economically important fish.

Construction of the first high-density SNP genetic map and identification of QTLs for the natural rubber content in Taraxacum kok-saghyz Rodin

BACKGROUND

Taraxacum kok-saghyz Rodin (TKS) is a promising commercial alternative natural rubber (NR) yielding plant. Cultivating TKS with a high NR content is an important breeding target, and developing molecular markers related to NR content can effectively accelerate the breeding process of TKS.

RESULTS

To construct a high-density SNP genetic map and uncover genomic regions related to the NR content in TKS, an F₁ mapping population of TKS was constructed by crossing two parents (l66 and X51) with significant differences in NR contents. The NR content of the F₁ plants ranged from 0.30 to 15.14% and was distributed normally with a coefficient of variation of 47.61%, indicating quantitative trait inheritance. Then, employing whole-genome resequencing (WGR), a TKS genetic linkage map of 12,680 bin markers comprising 322,439 SNPs was generated. Based on the genetic map and NR content of the F₁ population, six quantitative trait loci (QTLs) for NR content with LOD > 4.0 were identified on LG01/Chr01 and LG06/Chr06. Of them, the 2.17 Mb genomic region between qHRC-C6-1 and qHRC-C6-2 on ChrA06, with 65.62% PVE in total, was the major QTL region. In addition, the six QTLs have significant additive genetic effects on NR content and could be used to develop markers for marker-assisted selection (MAS) in TKS with a high NR content.

CONCLUSION

This work constructed the first high-density TKS genetic map and identified the QTLs and genomic regions controlling the NR content, which provides useful information for fine mapping, map-based cloning, and MAS in TKS.

Assessing the Role of Ancestral Fragments and Selection Signatures by Whole-Genome Scanning in Dehong Humped Cattle at the China-Myanmar Border

Dehong humped cattle are precious livestock resources of Yunnan Province, China; they have typical zebu traits. Here, we investigated their genetic characteristics using whole-genome resequencing data of Dehong humped animals (n = 18). When comparing our data with the publicly-available data, we found that Dehong humped cattle have high nucleotide diversity. Based on clustering models in a population structure analysis, Dehong humped cattle had a mutual genome ancestor with Chinese and Indian indicine cattle. While using the RFMix method, it is speculated that the body sizes of Dehong humped cattle were influenced by the Chinese indicine segments and that the immune systems of Dehong humped cattle were affected by additional ancestral segments (Indian indicine). Furthermore, we explored the position selection regions harboring genes in the Dehong humped cattle, which were related to heat tolerance (FILIP1L, ABHD6) and immune responses (GZMM, PRKCZ, STOML2, LRBA, PIK3CD). Notably, missense mutations were detected in the candidate gene ABHD6 (c.870C>A p.Asp290Glu; c.987C>A p.Ser329Arg). The missense mutations may have implications for Dehong humped cattle adaptation to hot environments. This study provides valuable genomic resource data at the genome-wide level and paves the way for future genetic breeding work in the Dehong humped cattle.

High Genetic Connectivity Inferred from Whole-Genome Resequencing Provides Insight into the Phylogeographic Pattern of Larimichthys polyactis

Compared with terrestrial biota, marine fishes usually present lower genetic differentiation among different geographical populations because of high-level gene flow and lack of physical barriers. Understanding the genetic structure of marine fishes is essential for dividing management unit and making reasonable protection measures. The small yellow croaker (Larimichthys polyactis) belongs to the family Sciaenidae, which is an economic fish and widely distributed in the Western Pacific. To delineate genetic diversity and phylogeographic pattern, whole-genome resequencing was used to evaluate genetic connectivity, genetic diversity, and spatial pattern of L. polyactis for the first time. We obtained 6,645,711 high-quality single nucleotide polymorphisms (SNPs) markers from 40 L. polyactis individuals. The phylogenetic analysis, STRUCTURE, principal component analysis, and F_st results all indicated that no genetic structure consistent with the distribution pattern was found. This result revealed high genetic connectivity of L. polyactis in different sampling sites. High genetic diversity was also detected, indicating that there was sufficient evolutionary potential to maintain its effective population size. Besides, a similar result of high genetic connectivity and genetic diversity was also detected by mitochondria DNA marker. Our study demonstrated the persistence of high levels of genetic connectivity and a lack of population structure across L. polyactis in different sea areas. This study aimed to analyze the division of population structure and the reason for the decline and not exhaustion of L. polyactis resource on a genetic level.

Whole-genome resequencing of wild and cultivated cannabis reveals the genetic structure and adaptive selection of important traits

BACKGROUND

Cannabis is an important industrial crop species whose fibre, seeds, flowers and leaves are widely used by humans. The study of cannabinoids extracted from plants has been popular research topic in recent years. China is one of the origins of cannabis and one of the few countries with wild cannabis plants. However, the genetic structure of Chinese cannabis and the degree of adaptive selection remain unclear.

RESULTS

The main morphological characteristics of wild cannabis in China were assessed. Based on whole-genome resequencing SNPs, Chinese cannabis could be divided into five groups in terms of geographical source and ecotype: wild accessions growing in the northwestern region; wild accessions growing in the northeastern region; cultivated accessions grown for fibre in the northeastern region; cultivated accessions grown for seed in northwestern region, and cultivated accessions in southwestern region. We further identified genes related to flowering time, seed germination, seed size, embryogenesis, growth, and stress responses selected during the process of cannabis domestication. The expression of flowering-related genes under long-day (LD) and short-day (SD) conditions showed that Chinese cultivated cannabis is adapted to different photoperiods through the regulation of Flowering locus T-like (FT-like) expression.

CONCLUSION

This study clarifies the genetic structure of Chinese cannabis and offers valuable genomic resources for cannabis breeding.

Deep whole-genome resequencing sheds light on the distribution and effect of amphioxus SNPs

Background

Amphioxus is a model organism for vertebrate evolutionary research. The significant contrast between morphological phenotypic similarity and high-level genetic polymorphism among amphioxus populations has aroused scientists' attention. Here we resequenced 21 amphioxus genomes to over 100X depth and mapped them to a haploid reference.

Results

More than 11.5 million common SNPs were detected in the amphioxus population, which mainly affect genes enriched in ion transport, signal transduction and cell adhesion, while protein structure analysis via AlphaFold2 revealed that these SNPs fail to bring effective structural variants.

Conclusions

Our work provides explanation for “amphioxus polymorphism paradox” in a micro view, and generates an enhanced genomic dataset for amphioxus research.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-022-01038-w.

Development of a versatile resource for post-genomic research through consolidating and characterizing 1500 diverse wild and cultivated soybean genomes

BACKGROUND

With advances in next-generation sequencing technologies, an unprecedented amount of soybean accessions has been sequenced by many individual studies and made available as raw sequencing reads for post-genomic research.

RESULTS

To develop a consolidated and user-friendly genomic resource for post-genomic research, we consolidated the raw resequencing data of 1465 soybean genomes available in the public and 91 highly diverse wild soybean genomes newly sequenced. These altogether provided a collection of 1556 sequenced genomes of 1501 diverse accessions (1.5 K). The collection comprises of wild, landraces and elite cultivars of soybean that were grown in East Asia or major soybean cultivating areas around the world. Our extensive sequence analysis discovered 32 million single nucleotide polymorphisms (32mSNPs) and revealed a SNP density of 30 SNPs/kb and 12 non-synonymous SNPs/gene reflecting a high structural and functional genomic diversity of the new collection. Each SNP was annotated with 30 categories of structural and/or functional information. We further identified paired accessions between the 1.5 K and 20,087 (20 K) accessions in US collection as genomic "equivalent" accessions sharing the highest genomic identity for minimizing the barriers in soybean germplasm exchange between countries. We also exemplified the utility of 32mSNPs in enhancing post-genomics research through in-silico genotyping, high-resolution GWAS, discovering and/or characterizing genes and alleles/mutations, identifying germplasms containing beneficial alleles that are potentially experiencing artificial selection.

CONCLUSION

The comprehensive analysis of publicly available large-scale genome sequencing data of diverse cultivated accessions and the newly in-house sequenced wild accessions greatly increased the soybean genome-wide variation resolution. This could facilitate a variety of genetic and molecular-level analyses in soybean. The 32mSNPs and 1.5 K accessions with their comprehensive annotation have been made available at the SoyBase and Ag Data Commons. The dataset could further serve as a versatile and expandable core resource for exploring the exponentially increasing genome sequencing data for a variety of post-genomic research.

Landscape and climatic variations shaped secondary contacts amid barn owls of the Western Palearctic

The combined actions of climatic variations and landscape barriers shape the history of natural populations. When organisms follow their shifting niches, obstacles in the landscape can lead to the splitting of populations, on which evolution will then act independently. When two such populations are reunited, secondary contact occurs in a broad range of admixture patterns, from narrow hybrid zones to the complete dissolution of lineages. A previous study suggested that barn owls colonized the Western Palearctic after the last glaciation in a ring-like fashion around the Mediterranean Sea, and conjectured an admixture zone in the Balkans. Here, we take advantage of whole-genome sequences of 94 individuals across the Western Palearctic to reveal the complex history of the species in the region using observational and modeling approaches. Even though our results confirm that two distinct lineages colonized the region, one in Europe and one in the Levant, they suggest that it predates the last glaciation and identify a secondary contact zone between the two in Anatolia. We also show that barn owls recolonized Europe after the glaciation from two distinct glacial refugia: a previously identified western one in Iberia and a new eastern one in Italy. Both glacial lineages now communicate via eastern Europe, in a wide and permeable contact zone. This complex history of populations enlightens the taxonomy of Tyto alba in the region, highlights the key role played by mountain ranges and large water bodies as barriers and illustrates the power of population genomics in uncovering intricate demographic patterns.

Genome-wide scan for selective footprints and genes related to cold tolerance in Chantecler chickens

The Chantecler chicken, a unique Canadian indigenous breed, is well adapted to extremely cold environments. However, its genetic characteristics have not been well studied. Here, we analyzed the whole genomes of 10 Chantecler chickens and 121 worldwide chickens, which indicated that Chantecler chickens were derived from commercial chickens and exhibit a high level of inbreeding. Based on a genome-wide scan, we identified two vital candidate regions containing ME3 and ZNF536, which are related to fat metabolism and nervous system in cold adaptation, respectively. We also found that the genetic mechanism of cold adaptation in Chantecler chickens differed from that of chickens from other cold regions, such as northern China. Our study indicated that specialized commercial chickens in the early 20th century contained sufficient genetic diversity to adapt to extreme cold environments over a very short time. These findings enrich our understanding of the adaptive potential of commercial species.

Whole-genome resequencing of Japanese whiting ( Sillago japonica) provide insights into local adaptations

The genetic adaptations of various organisms to heterogeneous environments in the northwestern Pacific remain poorly understood. Heterogeneous genomic divergence among populations may reflect environmental selection. Advancing our understanding of the mechanisms by which organisms adapt to different temperatures in response to climate change and predicting the adaptive potential and ecological consequences of anthropogenic global warming are critical. We sequenced the whole genomes of Japanese whiting ( Sillago japonica) specimens collected from different latitudinal locations along the coastal waters of China and Japan to detect possible thermal adaptations. Using population genomics, a total of 5.48 million single nucleotide polymorphisms (SNPs) from five populations revealed a complete genetic break between the Chinese and Japanese groups, which was attributed to both geographic distance and local adaptation. The shared natural selection genes between two isolated populations (i.e., Zhoushan and Ise Bay/Tokyo Bay) indicated possible parallel evolution at the genetic level induced by temperature. These genes also indicated that the process of temperature selection on isolated populations is repeatable. Moreover, we observed natural candidate genes related to membrane fluidity, possibly underlying adaptation to cold environmental stress. These findings advance our understanding of the genetic mechanisms underlying the rapid adaptations of fish species. Species distribution projection models suggested that the Chinese and Japanese groups may have different responses to future climate change, with the former expanding and the latter contracting. The findings of this study enhance our understanding of genetic differentiation and adaptation to changing environments.

Whole genome resequencing data sets of different species from Pistacia genus

OBJECTIVES

Pistacia genus belongs to the flowering plants in the cashew family and contains at least 11 species. The whole-genome resequencing data of different species from Pistacia genus are described herein. The data reported here will be useful for better understand the adaptive evolution, demographic history, genetic diversity, population structure, and domestication of pistachio.

DATA DESCRIPTION

Genomic DNA was isolated from fresh leaves and used to construct libraries with insert size of 350 bp. Sequence libraries were made and sequenced on the Illumina Hiseq 4000 platform to produce 150 bp paired-end reads. A total number of 4,851,118,730 billion reads (ranging from 33,305,900 to 34,990,618 reads per sample) were created across all samples. We produced a total of 727.67 Gbp data which have been deposited in the Genome Sequence Archive (GSA) database with the Accession of CRA000978. All of the data are also available as the sequence read archive (SRA) format in the National Center for Biotechnology Information (NCBI) with identifier of SRP189222, mirroring our deposited data in GSA.

Genomic analyses reveal distinct genetic architectures and selective pressures in Chinese donkeys

Donkey (Equus asinus) is an important livestock animal in China because of its draft and medicinal value. After a long period of natural and artificial selection, the variety and phenotype of donkeys have become abundant. We clarified the genetic and demographic characteristics of Chinese domestic donkeys and the selection pressures by analyzing 78 whole genomes from 12 breeds. According to population structure, most Chinese domestic donkeys showed a dominant ancestral type. However, the Chinese donkeys still represented a significant geographical distribution trend. In the selective sweep, gene annotation, functional enrichment, and differential expression analyses between large and small donkey groups, we identified selective signals, including NCAPG and LCORL, which are related to rapid growth and large body size. Our findings elucidate the evolutionary history and formation of different donkey breeds and provide theoretical insights into the genetic mechanism underlying breed characteristics and molecular breeding programs of donkey clades.

Identification of QTLs associated with potassium use efficiency and underlying candidate genes by whole-genome resequencing of two parental lines in Brassica napus

Breeding crops that acquire and/or utilize potassium (K) more effectively could reduce the use of K fertilizers. Sixteen traits affecting K use efficiency (KUE) at the seedling stage were investigated in a B. napus double haploid population grown at an optimal K supply (OK) and a low K supply (LK) in a hydroponic culture system. In total, 50 and 62 QTLs associated with these traits were identified at OK and LK, respectively. A total of 25 orthologues of 23 Arabidopsis genes regulating K transport were identified in the confidence intervals of nine QTLs impacting shoot dry weight at LK, and 22 of these showed variations in coding sequences and/or exhibited significant differences in mRNA abundances in roots at LK between the two parental lines. This study provided insights to the genetic basis of KUE in B. napus, which will accelerate the breeding of K-efficient rapeseed cultivars by marker-assisted selection.

A global analysis of CNVs in Chinese indigenous fine-wool sheep populations using whole-genome resequencing

Background

Copy number variation (CNV) is an important source of genetic variation that has a significant influence on phenotypic diversity, economically important traits and the evolution of livestock species. In this study, the genome-wide CNV distribution characteristics of 32 fine-wool sheep from three breeds were analyzed using resequencing.

Results

A total of 1,747,604 CNVs were detected in this study, and 7228 CNV regions (CNVR) were obtained after merging overlapping CNVs; these regions accounted for 2.17% of the sheep reference genome. The average length of the CNVRs was 4307.17 bp. “Deletion” events took place more frequently than “duplication” or “both” events. The CNVRs obtained overlapped with previously reported sheep CNVRs to variable extents (4.39–55.46%). Functional enrichment analysis showed that the CNVR-harboring genes were mainly involved in sensory perception systems, nutrient metabolism processes, and growth and development processes. Furthermore, 1855 of the CNVRs were associated with 166 quantitative trait loci (QTL), including milk QTLs, carcass QTLs, and health-related QTLs, among others. In addition, the 32 fine-wool sheep were divided into horned and polled groups to analyze for the selective sweep of CNVRs, and it was found that the relaxin family peptide receptor 2 (RXFP2) gene was strongly influenced by selection.

Conclusions

In summary, we constructed a genomic CNV map for Chinese indigenous fine-wool sheep using resequencing, thereby providing a valuable genetic variation resource for sheep genome research, which will contribute to the study of complex traits in sheep.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07387-7.

Whole-genome resequencing reveals diversity and selective signals in Longlin goat

The Longlin goat is one of the most valuable livestock species in Guangxi Autonomous Region of China, but its genomic diversity and selective signals are not clearly elucidated. Here we compared 20 genomes of Longlin goat to 66 genomes of other seven goat breeds worldwide to analyze patterns of Longlin goat genetic variation. We found the lowest linkage disequilibrium at the large distances between SNPs associated with the highest effective population size in the recent generations ago in Longlin goat. The eight goat breeds could be divided into Euro-African and East Asian goat population. Interestingly, like East Asian taurine, the same two migration phases might have occurred in the history of East Asian goat. More importantly, we identified selective signals implicated in immune resistance to disease, especially for skin disease, in Longlin goat. Our findings will not only help understand the evolutionary history and breed characteristic but can provide valuable resources for conservation of germplasm resources and implementation of crossbreeding programs.

Genome resequencing data for Iranian local dogs and wolves

OBJECTIVE

The data provided herein represent the whole-genome resequencing data related to three wolves and three Iranian local dogs. The understanding of genome evolution during animal domestication is an interesting subject in genome biology. Dog is an excellent model for understanding of domestication due to its considerable variety of behavioral and physical traits. The Zagros area of current day Iran has been identified as one of the initial centers of animal domestication. The availability of the complete genome sequences of Iranian local canids can be a valuable resource for researchers to address questions and testing hypotheses on the dog domestication process.

DATA DESCRIPTION

We collected blood samples from six Iranian local canids including two hunting dogs (Saluki breed), a mastiff dog (Qahderijani ecotype) and three wolves. We extracted genomic DNA from blood samples. Sequence data were produced using the Illumina HiSeq 2500 system. All sequence data are available in the National Genomics Data Center (NGDC), Genome Sequence Archive (GSA) database under the accession of CRA001324 and the National Center for Biotechnology Information (NCBI) under the accession of PRJNA639312. The short-read sequences with the mean depth of 16X were aligned to the dog reference genome (CanFam3.1) and achieved 99% coverage of the reference assembly. The obtained information from this experiment will be useful in evolutionary biology.

Dietary antioxidant vitamin C influences the evolutionary path of insecticide resistance in Drosophila melanogaster

Herbivorous insects encounter a variety of toxic environmental substances ranging from ingested plant defensive compounds to human-introduced insecticidal agents. Dietary antioxidants are known to reduce the negative physiological impacts of toxins in mammalian systems through amelioration of reactive oxygen-related cellular damage. The analogous impacts to insects caused by multigenerational exposure to pesticides and the effects on adaptive responses within insect populations, however, are currently unknown. To address these research gaps, we used Drosophila as a model system to explore adaptive phenotypic responses to acute dichlorodiphenyltrichloroethane (DDT) exposure in the presence of the dietary antioxidant vitamin C and to examine the structural genomic consequences of this exposure. DDT resistance increased significantly among four replicates exposed to a low concentration of DDT for 10 generations. In contrast, dietary intake of vitamin C significantly reduced DDT resistance after mutigenerational exposure to the same concentration of DDT. As to the genomic consequences, no significant differences were predicted in overall nucleotide substitution rates across the genome between any of the treatments. Despite this, replicates exposed to a low concentration of DDT without vitamin C showed the highest number of synonymous and non-synonymous variants (3196 in total), followed by the DDT plus vitamin C (1174 in total), and vitamin C alone (728 in total) treatments. This study demonstrates the potential role of diet (specifically, antioxidant intake) on adaptive genome responses, and thus on the evolution of pesticide resistance within insect populations.

Genomic variation between PRSV resistant transgenic SunUp and its progenitor cultivar Sunset

BACKGROUND

The safety of genetically transformed plants remains a subject of scrutiny. Genomic variants in PRSV resistant transgenic papaya will provide evidence to rationally address such concerns.

RESULTS

In this study, a total of more than 74 million Illumina reads for progenitor 'Sunset' were mapped onto transgenic papaya 'SunUp' reference genome. 310,364 single nucleotide polymorphisms (SNPs) and 34,071 small Inserts/deletions (InDels) were detected between 'Sunset' and 'SunUp'. Those variations have an uneven distribution across nine chromosomes in papaya. Only 0.27% of mutations were predicted to be high-impact mutations. ATP-related categories were highly enriched among these high-impact genes. The SNP mutation rate was about 8.4 × 10^- 4 per site, comparable with the rate induced by spontaneous mutation over numerous generations. The transition-to-transversion ratio was 1.439 and the predominant mutations were C/G to T/A transitions. A total of 3430 nuclear plastid DNA (NUPT) and 2764 nuclear mitochondrial DNA (NUMT) junction sites have been found in 'SunUp', which is proportionally higher than the predicted total NUPT and NUMT junction sites in 'Sunset' (3346 and 2745, respectively). Among all nuclear organelle DNA (norgDNA) junction sites, 96% of junction sites were shared by 'SunUp' and 'Sunset'. The average identity between 'SunUp' specific norgDNA and corresponding organelle genomes was higher than that of norgDNA shared by 'SunUp' and 'Sunset'. Six 'SunUp' organelle-like borders of transgenic insertions were nearly identical to corresponding sequences in organelle genomes (98.18 ~ 100%). None of the paired-end spans of mapped 'Sunset' reads were elongated by any 'SunUp' transformation plasmid derived inserts. Significant amounts of DNA were transferred from organelles to the nuclear genome during bombardment, including the six flanking sequences of the three transgenic insertions.

CONCLUSIONS

Comparative whole-genome analyses between 'SunUp' and 'Sunset' provide a reliable estimate of genome-wide variations and evidence of organelle-to-nucleus transfer of DNA associated with biolistic transformation.

Comparative population genomics reveals genetic divergence and selection in lotus, Nelumbo nucifera

Background

Lotus (Nelumbo nucifera) is an aquatic plant with important agronomic, horticulture, art and religion values. It was the basal eudicot species occupying a critical phylogenetic position in flowering plants. After the domestication for thousands of years, lotus has differentiated into three cultivated types -flower lotus, seed lotus and rhizome lotus. Although the phenotypic and genetic differentiations based on molecular markers have been reported, the variation on whole-genome level among the different lotus types is still ambiguous.

Results

In order to reveal the evolution and domestication characteristics of lotus, a total of 69 lotus accessions were selected, including 45 cultivated accessions, 22 wild sacred lotus accessions, and 2 wild American lotus accessions. With Illumina technology, the genomes of these lotus accessions were resequenced to > 13× raw data coverage. On the basis of these genomic data, 25 million single-nucleotide polymorphisms (SNPs) were identified in lotus. Population analysis showed that the rhizome and seed lotus were monophyletic and genetically homogeneous, whereas the flower lotus was biphyletic and genetically heterogeneous. Using population SNP data, we identified 1214 selected regions in seed lotus, 95 in rhizome lotus, and 37 in flower lotus. Some of the genes in these regions contributed to the essential domestication traits of lotus. The selected genes of seed lotus mainly affected lotus seed weight, size and nutritional quality. While the selected genes were responsible for insect resistance, antibacterial immunity and freezing and heat stress resistance in flower lotus, and improved the size of rhizome in rhizome lotus, respectively.

Conclusions

The genome differentiation and a set of domestication genes were identified from three types of cultivated lotus- flower lotus, seed lotus and rhizome lotus, respectively. Among cultivated lotus, flower lotus showed the greatest variation. The domestication genes may show agronomic importance via enhancing insect resistance, improving seed weight and size, or regulating lotus rhizome size. The domestication history of lotus enhances our knowledge of perennial aquatic crop evolution, and the obtained dataset provides a basis for future genomics-enabled breeding.

Genome-wide association study and protein network analysis for understanding candidate genes involved in root development at the rapeseed seedling stage

Root system is essential for plants to absorb water and nutrients. The root related traits are complex quantitative traits and regulated by genetic control. Here, we used two association mapping panels to perform a genome-wide association study (GWAS) on seven root related traits in Brassica napus at the seedling stage and obtained 27 SNP loci significantly associated with the phenotypes. We further conducted a genome-wide LD block analysis of the candidate peak regions and obtained 295 candidate genes with high association peaks across seven phenotypes in LD region. In addition, a protein interaction network using the candidate genes identified here was constructed, and 113 genes were associated. Seven genes, BnaA03g47330D, BnaC09g16810D, BnaA06g22840D, BnaA03g28390D, BnaA08g19920D, BnaA03g28930D and BnaA03g11440D were in a large cluster, and may play important roles in interacting with other related genes. Our data may provide resources for molecular breeding and functional analysis of root growth and development in rapeseed.

Whole-genome sequencing identifies potential candidate genes for reproductive traits in pigs

Reproductive performance is a complex quantitative trait, that is determined by multiple genes, regulatory pathways and environmental factors. A list of major genes with large effect have been detected, although multiple QTLs are identified. To identify candidate genes for pig prolificacy, whole genome variants from five high- and five low-prolificacy Yorkshire sows were collected using whole-genome resequencing. A total of 13,955,609 SNPs and 2,666,366 indels were detected across the genome. Common differential SNPs and indels were identified between the two groups of sows. Genes encoding components of the TGF-beta signaling pathway were enriched with the variations, including BMP5, BMP6, BMP7, ACVR1, INHBA, ZFYVE9, TGFBR2, DCN, ID4, BAMBI, and ACVR2A. Several differential variants within these genes related to reproductive traits were identified to be associated with litter size. A comparison of selective regions and published QTL data suggests that NEDD9, SLC39A11, SNCA, and UNC5D are candidate genes for reproduction traits.

Development of sequence-based markers for seed protein content in pigeonpea

Pigeonpea is an important source of dietary protein to over a billion people globally, but genetic enhancement of seed protein content (SPC) in the crop has received limited attention for a long time. Use of genomics-assisted breeding would facilitate accelerating genetic gain for SPC. However, neither genetic markers nor genes associated with this important trait have been identified in this crop. Therefore, the present study exploited whole genome re-sequencing (WGRS) data of four pigeonpea genotypes (~ 12X coverage) to identify sequence-based markers and associated candidate genes for SPC. By combining a common variant filtering strategy on available WGRS data with knowledge of gene functions in relation to SPC, 108 sequence variants from 57 genes were identified. These genes were assigned to 19 GO molecular function categories with 56% belonging to only two categories. Furthermore, Sanger sequencing confirmed presence of 75.4% of the variants in 37 genes. Out of 30 sequence variants converted into CAPS/dCAPS markers, 17 showed high level of polymorphism between low and high SPC genotypes. Assay of 16 of the polymorphic CAPS/dCAPS markers on an F₂ population of the cross ICP 5529 (high SPC) × ICP 11605 (low SPC), resulted in four of the CAPS/dCAPS markers significantly (P < 0.05) co-segregated with SPC. In summary, four markers derived from mutations in four genes will be useful for enhancing/regulating SPC in pigeonpea crop improvement programs.

Genomic resequencing combined with quantitative proteomic analyses elucidate the survival mechanisms of Lactobacillus plantarum P-8 in a long-term glucose-limited experiment

Lactobacillus plantarum, commonly isolated from plant material, is widely used to produce various types of fermented foods. However, nutrient-limiting conditions are often encountered during industrial applications. The present study aimed to investigate the response of L. plantarum P-8 to glucose-limited conditions in a long-term experiment. Genotypic and proteomic changes in L. plantarum P-8 were monitored over 3 years in glucose-limited and glucose-normal media using whole-genome resequencing and tandem mass tag-based quantitative proteomic analysis. Results showed that L. plantarum employed numerous survival mechanisms, including alteration of the cell envelope, activation of the PTS system, accumulation and consumption of amino acids, increase in the metabolism of carbohydrates (via glycolysis, citric acid cycle, and pyruvate metabolism), and increase in the production of ATP in response to glucose starvation. This study demonstrates the feasibility of experimental evolution of L. plantarum P-8, while whole-genome resequencing of adapted isolates provided clues toward bacterial functions involved and a deeper mechanistic understanding of the adaptive response of L. plantarum to glucose-limited conditions.

SIGNIFICANCE

We have conducted a 3-year experiment monitoring genotypic and proteomic changes in Lactobacillus plantarum P-8 in glucose-limited and glucose-normal media. Whole-genome resequencing and tandem mass tag-based quantitative proteomics were performed for analyzing genomic evolution of L. plantarum P-8 in glucose-limited and glucose-normal conditions. In addition, differential expressed proteins in all generations between these two conditions were identified and functions of these proteins specific to L group were predicted. L. plantarum employed numerous survival mechanisms, including alteration of the cell envelope, activation of the PTS system, accumulation and consumption of amino acids, increase in the metabolism of carbohydrates (glycolysis, citric acid cycle, and pyruvate metabolism), and increase in the production of ATP in response to glucose starvation.

Forward genetics screen coupled with whole-genome resequencing identifies novel gene targets for improving heterologous enzyme production in Aspergillus niger

Plant biomass, once reduced to its composite sugars, can be converted to fuel substitutes. One means of overcoming the recalcitrance of lignocellulose is pretreatment followed by enzymatic hydrolysis. However, currently available commercial enzyme cocktails are inhibited in the presence of residual pretreatment chemicals. Recent studies have identified a number of cellulolytic enzymes from bacteria that are tolerant to pretreatment chemicals such as ionic liquids. The challenge now is generation of these enzymes in copious amounts, an arena where fungal organisms such as Aspergillus niger have proven efficient. Fungal host strains still need to be engineered to increase production titers of heterologous protein over native enzymes, which has been a difficult task. Here, we developed a forward genetics screen coupled with whole-genome resequencing to identify specific lesions responsible for a protein hyper-production phenotype in A. niger. This strategy successfully identified novel targets, including a low-affinity glucose transporter, MstC, whose deletion significantly improved secretion of recombinant proteins driven by a glucoamylase promoter.

Genetic diversity and natural selection in wild fruit flies revealed by whole-genome resequencing

We characterized 26 wild fruit flies comparative population genomics from six different altitude and latitude locations by whole genome resequencing. Genetic diversity was relatively higher in Ganzi and Chongqing populations. We also found 13 genes showing selection signature between different altitude flies and variants related to hypoxia and temperature stimulus, were preferentially selected during the flies evolution. One of the most striking selective sweeps found in all high altitude flies occurred in the region harboring Hsp70Aa and Hsp70Ab on chromosome 3R. Interestingly, these two genes are involved in GO terms including response to hypoxia, unfolded protein, temperature stimulus, heat, oxygen levels. Mutation in HPH gene, a candidate gene in the hypoxia inducible factor pathway, might contributes to hypoxic high-altitude adaptation. Intriguingly, some of the selected genes, primarily utilized in humans, were involved in the response to hypoxia, which could imply a conserved molecular mechanisms underlying high-altitude adaptation between insects and humans.

Fine-mapping and candidate gene analysis of the Brassica juncea white-flowered mutant Bjpc2 using the whole-genome resequencing

Flower color in Brassica spp. is an important trait and considered a major visual signal for insect-pollinated plants. In previous study, we isolated and identified two genes (BjPC1 and BjPC2) that control the flower color in B. juncea, and mapped BjPC1 to a 0.13-cM region. In this study, we report the fine-mapping and candidate analysis of BjPC2. We conducted whole-genome resequencing, using bulked segregant analysis (BSA) to determine the BjPC2 candidate intervals. Crossing, allelism testing, and repeated full-sib mating were used to generate XG3, a near isogenic line (NIL) population that segregated on the BjPC2 locus. Through a genome-wide comparison of single nucleotide polymorphism (SNP) profiles between the yellow- and white-flowered bulks, a candidate interval for BjPC2 was identified on chromosome B04 (2.45 Mb). The BjPC2 linkage map was constructed with the newly developed simple sequence repeat (SSR) markers in the candidate interval to narrow the candidate BjPC2 region to 31-kb. Expression profiling and RNA-seq analysis partially confirmed that the AtPES2 homolog, BjuB027334 is the most promising candidate gene for BjPC2. Furthermore, analyses with high pressure liquid chromatography and transmission electron microscopy demonstrated that BjPC2 might be important in xanthophyll esterification, a process that limits xanthophyll degradation and increases sequestration. Overall, we mapped the BjPC2 to a 31-kb region on the B04 in B. juncea and identified BjuB027334 as a valuable candidate gene. Our results provide a basis for understanding the molecular mechanisms underlying the white-flowered trait and for molecular marker-assisted selection of flower color in B. juncea breeding.

A GJA9 frameshift variant is associated with polyneuropathy in Leonberger dogs

Background

Many inherited polyneuropathies (PN) observed in dogs have clinical similarities to the genetically heterogeneous group of Charcot-Marie-Tooth (CMT) peripheral neuropathies in humans. The canine disorders collectively show a variable expression of progressive clinical signs and ages of onset, and different breed prevalences. Previously in the Leonberger breed, a variant highly associated with a juvenile-onset PN was identified in the canine orthologue of a CMT-associated gene. As this deletion in ARHGEF10 (termed LPN1) does not explain all cases, PN in this breed may encompass variants in several genes with similar clinical and histopathological features.

Results

A genome-wide comparison of 173 k SNP genotypes of 176 cases, excluding dogs homozygous for the ARHGEF10 variant, and 138 controls, was carried out to detect further PN-associated variants. A single suggestive significant association signal on CFA15 was found. The genome of a PN-affected Leonberger homozygous for the associated haplotype was sequenced and variants in the 7.7 Mb sized critical interval were identified. These variants were filtered against a database of variants observed in 202 genomes of various dog breeds and 3 wolves, and 6 private variants in protein-coding genes, all in complete linkage disequilibrium, plus 92 non-coding variants were revealed. Five of the coding variants were predicted to have low or moderate effect on the encoded protein, whereas a 2 bp deletion in GJA9 results in a frameshift of high impact. GJA9 encodes connexin 59, a connexin gap junction family protein, and belongs to a group of CMT-associated genes that have emerged as important components of peripheral myelinated nerve fibers. The association between the GJA9 variant and PN was confirmed in an independent cohort of 296 cases and 312 controls. Population studies showed a dominant mode of inheritance, an average age of onset of approximately 6 years, and incomplete penetrance.

Conclusions

This GJA9 variant represents a highly probable candidate variant for another form of PN in Leonberger dogs, which we have designated LPN2, and a new candidate gene for CMT disease. To date, approximately every third PN-diagnosed Leonberger dog can be explained by the ARHGEF10 or GJA9 variants, and we assume that additional genetic heterogeneity in this condition exists in the breed.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4081-z) contains supplementary material, which is available to authorized users.

Genome-wide analysis reveals signatures of selection for important traits in domestic sheep from different ecoregions

Background

Throughout a long period of adaptation and selection, sheep have thrived in a diverse range of ecological environments. Mongolian sheep is the common ancestor of the Chinese short fat-tailed sheep. Migration to different ecoregions leads to changes in selection pressures and results in microevolution. Mongolian sheep and its subspecies differ in a number of important traits, especially reproductive traits. Genome-wide intraspecific variation is required to dissect the genetic basis of these traits.

Results

This research resequenced 3 short fat-tailed sheep breeds with a 43.2-fold coverage of the sheep genome. We report more than 17 million single nucleotide polymorphisms and 2.9 million indels and identify 143 genomic regions with reduced pooled heterozygosity or increased genetic distance to each other breed that represent likely targets for selection during the migration. These regions harbor genes related to developmental processes, cellular processes, multicellular organismal processes, biological regulation, metabolic processes, reproduction, localization, growth and various components of the stress responses. Furthermore, we examined the haplotype diversity of 3 genomic regions involved in reproduction and found significant differences in TSHR and PRL gene regions among 8 sheep breeds.

Conclusions

Our results provide useful genomic information for identifying genes or causal mutations associated with important economic traits in sheep and for understanding the genetic basis of adaptation to different ecological environments.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3212-2) contains supplementary material, which is available to authorized users.

Genome-wide resequencing of KRICE_CORE reveals their potential for future breeding, as well as functional and evolutionary studies in the post-genomic era

Background

Rice germplasm collections continue to grow in number and size around the world. Since maintaining and screening such massive resources remains challenging, it is important to establish practical methods to manage them. A core collection, by definition, refers to a subset of the entire population that preserves the majority of genetic diversity, enhancing the efficiency of germplasm utilization.

Results

Here, we report whole-genome resequencing of the 137 rice mini core collection or Korean rice core set (KRICE_CORE) that represents 25,604 rice germplasms deposited in the Korean genebank of the Rural Development Administration (RDA). We implemented the Illumina HiSeq 2000 and 2500 platform to produce short reads and then assembled those with 9.8 depths using Nipponbare as a reference. Comparisons of the sequences with the reference genome yielded more than 15 million (M) single nucleotide polymorphisms (SNPs) and 1.3 M INDELs. Phylogenetic and population analyses using 2,046,529 high-quality SNPs successfully assigned rice accessions to the relevant rice subgroups, suggesting that these SNPs capture evolutionary signatures that have accumulated in rice subpopulations. Furthermore, genome-wide association studies (GWAS) for four exemplary agronomic traits in the KRIC_CORE manifest the utility of KRICE_CORE; that is, identifying previously defined genes or novel genetic factors that potentially regulate important phenotypes.

Conclusion

This study provides strong evidence that the size of KRICE_CORE is small but contains high genetic and functional diversity across the genome. Thus, our resequencing results will be useful for future breeding, as well as functional and evolutionary studies, in the post-genomic era.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2734-y) contains supplementary material, which is available to authorized users.