Scans for Signatures of Selection in Genomes of Wagyu and Buryat Cattle Breeds Reveal Candidate Genes and Genetic Variants for Adaptive Phenotypes and Production Traits.

Past and ongoing selection shapes the genomes of livestock breeds. Identifying such signatures of selection allows for uncovering the genetic bases of affected phenotypes, including economically important traits and environmental adaptations, for the further improvement of breed genetics to respond to climate and economic challenges. Turano-Mongolian cattle are a group of taurine breeds known for their adaptation to extreme environmental conditions and outstanding production performance. Buryat Turano-Mongolian cattle are among the few breeds adapted to cold climates and poor forage. Wagyu, on the other hand, is famous for high productivity and unique top-quality marbled meat. We used hapFLK, the de-correlated composite of multiple signals (DCMS), PBS, and FST methods to search for signatures of selection in their genomes. The scans revealed signals in genes related to cold adaptation (e.g., STAT3, DOCK5, GSTM3, and CXCL8) and food digestibility (SI) in the Buryat breed, and growth and development traits (e.g., RBFOX2 and SHOX2) and marbling (e.g., DGAT1, IQGAP2, RSRC1, and DIP2B) in Wagyu. Several putatively selected genes associated with reproduction, immunity, and resistance to pathogens were found in both breed genomes. The results of our work could be used for creating new productive adapted breeds or improving the extant breeds.

The chromosome-scale genome assembly for the West Nile vector Culex quinquefasciatus uncovers patterns of genome evolution in mosquitoes.

BACKGROUND: Understanding genome organization and evolution is important for species involved in transmission of human diseases, such as mosquitoes. Anophelinae and Culicinae subfamilies of mosquitoes show striking differences in genome sizes, sex chromosome arrangements, behavior, and ability to transmit pathogens. However, the genomic basis of these differences is not fully understood. METHODS: In this study, we used a combination of advanced genome technologies such as Oxford Nanopore Technology sequencing, Hi-C scaffolding, Bionano, and cytogenetic mapping to develop an improved chromosome-scale genome assembly for the West Nile vector Culex quinquefasciatus. RESULTS: We then used this assembly to annotate odorant receptors, odorant binding proteins, and transposable elements. A genomic region containing male-specific sequences on chromosome 1 and a polymorphic inversion on chromosome 3 were identified in the Cx. quinquefasciatus genome. In addition, the genome of Cx. quinquefasciatus was compared with the genomes of other mosquitoes such as malaria vectors An. coluzzi and An. albimanus, and the vector of arboviruses Ae. aegypti. Our work confirms significant expansion of the two chemosensory gene families in Cx. quinquefasciatus, as well as a significant increase and relocation of the transposable elements in both Cx. quinquefasciatus and Ae. aegypti relative to the Anophelines. Phylogenetic analysis clarifies the divergence time between the mosquito species. Our study provides new insights into chromosomal evolution in mosquitoes and finds that the X chromosome of Anophelinae and the sex-determining chromosome 1 of Culicinae have a significantly higher rate of evolution than autosomes. CONCLUSION: The improved Cx. quinquefasciatus genome assembly uncovered new details of mosquito genome evolution and has the potential to speed up the development of novel vector control strategies.

Patterns of genetic differentiation imply distinct phylogeographic history of the mosquito species Anopheles messeae and Anopheles daciae in Eurasia.

Detailed knowledge of phylogeography is important for control of mosquito species involved in the transmission of human infectious diseases. Anopheles messeae is a geographically widespread and genetically diverse dominant vector of malaria in Eurasia. A closely related species, An. daciae, was originally distinguished from An. messeae based on five nucleotide substitutions in its ribosomal DNA (rDNA). However, the patterns of phylogeographic history of these species in Eurasia remain poorly understood. Here, using internal transcribed spacer 2 (ITS2) of rDNA and karyotyping for the species identification we determined the composition of five Anopheles species in 28 locations in Eurasia. Based on the frequencies of 11 polymorphic chromosomal inversions used as genetic markers, a large-scale population genetics analysis was performed of 1932 mosquitoes identified as An. messeae, An. daciae and their hybrids. The largest genetic differences between the species were detected in the X sex chromosome suggesting a potential involvement of this chromosome in speciation. The frequencies of autosomal inversions in the same locations differed by 13%-45% between the species demonstrating a restricted gene flow between the species. Overall, An. messeae was identified as a diverse species with a more complex population structure than An. daciae. The clinal gradients in frequencies of chromosomal inversions were determined in both species implicating their possible involvement in climate adaptations. The frequencies of hybrids were low ~1% in northern Europe but high up to 50% in south-eastern populations. Thus, our study revealed critical differences in patterns of phylogeographic history between An. messeae and An. daciae in Eurasia. This knowledge will help to predict the potential of the malaria transmission in the northern territories of the continent.

