Complete mitochondrial genomes of patients from Thailand with cardiovascular diseases.

Several previous studies have reported that both variation and haplogroups of mitochondrial (mt) DNA were associated with various kinds of diseases, including cardiovascular diseases, in different populations, but such studies have not been carried out in Thailand. Here, we sequenced complete mtDNA genomes from 82 patients diagnosed with three types of cardiovascular disease, i.e., Hypertrophic Cardiomyopathy (HCM) (n = 26), Long Q-T Syndrome (LQTS) (n = 7) and Brugada Syndrome (BrS) (n = 49) and compared these with 750 previously published mitogenome sequences from interviewed normal individuals as a control group. Both patient and control groups are from the same geographic region of northeastern Thailand. We found 9, 2, and 5 novel mutations that were not both damaging and deleterious in HCM, LQTS, and BrS patients, respectively. Haplogroup R9c was significantly associated with HCM (P = 0.0032; OR = 62.42; 95%CI = 6.892-903.4) while haplogroup M12b was significantly associated with LQTS (P = 0.0039; OR = 32.93; 95% CI = 5.784-199.6). None of the haplogroups was found to be significantly associated with BrS. A significantly higher density of mtDNA variants in the rRNA genes was found in patients with HCM and BrS (P < 0.001) than in those with LQTS or the control group. Effects of detected SNPs in either protein coding or tRNA genes of all the mitogenome sequences were also predicted. Interestingly, three SNPs in two tRNA genes (MT-TA m.5618T>C and m.5631G>A heteroplasmic variants in two BrS patients and MT-TQ m.4392C>T novel homoplasmic variant in a HCM patient) were predicted to alter tRNA secondary structure, possibly leading to abnormal tRNA function.

Association between HLA alleles and haplotypes with age at diagnosis of type 1 diabetes in an admixed Brazilian population: A nationwide study.

To investigate the potential relationship between HLA alleles and haplotypes and the age at diagnosis of type 1 diabetes (T1DAgeD) in an admixed Brazilian population. This nationwide study was conducted in public clinics across 12 Brazilian cities. We collected demographic and genetic data from 1,600 patients with T1D. DNA samples were utilised to determine genomic ancestry (GA) and perform HLA typings for DRB1, DQA1 and DQB1. We explored allele and haplotype frequencies and GA in patients grouped by T1DAgeD categories (<6 years, ≥6-<11 years, ≥11-<19 years and ≥19 years) through univariate and multivariate analyses and primary component analyses. Additionally, we considered self-reported colour-race and identified a familiar history of T1D in first-degree relatives. The homozygosity index for DRB1~DQA1~DQB1 haplotypes exhibited the highest variation among T1DAgeD groups, and the percentages of Sub-Saharan African and European ancestries showed opposite trends in principal component analysis (PCA) analyses. Regarding the association of alleles and haplotypes with T1DAgeD, risk alleles such as HLA-DQB1*03:02g, -DQA1*03:01g, -02:01g, DRB1*04:05g and -04:02g were more frequently observed in heterozygosity or homozygosity in T1D patients with an early disease onset. Conversely, alleles such as DRB1*07:01g, -13:03g, DQB1*06:02g and DQA1*02:01 were more prevalent in older T1D patients. The combination DR3/DR4.5 was significantly associated with early disease onset. However, gender, GA, familiar history of T1D and self-reported colour-race identity did not exhibit significant associations with the onset of T1D. It is worth noting that the very common risk haplotype DRB1*03:01g~DQA1*05:01g~DQB1*02:01g did not differentiate between T1DAgeD groups. In the admixed Brazilian population, the high-risk haplotype DRB1*04:05~DQA1*03:01~DQB1*03:02 was more prevalent in individuals diagnosed before 6 years of age. In contrast, the protective alleles DQA1*01:02g, DQB1*06:02g, DRB1*07:01g and DRB1*13:03g and haplotypes DRB1*13:03g~DQA1*05:01g~DQB1*03:01g and DRB1*16:02g~DQA1*01:02g~DQB1*05:02g were more frequently observed in patients diagnosed in adulthood. Notably, these associations were independent of factors such as sex, economic status, GA, familiar history of T1D and region of birth in Brazil. These alleles and haplotypes contribute to our understanding of the disease onset heterogeneity and may have implications for early interventions when detected in association with well-known genomic risk or protection factors for T1D.

[Genetic study of a rare Chinese pedigree with a recombination occurring between the HLA-A/C loci in both parents].

