A Chinese indicine pangenome reveals a wealth of novel structural variants introgressed from other Bos species.

Chinese indicine cattle harbor a much higher genetic diversity compared with other domestic cattle, but their genome architecture remains uninvestigated. Using PacBio HiFi sequencing data from 10 Chinese indicine cattle across southern China, we assembled 20 high-quality partially phased genomes and integrated them into a multiassembly graph containing 148.5 Mb (5.6%) of novel sequence. We identified 156,009 high-confidence nonredundant structural variants (SVs) and 206 SV hotspots spanning ∼195 Mb of gene-rich sequence. We detected 34,249 archaic introgressed fragments in Chinese indicine cattle covering 1.93 Gb (73.3%) of the genome. We inferred an average of 3.8%, 3.2%, 1.4%, and 0.5% of introgressed sequence originating, respectively, from banteng-like, kouprey-like, gayal-like, and gaur-like Bos species, as well as 0.6% of unknown origin. Introgression from multiple donors might have contributed to the genetic diversity of Chinese indicine cattle. Altogether, this study highlights the contribution of interspecies introgression to the genomic architecture of an important livestock population and shows how exotic genomic elements can contribute to the genetic variation available for selection.

A sheep pangenome reveals the spectrum of structural variations and their effects on tail phenotypes.

Structural variations (SVs) are a major contributor to genetic diversity and phenotypic variations, but their prevalence and functions in domestic animals are largely unexplored. Here we generated high-quality genome assemblies for 15 individuals from genetically diverse sheep breeds using Pacific Biosciences (PacBio) high-fidelity sequencing, discovering 130.3 Mb nonreference sequences, from which 588 genes were annotated. A total of 149,158 biallelic insertions/deletions, 6531 divergent alleles, and 14,707 multiallelic variations with precise breakpoints were discovered. The SV spectrum is characterized by an excess of derived insertions compared to deletions (94,422 vs. 33,571), suggesting recent active LINE expansions in sheep. Nearly half of the SVs display low to moderate linkage disequilibrium with surrounding single-nucleotide polymorphisms (SNPs) and most SVs cannot be tagged by SNP probes from the widely used ovine 50K SNP chip. We identified 865 population-stratified SVs including 122 SVs possibly derived in the domestication process among 690 individuals from sheep breeds worldwide. A novel 168-bp insertion in the 5' untranslated region (5' UTR) of HOXB13 is found at high frequency in long-tailed sheep. Further genome-wide association study and gene expression analyses suggest that this mutation is causative for the long-tail trait. In summary, we have developed a panel of high-quality de novo assemblies and present a catalog of structural variations in sheep. Our data capture abundant candidate functional variations that were previously unexplored and provide a fundamental resource for understanding trait biology in sheep.

Long divergent haplotypes introgressed from wild sheep are associated with distinct morphological and adaptive characteristics in domestic sheep.

The worldwide sheep population comprises more than 1000 breeds. Together, these exhibit a considerable morphological diversity, which has not been extensively investigated at the molecular level. Here, we analyze whole-genome sequencing individuals of 1,098 domestic sheep from 154 breeds, and 69 wild sheep from seven Ovis species. On average, we detected 6.8%, 1.0% and 0.2% introgressed sequence in domestic sheep originating from Iranian mouflon, urial and argali, respectively, with rare introgressions from other wild species. Interestingly, several introgressed haplotypes contributed to the morphological differentiations across sheep breeds, such as a RXFP2 haplotype from Iranian mouflon conferring the spiral horn trait, a MSRB3 haplotype from argali strongly associated with ear morphology, and a VPS13B haplotype probably originating from urial and mouflon possibly associated with facial traits. Our results reveal that introgression events from wild Ovis species contributed to the high rate of morphological differentiation in sheep breeds, but also to individual variation within breeds. We propose that long divergent haplotypes are a ubiquitous source of phenotypic variation that allows adaptation to a variable environment, and that these remain intact in the receiving population probably due to reduced recombination.

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

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