Early-onset gynecological tumors in DNA repair-deficient xeroderma pigmentosum group C patients: a case series.

BACKGROUND: Xeroderma pigmentosum (XP) is a group of rare hereditary disorders with highly increased risk of skin tumors due to defective DNA repair. Recently we reported 34-fold increased risk of internal tumors in XP patients in comparison with general population. The molecular data and clinical practice on the internal tumors treatment in XP patients is limited and scarcely represented in the medical literature. In this work, we describe young patients with constitutive biallelic deactivation of the XPC gene developing gynecological tumors with somatic DICER1 mutations. METHODS: Whole genome sequencing was used to analyze in detail somatic mutational landscape and driver events of these rare tumors. RESULTS: We describe five early-onset gynecological tumors in four xeroderma pigmentosum group C (XP-C) young patients (11 to 19 years old) including vaginal embryonal rhabdomyosarcomas in monozygotic twin sisters, juvenile granulosa-cell tumor of the ovary and poorly differentiated stage IA Sertoli-Leydig cell tumor in 19-years old patient, and FIGO stage IC1 tumor of ovary in 13-years old patient. XP-C ovarian tumors harbor 4.4 times more single base substitutions than sporadic tissue-matched cancers and demonstrate XP-C specific mutation signature with strong transcriptional bias indicating inability of the cells to repair bulky DNA lesions of unknown etiology. A special mode of treatment was applied to avoid usage of chemotherapy which is toxic for XP patients. CONCLUSIONS: XP-C status should be accounted for prevention and specific treatment of gynecological tumors in young DNA repair-deficient XP patients.

Xeroderma pigmentosum group C (XP-C) is a rare inherited disorder resulting in a highly increased risk of skin and internal cancers due to the inability to efficiently repair DNA. In this study, we described four young XP-C patients who developed early-onset tumors affecting the female reproductive organs. We describe how we cared for these patients in the clinic. We looked at the genetic material within the tumors to better understand the mechanisms through which these tumors developed. We observed high numbers of specific types of changes in DNA, which are not typical for sporadic (non-inherited) gynecological tumors, but are characteristic of internal XP-C tumors. Further studies are needed to better understand the nature of these changes. Our findings highlight the important role of DNA repair in human tissues and cancer risk, and might inform future strategies for tumor prevention in XP-C patients.

Brazilian XP-E siblings carrying a novel DDB2 variant developed early-onset melanoma: a case report.

BACKGROUND: Xeroderma pigmentosum group E (XP-E) is one of the least common forms of XP, a rare syndrome where patients are prone to develop skin cancer in exposed sunlight areas. XP-E patients are generally not diagnosed until they are adults due to the mild phenotype. CASE PRESENTATION: two XP-E siblings, female, 23 years, and male, 25 years, from a Brazilian consanguineous family carrying the novel missense pathogenic variant in DDB2 gene, NM_000107.3:c.1027G > C, associated with skin cancer early-onset and severe phenotype, as nodular melanoma in the cornea and in the ear. CONCLUSION: The assessment of genomic variant pathogenicity was a challenge since this family belongs to an underrepresented population in genomic databases. Given the scarcity of literature documenting XP-E cases and the challenges encountered in achieving an early diagnosis, this report emphasizes the imperative of sun protection measures in XP-E patients. Additionally, it highlights the detrimental impact of the COVID-19 pandemic on cancer diagnosis, leading to the manifestation of a severe phenotype in affected individuals.

Genomic mutation landscape of skin cancers from DNA repair-deficient xeroderma pigmentosum patients.