OBJECTIVE: To analyze a Chinese pedigree with a recombination occurring between the HLA-A/C loci in both parents. METHODS: A patient who was planning to undergo hematopoietic stem cell transplantation due to "aplastic anemia" in February 2022 was selected as the study subject. Peripheral blood samples were collected from the patient, his parents and brother. HLA-A/C/B/DRB1/DQB1 high-resolution typing was carried out by using sequence-based typing and sequence-specific oligonucleotides. The recombination was identified by pedigree analysis. The HLA haplotype of each individual was identified by genealogical analysis. The parentage possibility was determined by short tandem repeat analysis. HLA-A/C/B/DRB1/DRB345/DQA1/DQB1/DPA1/DPB1 were determined with next-generation high-throughput sequence-based typing. The recombination sites were analyzed by family study. RESULTS: The high parentage possibilities of the family was confirmed by short tandem repeat analysis. Recombination was found between the HLA-A*24:02 A*33:03/C*14:03 in the paternally transmitted haplotype, whilst HLA-A*01:01 A*03:01/C*08:02 was found in the maternally transmitted haplotype, which had resulted in two novel HLA haplotypes in the proband. CONCLUSION: A rare case with simultaneous recombination of the paternal and maternal HLA-A/C loci has been discovered, which may facilitate further study of the mechanisms of the HLA recombination.

Genome-wide association study for cashmere traits in Inner Mongolia cashmere goat population reveals new candidate genes and haplotypes.

BACKGROUND: The cashmere goat industry is one of the main pillars of animal husbandry in Inner Mongolia Autonomous Region, and plays an irreplaceable role in local economic development. With the change in feeding methods and environment, the cashmere produced by Inner Mongolia cashmere goats shows a tendency of coarser, and the cashmere yield can not meet the consumption demand of people. However, the genetic basis behind these changes is not fully understood. We measured cashmere traits, including cashmere yield (CY), cashmere diameter (CD), cashmere thickness (CT), and fleece length (FL) traits for four consecutive years, and utilized Genome-wide association study of four cashmere traits in Inner Mongolia cashmere goats was carried out using new genomics tools to infer genomic regions and functional loci associated with cashmere traits and to construct haplotypes that significantly affect cashmere traits. RESULTS: We estimated the genetic parameters of cashmere traits in Inner Mongolia cashmere goats. The heritability of cashmere yield, cashmere diameter, and fleece length traits of Inner Mongolia cashmere goats were 0.229, 0.359, and 0.250, which belonged to the medium heritability traits (0.2 ~ 0.4). The cashmere thickness trait has a low heritability of 0.053. We detected 151 genome-wide significantly associated SNPs with four cashmere traits on different chromosomes, which were very close to the chromosomes of 392 genes (located within the gene or within ± 500 kb). Notch3, BMPR1B, and CCNA2 have direct functional associations with fibroblasts and follicle stem cells, which play important roles in hair follicle growth and development. Based on GO functional annotation and KEGG enrichment analysis, potential candidate genes were associated with pathways of hair follicle genesis and development (Notch, P13K-Akt, TGF-beta, Cell cycle, Wnt, MAPK). We calculated the effective allele number of the Inner Mongolia cashmere goat population to be 1.109-1.998, the dominant genotypes of most SNPs were wild-type, the polymorphic information content of 57 SNPs were low polymorphism (0 < PIC < 0.25), and the polymorphic information content of 79 SNPs were moderate polymorphism (0.25 < PIC < 0.50). We analyzed the association of SNPs with phenotypes and found that the homozygous mutant type of SNP1 and SNP3 was associated with the highest cashmere yield, the heterozygous mutant type of SNP30 was associated with the lowest cashmere thickness, the wild type of SNP76, SNP77, SNP78, SNP80, and SNP81 was associated with the highest cashmere thickness, and the wild type type of SNP137 was associated with the highest fleece length. 21 haplotype blocks and 68 haplotype combinations were constructed. Haplotypes A2A2, B2B2, C2C2, and D4D4 were associated with increased cashmere yield, haplotypes E2E2, F1F1, G5G5, and G1G5 were associated with decreased cashmere fineness, haplotypes H2H2 was associated with increased cashmere thickness, haplotypes I1I1, I1I2, J1J4, L5L3, N3N2, N3N3, O2O1, P2P2, and Q3Q3 were associated with increased cashmere length. We verified the polymorphism of 8 SNPs by KASP, and found that chr7_g.102631194A > G, chr10_g.82715068 T > C, chr1_g.124483769C > T, chr24_g.12811352C > T, chr6_g.114111249A > G, and chr6_g.115606026 T > C were significantly genotyped in verified populations (P < 0.05). CONCLUSIONS: In conclusion, the genetic effect of single SNP on phenotypes is small, and SNPs are more inclined to be inherited as a whole. By constructing haplotypes from SNPs that are significantly associated with cashmere traits, it will help to reveal the complex and potential causal variations in cashmere traits of Inner Mongolia cashmere goats. This will be a valuable resource for genomics and breeding of the cashmere goat.

