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.

[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.

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.

[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.

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.

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.

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.

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.

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.

Haplotype-based association mapping of genomic regions associated with Zymoseptoria tritici resistance using 217 diverse wheat genotypes.

BACKGROUND: Septoria tritici blotch (STB) is considered to be one of the most destructive foliar wheat diseases and is caused by Zymoseptoria tritici. The yield losses are severe and in Northwestern Europe can reach up to 50%. The efficacy of fungicides is diminishing due to changes in the genetic structure of the pathogen. Therefore, resistance breeding is the most effective strategy of disease management. Recently, genome-wide association studies (GWAS) have become more popular due to their robustness in dissecting complex traits, including STB resistance in wheat. This was made possible by the use of large mapping populations and new sequencing technologies. High-resolution mapping benefits from historical recombination and greater allele numbers in GWAS. RESULTS: In our study, 217 wheat genotypes of diverse origin were phenotyped against five Z. tritici isolates (IPO323, IPO88004, IPO92004, IPO86036 and St1-03) and genotyped on the DArTseq platform. In polytunnel tests two disease parameters were evaluated: the percentage of leaf area covered by necrotic lesions (NEC) and the percentage of leaf area covered by lesions bearing pycnidia (PYC). The disease escape parameters heading date (Hd) and plant height (Ht) were also measured. Pearson's correlation showed a positive effect between disease parameters, providing additional information. The Structure analysis indicated four subpopulations which included from 28 (subpopulation 2) to 79 genotypes (subpopulation 3). All of the subpopulations showed a relatively high degree of admixture, which ranged from 60% of genotypes with less than 80% of proportions of the genome attributed to assigned subpopulation for group 2 to 85% for group 4. Haplotype-based GWAS analysis allowed us to identify 27 haploblocks (HBs) significantly associated with analysed traits with a p-value above the genome-wide significance threshold (5%, which was -log10(p) > 3.64) and spread across the wheat genome. The explained phenotypic variation of identified significant HBs ranged from 0.2% to 21.5%. The results of the analysis showed that four haplotypes (HTs) associated with disease parameters cause a reduction in the level of leaf coverage by necrosis and pycnidia, namely: Chr3A_HB98_HT2, Chr5B_HB47_HT1, Chr7B_HB36_HT1 and Chr5D_HB10_HT3. CONCLUSIONS: GWAS analysis enabled us to identify four significant chromosomal regions associated with a reduction in STB disease parameters. The list of valuable HBs and wheat varieties possessing them provides promising material for further molecular analysis of resistance loci and development of breeding programmes.

Population genetic structure and phylogenetic analysis of Anopheles hyrcanus (Diptera: Culicidae) inferred from DNA sequences of nuclear ITS2 and the mitochondrial COI gene in the northern part of Iran.

BACKGROUND: The Anopheles hyrcanus group is distributed throughout the Oriental and Palaearctic regions and can transmit diseases such as malaria, Japanese encephalitis virus, and filariasis. This investigation marks the inaugural comprehensive study to undertake a phylogenetic analysis of the constituents of this malaria vector group in the northeastern region of Iran, juxtaposed with documented occurrences from different areas within Iran and worldwide. METHODS: Mosquitoes were collected using various methods from nine different locations in Golestan province from April to December 2023. The collected mosquitoes were identified morphologically using valid taxonomic keys. DNA was isolated using the Sambio™ Kit. COI and ITS2 primers were designed using Oligo7 and GeneRunner. PCR and purification were performed with the Qiagen kit. Subsequently, sequencing was carried out at the Mehr Mam GENE Center using an Applied Biosystems 3730XL sequencer. The nucleotide sequences were then analyzed and aligned with GenBank data using BioEdit. Kimura 2-parameter was Utilized for base substitutions. DNA models were selected based on AIC and BIC criteria. Bayesian and Maximum Likelihood trees were constructed, along with a haplotype network. Molecular diversity statistics computed using DnaSP software. RESULTS: In this study, a total of 819 adult mosquitoes were collected. An. hyrcanus was the second most abundant species, predominantly found in Kalaleh and Turkman counties. The sequenced and edited COI and ITS2 sequences were deposited in GenBank under specific accession numbers. Phylogenetic analyses using ML, BI, and NJ methods confirmed a monophyletic lineage for An. hyrcanus with strong support. Molecular analysis of Iranian An. hyrcanus found 11 diverse haplotypes, with the COI gene displaying low diversity. The ITS2 gene revealed two clades - one associating with Iran, Europe, and Asia; the other originating from southwestern Iran. The haplotype network showed two main groups - one from southwest Iran and the other from north Iran. Iran exhibited six distinct haplotypes, while Turkey showcased the highest diversity. CONCLUSIONS: An. hyrcanus in southwestern Iran exhibits a distinct haplogroup, suggesting possible subspecies differentiation. Additional studies are required to validate this phenomenon.

Association between CYP1A2 gene variants -163 C/A (rs762551) and -3860 G/A (rs2069514) and bladder cancer susceptibility.

BACKGROUND: Bladder cancer (BLCA) poses a significant global health challenge due to its high incidence, poor prognosis, and limited treatment options. AIMS AND OBJECTIVES: This study aims to investigate the association between two specific polymorphisms, CYP1A2-163 C/A and CYP1A2-3860G/A, within the Cytochrome P450 1A2 (CYP1A2) gene and susceptibility to BLCA. METHODS: The study employed a case-control design, genotyping 340 individuals using Polymerase Chain Reaction-High-Resolution Melting Curve (PCR-HRM). Various genetic models were applied to evaluate allele and genotype frequencies. Genetic linkage analysis was facilitated using R packages. RESULTS: The study reveals a significant association with the - 163 C/A allele, particularly in the additive model. Odds ratio (OR) analysis links CYP1A2-163 C/A (rs762551) and CYP1A2-3860G/A(rs2069514) polymorphisms to BLCA susceptibility. The rs762551 C/A genotype is prevalent in 55% of BLCA cases and exhibits an OR of 2.21. The A/A genotype has an OR of 1.54. Regarding CYP1A2-3860G/A, the G/A genotype has an OR of 1.54, and the A/A genotype has an OR of 2.08. Haplotype analysis shows a predominant C-C haplotype at 38.2%, followed by a C-A haplotype at 54.7%, and a less frequent A-A haplotype at 7.1%. This study underscores associations between CYP1A2 gene variants, particularly rs762551 (CYP1A2-163 C/A), and an increased susceptibility to BLCA. Haplotype analysis of 340 individuals reveals a predominant C-C haplotype at 38.2%, followed by a C-A haplotype at 54.7%, and a less frequent A-A haplotype at 7.1%. CONCLUSION: In conclusion, the - 163 C/A allele, C/A genotype of rs762551, and G/A genotype of rs2069514 emerge as potential genetic markers associated with elevated BLCA risk.

Optimisation of global stem cell donor recruitment based on analysis of unsuccessful donor searches.

Despite over 41 million registered potential volunteer stem cell donors worldwide, many patients in need of a transplant do not find an HLA-matched unrelated donor or cord blood units, with the respective odds differing significantly between various populations. In this study, we analysed data of 2205 unsuccessful real-life donor searches sent to the DKMS Registry to identify populations in which further donor recruitment would be associated with particularly large patient benefits. For that purpose, we estimated haplotype frequencies of 67 donor populations at various sample sizes and entered them into two different mathematical models. These models assessed patient benefits from population-specific donor recruitment, operationalised by the number of originally unsuccessful searches that may become successful due to new donors. Consistently, across the different mathematical models and sample sizes, we obtained several countries from East and Southeast Asia (Thailand, Vietnam, China, and the Philippines) and the population of Asians in the USA as countries/populations where donor recruitment activities would be particularly beneficial for patients. We also identified various countries in Southeast and Central Europe as possible target regions for donor recruitment with above-average patient benefits. The results presented are registry-specific in the sense that they were obtained by optimising unsuccessful searches that had been sent to the DKMS Registry. Therefore, it would be desirable to apply the presented methods to a global data set that includes all unsuccessful stem cell donor searches worldwide and uses population-specific haplotype frequencies based on all donors available in the WMDA Search & Match Service.

HLA-DQA1*01:03 and DQB1*06:01 are risk factors for severe COVID-19 pneumonia.

The clinical spectrum of COVID-19 includes a wide range of manifestations, from mild symptoms to severe pneumonia. HLA system plays a pivotal role in immune responses to infectious diseases. The purpose of our study was to investigate the association between HLA and COVID-19 severity in a Japanese population. The study included 209 Japanese COVID-19 patients aged ≥20 years. Saliva samples were collected and used to determine the HLA genotype by HLA imputation through genome-wide association analyses. The association between HLA genotype and COVID-19 severity was then evaluated. The allele frequency was compared between patients with respiratory failure (severe group: 91 cases) and those without respiratory failure (non-severe group: 118 cases), categorising the data into three time periods: pre-Omicron epidemic period, Omicron epidemic period, and total period of this study (from January 2021 to May 2023). In comparing the severe and non-severe groups, the frequencies of the HLA-DQA1*01:03 (35.1% vs. 10.5%, odds ratio [OR] = 4.57, corrected p [pc] = 0.041) and -DQB1*06:01 (32.4% vs. 7.9%, OR = 5.54, pc = 0.030) alleles were significantly higher in the severe group during the pre-Omicron epidemic period. During the Omicron epidemic period, HLA-DQB1*06 (32.4% vs. 7.9%, OR = 5.54, pc = 0.030) was significantly higher in the severe group. During total period of this study, HLA-DQA1*01:03 (30.2% vs. 14.4%, OR = 2.57, corrected pc = 0.0013) and -DQB1*06:01 (44.5% vs. 26.7%, OR = 2.20, pc = 0.013) alleles were significantly higher in the severe group. HLA-DQB1*06:01 and -DQA1*01:03 were in strong linkage disequilibrium with each other (r2 = 0.91) during total period of this study, indicating that these two alleles form a haplotype. The frequency of the HLA-DQA1*01:03-DQB1*06:01 in the severe group was significantly higher than in the non-severe group during pre-Omicron epidemic period (32.4% vs. 7.9%, OR = 5.59, pc = 0.00072), and total period of this study (28.6% vs. 13.1%, OR = 2.63, pc = 0.0013). During Omicron epidemic period, the haplotype did not demonstrate statistical significance, although the odds ratio indicated a value greater 1. Frequencies of the HLA-DQA1*01:03 and -DQB1*06:01 alleles were significantly higher in severe COVID-19 patients, suggesting that these alleles are risk factors for severe COVID-19 pneumonia in the Japanese population.

Integrating wing morphometrics and mitochondrial DNA analysis to assess the filaria vector Mansonia uniformis (Diptera: Culicidae) populations in Thailand.

Mansonia uniformis (Diptera: Culicidae) is recognized as a vector of Brugia malayi and has been reported to transmit Wuchereria bancrofti, both causing lymphatic filariasis in humans. This study employed geometric morphometrics (GM) to investigate wing shape variation and analyzed genetic diversity through cytochrome c oxidase subunit 1 (COI) gene analyses in Ma. uniformis populations across Thailand. Wing GM analyses indicated significant differences in wing shape based on Mahalanobis distances among nearly all population pairs (p < 0.05), with no significant correlation between wing shape and geographic distance (r = 0.210, p > 0.05). Genetic analyses identified 63 haplotypes and 49 polymorphic sites, with the overall population exhibiting a nucleotide diversity of 0.006 (± 0.001) and a haplotype diversity of 0.912 (± 0.017). Deviations from neutrality, as indicated by Tajima's D and Fu's FS tests for the overall Ma. uniformis populations in Thailand, were statistically significant and negative, suggesting population expansion (both p < 0.05). Analysis of molecular variance revealed no significant genetic structure when all populations were categorized based on collection sites and geographic regions. However, significant differences in FST values were observed between some populations. These findings enhance our understanding of the geographical and genetic factors influencing Ma. uniformis populations, which are crucial for developing effective control strategies in Thailand.

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.

A systematic strategy for identifying causal single nucleotide polymorphisms and their target genes on Juvenile arthritis risk haplotypes.

BACKGROUND: Although genome-wide association studies (GWAS) have identified multiple regions conferring genetic risk for juvenile idiopathic arthritis (JIA), we are still faced with the task of identifying the single nucleotide polymorphisms (SNPs) on the disease haplotypes that exert the biological effects that confer risk. Until we identify the risk-driving variants, identifying the genes influenced by these variants, and therefore translating genetic information to improved clinical care, will remain an insurmountable task. We used a function-based approach for identifying causal variant candidates and the target genes on JIA risk haplotypes. METHODS: We used a massively parallel reporter assay (MPRA) in myeloid K562 cells to query the effects of 5,226 SNPs in non-coding regions on JIA risk haplotypes for their ability to alter gene expression when compared to the common allele. The assay relies on 180 bp oligonucleotide reporters ("oligos") in which the allele of interest is flanked by its cognate genomic sequence. Barcodes were added randomly by PCR to each oligo to achieve > 20 barcodes per oligo to provide a quantitative read-out of gene expression for each allele. Assays were performed in both unstimulated K562 cells and cells stimulated overnight with interferon gamma (IFNg). As proof of concept, we then used CRISPRi to demonstrate the feasibility of identifying the genes regulated by enhancers harboring expression-altering SNPs. RESULTS: We identified 553 expression-altering SNPs in unstimulated K562 cells and an additional 490 in cells stimulated with IFNg. We further filtered the SNPs to identify those plausibly situated within functional chromatin, using open chromatin and H3K27ac ChIPseq peaks in unstimulated cells and open chromatin plus H3K4me1 in stimulated cells. These procedures yielded 42 unique SNPs (total = 84) for each set. Using CRISPRi, we demonstrated that enhancers harboring MPRA-screened variants in the TRAF1 and LNPEP/ERAP2 loci regulated multiple genes, suggesting complex influences of disease-driving variants. CONCLUSION: Using MPRA and CRISPRi, JIA risk haplotypes can be queried to identify plausible candidates for disease-driving variants. Once these candidate variants are identified, target genes can be identified using CRISPRi informed by the 3D chromatin structures that encompass the risk haplotypes.

Candidate Genes and Favorable Haplotypes Associated with Iron Toxicity Tolerance in Rice.

Iron (Fe) toxicity is a major issue adversely affecting rice production worldwide. Unfortunately, the physiological and genetic mechanisms underlying Fe toxicity tolerance in rice remain relatively unknown. In this study, we conducted a genome-wide association study using a diverse panel consisting of 551 rice accessions to identify genetic mechanisms and candidate genes associated with Fe toxicity tolerance. Of the 29 quantitative trait loci (QTL) for Fe toxicity tolerance detected on chromosomes 1, 2, 5, and 12, five (qSH_Fe5, qSFW_Fe2.3, qRRL5.1, qRSFW1.1, and qRSFW12) were selected to identify candidate genes according to haplotype and bioinformatics analyses. The following five genes were revealed as promising candidates: LOC_Os05g40160, LOC_Os05g40180, LOC_Os12g36890, LOC_Os12g36900, and LOC_Os12g36940. The physiological characteristics of rice accessions with contrasting Fe toxicity tolerance reflected the importance of reactive oxygen species-scavenging antioxidant enzymes and Fe homeostasis for mitigating the negative effects of Fe toxicity on rice. Our findings have clarified the genetic and physiological mechanisms underlying Fe toxicity tolerance in rice. Furthermore, we identified valuable genetic resources for future functional analyses and the development of Fe toxicity-tolerant rice varieties via marker-assisted selection.

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.