Whole-genome Resequencing Data Research Articles

Revealing genes related teat number traits via genetic variation in Yorkshire pigs based on whole-genome sequencing

BackgroundTeat number is one of the most important indicators to evaluate the lactation performance of sows, and increasing the teat number has become an important method to improve the economic efficiency of farms. Therefore, it is particularly important to deeply analyze the genetic mechanism of teat number traits in pigs. In this study, we detected Single Nucleotide Ploymorphism (SNP), Insertion-Deletion (InDel) and Structural variant (SV) by high-coverage whole-genome resequencing data, and selected teat number at birth and functional teat number as two types of teat number traits for genome-wide association study (GWAS) to reveal candidate genes associated with pig teat number traits.ResultsIn this study, we used whole genome resequencing data from 560 Yorkshire sows to detect SNPs, InDels and SVs, and performed GWAS for the traits of born teat number and functional teat number, and detected a total of 85 significant variants and screened 214 candidate genes, including HEG1, XYLT1, SULF1, MUC13, VRTN, RAP1A and NPVF. Among them, HEG1 and XYLT1 were the new candidate genes in this study. The co-screening and population validation of multiple traits suggested that HEG1 may have a critical effect on the born teat number.ConclusionOur study shows that more candidate genes associated with pig teat number traits can be identified by GWAS with different variant types. Through large population validation, we found that HEG1 may be a new key candidate gene affecting pig teat number traits. In conclusion, the results of this study provide new information for exploring the genetic mechanisms affecting pig teat number traits and genetic improvement of pigs.

Whole-genome resequencing to investigate the genetic diversity and mechanisms of plateau adaptation in Tibetan sheep

IntroductionTibetan sheep, economically important animals on the Qinghai–Tibet Plateau, have diversified into numerous local breeds with unique characteristics through prolonged environmental adaptation and selective breeding. However, most current research focuses on one or two breeds, and lacks a comprehensive representation of the genetic diversity across multiple Tibetan sheep breeds. This study aims to fill this gap by investigating the genetic structure, diversity and high-altitude adaptation of 6 Tibetan sheep breeds using whole-genome resequencing data.ResultsSix Tibetan sheep breeds were investigated in this study, and whole-genome resequencing data were used to investigate their genetic structure and population diversity. The results showed that the 6 Tibetan sheep breeds exhibited distinct separation in the phylogenetic tree; however, the levels of differentiation among the breeds were minimal, with extensive gene flow observed. Population structure analysis broadly categorized the 6 breeds into 3 distinct ecological types: plateau-type, valley-type and Euler-type. Analysis of unique single-nucleotide polymorphisms (SNPs) and selective sweeps between Argali and Tibetan sheep revealed that Tibetan sheep domestication was associated primarily with sensory and signal transduction, nutrient absorption and metabolism, and growth and reproductive characteristics. Finally, comprehensive analysis of selective sweep and transcriptome data suggested that Tibetan sheep breeds inhabiting different altitudes on the Qinghai–Tibet Plateau adapt by enhancing cardiopulmonary function, regulating body fluid balance through renal reabsorption, and modifying nutrient digestion and absorption pathways.ConclusionIn this study, we investigated the genetic diversity and population structure of 6 Tibetan sheep breeds in Qinghai Province, China. Additionally, we analyzed the domestication traits and investigated the unique adaptation mechanisms residing varying altitudes in the plateau region of Tibetan sheep. This study provides valuable insights into the evolutionary processes of Tibetan sheep in extreme environments. These findings will also contribute to the preservation of genetic diversity and offer a foundation for Tibetan sheep diversity preservation and plateau animal environmental adaptation mechanisms.

Genome-Wide Association Study of Seed Quality and Yield Traits in a Soybean Collection from Southeast Kazakhstan