Markhor-derived Introgression of a Genomic Region Encompassing PAPSS2 Confers High-altitude Adaptability in Tibetan Goats.

Understanding the genetic mechanism of how animals adapt to extreme conditions is fundamental to determine the relationship between molecular evolution and changing environments. Goat is one of the first domesticated species and has evolved rapidly to adapt to diverse environments, including harsh high-altitude conditions with low temperature and poor oxygen supply but strong ultraviolet radiation. Here, we analyzed 331 genomes of domestic goats and wild caprid species living at varying altitudes (high > 3000 m above sea level and low < 1200 m), along with a reference-guided chromosome-scale assembly (contig-N50: 90.4 Mb) of a female Tibetan goat genome based on PacBio HiFi long reads, to dissect the genetic determinants underlying their adaptation to harsh conditions on the Qinghai-Tibetan Plateau (QTP). Population genomic analyses combined with genome-wide association studies (GWAS) revealed a genomic region harboring the 3'-phosphoadenosine 5'-phosphosulfate synthase 2 (PAPSS2) gene showing strong association with high-altitude adaptability (PGWAS = 3.62 × 10-25) in Tibetan goats. Transcriptomic data from 13 tissues revealed that PAPSS2 was implicated in hypoxia-related pathways in Tibetan goats. We further verified potential functional role of PAPSS2 in response to hypoxia in PAPSS2-deficient cells. Introgression analyses suggested that the PAPSS2 haplotype conferring the high-altitude adaptability in Tibetan goats originated from a recent hybridization between goats and a wild caprid species, the markhor (Capra falconeri). In conclusion, our results uncover a hitherto unknown contribution of PAPSS2 to high-altitude adaptability in Tibetan goats on QTP, following interspecific introgression and natural selection.

Genome-wide association study for wattles trait in the dairy goat breed.

Dairy goats are significant livestock that provide high-quality milk sources in the world. The wattles trait is an evident phenotypic character on the neck of a dairy goat, which is considered to be under genetic control. We collected samples of 189 dairy goats, including 94 with wattles and 95 without wattles, from four different farms and multiple dairy goat breeds. The samples were genotyped with the GeneSeek Genomic Profiler Goat 70 K SNP chip. Genome-wide association studies (GWAS) in wattles have identified associations with single nucleotide polymorphisms (SNPs) at chromosome 10. In this area, an extremely strong association locus was assigned to FMN1 (Formin 1) belongs to the formin homology family and is associated with limb deformity, other candidate genes of interest confirmed for wattles were ARHGAP11A (Rho GTPase Activating Protein 11 A) and GJD2 (Gap Junction Protein Delta 2). Meanwhile, we found the presence or absence of wattles had no significant effect on milk yield. This research will provide genetic resources useful to explore genetic factors affecting the trait.

Ancient Genomes Reveal the Evolutionary History and Origin of Cashmere-Producing Goats in China.

Goats are one of the most widespread farmed animals across the world; however, their migration route to East Asia and local evolutionary history remain poorly understood. Here, we sequenced 27 ancient Chinese goat genomes dating from the Late Neolithic period to the Iron Age. We found close genetic affinities between ancient and modern Chinese goats, demonstrating their genetic continuity. We found that Chinese goats originated from the eastern regions around the Fertile Crescent, and we estimated that the ancestors of Chinese goats diverged from this population in the Chalcolithic period. Modern Chinese goats were divided into a northern and a southern group, coinciding with the most prominent climatic division in China, and two genes related to hair follicle development, FGF5 and EDA2R, were highly divergent between these populations. We identified a likely causal de novo deletion near FGF5 in northern Chinese goats that increased to high frequency over time, whereas EDA2R harbored standing variation dating to the Neolithic. Our findings add to our understanding of the genetic composition and local evolutionary process of Chinese goats.

A near complete genome for goat genetic and genomic research.