Cold tolerance SNPs and candidate gene mining in the soybean germination stage based on genome-wide association analysis.

KEY MESSAGE: Three QTLs associated with low-temperature tolerance were identified by genome-wide association analysis, and 15 candidate genes were identified by haplotype analysis and gene expression analyses. Low temperature is a critical factor affecting the geographical distribution, growth, development, and yield of soybeans, with cold stress during seed germination leading to substantial productivity loss. In this study, an association panel comprising 260 soybean accessions was evaluated for four germination traits and four cold tolerance index traits, revealing extensive variation in cold tolerance. Genome-wide association study (GWAS) identified 10 quantitative trait nucleotides (QTNs) associated with cold tolerance, utilizing 30,799 single nucleotide polymorphisms (SNPs) and four GWAS models. Linkage disequilibrium (LD) analysis positioned these QTNs within three cold-tolerance quantitative trait loci (QTL) and, with QTL19-1, was positioned by three multi-locus models, underscoring its importance as a key QTL. Integrative haplotype analysis, supplemented by transcriptome analysis, uncovered 15 candidate genes. The haplotypes within the genes Glyma.18G044200, Glyma.18G044300, Glyma.18G044900, Glyma.18G045100, Glyma.19G222500, and Glyma.19G222600 exhibited significant phenotypic variations, with differential expression in materials with varying cold tolerance. The QTNs and candidate genes identified in this study offer substantial potential for marker-assisted selection and gene editing in breeding cold-tolerant soybeans, providing valuable insights into the genetic mechanisms underlying cold tolerance during soybean germination.

[Molecular tracing of Biomphalaria straminea in China].

OBJECTIVE: To investigate the origin of Biomphalaria straminea in China, so as to provide insights into assessment of schistosomiasis mansoni transmission risk and B. straminea control. METHODS: Guanlan River, Dasha River, Shenzhen Reservoir, upper and lower reaches of Kuiyong River, and Xinzhen River in Shenzhen, China, were selected as sampling sites. Ten Biomphalaria samples were collected from each site, and genomic DNA was extracted from Biomphalaria samples. DNA samples were obtained from 15 B. straminea sampled from 5 sampling sites in Minas Gerais State, Pará State, Federal District, Pernambuco State, and Sao Paulo State in Brazil, South America. Cytochrome c oxidase I (COI) and mitochondrial 16S ribosomal RNA (16S rRNA) genes were sampled using the above DNA templates, and the amplified products were sequenced. The COI and 16S rRNA gene sequences were downloaded from GenBank, and the sampling sites were acquired. All COI and 16S rRNA gene sequences were aligned and evolutionary trees of B. straminea were created based on COI and 16S rRNA gene sequences to identify the genetic similarity and evolutionary relationship between B. straminea samples from China and South America. RESULTS: A total of 60 COI gene sequences with a length of 529 bp and 3 haplotypes were obtained from B. straminea sampled from China. There were 165 COI gene sequences of B. straminea retrieved from GenBank, and following alignment with the above 60 gene sequences, a total of 33 haplotypes were obtained. Phylogenetic analysis showed that the three haplotypes of B. straminea from China were clustered into one clade, among which the haplotype China11 and three B. straminea samples from Brazil retrieved from GenBank belonged to the same haplotype. Geographical evolution analysis showed that the B. straminea samples from three sampling sites along eastern coasts of Brazil had the same haplotype with China11, and B. straminea samples from other two sampling sites were closely, genetically related to China11. A total of 60 16S rDNA gene sequences with approximately 322 bp in length were amplified from B. straminea in China, with 2 haplotypes identified. A total of 70 16S rDNA gene sequences of B. straminea were captured from GenBank. Phylogenetic analysis showed that Biomphalaria snails collected from China were clustered into a clade, and the haplotype China64 and the haplotype 229BS from Brazil shared the same haplotype. The 49 16S rDNA gene sequences of B. straminea from 25 sampling sites in southern Brazil, which were captured from GenBank, were included in the present analysis, and the B. straminea from 3 sampling sites shared the same haplotype with China64 in China. Geographical evolution analysis based on COI and 16S rRNA gene sequences showed that B. straminea sampled from eastern coastal areas of Brazil shared the same haplotypes in two gene fragment sequences with Biomphalaria snails collected from China. CONCLUSIONS: The Biomphalaria snails in China are characterized as B. straminea, which have a low genetic diversity. The Biomphalaria snails in China have a high genetic similarity with B. straminea sampled from eastern coastal areas of Brazil, which may have originated from the eastern coastal areas of Brazil.

