Using bioinformatics to investigate functional diversity: a case study of MHC diversity in koalas.

Conservation genomics can greatly improve conservation outcomes of threatened populations, including those impacted by disease. Understanding diversity within immune gene families, including the major histocompatibility complex (MHC) and toll-like receptors (TLR), is important due to the role they play in disease resilience and susceptibility. With recent advancements in sequencing technologies and bioinformatic tools, the cost of generating high-quality sequence data has significantly decreased and made it possible to investigate diversity across entire gene families in large numbers of individuals compared to investigating only a few genes or a few populations previously. Here, we use the koala as a case study for investigating functional diversity across populations. We utilised previous target enrichment data and 438 whole genomes to firstly, determine the level of sequencing depth required to investigate MHC diversity and, secondly, determine the current level of diversity in MHC genes in koala populations. We determined for low complexity, conserved genes such as TLR genes 10 × sequencing depth is sufficient to reliably genotype more than 90% of variants, whereas for complex genes such as the MHC greater than 20 × and preferably 30 × sequencing depth is required. We used whole genome data to identify 270 biallelic SNPs across 24 MHC genes as well as copy number variation (CNV) within class I and class II genes and conduct supertype analysis. Overall, we have provided a bioinformatic workflow for investigating variation in a complex immune gene family from whole genome sequencing data and determined current levels of diversity within koala MHC genes.

Introgression drives adaptation to the plateau environment in a subterranean rodent.

BACKGROUND: Introgression has repeatedly been shown to play an important role in the adaptation of species to extreme environments, yet how introgression enables rodents with specialized subterranean lifestyle to acclimatize to high altitudes is still unclear. Myospalacinae is a group of subterranean rodents, among which the high-altitude plateau zokors (Eospalax baileyi) and the low-altitude Gansu zokors (E. cansus) are sympatrically distributed in the grassland ecosystems of the Qinghai-Tibet Plateau (QTP). Together, they provide a model for the study of the role of introgression in the adaptation of low-altitude subterranean rodents to high altitudes. RESULTS: Applying low-coverage whole-genome resequencing and population genetics analyses, we identified evidence of adaptive introgression from plateau zokors into Gansu zokors, which likely facilitated the adaptation of the latter to the high-altitude environment of the QTP. We identified positively selected genes with functions related to energy metabolism, cardiovascular system development, calcium ion transport, and response to hypoxia which likely made critical contributions to adaptation to the plateau environment in both plateau zokors and high-altitude populations of Gansu zokors. CONCLUSIONS: Introgression of genes associated with hypoxia adaptation from plateau zokors may have played a role in the adaptation of Gansu zokors to the plateau environment. Our study provides new insights into the understanding of adaptive evolution of species on the QTP and the importance of introgression in the adaptation of species to high-altitude environments.

Whole genome re-sequencing reveals high altitude adaptation signatures and admixture in Ladakhi cattle.

Ladakhi cattle, native to the high-altitude region of Ladakh in northern India (ranging from 3,000 to 5,000 m above sea level), have evolved unique genetic adaptations to thrive in harsh environmental conditions, such as hypoxia, extreme cold, and low humidity. This study explored the genome of Ladakhi cattle to investigate genetic structure, selection signatures, and adaptive mechanisms. Whole genome sequencing reads, generated on Illumina NovaSeq 6000 platform, were aligned to the Bos taurus reference genome with BWA-MEM. SNPs were identified and filtered using GATK and bcftools, and functionally annotated with SnpEff. For population genomic analysis, PCA and admixture modeling assessed genetic structure, while Neighbor-Joining trees, LD decay, nucleotide diversity (π), and FST evaluated phylogenetic relationships and genetic variation. Selective sweeps were detected using RAiSD, and gene-set enrichment and protein-protein interaction analyses were conducted to explore functional pathways related to adaptation. The study revealed 3,759,279 unique SNPs and demonstrated that Ladakhi cattle form a distinct genetic cluster with an estimated admixture of 68 % Bos indicus and 32 % Bos taurus ancestry. Key findings include rapid linkage disequilibrium decay, low inbreeding level, and the identification of selection signatures and genes associated with hypoxia response, energy metabolism, and cold adaptation. Mean nucleotide diversity (π, 0.0037) and FST values indicated moderate genetic differentiation from other breeds. The analysis highlighted selection signatures for genes like HIF1A, ENO4, ANGPT1, EPO, NOS3, MAPK3, HMOX1, BCL2,CAMK2D, MTOR, AKT2,PIK3CB, and MAP2K1, among others, including various keratin and heat shock proteins. The interaction between genes associated with hypoxia signaling (HIF-1) and other enriched pathways such as PI3K, mTOR, NFκB, ERK, and ER stress, reveals a complex mechanism for managing hypoxic stress in Ladakhi cattle. These findings offer valuable insights for breeding programs aimed at enhancing livestock resilience in extreme environments and enhance understanding of mammalian adaptation to high-altitude conditions.