BACKGROUND: Goat, one of the first domesticated livestock, is a worldwide important species both culturally and economically. The current goat reference genome, known as ARS1, is reported as the first nonhuman genome assembly using 69× PacBio sequencing. However, ARS1 suffers from incomplete X chromosome and highly fragmented Y chromosome scaffolds. RESULTS: Here, we present a very high-quality de novo genome assembly, Saanen_v1, from a male Saanen dairy goat, with the first goat Y chromosome scaffold based on 117× PacBio long-read sequencing and 118× Hi-C data. Saanen_v1 displays a high level of completeness thanks to the presence of centromeric and telomeric repeats at the proximal and distal ends of two-thirds of the autosomes, and a much reduced number of gaps (169 vs. 773). The completeness and accuracy of the Saanen_v1 genome assembly are also evidenced by more assembled sequences on the chromosomes (2.63 Gb for Saanen_v1 vs. 2.58 Gb for ARS1), a slightly increased mapping ratio for transcriptomic data, and more genes anchored to chromosomes. The eight putative large assembly errors (1 to ~ 7 Mb each) found in ARS1 were amended, and for the first time, the substitution rate of this ruminant Y chromosome was estimated. Furthermore, sequence improvement in Saanen_v1, compared with ARS1, enables us to assign the likely correct positions for 4.4% of the single nucleotide polymorphism (SNP) probes in the widely used GoatSNP50 chip. CONCLUSIONS: The updated goat genome assembly including both sex chromosomes (X and Y) and the autosomes with high-resolution quality will serve as a valuable resource for goat genetic research and applications.

Enantioselective Synthesis of Multisubstituted Allenes by Cooperative Cu/Pd-Catalyzed 1,4-Arylboration of 1,3-Enynes.

A cooperative Cu/Pd-catalyzed enantioselective synthesis of multisubstituted allenes is established. By employing chiral sulfoxide phosphine (SOP)/Cu and PdCl2 (dppf) complexes as catalysts, the 1,4-arylboration of 1,3-enynes provides an efficient approach to trisubstituted chiral allenes with up to 92 % yield and 97:3 er. Furthermore, by using 2-substituted 1,3-enynes as substrates, the tetrasubstituted chiral allenes were successfully generated using this strategy. Finally, theoretical calculations indicate that the transmetallation of the allenylcopper species is the rate-limiting step of this transformation.

Recovery of non-reference sequences missing from the human reference genome.

BACKGROUND: The non-reference sequences (NRS) represent structure variations in human genome with potential functional significance. However, besides the known insertions, it is currently unknown whether other types of structure variations with NRS exist. RESULTS: Here, we compared 31 human de novo assemblies with the current reference genome to identify the NRS and their location. We resolved the precise location of 6113 NRS adding up to 12.8 Mb. Besides 1571 insertions, we detected 3041 alternate alleles, which were defined as having less than 90% (or none) identity with the reference alleles. These alternate alleles overlapped with 1143 protein-coding genes including a putative novel MHC haplotype. Further, we demonstrated that the alternate alleles and their flanking regions had high content of tandem repeats, indicating that their origin was associated with tandem repeats. CONCLUSIONS: Our study detected a large number of NRS including many alternate alleles which are previously uncharacterized. We suggested that the origin of alternate alleles was associated with tandem repeats. Our results enriched the spectrum of genetic variations in human genome.

Palladium-Catalyzed Enantioselective Thiocarbonylation of Styrenes.

A highly enantioselective thiocarbonylation of styrenes with CO and thiols has been achieved by Pd catalysis, providing highly enantioenriched thioesters in good to excellent yields. Key to the successful execution of this reaction is the use of a chiral sulfoxide-(P-dialkyl)-phosphine (SOP) ligands. This thiocarbonylation proceeds smoothly under mild reaction conditions (1 atm CO and 0 °C) and displays broad substrate scope. Also demonstrated is that this transformation can be conducted using surrogates of CO, greatly increasing the safety aspects of running the reaction. The generality and utility of the method is manifested by its application to the synthetic transformations of thioester products and the direct acylation of cysteine-containing dipeptides. A primary mechanism was investigated and a plausible catalytic cycle was proposed.