Drug resistance markers in Plasmodium vivax isolates from a Kanchanaburi province, Thailand between January to May 2023.

BACKGROUND: Plasmodium vivax has become the predominant species in the border regions of Thailand. The emergence and spread of antimalarial drug resistance in P. vivax is one of the significant challenges for malaria control. Continuous surveillance of drug resistance is therefore necessary for monitoring the development of drug resistance in the region. This study aims to investigate the prevalence of the mutation in the P. vivax multidrug resistant 1 (Pvmdr1), dihydrofolate reductase (Pvdhfr), and dihydropteroate synthetase (Pvdhps) genes conferred resistance to chloroquine (CQ), pyrimethamine (P) and sulfadoxine (S), respectively. METHOD: 100 P. vivax isolates were obtained between January to May 2023 from a Kanchanaburi province, western Thailand. Nucleotide sequences of Pvmdr1, Pvdhfr, and Pvdhps genes were amplified and sequenced. The frequency of single nucleotide polymorphisms (SNPs)-haplotypes of drug-resistant alleles was assessed. The linkage disequilibrium (LD) tests were also analyzed. RESULTS: In Pvmdr1, T958M, Y976F, and F1076L, mutations were detected in 100%, 21%, and 23% of the isolates, respectively. In Pvdhfr, the quadruple mutant allele (I57R58M61T117) prevailed in 84% of the samples, followed by (L57R58M61T117) in 11%. For Pvdhps, the double mutant allele (G383G553) was detected (48%), followed by the triple mutant allele (G383M512G553) (47%) of the isolates. The most prevalent combination of Pvdhfr (I57R58M61T117) and Pvdhps (G383G553) alleles was sextuple mutated haplotypes (48%). For LD analysis, the association in the SNPs pairs was found between the intragenic and intergenic regions of the Pvdhfr and Pvdhps genes. CONCLUSION: The study has recently updated the high prevalence of three gene mutations associated with CQ and SP resistance. Genetic monitoring is therefore important to intensify in the regions to further assess the spread of drug resistant. Our data also provide evidence on the distribution of drug resistance for the early warning system, thereby threatening P. vivax malaria treatment policy decisions at the national level.

Detecting haplotype-specific transcript variation in long reads with FLAIR2.

BACKGROUND: RNA-seq has brought forth significant discoveries regarding aberrations in RNA processing, implicating these RNA variants in a variety of diseases. Aberrant splicing and single nucleotide variants (SNVs) in RNA have been demonstrated to alter transcript stability, localization, and function. In particular, the upregulation of ADAR, an enzyme that mediates adenosine-to-inosine editing, has been previously linked to an increase in the invasiveness of lung adenocarcinoma cells and associated with splicing regulation. Despite the functional importance of studying splicing and SNVs, the use of short-read RNA-seq has limited the community's ability to interrogate both forms of RNA variation simultaneously. RESULTS: We employ long-read sequencing technology to obtain full-length transcript sequences, elucidating cis-effects of variants on splicing changes at a single molecule level. We develop a computational workflow that augments FLAIR, a tool that calls isoform models expressed in long-read data, to integrate RNA variant calls with the associated isoforms that bear them. We generate nanopore data with high sequence accuracy from H1975 lung adenocarcinoma cells with and without knockdown of ADAR. We apply our workflow to identify key inosine isoform associations to help clarify the prominence of ADAR in tumorigenesis. CONCLUSIONS: Ultimately, we find that a long-read approach provides valuable insight toward characterizing the relationship between RNA variants and splicing patterns.

The first complete genome of the extinct European wild ass (Equus hemionus hydruntinus).