An ultra-dense linkage map identified quantitative trait loci corresponding to fruit quality- and size-related traits in red goji berry.

Goji berries are a small-fruited shrub with industrial importance whose fruit considered beneficial in both fresh and dried forms. Current germplasms of goji berries include small fruits with a short shelf life, less sweet and bitter taste, and a lack of appropriate genetic information. This study aimed to employ whole genome resequencing to generate an ultra-dense bin linkage map and to elucidate the genetic basis of goji fruit quality and size using quantitative trait loci (QTL) mapping analysis in a cross-pollinated hybrid population. To achieve this goal, human sensory tests were carried out to determine the bitter taste (BT) and sweet taste (ST), and to quantify the soluble solid content (SSC), fruit firmness (FF), and fruit size-related traits of fresh goji fruits over three or four years. The results revealed that the goji bin linkage map based on resequencing spanned a total length of 966.42 cM and an average bin interval of 0.03 cM. Subsequent variant calling and ordering resulted in 3,058 bins containing 35,331 polymorphic markers across 12 chromosomes. A total of 99 QTLs, with individual loci in different environments explaining a phenotypic variance of 1.21-16.95% were identified for the studied traits. Ten major effects, including colocalized QTLs corresponding to different traits, were identified on chromosomes 1, 3, 5, 6, 7, and 8, with a maximum Logarithm of Odds (LOD) of 29.25 and 16.95% of explained phenotypic variance (PVE). In addition, four stable loci, one for FF, one for fruit weight (FW), and two for fruit shape index (FSI), were mainly mapped on chromosomes 5, 6, and 7, elucidating 2.10-16.95% PVE. These findings offer valuable insights into the genetic architecture of goji fruit traits along with identified specific loci and markers to further improve and develop sweeter, less bitter and larger fruited goji berry cultivars with extended shelf life.

Whole-Genome Resequencing Identifies SNPs in Sucrose Synthase and Sugar Transporter Genes Associated with Sweetness in Coconut.

Coconut (Cocos nucifera L.) is an important agricultural commodity with substantial economic and nutritional value, widely used for various products, including coconut water. The sweetness is an important quality trait of coconut water, which is influenced by genetic and environmental factors. In this study, we utilized next-generation sequencing to identify genetic variations in the coconut genome associated with the sweetness of coconut water. Whole-genome resequencing of 49 coconut accessions, including diverse germplasm and an F2 population of 81 individuals, revealed ~27 M SNPs and ~1.5 M InDels. Sugar content measured by °Bx was highly variable across all accessions tested, with dwarf varieties generally sweeter. A comprehensive analysis of the sugar profiles revealed that sucrose was the major sugar contributing to sweetness. Allele mining of the 148 genes involved in sugar metabolism and transport and genotype-phenotype association tests revealed two significant SNPs in the hexose carrier protein (Cnu01G018720) and sucrose synthase (Cnu09G011120) genes associated with the higher sugar content in both the germplasm and F2 populations. This research provides valuable insights into the genetic basis of coconut sweetness and offers molecular markers for breeding programs aimed at improving coconut water quality. The identified variants can improve the selection process in breeding high-quality sweet coconut varieties and thus support the economic sustainability of coconut cultivation.

Whole-genome resequencing improves the utility of otoliths as a critical source of DNA for fish stock research and monitoring.