Xeroderma pigmentosum (XP) is a genetic disorder caused by mutations in genes of the Nucleotide Excision Repair (NER) pathway (groups A-G) or in Translesion Synthesis DNA polymerase η (V). XP is associated with an increased skin cancer risk, reaching, for some groups, several thousand-fold compared to the general population. Here, we analyze 38 skin cancer genomes from five XP groups. We find that the activity of NER determines heterogeneity of the mutation rates across skin cancer genomes and that transcription-coupled NER extends beyond the gene boundaries reducing the intergenic mutation rate. Mutational profile in XP-V tumors and experiments with POLH knockout cell line reveal the role of polymerase η in the error-free bypass of (i) rare TpG and TpA DNA lesions, (ii) 3' nucleotides in pyrimidine dimers, and (iii) TpT photodimers. Our study unravels the genetic basis of skin cancer risk in XP and provides insights into the mechanisms reducing UV-induced mutagenesis in the general population.

Phylogenomics revealed migration routes and adaptive radiation timing of Holarctic malaria mosquito species of the Maculipennis Group.

BACKGROUND: Phylogenetic analyses of closely related species of mosquitoes are important for better understanding the evolution of traits contributing to transmission of vector-borne diseases. Six out of 41 dominant malaria vectors of the genus Anopheles in the world belong to the Maculipennis Group, which is subdivided into two Nearctic subgroups (Freeborni and Quadrimaculatus) and one Palearctic (Maculipennis) subgroup. Although previous studies considered the Nearctic subgroups as ancestral, details about their relationship with the Palearctic subgroup, and their migration times and routes from North America to Eurasia remain controversial. The Palearctic species An. beklemishevi is currently included in the Nearctic Quadrimaculatus subgroup adding to the uncertainties in mosquito systematics. RESULTS: To reconstruct historic relationships in the Maculipennis Group, we conducted a phylogenomic analysis of 11 Palearctic and 2 Nearctic species based on sequences of 1271 orthologous genes. The analysis indicated that the Palearctic species An. beklemishevi clusters together with other Eurasian species and represents a basal lineage among them. Also, An. beklemishevi is related more closely to An. freeborni, which inhabits the Western United States, rather than to An. quadrimaculatus, a species from the Eastern United States. The time-calibrated tree suggests a migration of mosquitoes in the Maculipennis Group from North America to Eurasia about 20-25 million years ago through the Bering Land Bridge. A Hybridcheck analysis demonstrated highly significant signatures of introgression events between allopatric species An. labranchiae and An. beklemishevi. The analysis also identified ancestral introgression events between An. sacharovi and its Nearctic relative An. freeborni despite their current geographic isolation. The reconstructed phylogeny suggests that vector competence and the ability to enter complete diapause during winter evolved independently in different lineages of the Maculipennis Group. CONCLUSIONS: Our phylogenomic analyses reveal migration routes and adaptive radiation timing of Holarctic malaria vectors and strongly support the inclusion of An. beklemishevi into the Maculipennis Subgroup. Detailed knowledge of the evolutionary history of the Maculipennis Subgroup provides a framework for examining the genomic changes related to ecological adaptation and susceptibility to human pathogens. These genomic variations may inform researchers about similar changes in the future providing insights into the patterns of disease transmission in Eurasia.

Investigating the genetic basis of vertebrate dispersal combining RNA-seq, RAD-seq and quantitative genetics.

Although animal dispersal is known to play key roles in ecological and evolutionary processes such as colonization, population extinction and local adaptation, little is known about its genetic basis, particularly in vertebrates. Untapping the genetic basis of dispersal should deepen our understanding of how dispersal behaviour evolves, the molecular mechanisms that regulate it and link it to other phenotypic aspects in order to form the so-called dispersal syndromes. Here, we comprehensively combined quantitative genetics, genome-wide sequencing and transcriptome sequencing to investigate the genetic basis of natal dispersal in a known ecological and evolutionary model of vertebrate dispersal: the common lizard, Zootoca vivipara. Our study supports the heritability of dispersal in semi-natural populations, with less variation attributable to maternal and natal environment effects. In addition, we found an association between natal dispersal and both variation in the carbonic anhydrase (CA10) gene, and in the expression of several genes (TGFB2, SLC6A4, NOS1) involved in central nervous system functioning. These findings suggest that neurotransmitters (serotonin and nitric oxide) are involved in the regulation of dispersal and shaping dispersal syndromes. Several genes from the circadian clock (CRY2, KCTD21) were also differentially expressed between disperser and resident lizards, supporting that the circadian rhythm, known to be involved in long-distance migration in other taxa, might affect dispersal as well. Since neuronal and circadian pathways are relatively well conserved across vertebrates, our results are likely to be generalisable, and we therefore encourage future studies to further investigate the role of these pathways in shaping dispersal in vertebrates.