We present palaeogenomes of three morphologically unidentified Anatolian equids dating to the first millennium BCE, sequenced to a coverage of 0.6-6.4×. Mitochondrial DNA haplotypes of the Anatolian individuals clustered with those of Equus hydruntinus (or Equus hemionus hydruntinus), the extinct European wild ass, secular name 'hydruntine'. Further, the Anatolian wild ass whole genome profiles fell outside the genomic diversity of other extant and past Asiatic wild ass (E. hemionus) lineages. These observations suggest that the three Anatolian wild asses represent hydruntines, making them the latest recorded survivors of this lineage, about a millennium later than the latest observations in the zooarchaeological record. Our mitogenomic and genomic analyses indicate that E. h. hydruntinus was a clade belonging to ancient and present-day E. hemionus lineages that radiated possibly between 0.6 and 0.8 Mya. We also find evidence consistent with recent gene flow between hydruntines and Middle Eastern wild asses. Analyses of genome-wide heterozygosity and runs of homozygosity suggest that the Anatolian wild ass population may have lost genetic diversity by the mid-first millennium BCE, a possible sign of its eventual demise.

The association of rs17713054 with Neanderthal origin at 3p21.31 locus with the severity of COVID-19 in Iranian patients.

Since the COVID-19 pandemic, the diversity of clinical manifestations in patients has been a tremendous challenge. It seems that genetic variations, as one of the players, contribute to the variety of symptoms. Genome-wide association studies have demonstrated the influence of certain genomic regions on the disease prognosis. Particularly, a haplotype at 3p21.31 locus, inherited from Neanderthals, showed an association with COVID-19 severity. Despite several studies regarding this haplotype, some key variants are not sufficiently addressed. In the present study, we investigated the association of rs17713054 at 3p21.31 with COVID-19 severity. We analyzed the genotype of 251 Iranian COVID-19 patients (151 patients with asymptomatic to mild form as control and 100 patients with severe to critical symptoms without any comorbidities as case group) using the ARMS-PCR method. Results demonstrated that the A allele confers an almost twofold increased risk for COVID-19 severity (P value = 0.008). The AA genotype also raises the risk by more than 11 times following the recessive model (P value = 0.013). In conclusion, the A allele in rs17713054 was a risk allele in Iranian patients and was independently associated with COVID-19 severity. More studies are beneficial to confirm these findings in other populations and to develop strategies for risk assessment, prevention, and personalized medicine.

High prevalence of m.1555A > G in patients with hearing loss in the Baikal Lake region of Russia as a result of founder effect.

Mitochondrial forms account approximately 1-2% of all nonsyndromic cases of hearing loss (HL). One of the most common causative variants of mtDNA is the m.1555A > G variant of the MT-RNR1 gene (OMIM 561000). Currently the detection of the m.1555A > G variant of the MT-RNR1 gene is not included in all research protocols. In this study this variant was screened among 165 patients with HL from the Republic of Buryatia, located in the Baikal Lake region of Russia. In our study, the total contribution of the m.1555A > G variant to the etiology of HL was 12.7% (21/165), while the update global prevalence of this variant is 1.8% (863/47,328). The m.1555A > G variant was notably more prevalent in Buryat (20.2%) than in Russian patients (1.3%). Mitogenome analysis in 14 unrelated Buryat families carrying the m.1555A > G variant revealed a predominant lineage: in 13 families, a cluster affiliated with sub-haplogroup A5b (92.9%) was identified, while one family had the D5a2a1 lineage (7.1%). In a Russian family with the m.1555A > G variant the lineage affiliated with sub-haplogroup F1a1d was found. Considering that more than 90% of Buryat families with the m.1555A > G variant belong to the single maternal lineage cluster we conclude that high prevalence of this variant in patients with HL in the Baikal Lake region can be attributed to a founder effect.

Haplotype-resolved assembly of the mule duck genome using high-fidelity sequencing technology.

Mule duck is vitally important to the production of global duck meat. Here, we present two high-quality haplotypes of a female mule duck (haplotype 1 (H1):1.28 Gb, haplotype 2 (H2): 1.40 Gb). The continuity (H1: contig N50 = 14.90 Mb, H2: contig N50 = 15.70 Mb) and completeness (BUSCO: H1 = 96.9%, H2 = 97.3%) are substantially better than those of other duck genomes. We detected the structural variations (SVs) in H1 and H2. We observed a positive correlation between autosome length and the number of SVs. Z chromosome was some deficient in deletions and insertions, but W chromosome was some excessive. A total of 1,451 genes were haplotype specific expression (HSEs). Among them, 737 specifically expressed in H1, and 714 specifically expressed in H2. We found that H1 and H2 HSEs tended to be involved in similar biological processes, such as myometrial relaxation and contraction pathways, muscle structure development and phosphorylation. Our haplotype-resolved genome assembly provides a powerful platform for future functional genomics, molecular breeding, and genome editing in mule duck.