Size-Dependent Effect of Cu2O Nanocubes in Electrochemical and Photocatalytic Properties.

Cu2O nanocubes with different size (ranging from 20 nm to 400 nm) were prepared by a seed-mediated method to systematically explore the strong size-dependent properties in photocatalytic degradation of methyl orange (MO). Cu2O nanotubes were characterized by TEM, XRD, UV-Vis measurements. The size-dependent photocatalytic efficiency of the Cu2O nanocubes was evaluated by degradation of methyl orange (MO) in water under visible light (λ > 420 nm) irradiation. Furthermore, the photocurrent, linear sweep voltammetry (LSV) and electrochemical impedance spectra (EIS) measurements were applied to elucidate the size-dependent properties of Cu2O nanocubes, which demonstrated that smaller Cu2O nanocubes with certain length (30 nm) showed higher current density, faster electron transfer and lower rate of charge recombination in their exposed (100) facet. Therefore, 30 nm Cu2O nanocubes showed stronger visible light absorption capacity and higher photocatalytic activity in MO degradation among a series of nanocubes (20, 30, 100, 130, 200 and 400 nm) and their corresponding photocatalytic activities decreased with increasing the particles sizes.

PtdIns(3,4,5)P3 is constitutively synthesized and required for spindle translocation during meiosis in mouse oocytes.

Prior to ovulation, mammalian oocytes complete their first meiotic division and arrest at metaphase II. During this marked asymmetric cell division, the meiotic spindle moves dramatically from the center of the oocyte to the cortex to facilitate segregation of half of its chromosomal content into the diminutive first polar body. Recent investigations have documented crucial roles for filamentous actin (F-actin) in meiotic spindle translocation. However, the identity of the upstream regulators responsible for these carefully orchestrated movements has remained elusive. Utilizing fluorescently tagged probes and time-lapse confocal microscopy, we document that phosphatidylinositol 3,4,5-trisphosphate [PtdIns(3,4,5)P3] is constitutively synthesized with spatial and temporal dynamics similar to that of F-actin and Formin 2 (Fmn2). Blockage of PtdIns(3,4,5)P3 synthesis by LY294002, a specific inhibitor of phosphoinositide 3-kinase (PI3K), disrupts cytoplasmic F-actin organization and meiotic spindle migration to the cortex. F-actin nucleator Fmn2 and Rho GTPase Cdc42 play roles in mediating the effect of PtdIns(3,4,5)P3 on F-actin assembly. Moreover, the spatial and temporal dynamics of PtdIns(3,4,5)P3 is impaired by depletion of MATER or Filia, two oocyte proteins encoded by maternal effect genes. Thus, PtdIns(3,4,5)P3 is synthesized during meiotic maturation and acts upstream of Cdc42 and Fmn2, but downstream of MATER/Filia proteins to regulate the F-actin organization and spindle translocation to the cortex during mouse oocyte meiosis.

The effect of chronic lymphocytic thyroiditis on patients with thyroid cancer.