Soybean (Glycine max (L.) Merr.) is a vital agricultural crop and a key source of protein and oil for food and feed production. The search for new genetic factors affecting the main agronomic traits of soybean is a significant step for efficient breeding strategies. This study aimed to identify marker–trait associations (MTAs) for seed protein and oil content and yield by conducting a genome-wide association study (GWAS). The collection of 252 soybean accessions of five different origins was analyzed over a period of five years. The GWAS was conducted using 44,385 SNP markers extracted from whole-genome resequencing data using Illumina HiSeq X Ten. The multiple-locus mixed linear model (MLMM) facilitated the identification of 38 stable MTAs: nine for protein content, nine for oil content, seven for the number of fertile nodes, six for the number of seeds per plant, four for thousand seeds weight, and three for yield per plant. Fifteen of these MTAs are presumed to be novel, with one linked to seed protein content, three linked to seed oil content, and the remaining MTAs linked to yield-related traits. These findings offer valuable insights for soybean breeding programs aimed at developing new, competitive cultivars with improved seed quality and yield characteristics.

An integrated microbiome- and metabolome-genome-wide association study reveals the role of heritable ruminal microbial carbohydrate metabolism in lactation performance in Holstein dairy cows

BackgroundDespite the growing number of studies investigating the connection between host genetics and the rumen microbiota, there remains a dearth of systematic research exploring the composition, function, and metabolic traits of highly heritable rumen microbiota influenced by host genetics. Furthermore, the impact of these highly heritable subsets on lactation performance in cows remains unknown. To address this gap, we collected and analyzed whole-genome resequencing data, rumen metagenomes, rumen metabolomes and short-chain fatty acids (SCFAs) content, and lactation performance phenotypes from a cohort of 304 dairy cows.ResultsThe results indicated that the proportions of highly heritable subsets (h2 ≥ 0.2) of the rumen microbial composition (55%), function (39% KEGG and 28% CAZy), and metabolites (18%) decreased sequentially. Moreover, the highly heritable microbes can increase energy-corrected milk (ECM) production by reducing the rumen acetate/propionate ratio, according to the structural equation model (SEM) analysis (CFI = 0.898). Furthermore, the highly heritable enzymes involved in the SCFA synthesis metabolic pathway can promote the synthesis of propionate and inhibit the acetate synthesis. Next, the same significant SNP variants were used to integrate information from genome-wide association studies (GWASs), microbiome-GWASs, metabolome-GWASs, and microbiome-wide association studies (mWASs). The identified single nucleotide polymorphisms (SNPs) of rs43470227 and rs43472732 on SLC30A9 (Zn2+ transport) (P < 0.05/nSNPs) can affect the abundance of rumen microbes such as Prevotella_sp., Prevotella_sp._E15-22, Prevotella_sp._E13-27, which have the oligosaccharide-degradation enzymes genes, including the GH10, GH13, GH43, GH95, and GH115 families. The identified SNPs of chr25:11,177 on 5s_rRNA (small ribosomal RNA) (P < 0.05/nSNPs) were linked to ECM, the abundance alteration of Pseudobutyrivibrio_sp. (a genus that was also showed to be linked to the ECM production via the mWASs analysis), GH24 (lysozyme), and 9,10,13-TriHOME (linoleic acid metabolism). Moreover, ECM, and the abundances of Pseudobutyrivibrio sp., GH24, and 9,10,13-TRIHOME were significantly greater in the GG genotype than in the AG genotype at chr25:11,177 (P < 0.05). By further the SEM analysis, GH24 was positively correlated with Pseudobutyrivibrio sp., which was positively correlated with 9,10,13-triHOME and subsequently positively correlated with ECM (CFI = 0.942).ConclusionOur comprehensive study revealed the distinct heritability patterns of rumen microbial composition, function, and metabolism. Additionally, we shed light on the influence of host SNP variants on the rumen microbes with carbohydrate metabolism and their subsequent effects on lactation performance. Collectively, these findings offer compelling evidence for the host-microbe interactions, wherein cows actively modulate their rumen microbiota through SNP variants to regulate their own lactation performance.C1VA69ck9huqCRa2wNrTKWVideo

Genome-Wide Selection Signals Reveal Candidate Genes Associated with Plateau Adaptation in Tibetan Sheep.