Fish ear bones, known as otoliths, are often collected in fisheries to assist in management, and are a common sample type in museum and national archives. Beyond their utility for ageing, morphological and trace element analysis, otoliths are a repository of valuable genomic information. Previous work has shown that DNA can be extracted from the trace quantities of tissue remaining on the surface of otoliths, despite the fact that they are often stored dry at room temperature. However, much of this work has used reduced representation sequencing methods in clean lab conditions, to achieve adequate yields of DNA, libraries and ultimately single-nucleotide polymorphisms (SNPs). Here, we pioneer the use of small-scale (spike-in) sequencing to screen contemporary otolith samples prepared in regular molecular biology (in contrast to clean) laboratories for contamination and quality levels, submitting for whole-genome resequencing only samples above a defined endogenous DNA threshold. Despite the typically low quality and quantity of DNA extracted from otoliths, we are able to produce whole-genome libraries and ultimately sets of filtered, unlinked and even putatively adaptive SNPs of ample numbers for downstream uses in population, climate and conservation genomics. By comparing with a set of tissue samples from the same species, we are able to highlight the quality and efficacy of otolith samples from DNA extraction and library preparation, to bioinformatic preprocessing and SNP calling. We provide detailed schematics, protocols and scripts of our approach, such that it can be adopted widely by the community, improving the use of otoliths as a source of valuable genomic data.

Microevolutionary change in wild stickleback: Using integrative time-series data to infer responses to selection.

A central goal in evolutionary biology is to understand how different evolutionary processes cause trait change in wild populations. However, quantifying evolutionary change in the wild requires linking trait change to shifts in allele frequencies at causal loci. Nevertheless, datasets that allow for such tests are extremely rare and existing theoretical approaches poorly account for the evolutionary dynamics that likely occur in ecological settings. Using a decade-long integrative phenome-to-genome time-series dataset on wild threespine stickleback (Gasterosteus aculeatus), we identified how different modes of selection (directional, episodic, and balancing) drive microevolutionary change in correlated traits over time. Most strikingly, we show that feeding traits changed by as much 25% across 10 generations which was driven by changes in the genetic architecture (i.e., in both genomic breeding values and allele frequencies at genetic loci for feeding traits). Importantly, allele frequencies at genetic loci related to feeding traits changed at a rate greater than expected under drift, suggesting that the observed change was a result of directional selection. Allele frequency dynamics of loci related to swimming traits appeared to be under fluctuating selection evident in periodic population crashes in this system. Our results show that microevolutionary change in a wild population is characterized by different modes of selection acting simultaneously on different traits, which likely has important consequences for the evolution of correlated traits. Our study provides one of the most thorough descriptions to date of how microevolutionary processes result in trait change in a natural population.

Whole-genome resequencing and RNA-seq analysis implicates GPR75 as a potential genetic basis related to retarded growth in South China carp (Cyprinus carpio rubrofuscus).

The south China carp (Cyprinus carpio rubrofuscus) is an indigenous and important fish species, widely cultured in south China. However, part of individuals experienced retarded growth, the genetic basis of which has yet to be elucidated. In this study, whole-genome resequencing of 35 fast-growing and 35 retarded-growing south China carp were conducted to identify promising genes associated with retarded growth. Twelve candidate SNPs were detected and annotated to the Gpr75 gene, which has been reported to be related with body weight through regulating insulin homeostasis. RNA-seq analysis of muscle suggested that differentially expressed genes were significantly enriched in the insulin signaling pathway. Additionally, the fasting serum insulin level was significantly lower while the blood glucose level was significantly higher in the retarded-growing group. Our preliminary study provides insights into the genetic basis underlying the retarded growth and may facilitate further genetic improvement of south China carp.

Genetic structure analysis of yak breeds and their response to adaptive evolution.