BACKGROUND: The purpose of this study was to investigate the association between chronic lymphocytic thyroiditis (CLT) and malignant tumors of the thyroid. METHODS: A retrospective review of 647 patients who underwent thyroid surgery at the Department of Breast and Thyroid Surgery in Anhui Provincial Hospital, China in 2012 was performed. The clinicopathological characteristics of patients with thyroid malignancies and CLT were collected. CLT was diagnosed by histopathological method. RESULTS: Among 647 patients, 144 patients had thyroid malignancies and 108 patients had been diagnosed with CLT. Moreover, in total, 44 patients had thyroid malignancies coexistent with CLT: forty-one (93.2%) patients had been diagnosed with the papillary thyroid cancer (PTC); two (4.5%) patients suffered from medullary carcinoma; and one (2.3%) patient suffered from lymphoma. The morbidity of thyroid malignancies in patients with CLT was significantly higher than that in patients without CLT (40.7% versus 18.6%; P <0.001). A female preponderance was observed in the patients with CLT compared with those without CLT (P <0.001). There was no statistically significant difference in the tumor size (P = 0.073), multifocality (P = 0.0871), neck lymph node metastasis (P = 0.350), age (P = 0.316), microcarcinoma (P = 0.983) and tumor-node-metastasis (TNM) stage (P = 0.949) between the patients of thyroid malignancies with CLT and without CLT. CONCLUSIONS: Female predominance was observed in patients with CLT. CLT may have no effect on the progression of thyroid malignant tumor. Nevertheless, the influences of CLT on the prognosis of the thyroid carcinoma still need to be investigated with a larger sample size.

Evaluation of compliance of CONSORT-CHM formula 2017 in randomized controlled trials of Chinese herbal medicine formulas: protocol of a five-year review.