Tibetan sheep have developed unique adaptations for survival in the Qinghai-Tibet Plateau environment. However, the functional genes and molecular mechanisms that regulate hypoxia adaptation have not been fully characterized. In this study, based on the whole-genome resequencing data for Tibetan sheep at different altitudes, the population differentiation index (FST) and nucleotide diversity ratio (θπ ratio) were evaluated in populations of 20 Oula sheep (3501 m altitude, OL), 20 Zashijia sheep (4369 m altitude, ZSJ), and 20 Awang sheep (4643 m altitude, AW) to reveal candidate loci related to high-altitude hypoxia. We screened 728 and 524 candidate genes in the AW vs. OL and ZSJ vs. OL groups, respectively, of which 134 genes were jointly screened. Candidate genes were mainly enriched in the Ras, melanoma, melanogenesis, VEGF, and PPAR signaling pathways. HIF1AN, PDGFA, PDGFD, ANXA2, SOCS2, NOXA1, WNT7B, MMP14, GNG2, ATF6, PGAM2, PPP3R1, GSTCD, and PPARA may play important roles in the high-altitude adaptation of Tibetan sheep. In conclusion, this study provides valuable insights into the genes and molecular mechanisms that underlie high-altitude hypoxia adaptation in Tibetan sheep.

Insights into Adaption and Growth Evolution: Genome-Wide Copy Number Variation Analysis in Chinese Hainan Yellow Cattle Using Whole-Genome Re-Sequencing Data.

Copy number variation (CNV) serves as a crucial source of genomic variation and significantly aids in the mining of genomic information in cattle. This study aims to analyze re-sequencing data from Chinese Hainan yellow cattle, to uncover breed CNV information, and to elucidate the resources of population genetic variation. We conducted whole-genome sequencing on 30 Chinese Hainan yellow cattle, thus generating 814.50 Gb of raw data. CNVs were called using CNVnator software, and subsequent filtering with Plink and HandyCNV yielded 197,434 high-quality CNVs and 5852 CNV regions (CNVRs). Notably, the proportion of deleted sequences (81.98%) exceeded that of duplicated sequences (18.02%), with the lengths of CNVs predominantly ranging between 20 and 500 Kb This distribution demonstrated a decrease in CNVR count with increasing fragment length. Furthermore, an analysis of the population genetic structure using CNVR databases from Chinese, Indian, and European commercial cattle breeds revealed differences between Chinese Bos indicus and Indian Bos indicus. Significant differences were also observed between Hainan yellow cattle and European commercial breeds. We conducted gene annotation for both Hainan yellow cattle and European commercial cattle, as well as for Chinese Bos indicus and Indian Bos indicus, identifying 206 genes that are expressed in both Chinese and Indian Bos indicus. These findings may provide valuable references for future research on Bos indicus. Additionally, selection signatures analysis based on Hainan yellow cattle and three European commercial cattle breeds identified putative pathways related to heat tolerance, disease resistance, fat metabolism, environmental adaptation, candidate genes associated with reproduction and the development of sperm and oocytes (CABS1, DLD, FSHR, HSD17B2, KDM2A), environmental adaptation (CNGB3, FAM161A, DIAPH3, EYA4, AAK1, ERBB4, ERC2), oxidative stress anti-inflammatory response (COMMD1, OXR1), disease resistance (CNTN5, HRH4, NAALADL2), and meat quality (EHHADH, RHOD, GFPT1, SULT1B1). This study provides a comprehensive exploration of CNVs at the molecular level in Chinese Hainan yellow cattle, offering theoretical support for future breeding and selection programs aimed at enhancing qualities of this breed.

Identifying low-density, ancestry-informative SNP markers through whole genome resequencing in Indian, Chinese, and wild yak