Yaks are crucial genetic resources in the Tibetan Plateau and surrounding regions. Throughout the long process of domestication, natural and artificial selection pressures have enabled yaks to demonstrate adaptive characteristics to the environment in terms of physiological structure and genetic molecules, but no systematic cell analysis has been carried out on this phenomenon of yaks. Here, the population structure and genetic diversity of yak were studied by WGRS, and the genes related to yak adaptability were excavated. Combined with scRNA-seq method, the transcription map of yak lung tissue and skin tissue was constructed, which provided a new comprehensive insight into yak adaptability. The analysis of yak population structure showed that there was obvious genetic differentiation between TZ _ yak and other seven yak populations, while there was significant genetic exchange between PL _ yak and SB _ yak at high altitude. WGRS and scRNA-seq analysis revealed that the gene HIF1A related to high altitude adaptation was expressed in various cell types, while EPAS1 was predominantly expressed in epithelial and endothelial cells of yak lung tissue. Endothelial cells play a critical role in hypoxia-adapted VEGF signaling, which correlates closely with the high expression of KDR and VEGFA genes in endothelial cells and monocytes. Furthermore, in the selection signal of High _ yak vs Low _ yak, 19.8 % of the genes overlapped with the genes screened by skin scRNA-seq, including genes related to coat color such as RORA, BNC2, and KIT. Notably, BNC2 is a gene associated with melanin deposition and shows high expression levels in HS cells. Additionally, GRN in melanocytes and SORT1 in IRS play an important role in cell communication between melanocytes and IRS. These findings offer new insights into the natural polymorphism of yaks and provide a valuable reference for future research on high-altitude mammals, and potentially even human genetics.

Whole-Genome Resequencing Reveals Signatures of Adaptive Evolution in Acanthopagrus latus and Rhabdosargus sarba.

Acanthopagrus latus and Rhabdosargus sarba are economically important marine species along the coast of China, with similar external morphological characteristics and living habits, with wide distribution and strong adaptability. To investigate the molecular mechanisms underlying the adaptive evolution of these two species, we conducted whole-genome resequencing of 10 individuals of both species from the coastal waters of Wuyu Island, Fujian, China, using high-throughput sequencing technology. We obtained SNP, InDel, CNV, and SV variation information and annotated these variations, constructing a genomic variation database for both species. By comparing the resequencing data with reference genomes, we identified 9,829,511 SNP loci in the population of A. latus and 34,051,056 SNP loci in the population of R. sarba. Using whole-genome SNP data, we employed Fst and ROD methods to identify candidate genomic regions under selection. Functional annotation and enrichment analysis using GO and KEGG databases revealed potential adaptive evolution in R. sarba associated with immune response, feeding, growth and development, and locomotion, while A. latus showed potential adaptive evolution associated with immune response, nervous system, growth and development, and metabolism.

Furfural tolerance of mutant Saccharomyces cerevisiae selected via ionizing radiation combined with adaptive laboratory evolution.

BACKGROUND: Lignocellulose is a renewable and sustainable resource used to produce second-generation biofuel ethanol to cope with the resource and energy crisis. Furfural is the most toxic inhibitor of Saccharomyces cerevisiae cells produced during lignocellulose treatment, and can reduce the ability of S. cerevisiae to utilize lignocellulose, resulting in low bioethanol yield. In this study, multiple rounds of progressive ionizing radiation was combined with adaptive laboratory evolution to improve the furfural tolerance of S. cerevisiae and increase the yield of ethanol. RESULTS: In this study, the strategy of multiple rounds of progressive X-ray radiation combined with adaptive laboratory evolution significantly improved the furfural tolerance of brewing yeast. After four rounds of experiments, four mutant strains resistant to high concentrations of furfural were obtained (SCF-R1, SCF-R2, SCF-R3, and SCF-R4), with furfural tolerance concentrations of 4.0, 4.2, 4.4, and 4.5 g/L, respectively. Among them, the mutant strain SCF-R4 obtained in the fourth round of radiation had a cellular malondialdehyde content of 49.11 nmol/mg after 3 h of furfural stress, a weakening trend in mitochondrial membrane potential collapse, a decrease in accumulated reactive oxygen species, and a cell death rate of 12.60%, showing better cell membrane integrity, stable mitochondrial function, and an improved ability to limit reactive oxygen species production compared to the other mutant strains and the wild-type strain. In a fermentation medium containing 3.5 g/L furfural, the growth lag phase of the SCF-R4 mutant strain was shortened, and its growth ability significantly improved. After 96 h of fermentation, the ethanol production of the mutant strain SCF-R4 was 1.86 times that of the wild-type, indicating that with an increase in the number of irradiation rounds, the furfural tolerance of the mutant strain SCF-R4 was effectively enhanced. In addition, through genome-transcriptome analysis, potential sites related to furfural detoxification were identified, including GAL7, MAE1, PDC6, HXT1, AUS1, and TPK3. CONCLUSIONS: These results indicate that multiple rounds of progressive X-ray radiation combined with adaptive laboratory evolution is an effective mutagenic strategy for obtaining furfural-tolerant mutants and that it has the potential to tap genes related to the furfural detoxification mechanism.