Background: The CONSORT Extension for Chinese Herbal Medicine Formula 2017 (CONSORT-CHM Formula 2017) has established a reporting standard for randomized controlled trials (RCTs) of Chinese Herbal Medicine Formula (CHMF) interventions; however, its adherence and implications for the design and execution of study design remain ambiguous. It is necessary to evaluate the level of compliance with the CONSORT-CHM Formula 2017 in RCTs conducted over the past 5 years, and to determine the reporting quality of clinical trials in this field. Methods: First, a systematic search is conducted for RCTs on CHMF in EBM Reviews, Allied and Complementary Medicine (AMED), Embase, Ovid-MEDLINE(R), Wanfang data, China National Knowledge Infrastructure (CNKI), VIP Chinese Medical Journal Database (VIP) and Chinese Biomedical Literature (CBM) database, that encompassed CHMF interventional RCTs published from 1 January 2018 to 8 June 2022, with language restriction to English or Chinese. Second, a descriptive analysis will be performed regarding the study design and general characteristics of the included trials. Third, for the quality assessment, we have subdivided the CONSORT-CHM Formula 2017 checklist (consisting of 22 extended items) into a total of 42 sub-questions to facilitate scoring, with a specific focus on the description, quality control, and safety assessment of CHMF interventions. Professional training and a pilot test on 100 randomly selected articles will be provided for all reviewers. Throughout this process, a standard operating procedure (SOP) for quality assessment will be developed to ensure consistency. Each item will be assessed by two reviewers in a paired back-to-back manner, and the compliance rate will be calculated to assess inter-rater agreement. Discussion: This review will identify the current reporting characteristics and quality of CHMF interventional studies and further evaluate the impact of CONSORT-CHM Formula 2017. The results may provide suggestions for future application or promotion of the guideline. Registration: The study has been registered on Open Science Framework (https://osf.io/xpn7f).

Xiaoyao Pills, a Chinese patent medicine, treats mild and moderate depression: A randomized clinical trial combined with DNA methylation analysis.

BACKGROUND: Xiaoyao pills (XYP) is a commercial Chinese patent medicine used in the treatment of depression. However, the mechanisms underlying its therapeutic effects, as well as the patients who can benefit from XYP, have not been evaluated so far. OBJECTIVES: To this end, we conducted a double-blinded, random, and placebo-controlled clinical trial of orally administered XYP in patients with depression. METHODS: The 17-item Hamilton Depression Rating Scale (HAMD-17) scores were recorded at baseline, and every 2 weeks after the start of treatment. To further elucidate the epigenetic mechanism of XYP, we performed mRNA sequencing and genome-wide DNA methylation sequencing using peripheral blood leukocytes of patients and healthy. RESULTS: XYP effectively alleviated the symptoms in patients with mild or moderate depressive disorders, particularly that of psychomotor retardation. XYP restored aberrant gene expression and DNA methylation patterns associated with depression, and the normalization of DNA methylation correlated with downregulation of several genes. In addition, altered DNA methylation levels in the XYP-treated samples were attributed to increased expression of the DNA methyltransferase DNMT1. CONCLUSIONS: Our study provides new insights into the epigenetic mechanism underlying depression and the therapeutic effects of XYP, along with an experimental basis for using XYP in the treatment of depression. TRIAL REGISTRATION INFORMATION: The name of the registry and number: U.S. CLINICAL TRIALS REGISTRY: The link to the registration: ClinicalTrials.gov ISRCTN12746343 (https://www.isrctn.com/ISRCTN12746343). The name of the trial register is "Efficacy and safety of the Xiaoyao pill for improving the clinical symptoms of stagnation of liver qi (chi) and spleen deficiency". The clinical trial registration number is ISRCTN12746343.

Stereospecific synthesis of silicon-stereogenic optically active silylboranes and general synthesis of chiral silyl Anions.

Silicon-stereogenic optically active silylboranes could potentially allow the formation of chiral silyl nucleophiles as well as the synthesis of various chiral silicon compounds. However, the synthesis of such silicon-stereogenic silylboranes has not been achieved so far. Here, we report the synthesis of silicon-stereogenic optically active silylboranes via a stereospecific Pt(PPh3)4-catalyzed Si-H borylation of chiral hydrosilanes, which are synthesized by stoichiometric and catalytic asymmetric synthesis, in high yield and very high or perfect enantiospecificity (99% es in one case, and >99% es in the others) with retention of the configuration. Furthermore, we report a practical approach to generate silicon-stereogenic silyl nucleophiles with high enantiopurity and configurational stability using MeLi activation. This protocol is suitable for the stereospecific and general synthesis of silicon-stereogenic trialkyl-, dialkylbenzyl-, dialkylaryl-, diarylalkyl-, and alkylary benzyloxy-substituted silylboranes and their corresponding silyl nucleophiles with excellent enantiospecificity (>99% es except one case of 99% es). Transition-metal-catalyzed C-Si bond-forming cross-coupling reactions and conjugate-addition reactions are also demonstrated. The mechanisms underlying the stability and reactivity of such chiral silyl anion were investigated by combining NMR spectroscopy and DFT calculations.

Diagnostic value of radiomics in predicting Ki-67 and cytokeratin 19 expression in hepatocellular carcinoma: a systematic review and meta-analysis.

Background: Radiomics have been increasingly used in the clinical management of hepatocellular carcinoma (HCC), such as markers prediction. Ki-67 and cytokeratin 19 (CK-19) are important prognostic markers of HCC. Radiomics has been introduced by many researchers in the prediction of these markers expression, but its diagnostic value remains controversial. Therefore, this review aims to assess the diagnostic value of radiomics in predicting Ki-67 and CK-19 expression in HCC. Methods: Original studies were systematically searched in PubMed, EMBASE, Cochrane Library, and Web of Science from inception to May 2023. All included studies were evaluated by the radiomics quality score. The C-index was used as the effect size of the performance of radiomics in predicting Ki-67and CK-19 expression, and the positive cutoff values of Ki-67 label index (LI) were determined by subgroup analysis and meta-regression. Results: We identified 34 eligible studies for Ki-67 (18 studies) and CK-19 (16 studies). The most common radiomics source was magnetic resonance imaging (MRI; 25/34). The pooled C-index of MRI-based models in predicting Ki-67 was 0.89 (95% CI:0.86-0.92) in the training set, and 0.87 (95% CI: 0.82-0.92) in the validation set. The pooled C-index of MRI-based models in predicting CK-19 was 0.86 (95% CI:0.81-0.90) in the training set, and 0.79 (95% CI: 0.73-0.84) in the validation set. Subgroup analysis suggested Ki-67 LI cutoff was a significant source of heterogeneity (I 2 = 0.0% P>0.05), and meta-regression showed that the C-index increased as Ki-67 LI increased. Conclusion: Radiomics shows promising diagnostic value in predicting positive Ki-67 or CK-19 expression. But lacks standardized guidelines, which makes the model and variables selection dependent on researcher experience, leading to study heterogeneity. Therefore, standardized guidelines are warranted for future research. Systematic Review Registration: https://www.crd.york.ac.uk/PROSPERO/, identifier CRD42023427953.

Design and characterization of a high-resolution multiple-SNP capture array by target sequencing for sheep.

The efficiency of molecular breeding largely depends on inexpensive genotyping arrays. In this study, we aimed to develop an ovine high-resolution multiple-single-nucleotide polymorphism (SNP) capture array, based on genotyping by target sequencing (GBTS) system with capture-in-solution (liquid chip) technology. All the markers were from 40K captured regions, including genes located within selective sweep regions, breed-specific regions, quantitative trait loci (QTL), and the potential functional SNPs on the sheep genome. The results showed that a total of 210K high-quality SNPs were identified in the 40K regions, indicating a high average capture ratio (99.7%) for the target genomic regions. Using genotyped data (n = 317) from liquid chip technology, we further performed genome-wide association studies (GWAS) to detect the genetic loci affecting sheep hair types and teat number. A single significant association signal for hair types was identified on 6.7-7.1 Mb of chromosome 25. The IRF2BP2 gene (chr25: 7,067,974-7,071,785), which is located within this genomic region, has been previously known to be involved in hair/wool traits in sheep. The results further showed a new candidate region around 26.4 Mb of chromosome 13, between the ARHGAP21 and KIAA1217 genes, that was significantly related to teat number in sheep. The haplotype patterns of this region also showed differences in animals with 2, 3, or 4 teats. Advances in using the high-accuracy and low-cost liquid chip are expected to accelerate sheep genomic and breeding studies in the coming years.

Large-scale genotyping platforms are valuable tools for animal selection and breeding programs. The bead chip has been widely used in both research and commercial applications for a long time. A highly efficient and economical genotyping platform has been developed recently. In the present study, by combining the advantages of resequencing and bead chips, we developed a high-resolution capture array based on target sequencing with capture-in-solution technology (liquid chip), including updated functional probes according to the latest research. We further evaluated this approach by using 317 individuals and found that 210K single-nucleotide polymorphisms can be accurately genotyped, confirming the ratio of the captured regions compared with the designed rations is around 99.7%. Genome-wide association studies conducted using this chip suggested IRF2BP2 gene may be involved in hair types and ARHGAP21-KIAA1217 locus may be related to teats number. The liquid chip with high accuracy and low cost can be widely used in genome-wide association studies and genome selection, supporting efforts in molecular breeding and genetic improvement of sheep.

Genomic analysis of indigenous goats in Southwest Asia reveals evidence of ancient adaptive introgression related to desert climate.

Understanding how evolutionary pressures related to climate change have shaped the current genetic background of domestic animals is a fundamental pursuit of biology. Here, we generated whole-genome sequencing data from native goat populations in Iraq and Pakistan. Combined with previously published data on modern, ancient (Late Neolithic to Medieval periods), and wild Capra species worldwide, we explored the genetic population structure, ancestry components, and signatures of natural positive selection in native goat populations in Southwest Asia (SWA). Results revealed that the genetic structure of SWA goats was deeply influenced by gene flow from the eastern Mediterranean during the Chalcolithic period, which may reflect adaptation to gradual warming and aridity in the region. Furthermore, comparative genomic analysis revealed adaptive introgression of the KITLG locus from the Nubian ibex ( C. nubiana) into African and SWA goats. The frequency of the selected allele at this locus was significantly higher among goat populations located near northeastern Africa. These results provide new insights into the genetic composition and history of goat populations in the SWA region.