The current investigation was undertaken to elucidate the population-stratifying and ancestry-informative markers in Indian, Chinese, and wild yak populations using whole genome resequencing (WGS) analysis while employing various selection strategies (Delta, Pairwise Wright’s Fixation Index-FST, and Informativeness of Assignment) and marker densities (5–25 thousand). The study used WGS data on 105 individuals from three separate yak cohorts i.e., Indian yak (n = 29), Chinese yak (n = 61), and wild yak (n = 15). Variant calling in the GATK program with strict quality control resulted in 1,002,970 high-quality and independent (LD-pruned) SNP markers across the yak autosomes. Analysis was undertaken in toolbox for ranking and evaluation of SNPs (TRES) program wherein three different criteria i.e., Delta, Pairwise Wright’s Fixation Index-FST, and Informativeness of Assignment were employed to identify population-stratifying and ancestry-informative markers across various datasets. The top-ranked 5,000 (5K), 10,000 (10K), 15,000 (15K), 20,000 (20K), and 25,000 (25K) SNPs were identified from each dataset while their composition and performance was assessed using different criteria. The average genomic breed clustering of Indian, Chinese, and wild yak cohorts with full density dataset (105 individuals with 1,002,970 markers) was 81.74%, 80.02%, and 83.62%, respectively. Informativeness of Assignment criterion with 10K density emerged as the best combination for three yak cohorts with 86.94%, 96.46%, and 98.20% clustering for Indian, Chinese, and wild yak, respectively. There was an average increase of 7.56%, 22.72%, and 30.35% in genomic breed clustering scores of Indian, Chinese, and wild yak cohorts over the estimates of the original dataset. The selected markers showed overlap multiple protein-coding genes within a 10 kb window including ADGRB3, ANK1, CACNG7, CALN1, CHCHD2, CREBBP, GLI3, KHDRBS2, and OSBPL10. This is the first report ever on elucidating low-density SNP marker sets with population-stratifying and ancestry-informative properties in three yak groups using WGS data. The results gain significance for application of genomic selection using cost-effective low-density SNP panels in global yak species.

Lineage Differentiation and Genomic Vulnerability in a Relict Tree From Subtropical Forests.

The subtropical forests of East Asia are renowned for their high plant diversity, particularly the abundance of ancient relict species. However, both the evolutionary history of these relict species and their capacity for resilience in the face of impending climatic changes remain unclear. Using whole-genome resequencing data, we investigated the lineage differentiation and demographic history of the relict and endangered tree, Bretschneidera sinensis (Akaniaceae). We employed a combination of population genomic and landscape genomic approaches to evaluate variation in mutation load and genomic offset, aiming to predict how different populations may respond to climate change. Our analysis revealed a profound genomic divergence between the East and West lineages, likely as the result of recurrent bottlenecks due to climatic fluctuations during the glacial period. Furthermore, we identified several genes potentially linked to growth characteristics and hypoxia response that had been subjected to positive selection during the lineage differentiation. Our assessment of genomic vulnerability uncovered a significantly higher mutation load and genomic offset in the edge populations of B. sinensis compared to their core counterparts. This implies that the edge populations are likely to experience the most significant impact from the predicted climate conditions. Overall, our research sheds light on the historical lineage differentiation and contemporary genomic vulnerability of B. sinensis. Broadening our understanding of the speciation history and future resilience of relict and endangered species such as B. sinensis, is crucial in developing effective conservation strategies in anticipation of future climatic changes.

Genomic insights into the specialisation and selection of the Jinding duck

The domestication of ducks represents a pivotal evolutionary shift from the unguided propagation of wild species to deliberate human-mediated selection, culminating in distinct behavioural, morphological, and physiological traits that differentiate domesticated ducks from their wild counterparts. This transition has yielded breeds with traits fine-tuned to specific economic roles, such as egg production, meat yield, or dual-purpose functionality. Duck domestication plays a significant role in poultry production globally, meeting the growing demand for eggs and meat in various regions. Here, we focus on the Jinding Duck (JDD), a breed renowned for its prolific egg-laying traits. Employing whole-genome resequencing data from 325 individuals across five Chinese indigenous duck breeds, we aimed to dissect the unique population structure and assess the genetic diversity within the JDD cohort. The findings reveal the distinct genetic heritage of JDD, diverged from other domesticated breeds, and show a relative paucity of genetic diversity. A salient discovery was a 200 kb genomic interval containing three genes (NCF2, SMG7, and ARPC5) with almost exclusive haplotypes, which were inherited from Anas platyrhynchos or Anas zonorhyncha, impacting the morphological attributes of JDD. The study highlights a c.28G>A non−synonymous mutation in the first exon of the LAMC1 gene, which is potentially influencing feather morphology in JDD. Our findings suggest that unique blue eggshell colouration in JDD is likely attributable to variations within the promoter element of the ABCG2 gene, distinguishing it from other breeds. Moreover, the MAP7 and FHL1 genes emerge as significant factors in the laying performance of JDD. These genetic insights are not only crucial for improving the JDD breed but also provide valuable information that could be applied to duck breeding programmes worldwide, helping enhance productivity and meet international demands for high-efficiency poultry breeds.