Development and validation of sex-linked molecular markers for rapid and accurate identification of male and female Hippophae tibetana plants.

Hippophae tibetana, one of the highest-altitude woody plants endemic to the Qinghai-Tibet Plateau, primarily thrives on riverbanks formed by glacial meltwater. As a dioecious species, it demonstrates significant ecological and economic value in extreme alpine environments. However, the lack of sex identification techniques outside of the flowering period severely limits research on sex ratio, differentiation, and breeding. There is an urgent need to develop effective sex-linked molecular markers that are independent of developmental stages, but current research in this area remains limited. This study developed a set of accurate sex-linked molecular markers for the rapid identification of male and female individuals of H. tibetana. Through whole-genome resequencing of 32 sexually differentiated H. tibetana samples, this study offers strong evidence supporting chromosome 2 as the sex chromosome and successfully identified key loci related to sex determination on this chromosome. Utilizing these loci, we, for the first time, developed three reliable pairs of sex-specific molecular markers, which exhibited high accuracy during validation across various geographic populations, offering an effective tool for the sex identification of H. tibetana. Additionally, this study lays the groundwork for further research into the mechanisms of sex determination and the evolution of sex chromosomes in H. tibetana.

Improvement of Saccharomyces cerevisiae strain tolerance to vanillin through heavy ion radiation combined with adaptive laboratory evolution.

Vanillin is an inhibitor of lignocellulose hydrolysate, which can reduce the ability of Saccharomyces cerevisiae to utilize lignocellulose, which is an important factor limiting the development of the ethanol fermentation industry. In this study, mutants of vanillin-tolerant yeast named H6, H7, X3, and X8 were bred by heavy ion irradiation (HIR) combined with adaptive laboratory evolution (ALE). Phenotypic tests revealed that the mutants outperformed the original strain WT in tolerance, growth rate, genetic stability and fermentation ability. At 1.6 g/L vanillin concentration, the average OD600 value obtained for mutant strains was 0.95 and thus about 3.4-fold higher than for the wild-type. When the concentration of vanillin was 2.0 g/L, the glucose utilization rate of the mutant was 86.3 % within 96 h, while that of the original strain was only 70.0 %. At this concentration of vanillin, the mitochondrial membrane potential of the mutant strain recovered faster than that of the original strain, and the ROS scavenging ability was stronger. We analyzed the whole transcriptome sequencing map and the whole genome resequencing of the mutant, and found that DEGs such as FLO9, GRC3, PSP2 and SWF1, which have large differential expression multiples and obvious mutation characteristics, play an important role in cell flocculation, rDNA transcription, inhibition of DNA polymerase mutation and protein palmitoylation. These functions can help cells resist vanillin stress. The results show that combining HIR with ALE is an effective mutagenesis strategy. This approach can efficiently obtain Saccharomyces cerevisiae mutants with improved vanillin tolerance, and provide reference for obtaining robust yeast strains with lignocellulose inhibitor tolerance.

Genome-wide association study reveals 2 copy number variations associated with the variation of plumage color in the white duck hybrid population.

Plumage color is an intuitive external poultry characteristic with rich manifestations and complex genetic mechanisms. In our previous study, we observed that there were more dark variations in plumage color in the F2 population derived from the hybridization of 2 white duck varieties. Therefore, based on the statistics of plumage color of 308 F2 populations, we further used the resequencing data of these individuals to detect copy number variations (CNVs) in the whole genome and conducted genome-wide association studies (GWAS) to determine the genetic basis related to plumage color traits. The CNV detection revealed 9,337 CNVs, with an average length of 15,950 bp and a total length of 142.02 MB, accounting for approximately 12.91% of the reference genome. The CNV distribution on the chromosomes was relatively uniform, and the number of CNVs on each chromosome positively correlated with the length of the chromosome. In the pure black plumage group, 2,101 CNVs were only identified, and 1,714 were specifically identified in the pure white plumage group. Ten CNVs were randomly selected for validation using quantitative real-time PCR, and 9 CNVs had the same CNV types as predicted, with an accuracy of 90%. Based on GWAS, we identified 2 CNVs potentially associated with plumage color variations, with the associated CNV regions covering 9 genes. Enrichment analysis of these 9 candidate genes showed significant enrichment of 3 pathways (ribosome biogenesis in eukaryotes, RNA transport, and protein export) and 17 gene ontology terms. Among these, VWA5A can downregulate MITF by binding to the regulatory factors SOX10. The occurrence of CNV may indirectly contribute to duck plumage color variation by affecting the regulatory factors of the switch gene MITF in the melanogenesis pathway. These findings have improved the understanding of the genetic basis of duck plumage color variation and have been beneficial for developing and using plumage color traits in subsequent poultry breeding.