Integrative Pan-Cancer Genomic and Transcriptomic Analyses of Refractory Metastatic Cancer.

Metastatic relapse after treatment is the leading cause of cancer mortality, and known resistance mechanisms are missing for most treatments administered to patients. To bridge this gap, we analyze a pan-cancer cohort (META-PRISM) of 1,031 refractory metastatic tumors profiled via whole-exome and transcriptome sequencing. META-PRISM tumors, particularly prostate, bladder, and pancreatic types, displayed the most transformed genomes compared with primary untreated tumors. Standard-of-care resistance biomarkers were identified only in lung and colon cancers-9.6% of META-PRISM tumors, indicating that too few resistance mechanisms have received clinical validation. In contrast, we verified the enrichment of multiple investigational and hypothetical resistance mechanisms in treated compared with nontreated patients, thereby confirming their putative role in treatment resistance. Additionally, we demonstrated that molecular markers improve 6-month survival prediction, particularly in patients with advanced breast cancer. Our analysis establishes the utility of the META-PRISM cohort for investigating resistance mechanisms and performing predictive analyses in cancer. SIGNIFICANCE: This study highlights the paucity of standard-of-care markers that explain treatment resistance and the promise of investigational and hypothetical markers awaiting further validation. It also demonstrates the utility of molecular profiling in advanced-stage cancers, particularly breast cancer, to improve the survival prediction and assess eligibility to phase I clinical trials. This article is highlighted in the In This Issue feature, p. 1027.

Genomic Instability and Protumoral Inflammation Are Associated with Primary Resistance to Anti-PD-1 + Antiangiogenesis in Malignant Pleural Mesothelioma.

Cancer immunotherapy combinations have recently been shown to improve the overall survival of advanced mesotheliomas, especially for patients responding to those treatments. We aimed to characterize the biological correlates of malignant pleural mesotheliomas' primary resistance to immunotherapy and antiangiogenics by testing the combination of pembrolizumab, an anti-PD-1 antibody, and nintedanib, a pan-antiangiogenic tyrosine kinase inhibitor, in the multicenter PEMBIB trial (NCT02856425). Thirty patients with advanced malignant pleural mesothelioma were treated and explored. Unexpectedly, we found that refractory patients were actively recruiting CD3+CD8+ cytotoxic T cells in their tumors through CXCL9 tumor release upon treatment. However, these patients displayed high levels of somatic copy-number alterations in their tumors that correlated with high blood and tumor levels of IL6 and CXCL8. Those proinflammatory cytokines resulted in higher tumor secretion of VEGF and tumor enrichment in regulatory T cells. Advanced mesothelioma should further benefit from stratified combination therapies adapted to their tumor biology. SIGNIFICANCE: Sequential explorations of fresh tumor biopsies demonstrated that mesothelioma resistance to anti-PD-1 + antiangiogenics is not due to a lack of tumor T-cell infiltration but rather due to adaptive immunosuppressive pathways by tumors, involving molecules (e.g., IL6, CXCL8, VEGF, and CTLA4) that are amenable to targeted therapies. This article is highlighted in the In This Issue feature, p. 799.

A mitochondria-specific mutational signature of aging: increased rate of A > G substitutions on the heavy strand.

The mutational spectrum of the mitochondrial DNA (mtDNA) does not resemble any of the known mutational signatures of the nuclear genome and variation in mtDNA mutational spectra between different organisms is still incomprehensible. Since mitochondria are responsible for aerobic respiration, it is expected that mtDNA mutational spectrum is affected by oxidative damage. Assuming that oxidative damage increases with age, we analyse mtDNA mutagenesis of different species in regards to their generation length. Analysing, (i) dozens of thousands of somatic mtDNA mutations in samples of different ages (ii) 70053 polymorphic synonymous mtDNA substitutions reconstructed in 424 mammalian species with different generation lengths and (iii) synonymous nucleotide content of 650 complete mitochondrial genomes of mammalian species we observed that the frequency of AH > GH substitutions (H: heavy strand notation) is twice bigger in species with high versus low generation length making their mtDNA more AH poor and GH rich. Considering that AH > GH substitutions are also sensitive to the time spent single-stranded (TSSS) during asynchronous mtDNA replication we demonstrated that AH > GH substitution rate is a function of both species-specific generation length and position-specific TSSS. We propose that AH > GH is a mitochondria-specific signature of oxidative damage associated with both aging and TSSS.