Microevolution of longan flowering behaviors from primitive to KClO3 induction is driven by artificial selection for agricultural benefit

Dimocarpus longan is an economically important fruit tree primarily cultivated in various Asian countries and the Indochina region encompassing Thailand, China, Vietnam, Laos and Cambodia. Here, we report an extensive dataset of 606.72Gb of the whole-genome resequencing data obtained from 31 Dimocarpus longan and 1 Dimocarpus obtusus samples performed on Illumina HiSeq PE150, providing the coverage over 30-fold of the reference genome (averaging 15Gb per sample) with an addition of 12 longan accessions from China. A phylogenomic tree inferred from this analyzed dataset (19,270,513 SNPs) was constructed and results indicated that there are three separate monophyletic clades of China-USA (1), Thai (2) and Vietnam longan (3). Interestingly, however, Thai and Vietnam clades appear to be closely linked in genetic relationship based on analyses of karyotype, genomic data and certain morphological characteristics. De novo transcriptome assembly further revealed variations in gene sequences, divergent gene expression and candidate genes associated with different flowering behaviors. The greatest change in differential gene expression (51.87%) observed between natural independent and seasonal flowering implied the microevolutionary shift observed in longan is influenced by artificial selection, resulting in gradual changes in flowering behaviors from natural independent flowering to artificial flowering induction using potassium chlorate (KClO3). In summary, the data obtained from this study serves as the essential evidence for elucidating the microevolution of longan and shedding new light on an agronomical application of artificial flowering induction via modulation of KClO3 responsive genes in longan or other fruits in the future.

Population Genetic Characteristics of the Cultivated Coffea arabica with Whole-Genome Resequencing

Coffea arabica L. (C. arabica) is an economically important agricultural crop and the most popular beverage worldwide. To analyze genetic diversity and provide genetic resources for the selection and breeding of superior varieties of C. arabica, 61 cultivated Arabica coffee accessions were analyzed in the study, including 12 resequencing accessions from previous research and 49 accessions that were resequenced in this study. Single nucleotide polymorphisms (SNPs) and insertion–deletions (InDels) were statistically analyzed. Based on SNP variations, a genetic structure analysis, phylogenetic tree construction, and principal component analysis were performed for the 61 coffee accessions. The results showed that a total of 805.46 Gb of raw whole-genome resequencing data was obtained from the 61 coffee accessions, with 781.29 Gb of high-quality sequencing data after filtering. In total, 7,013,820 SNP sites and 1,074,329 InDel sites were detected. The average sequencing depth ranged from 6.69× to 19.35×, and the coverage ranged from 85.49% to 96.43%. The population genetic structure and phylogenetic analysis of the 61 coffee accessions revealed four lineages, suggesting that they had at least four ancestral genetic components. Catimor exhibited the highest genetic diversity, while Geisha had the lowest genetic diversity. The selective sweep analysis indicated that among the selected genes in Catimor, disease-resistance genes were significantly more numerous than in other coffee varieties. The genome resequencing data and genetic markers identified from the 61 cultivated Arabica coffee materials provided insights into the genetic variation in Arabica coffee germplasm and facilitated extensive genetic research.

Genome-wide association mapping identifies novel SNPs for root nodulation and agronomic traits in chickpea.