Drug resistance and genetic characteristics of one Eimeria tenella isolate from Xiantao, Hubei Province, China.

Chicken coccidiosis causes retarded growth and low production performance in poultry, resulting in huge economic losses to the poultry industry. In order to prevent and control chicken coccidiosis, great efforts have been made to develop new drugs and vaccines, which require pure isolates of Eimeria spp. In this study, we obtained the Eimeira tenella Xiantao isolate by single oocyst isolation technology and compared its genome with the reference genome GCF_000499545.2_ETH001 of the Houghton strain. The results of the comparative genomic analysis indicated that the genome of this isolate contained 46,888 single nucleotide polymorphisms (SNPs). There were 15,107 small insertion and deletion variations (indels), 1693 structural variations (SV), and 3578 copy number variations (CNV). In addition, 64 broilers were used to determine the resistance profile of Xiantao strain. Drug susceptibility testing revealed that this isolate was completely resistant to monensin, diclazuril, halofuginone, sulfachlorpyrazine sodium, and toltrazuril, but sensitive to decoquinate. These data improve our understanding of drug resistance in avian coccidia.

Comparison of Mutations Induced by Different Doses of Fast-Neutron Irradiation in the M1 Generation of Sorghum (Sorghum bicolor).

Sorghum is an important C4 crop with various food and nonfood uses. Although improvements through hybridization and selection have been exploited, the introduction of genetic variation and the development of new genotypes in sorghum are still limited. Fast-neutron (FN) mutagenesis is a very effective method for gene functional studies and to create genetic variability. However, the full spectrum of FN-induced mutations in sorghum is poorly understood. To address this, we generated an FN-induced mutant population from the inbred line 'BTx623' and sequenced 40 M1 seedlings to evaluate the mutagenic effects of FNs on sorghum. The results show that each line had an average of 43.7 single-base substitutions (SBSs), 3.7 InDels and 35.15 structural variations (SVs). SBSs accounted for approximately 90.0% of the total number of small mutations. Among the eight treatment groups, FN irradiation at a dose of 19 Gy generated the highest number of mutations. The ratio of transition/transversion ranged from 1.77 to 2.21, and the G/C to A/T transition was the most common substitution in all mutant lines. The distributions of the identified SBSs and InDels were similar and uneven across the genome. An average of 3.63 genes were mutated in each mutant line, indicating that FN irradiation resulted in a suitable density of mutated genes, which can be advantageous for improving elite material for one specific or a few traits. These results provide a basis for the selection of the suitable dose of mutagen and new genetic resources for sorghum breeding.

Loss of genetic variation and ancestral sex determination system in North American northern pike characterized by whole-genome resequencing.