Concurrent of compound heterozygous variant of a novel in-frame deletion and the common hypomorphic haplotype in TBX6 and inherited 17q12 microdeletion in a fetus.

BACKGROUND: TBX6, a member of the T-box gene family, encodes the transcription factor box 6 that is critical for somite segmentation in vertebrates. It is known that the compound heterozygosity of disruptive variants in trans with a common hypomorphic risk haplotype (T-C-A) in the TBX6 gene contribute to 10% of congenital scoliosis (CS) cases. The deletion of chromosome 17q12 is a rare cytogenetic abnormality, which often leads to renal cysts and diabetes mellitus. However, the affected individuals often exhibit clinical heterogeneity and incomplete penetrance. METHODS: We here present a Chinese fetus who was shown to have CS by ultrasound examination at 17 weeks of gestation. Trio whole-exome sequencing (WES) was performed to investigate the underlying genetic defects of the fetus. In vitro functional experiments, including western-blotting and luciferase transactivation assay, were performed to determine the pathogenicity of the novel variant of TBX6. RESULTS: WES revealed the fetus harbored a compound heterozygous variant of c.338_340del (p.Ile113del) and the common hypomorphic risk haplotype of the TBX6 gene. In vitro functional study showed the p.Ile113del variant had no impact on TBX6 expression, but almost led to complete loss of its transcriptional activity. In addition, we identified a 1.85 Mb deletion on 17q12 region in the fetus and the mother. Though there is currently no clinical phenotype associated with this copy number variation in the fetus, it can explain multiple renal cysts in the pregnant woman. CONCLUSIONS: This study is the first to report a Chinese fetus with a single amino acid deletion variant and a T-C-A haplotype of TBX6. The clinical heterogeneity of 17q12 microdeletion poses significant challenges for prenatal genetic counseling. Our results once again suggest the complexity of prenatal genetic diagnosis.

Population structure and haplotype network analyses of Hyalomma anatolicum based on the large subunit ribosomal RNA (16S rRNA) gene.

Hyalomma anatolicum, an Anatolian hard tick is a well-recognized vector involved in the transmission of various pathogens to animals and humans. The present study elucidated the population structure and haplotype network of H. anatolicum based on the mitochondrial large subunit ribosomal RNA (16S rRNA) gene sequence. The population structure and haplotype network analysis of 75 sequences archived in the GenBank, including the 15 sequences generated herein, yielded 24 haplotypes. Haplotype 1 (Hap_1) was the predominant haplotype consisting of 45 sequences from India, China, Pakistan, Turkey, Egypt, Iraq, and Tajikistan. The complete haplotype network exhibited a stellate conformation, highlighting a recent population expansion. The overall dataset, together with the sequences corresponding to India, China, and Pakistan, showed a high haplotype (0.638 ± 0.065, 0.671 ± 0.103, 0.753 ± 0.099, and 0.854 ± 0.061, respectively) and low nucleotide (0.00407 ± 0.00090, 0.00525 ± 0.00196, 0.00680 ± 0.00233, and 0.00453 ± 0.00056, respectively) diversity, further emphasized a recent population expansion. The neutrality indices including Tajima's D, Fu and Li's D, and Fu and Li's F for the complete dataset (- 2.661, - 6.008, and - 5.649, respectively) as well as for the sequences from India (- 2.223, - 3.414, and - 3.567, respectively) were negative, suggesting deviation from neutrality and a recent population expansion. The present study provided novel insights into the population structure and haplotype networks of H. anatolicum based on the mitochondrial 16S rRNA gene, and the different tests inferred a low genetic differentiation and suggested a recent population expansion of this economically important tick species.

Population structure and selection signal analysis of indigenous sheep from the southern edge of the Taklamakan Desert.