The chickpea (Cicer arietinum L.) is well-known for having climate resilience and atmospheric nitrogen fixation ability. Global demand for nitrogenous fertilizer is predicted to increase by 1.4% annually, and the loss of billions of dollars in farm profit has drawn attention to the need for alternative sources of nitrogen. The ability of chickpea to obtain sufficient nitrogen via its symbiotic relationship with Mesorhizobium ciceri is of critical importance in determining the growth and production of chickpea. To support findings on nodule formation in chickpea and to map the genomic regions for nodulation, an association panel consisting of 271 genotypes, selected from the global chickpea germplasm including four checks at four locations, was evaluated, and data were recorded for nodulation and 12 yield-related traits. A genome-wide association study (GWAS) was conducted using phenotypic data and genotypic data was extracted from whole-genome resequencing data of chickpea by creating a hap map file consisting of 602,344 single-nucleotide polymorphisms (SNPs) in the working set with best-fit models of association mapping. The GWAS panel was found to be structured with sufficient diversity among the genotypes. Linkage disequilibrium (LD) analysis showed an LD decay value of 37.3 MB, indicating that SNPs within this distance behave as inheritance blocks. A total of 450 and 632 stringent marker-trait associations (MTAs) were identified from the BLINK and FarmCPU models, respectively, for all the traits under study. The 75 novel MTAs identified for nodulation traits were found to be stable. SNP annotations of associated markers were found to be related to various genes including a few auxins encoding as well as nod factor transporter genes. The identified significant MTAs, candidate genes, and associated markers have the potential for use in marker-assisted selection for developing high-nodulation cultivars after validation in the breeding populations.

Unveiling the Genome-Wide Consequences of Range Expansion and Mating System Transitions in Primula vulgaris.

Genetic diversity is heterogeneously distributed among populations of the same species, due to the joint effects of multiple demographic processes, including range contractions and expansions, and mating systems shifts. Here, we ask how both processes shape genomic diversity in space and time in the classical Primula vulgaris model. This perennial herb originated in the Caucasus region and was hypothesized to have expanded westward following glacial retreat in the Quaternary. Moreover, this species is a long-standing model for mating system transitions, exemplified by shifts from heterostyly to homostyly. Leveraging a high-quality reference genome of the closely related Primula veris and whole-genome resequencing data from both heterostylous and homostylous individuals from populations encompassing a wide distribution of P. vulgaris, we reconstructed the demographic history of P. vulgaris. Results are compatible with the previously proposed hypothesis of range expansion from the Caucasus region approximately 79,000 years ago and suggest later shifts to homostyly following rather than preceding postglacial colonization of England. Furthermore, in accordance with population genetic theoretical predictions, both processes are associated with reduced genetic diversity, increased linkage disequilibrium, and reduced efficacy of purifying selection. A novel result concerns the contrasting effects of range expansion versus shift to homostyly on transposable elements, for the former, process is associated with changes in transposable element genomic content, while the latter is not. Jointly, our results elucidate how the interactions among range expansion, transitions to selfing, and Quaternary climatic oscillations shape plant evolution.

Chromosome-level genome assembly and population genomic analysis provide insights into the genetic diversity and adaption of Schizopygopsis younghusbandi on the Tibetan Plateau.

The Yarlung Tsangpo River on the Tibetan Plateau provides a unique natural environment for studying fish evolution and ecology. However, the genomes and genetic diversity of plateau fish species have been rarely reported. Schizopygopsis younghusbandi, a highly specialized Schizothoracine species and economically important fish inhabiting the Yarlung Tsangpo River, is threatened by overfishing and biological invasion. Herein, we generated a chromosome-level genome of S. younghusbandi and whole-genome resequencing data for 59 individuals from six locations of the river. The results showed that the divergence time between S. younghusbandi and other primitive Schizothoracine species was ∼4.2 Mya, coinciding with the major phase of the Neogene Tibetan uplift. The expanded gene families enriched in DNA integration and replication, ion binding and transport, energy storage, and metabolism likely contribute to the adaption of this species. The S. younghusbandi may have diverged from other highly specialized Schizothoracine species in the Zanda basin during the Pliocene epoch, which underwent major population reduction possibly due to the drastic climate change during the last glacial period. Population analysis indicated that the ancient population might have originated upstream before gradually adapting to evolve into the populations inhabiting the mid-stream and downstream regions of the Yarlung Tsangpo River. In conclusion, the chromosome-level genome and population diversity of S. younghusbandi provide valuable genetic resources for the evolution, ecology, and conservation studies of endemic fishes on the Tibetan Plateau.

Genomic Diversity as a Key Conservation Criterion: Proof-of-Concept From Mammalian Whole-Genome Resequencing Data.