Overcoming resistance to αPD-1 of MMR-deficient tumors with high tumor-induced neutrophils levels by combination of αCTLA-4 and αPD-1 blockers.

BACKGROUND: Clinical studies have highlighted the efficacy of anti-programmed death 1 (αPD-1) monoclonal antibodies in patients with DNA mismatch repair-deficient (MMRD) tumors. However, the responsiveness of MMRD cancers to αPD-1 therapy is highly heterogeneous, and the origins of this variability remain not fully understood. METHODS: 4T1 and CT26 mouse tumor cell lines were inactivated for the MMRD gene Msh2, leading to a massive accumulation of mutations after serial passages of cells. Insertions/deletion events and mutation load were evaluated by whole exome sequencing. Mice bearing highly mutated MMRD tumor or parental tumors were treated with αPD-1 and tumor volume was monitored. Immune cell type abundance was dynamically assessed in the tumor microenvironment and the blood by flow cytometry. Neutrophils were depleted in mice using αLY6G antibody, and regulatory T (Treg) cell population was reduced with αCD25 or anti-cytotoxic T-lymphocytes-associated protein 4 (αCTLA-4) antibodies. Patients with MMRD tumors treated with immune checkpoint blockade-based therapy were retrospectively identified and neutrophil-to-lymphocyte ratio (NLR) was evaluated and examined for correlation with clinical benefit. RESULTS: By recapitulating mismatch repair deficiency in different mouse tumor models, we revealed that elevated circulating tumor-induced neutrophils (TIN) in hypermutated MMRD tumors hampered response to αPD-1 monotherapy. Importantly, depletion of TIN using αLy-6G antibody reduced Treg cells and restored αPD-1 response. Conversely, targeting Treg cells by αCD25 or αCTLA-4 antibodies limited peripheral TIN accumulation and elicited response in αPD-1-resistant MMRD tumors, highlighting a crosstalk between TIN and Treg cells. Thus, αPD-1+αCTLA-4 combination overcomes TIN-induced resistance to αPD-1 in mice bearing MMRD tumors. Finally, in a cohort of human (high microsatellite instability)/MMRD tumors we revealed that early on-treatment change in the NLR ratio may predict resistance to αPD-1 therapy. CONCLUSIONS: TIN countered αPD-1 efficacy in MMRD tumors. Since αCTLA-4 could restrict TIN accumulation, αPD-1+αCTLA-4 combination overcomes αPD-1 resistance in hosts with hypermutated MMRD tumors displaying abnormal neutrophil accumulation.

Increased risk of internal tumors in DNA repair-deficient xeroderma pigmentosum patients: analysis of four international cohorts.

BACKGROUND: Xeroderma pigmentosum (XP) is a rare, autosomal, recessive DNA repair-deficiency disorder with a frequency of 1-3 per million livebirths in Europe and USA but with higher frequencies in isolated islands or in countries with a high level of consanguinity. XP is characterized by high incidence of skin cancers on sun-exposed sites. Recent improvement in life expectancy of XP patients suggests an increased risk of frequently aggressive and lethal internal tumors. Our purpose was to quantify relative risks of internal tumor development for XP patients by tumor type, XP-subtype, patients' ages and ethnicity through comparison with the US general population. METHODS: We analyzed four independent international well-characterized XP cohorts (from USA, UK, France and Brazil) with a total of 434 patients, where 11.3% developed internal tumors and compared them to the American general population. We also compiled, through PubMed/Medline, a dataset of 89 internal tumors in XP patients published between 1958 and 2020. RESULTS: In the combined 4-XP cohort, relative risk of internal tumors was 34 (95% confidence interval (CI) 25-47) times higher than in the general population (p-value = 1.0E-47) and tumor arose 50 years earlier. The XP-C group was at the highest risk for the 0-20 years old-patients (OR = 665; 95% CI 368-1200; p-value = 4.3E-30). The highest risks were observed for tumors of central nervous system (OR = 331; 95% CI 171-641; p-value = 2.4E-20), hematological malignancies (OR = 120; 95% CI 77-186; p-value = 3.7E-36), thyroid (OR = 74; 95% CI 31-179; p-value = 1.2E-8) and gynecological tumors (OR = 91; 95% CI 42-193; p-value = 3.5E-12). The type of mutation on the XPC gene is associated with different classes of internal tumors. The majority of French XP-C patients (80%) are originated from North Africa and carried the XPC delTG founder mutation specific from the South Mediterranean area. The OR is extremely high for young (0-20 years) patients with more than 1300-fold increase for the French XPs carrying the founder mutation. CONCLUSION: Because the age of XP population is increasing due to better sun-protection and knowledge of the disease, these results are of particular importance for the physicians to help in early prevention and detection of internal tumors in their XP patients. Few preventive blood analyses or simple medical imaging may help to better detect early cancer appearance in this population.