Analyzing the genetic diversity and selection characteristics of sheep (Ovis aries) holds significant value in understanding their environmental adaptability, enhancing breeding efficiency, and achieving effective conservation and rational utilization of genetic resources. In this study, we utilized Illumina Ovine SNP 50 K BeadChip data from four indigenous sheep breeds from the southern margin of the Taklamakan Desert (Duolang sheep: n = 36, Hetian sheep: n = 74, Kunlun sheep: n = 27, Qira black sheep: n = 178) and three foreign meat sheep breeds (Poll Dorset sheep: n = 105, Suffolk sheep: n = 153, Texel sheep: n = 150) to investigate the population structure, genetic diversity, and genomic signals of positive selection within the indigenous sheep. According to the Principal component analysis (PCA), the Neighbor-Joining tree (NJ tree), and Admixture, we revealed distinct clustering patterns of these seven sheep breeds based on their geographical distribution. Then used Cross Population Extended Haplotype Homozygosity (XP-EHH), Fixation Index (FST), and Integrated Haplotype Score (iHS), we identified a collective set of 32 overlapping genes under positive selection across four indigenous sheep breeds. These genes are associated with wool follicle development and wool traits, desert environmental adaptability, disease resistance, reproduction, and high-altitude adaptability. This study reveals the population structure and genomic selection characteristics in the extreme desert environments of native sheep breeds from the southern edge of the Taklimakan Desert, providing new insights into the conservation and sustainable use of indigenous sheep genetic resources in extreme environments. Additionally, these findings offer valuable genetic resources for sheep and other mammals to adapt to global climate change.

Association between OX40L polymorphism and type 2 diabetes mellitus in Iranians.

INTRODUCTION: Diabetes mellitus (DM) is one of the leading causes of morbidity and mortality worldwide. It is a multifactorial disease that genetic and environmental factors contribute to its development. The aim of the study was to investigate the association of OX40L promoter gene polymorphisms with type 2 diabetes mellitus (T2DM) in Iranians. MATERIALS AND METHODS: Three hundred and sixty-eight subjects including 184 healthy subjects and 184 T2DM patients were enrolled in our study. Polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) was applied to detect genotype and allele frequencies of rs3850641, rs1234313 and rs10912580. In addition, SNPStats web tool was applied to estimate haplotype frequency and linkage disequilibrium (LD). RESULTS: The distribution of tested polymorphisms was statistically different between the T2DM patients and healthy subjects (P < 0.01). rs1234313 AG (OR = 0.375, 95% CI = 0.193-0.727, P = 0.004) and rs10912580 AG (OR = 0.351, 95% CI = 0.162-0.758, P = 0.008) genotypes were associated with the decreased risk of T2DM in Iranians. Moreover, our prediction revealed that AAG (OR = 0.46, 95% CI= (0.28-0.76), P = 0.0028) and GAG (OR = 0.24, 95% CI= (0.13-0.45), P < 0.0001) haplotypes were related to the reduced risk of the disease. However, the tested polymorphisms had no effect on biochemical parameters and body mass index (BMI) in the patient group (P > 0.05). CONCLUSION: Our findings revealed that OX40L promoter gene polymorphisms are associated with T2DM. Moreover, genotype and allelic variations were related to the decreased risk of T2DM in Iranians. Further studies are recommended to show whether these polymorphic variations could affect OX40/OX40L interaction or OX40L phenotype.

Two decades of molecular surveillance in Senegal reveal rapid changes in known drug resistance mutations over time.

BACKGROUND: Drug resistance in Plasmodium falciparum is a major threat to malaria control efforts. Pathogen genomic surveillance could be invaluable for monitoring current and emerging parasite drug resistance. METHODS: Data from two decades (2000-2020) of continuous molecular surveillance of P. falciparum parasites from Senegal were retrospectively examined to assess historical changes in malaria drug resistance mutations. Several known drug resistance markers and their surrounding haplotypes were profiled using a combination of single nucleotide polymorphism (SNP) molecular surveillance and whole genome sequence based population genomics. RESULTS: This dataset was used to track temporal changes in drug resistance markers whose timing correspond to historically significant events such as the withdrawal of chloroquine (CQ) and the introduction of sulfadoxine-pyrimethamine (SP) in 2003. Changes in the mutation frequency at Pfcrt K76T and Pfdhps A437G coinciding with the 2014 introduction of seasonal malaria chemoprevention (SMC) in Senegal were observed. In 2014, the frequency of Pfcrt K76T increased while the frequency of Pfdhps A437G declined. Haplotype-based analyses of Pfcrt K76T showed that this rapid increase was due to a recent selective sweep that started after 2014. DISCUSSION (CONCLUSION): The rapid increase in Pfcrt K76T is troubling and could be a sign of emerging amodiaquine (AQ) resistance in Senegal. Emerging AQ resistance may threaten the future clinical efficacy of artesunate-amodiaquine (ASAQ) and AQ-dependent SMC chemoprevention. These results highlight the potential of molecular surveillance for detecting rapid changes in parasite populations and stress the need to monitor the effectiveness of AQ as a partner drug for artemisinin-based combination therapy (ACT) and for chemoprevention.