Many international, national, state, and local organizations prioritize the ranking of threatened and endangered species to help direct conservation efforts. For example, the International Union for Conservation of Nature (IUCN) assesses the Green Status of species and publishes the influential Red List of threatened species. Unfortunately, such conservation yardsticks do not explicitly consider genetic or genomic diversity (GD), even though GD is positively associated with contemporary evolutionary fitness, individual viability, and with future evolutionary potential. To test whether populations of genome sequences could help improve conservation assessments, we estimated GD metrics from 82 publicly available mammalian datasets and examined their statistical association with attributes related to conservation. We also considered intrinsic biological factors, including trophic level and body mass, that could impact GD and quantified their relative influences. Our results identify key population GD metrics that are both reflective and predictive of IUCN conservation categories. Specifically, our analyses revealed that Watterson's theta (the population mutation rate) and autozygosity (a product of inbreeding) are associated with the current Red List categorization, likely because demographic declines that lead to "listing" decisions also reduce levels of standing genetic variation. We argue that by virtue of this relationship, conservation organizations like IUCN could leverage emerging genome sequence data to help categorize Red List threat rankings (especially in otherwise data-deficient species) and/or enhance Green Status assessments to establish a baseline for future population monitoring. Thus, our paper (1) outlines the theoretical and empirical justification for a new GD-based assessment criterion, (2) provides a bioinformatic pipeline for estimating GD from population genomic data, and (3) suggests an analytical framework that can be used to measure baseline GD while providing quantitative GD context for consideration by conservation authorities.

Contribution of homoeologous exchange to domestication of polyploid Brassica

BackgroundPolyploidy is widely recognized as a significant evolutionary force in the plant kingdom, contributing to the diversification of plants. One of the notable features of allopolyploidy is the occurrence of homoeologous exchange (HE) events between the subgenomes, causing changes in genomic composition, gene expression, and phenotypic variations. However, the role of HE in plant adaptation and domestication remains unclear.ResultsHere we analyze the whole-genome resequencing data from Brassica napus accessions representing the different morphotypes and ecotypes, to investigate the role of HE in domestication. Our findings demonstrate frequent occurrence of HEs in Brassica napus, with substantial HE patterns shared across populations, indicating their potential role in promoting crop domestication. HE events are asymmetric, with the A genome more frequently replacing C genome segments. These events show a preference for specific genomic regions and vary among populations. We also identify candidate genes in HE regions specific to certain populations, which likely contribute to flowering-time diversification across diverse morphotypes and ecotypes. In addition, we assemble a new genome of a swede accession, confirming the HE signals on the genome and their potential involvement in root tuber development. By analyzing HE in another allopolyploid species, Brassica juncea, we characterize a potential broader role of HE in allopolyploid crop domestication.ConclusionsOur results provide novel insights into the domestication of polyploid Brassica species and highlight homoeologous exchange as a crucial mechanism for generating variations that are selected for crop improvement in polyploid species.

Genome wide detection of CNV and their association with body size in Danzhou chickens

Copy number variation (CNV) is a crucial component of genetic diversity in the genome, serving as the foundation for the genetic architecture and phenotypic variability of complex traits. In this study, we examined CNVs in the Danzhou (DZ) chicken, an indigenous breed exclusive to Hainan Province, China. By employing whole-genome resequencing data from 200 DZ chickens, we conducted a comprehensive genome-wide analysis of CNVs using CNVpytor and performed CNV-based genome-wide association studies (GWAS) on six body size traits, including body slope length (BSL), keel length (KeL), tibial length (TiL), tibial circumference (TiC), chest width (ChW), and chest depth (ChD) utilizing linear mixed model methods considering a genomic relationship matrix. We identified a total of 144,265 autosomal CNVs among the 200 individuals, comprising 67,818 deletions and 76,447 duplications. After merging these variants together, we obtained 4,824 distinct copy number variant regions, which accounted for approximately 20% of the chicken autosomal genome. Furthermore, we discovered several significantly associated CNV segments with body size traits located proximal to genes such as IHH, WNT6, WNT10A, LPR4, FZD2, WNT7B, and GNAS that have been extensively implicated in skeletal development and growth processes. These findings enhance our understanding of CNVs in chickens and their potential impact on body size traits by revealing candidate genes involved in the regulation of these traits. This establishes a solid framework for future studies and may prove particularly beneficial for exploring genetic structural variation in chickens.