Frequency and Genomic Aspects of Intrinsic Resistance to Vismodegib in Locally Advanced Basal Cell Carcinoma.

PURPOSE: Vismodegib is approved for the treatment of locally advanced basal cell carcinoma (laBCC), but some cases demonstrate intrinsic resistance (IR) to the drug. We sought to assess the frequency of IR to vismodegib in laBCC and its underlying genomic mechanisms. EXPERIMENTAL DESIGN: Response to vismodegib was evaluated in a cohort of 148 laBCC patients. Comprehensive genomic and transcriptomic profiling was performed in a subset of five intrinsically resistant BCC (IR-BCC). RESULTS: We identified that IR-BCC represents 6.1% of laBCC in the studied cohort. Prior treatment with chemotherapy was associated with IR. Genetic events that were previously associated with acquired resistance (AR) in BCC or medulloblastoma were observed in three out of five IR-BCC. However, IR-BCCs were distinct by highly rearranged polyploid genomes. Functional analyses identified hyperactivation of the HIPPO-YAP and WNT pathways at RNA and protein levels in IR-BCC. In vitro assay on the BCC cell line further confirmed that YAP1 overexpression increases the cell proliferation rate. CONCLUSIONS: IR to vismodegib is a rare event in laBCC. IR-BCCs frequently harbor resistance mutations in the Hh pathway, but also are characterized by hyperactivation of the HIPPO-YAP and WNT pathways.

The functional genetic architecture of egg-laying and live-bearing reproduction in common lizards.

All amniotes reproduce either by egg-laying (oviparity), which is ancestral to vertebrates or by live-bearing (viviparity), which has evolved many times independently. However, the genetic basis of these parity modes has never been resolved and, consequently, its convergence across evolutionary scales is currently unknown. Here, we leveraged natural hybridizations between oviparous and viviparous common lizards (Zootoca vivipara) to describe the functional genes and genetic architecture of parity mode and its key traits, eggshell and gestation length, and compared our findings across vertebrates. In these lizards, parity trait genes were associated with progesterone-binding functions and enriched for tissue remodelling and immune system pathways. Viviparity involved more genes and complex gene networks than did oviparity. Angiogenesis, vascular endothelial growth and adrenoreceptor pathways were enriched in the viviparous female reproductive tissue, while pathways for transforming growth factor were enriched in the oviparous. Natural selection on these parity mode genes was evident genome-wide. Our comparison to seven independent origins of viviparity in mammals, squamates and fish showed that genes active in pregnancy were related to immunity, tissue remodelling and blood vessel generation. Therefore, our results suggest that pre-established regulatory networks are repeatedly recruited for viviparity and that these are shared at deep evolutionary scales.

Whole-Genome Resequencing Points to Candidate DNA Loci Affecting Body Temperature under Cold Stress in Siberian Cattle Populations.

Despite the economic importance of creating cold resilient cattle breeds, our knowledge of the genetic basis of adaptation to cold environments in cattle is still scarce compared to information on other economically important traits. Herein, using whole-genome resequencing of animals showing contrasting phenotypes on temperature maintenance under acute cold stress combined with the existing SNP (single nucleotide polymorphism) functional annotations, we report chromosomal regions and candidate SNPs controlling body temperature in the Siberian cattle populations. The SNP ranking procedure based on regional FST calculations, functional annotations, and the allele frequency difference between cold-tolerant and cold-sensitive groups of animals pointed to multiple candidate genes. Among these, GRIA4, COX17, MAATS1, UPK1B, IFNGR1, DDX23, PPT1, THBS1, CCL5, ATF1, PLA1A, PRKAG1, and NR1I2 were previously related to thermal adaptations in cattle. Other genes, for example KMT2D and SNRPA1, are known to be related to thermogenesis in mice and cold adaptation in common carp, respectively. This work could be useful for cattle breeding strategies in countries with harsh climates, including the Russian Federation.

Genotyping and Whole-Genome Resequencing of Welsh Sheep Breeds Reveal Candidate Genes and Variants for Adaptation to Local Environment and Socioeconomic Traits.