Climate change will lead to range shifts and genetic diversity losses of dung beetles in the Gobi Desert and Mongolian Steppe.

Desertification is known to be a major threat to biodiversity, yet our understanding of the consequent decline in biodiversity remains insufficient. Here, we predicted climate change-induced range shifts and genetic diversity losses in three model dung beetles: Colobopterus erraticus, Cheironitis eumenes, and Gymnopleurus mopsus, distributed across the Gobi Desert and Mongolian Steppe, areas known for desertification. Phylogeographic analyses of mitochondrial COI sequences and species distribution modeling, based on extensive field investigations spanning 14 years, were performed. Species confined to a single biome were predicted to contract and shift their distribution in response to climate change, whereas widespread species was predicted to expand even if affected by range shifts. We indicated that all species are expected to experience significant haplotype losses, yet the presence of high singleton frequencies and low genetic divergence across geographic configurations and lineages mitigate loss of genetic diversity. Notably, Cheironitis eumenes, a desert species with low genetic diversity, appears to be the most vulnerable to climate change due to the extensive degradation in the Gobi Desert. This is the first study to predict the response of insects to desertification in the Gobi Desert. Our findings highlight that dung beetles in the Gobi Desert and Mongolian Steppe might experience high rates of occupancy turnover and genetic loss, which could reshuffle the species composition.

Impact of patrilocality on contrasting patterns of paternal and maternal heritage in Central-West Africa.

Despite their ancient past and high diversity, African populations are the least represented in human population genetic studies. In this study, uniparental markers (mtDNA and Y chromosome) were used to investigate the impact of sociocultural factors on the genetic diversity and inter-ethnolinguistic gene flow in the three major Nigerian groups: Hausa (n = 89), Yoruba (n = 135) and Igbo (n = 134). The results show a distinct history from the maternal and paternal perspectives. The three Nigerian groups present a similar substrate for mtDNA, but not for the Y chromosome. The two Niger-Congo groups, Yoruba and Igbo, are paternally genetically correlated with populations from the same ethnolinguistic affiliation. Meanwhile, the Hausa is paternally closer to other Afro-Asiatic populations and presented a high diversity of lineages from across Africa. When expanding the analyses to other African populations, it is observed that language did not act as a major barrier to female-mediated gene flow and that the differentiation of paternal lineages is better correlated with linguistic than geographic distances. The results obtained demonstrate the impact of patrilocality, a common and well-established practice in populations from Central-West Africa, in the preservation of the patrilineage gene pool and in the affirmation of identity between groups.

Natural variation of domestication-related genes contributed to latitudinal expansion and adaptation in soybean.

Soybean is a major source of protein and edible oil worldwide. Originating from the Huang-Huai-Hai region, which has a temperate climate, soybean has adapted to a wide latitudinal gradient across China. However, the genetic mechanisms responsible for the widespread latitudinal adaptation in soybean, as well as the genetic basis, adaptive differentiation, and evolutionary implications of theses natural alleles, are currently lacking in comprehensive understanding. In this study, we examined the genetic variations of fourteen major gene loci controlling flowering and maturity in 103 wild species, 1048 landraces, and 1747 cultivated species. We found that E1, E3, FT2a, J, Tof11, Tof16, and Tof18 were favoured during soybean improvement and selection, which explained 75.5% of the flowering time phenotypic variation. These genetic variation was significantly associated with differences in latitude via the LFMM algorithm. Haplotype network and geographic distribution analysis suggested that gene combinations were associated with flowering time diversity contributed to the expansion of soybean, with more HapA clustering together when soybean moved to latitudes beyond 35°N. The geographical evolution model was developed to accurately predict the suitable planting zone for soybean varieties. Collectively, by integrating knowledge from genomics and haplotype classification, it was revealed that distinct gene combinations improve the adaptation of cultivated soybeans to different latitudes. This study provides insight into the genetic basis underlying the environmental adaptation of soybean accessions, which could contribute to a better understanding of the domestication history of soybean and facilitate soybean climate-smart molecular breeding for various environments.