A Whole-Genome Scan Revealed Genomic Features and Selection Footprints of Mengshan Cattle.

(1) Background: Mengshan cattle from the Yimeng mountainous region in China stand out as a unique genetic resource, known for their adaptive traits and environmental resilience. However, these cattle are currently endangered and comprehensive genomic characterization remains largely unexplored. This study aims to address this gap by investigating the genomic features and selection signals in Mengshan cattle. (2) Methods: Utilizing whole-genome resequencing data from 122 cattle, including 37 newly sequenced Mengshan cattle, we investigated population structure, genetic diversity, and selection signals. (3) Results: Our analyses revealed that current Mengshan cattle primarily exhibit European taurine cattle ancestry, with distinct genetic characteristics indicative of adaptive traits. We identified candidate genes associated with immune response, growth traits, meat quality, and neurodevelopment, shedding light on the genomic features underlying the unique attributes of Mengshan cattle. Enrichment analysis highlighted pathways related to insulin secretion, calcium signaling, and dopamine synapse, further elucidating the genetic basis of their phenotypic traits. (4) Conclusions: Our results provide valuable insights for further research and conservation efforts aimed at preserving this endangered genetic resource. This study enhances the understanding of population genetics and underscores the importance of genomic research in informing genetic resources and conservation initiatives for indigenous cattle breeds.

Recurrent selection shapes the genomic landscape of differentiation between a pair of host-specialized haplodiploids that diverged with gene flow.

Understanding the genetics of adaptation and speciation is critical for a complete picture of how biodiversity is generated and maintained. Heterogeneous genomic differentiation between diverging taxa is commonly documented, with genomic regions of high differentiation interpreted as resulting from differential gene flow, linked selection and reduced recombination rates. Disentangling the roles of each of these non-exclusive processes in shaping genome-wide patterns of divergence is challenging but will enhance our knowledge of the repeatability of genomic landscapes across taxa. Here, we combine whole-genome resequencing and genome feature data to investigate the processes shaping the genomic landscape of differentiation for a sister-species pair of haplodiploid pine sawflies, Neodiprion lecontei and Neodiprion pinetum. We find genome-wide correlations between genome features and summary statistics are consistent with pervasive linked selection, with patterns of diversity and divergence more consistently predicted by exon density and recombination rate than the neutral mutation rate (approximated by dS). We also find that both global and local patterns of FST, dXY and π provide strong support for recurrent selection as the primary selective process shaping variation across pine sawfly genomes, with some contribution from balancing selection and lineage-specific linked selection. Because inheritance patterns for haplodiploid genomes are analogous to those of sex chromosomes, we hypothesize that haplodiploids may be especially prone to recurrent selection, even if gene flow occurred throughout divergence. Overall, our study helps fill an important taxonomic gap in the genomic landscape literature and contributes to our understanding of the processes that shape genome-wide patterns of genetic variation.

Beyond gene flow: (non)-parallelism of secondary contact in a pair of highly differentiated sibling species.

Replicated secondary contact zones can provide insights into the barriers to gene flow that are important during speciation and can reveal to which degree secondary contact may result in similar evolutionary outcomes. Here, we studied two secondary contact zones between highly differentiated Alpine butterflies of the genus Erebia using whole-genome resequencing data. We assessed the genomic relationships between populations and species and found hybridization to be rare, with no to little current or historical introgression in either contact zone. There are large similarities between contact zones, consistent with an allopatric origin of interspecific differentiation, with no indications for ongoing reinforcing selection. Consistent with expected reduced effective population size, we further find that scaffolds related to the Z-chromosome show increased differentiation compared to the already high levels across the entire genome, which could also hint towards a contribution of the Z chromosome to species divergence in this system. Finally, we detected the presence of the endosymbiont Wolbachia, which can cause reproductive isolation between its hosts, in all E. cassioides, while it appears to be fully or largely absent in contact zone populations of E. tyndarus. We discuss how this rare pattern may have arisen and how it may have affected the dynamics of speciation upon secondary contact.