The northern pike Esox lucius is a freshwater fish with low genetic diversity but ecological success throughout the Northern Hemisphere. Here, we generate an annotated chromosome-level genome assembly of 941 Mbp in length with 25 chromosome-length scaffolds. We then genotype 47 northern pike from Alaska through New Jersey at a genome-wide scale and characterize a striking decrease in genetic diversity along the sampling range. Individuals west of the North American Continental Divide have substantially higher diversity than those to the east (e.g. Interior Alaska and St. Lawrence River have on average 181 and 64K heterozygous SNPs per individual, or a heterozygous SNP every 5.2 and 14.6 kbp, respectively). Individuals clustered within each population with strong support, with numerous private alleles observed within each population. Evidence for recent population expansion was observed for a Manitoba hatchery and the St. Lawrence population (Tajima's D = -1.07 and -1.30, respectively). Several chromosomes have large regions with elevated diversity, including LG24, which holds amhby, the ancestral sex determining gene. As expected amhby was largely male-specific in Alaska and the Yukon and absent southeast to these populations, but we document some amhby(-) males in Alaska and amhby(+) males in the Columbia River, providing evidence for a patchwork of presence of this system in the western region. These results support the theory that northern pike recolonized North America from refugia in Alaska and expanded following deglaciation from west to east, with probable founder effects resulting in loss of both neutral and functional diversity (e.g. amhby).

Whole-Genome Resequencing Reveals Genetic Diversity and Growth Trait-Related Genes in Pinan Cattle.

The breeding of high-quality beef cattle breeds is crucial for the development of animal husbandry, and whole-genome resequencing is widely applicated in the field of molecular breeding. Advantages in growth and reproductive traits exist in Pinan cattle compared with other cattle breeds, but there is limited research on its genomic mechanism. Using whole-genome resequencing, the genetic structure and genomic selection signatures in Pinan cattle were investigated in this study. Phylogenetic, cluster, and admixture analysis results indicated that Pinan cattle have a closer genetic relationship with Kholmogory cattle and China north cattle breeds. Through a selective sweep strategy, 207 and 54 candidate genes related to growth and reproduction and immunity, respectively, were identified in the Pinan cattle population. Given the crucial role of the glutamate-cysteine ligase catalytic (GCLC) gene in muscle antioxidative defense, the high frequency of allele T of the GCLC c.429 C>T locus in the Pinan cattle population might partially contribute to the advantages of Pinan cattle in growth performance. This study laid the foundation for the genetic improvement in Chinese local beef cattle and provide background for the studies on the growth and development of Pinan cattle.

Discordance between taxonomy and population genomic data: An avian example relevant to the United States Endangered Species Act.

Population genomics can reveal cryptic biological diversity that may impact fitness while simultaneously serving to delineate relevant conservation units. Here, we leverage the power of whole-genome resequencing for conservation by studying 433 individual lesser prairie-chicken (Tympanuchus pallidicinctus; LEPC, a federally endangered species of conservation concern in the United States) and greater prairie-chicken (Tympanuchus cupido; GRPC, a legally huntable species throughout much of its range). The genomic diversity of two formally recognized distinct population segments (DPSs) of LEPCs is similar, but they are genetically distinct. Neither DPS is depleted of its genomic diversity, neither is especially inbred, and temporal diversity is relatively stable in both conservation units. Interspecific differentiation between the two species was only slightly higher than that observed between LEPC DPSs, due largely to bidirectional introgression. The high resolution provided by our dataset identified a genomic continuum between the two species such that individuals sampled from the hybrid zone were imperfectly assigned to their presumptive species when considering only their physical characteristics. The admixture between the two species is reflected in the spectrum of individual ancestry coefficients, which has legal implications for the "take" of individuals under the Endangered Species Act. Overall, our data highlight the recurring dissonance between static policies and dynamic species boundaries that are increasingly obvious in the population genomic era.

The genetic characterization of germplasm and identification of the litter size trait associated candidate genes in Dexin mutton and fine-wool sheep.

Xinjiang is a major province of sheep breeding in China, which plays an important role in meeting people's needs for meat products, increasing farmers' income and sustainable development of animal husbandry. However, the genetic differentiation relationship between breeds was not clear, and most sheep had low fecundity, which seriously restricted the efficient development of sheep industry. Therefore, this study used the whole genome resequencing to detect the genetic variation of Dexin mutton and fine-wool sheep, explored the selected regions and important genes of the litter size traits, analyzed the genetic mechanism of reproductive traits, and provided new insights for the high fecundity breeding of sheep. A total of 5,236.338 G genome data and 35,884,037 SNPs were obtained. Furthermore, we identified 39 selection signals spanning candidate genes, 99 genes were significantly associated related to growth, reproduction and immunity, among which, BRIP1, BMPR1B, BMP4, NGF, etc. genes, and MAKP signaling pathway, Fanconi anemia pathway and Thyroid hormone signaling pathway and other signaling pathways were significantly correlated with litter size trait. Among them, we identified NGF, TrKA and BRIP1 genes was the important genes for sheep litter size traits and the mutation frequencies of 9 SNPs in BRIP1 gene were significantly different in domestic sheep in the world. The research provided new insights for the breeding of self-cultivated meat fine-wool sheep.