BACKGROUND: Advances in genetic tools applied to livestock breeding has prompted research into the previously neglected breeds adapted to harsh local environments. One such group is the Welsh mountain sheep breeds, which can be farmed at altitudes of 300 m above sea level but are considered to have a low productive value because of their poor wool quality and small carcass size. This is contrary to the lowland breeds which are more suited to wool and meat production qualities, but do not fare well on upland pasture. Herein, medium-density genotyping data from 317 individuals representing 15 Welsh sheep breeds were used alongside the whole-genome resequencing data of 14 breeds from the same set to scan for the signatures of selection and candidate genetic variants using haplotype- and SNP-based approaches. RESULTS: Haplotype-based selection scan performed on the genotyping data pointed to a strong selection in the regions of GBA3, PPARGC1A, APOB, and PPP1R16B genes in the upland breeds, and RNF24, PANK2, and MUC15 in the lowland breeds. SNP-based selection scan performed on the resequencing data pointed to the missense mutations under putative selection relating to a local adaptation in the upland breeds with functions such as angiogenesis (VASH1), anti-oxidation (RWDD1), cell stress (HSPA5), membrane transport (ABCA13 and SLC22A7), and insulin signaling (PTPN1 and GIGFY1). By contrast, genes containing candidate missense mutations in the lowland breeds are related to cell cycle (CDK5RAP2), cell adhesion (CDHR3), and coat color (MC1R). CONCLUSION: We found new variants in genes with potentially functional consequences to the adaptation of local sheep to their environments in Wales. Knowledge of these variations is important for improving the adaptative qualities of UK and world sheep breeds through a marker-assisted selection.

Demographic History, Adaptation, and NRAP Convergent Evolution at Amino Acid Residue 100 in the World Northernmost Cattle from Siberia.

Native cattle breeds represent an important cultural heritage. They are a reservoir of genetic variation useful for properly responding to agriculture needs in the light of ongoing climate changes. Evolutionary processes that occur in response to extreme environmental conditions could also be better understood using adapted local populations. Herein, different evolutionary histories of the world northernmost native cattle breeds from Russia were investigated. They highlighted Kholmogory as a typical taurine cattle, whereas Yakut cattle separated from European taurines approximately 5,000 years ago and contain numerous ancestral and some novel genetic variants allowing their adaptation to harsh conditions of living above the Polar Circle. Scans for selection signatures pointed to several common gene pathways related to adaptation to harsh climates in both breeds. But genes affected by selection from these pathways were mostly different. A Yakut cattle breed-specific missense mutation in a highly conserved NRAP gene represents a unique example of a young amino acid residue convergent change shared with at least 16 species of hibernating/cold-adapted mammals from six distinct phylogenetic orders. This suggests a convergent evolution event along the mammalian phylogenetic tree and fast fixation in a single isolated cattle population exposed to a harsh climate.

The new Internal Transcribed Spacer 2 diagnostic tool clarifies the taxonomic position and geographic distribution of the North American malaria vector Anopheles punctipennis.

BACKGROUND: The malaria mosquito Anopheles punctipennis, a widely distributed species in North America, is capable of transmitting human malaria and is actively involved in the transmission of the ungulate malaria parasite Plasmodium odocoilei. However, molecular diagnostic tools based on Internal Transcribed Spacer 2 (ITS2) of ribosomal DNA are lacking for this species. Anopheles punctipennis is a former member of the Anopheles maculipennis complex but its systematic position remains unclear. METHODS: In this study, ITS2 sequences were obtained from 276 An. punctipennis specimens collected in the eastern and midwestern United States and a simple and robust Restriction Fragment Length Polymorphism approach for species identification was developed. The maximum-likelihood phylogenetic tree was constructed based on ITS2 sequences available through this study and from GenBank for 20 species of Anopheles. RESULTS: The analysis demonstrated a consistent ITS2 sequence length and showed no indications of intragenomic variation among the samples based on ITS2, suggesting that An. punctipennis represents a single species in the studied geographic locations. In this study, An. punctipennis was found in urban, rural, and forest settings, suggesting its potential broad role in pathogen transmission. Phylogeny based on ITS2 sequence comparison demonstrated the close relationship of this species with other members of the Maculipennis group. CONCLUSIONS: This study developed molecular tools based on ITS2 sequences for the malaria vector An. punctipennis and clarified the phylogenetic position of the species within the Maculipennis group.