Far-East Asian Toxoplasma isolates share ancestry with North and South/Central American recombinant lineages.

Toxoplasma gondii is a global protozoan pathogen. Clonal lineages predominate in Europe, North America, Africa, and China, whereas highly recombinant parasites are endemic in South/Central America. Far East Asian T. gondii isolates are not included in current global population genetic structure analyses at WGS resolution. Here we report a genome-wide population study that compared eight Japanese and two Chinese isolates against representative worldwide T. gondii genomes using POPSICLE, a novel population structure analyzing software. Also included were 7 genomes resurrected from non-viable isolates by target enrichment sequencing. Visualization of the genome structure by POPSICLE shows a mixture of Chinese haplogroup (HG) 13 haploblocks introgressed within the genomes of Japanese HG2 and North American HG12. Furthermore, two ancestral lineages were identified in the Japanese strains; one lineage shares a common ancestor with HG11 found in both Japanese strains and North American HG12. The other ancestral lineage, found in T. gondii isolates from a small island in Japan, is admixed with genetically diversified South/Central American strains. Taken together, this study suggests multiple ancestral links between Far East Asian and American T. gondii strains and provides insight into the transmission history of this cosmopolitan organism.

The monsoon-associated equine South African pointy mosquito 'Aedes caballus'; the first comprehensive record from southeastern Iran with a description of ecological, morphological, and molecular aspects.

The equine South African pointy vector mosquito, Aedes caballus, poses a significant threat to human health due to its capacity for transmitting arboviruses. Despite favorable climate for its existence in southeast Iran, previous records of this species in the area have indicated very low abundance. This comprehensive field and laboratory study aimed to assess its current adult population status in this region, utilizing a combination of ecological, morphological and molecular techniques. Four distinct types of traps were strategically placed in three fixed and two variable mosquito sampling sites in the southern strip of Sistan and Baluchistan Province. Subsequently, DNA was extracted from trapped mosquitoes and subjected to PCR amplification using the molecular markers COI, ITS2, and ANT. In total, 1734 adult Ae. caballus specimens were collected from rural areas, with the majority being captured by CO2-baited bednet traps. A notable increase in the abundance of this species was observed following rainfall in February. The genetic analysis revealed multiple haplotypes based on COI and ITS2 sequences, with COI gene divergence at 0.89%, and ITS2 sequence divergence at 1.6%. This suggests that previous challenges in morphological identification may have led to misidentifications, with many adults previously classified as Ae. vexans potentially being Ae. caballus. The findings of this study hold significant implications for public health authorities, providing valuable insights for integrated and targeted vector control and disease management efforts.

Defining a Haplotype Encompassing the LCORL-NCAPG Locus Associated with Increased Lean Growth in Beef Cattle.

Numerous studies have shown genetic variation at the LCORL-NCAPG locus is strongly associated with growth traits in beef cattle. However, a causative molecular variant has yet to be identified. To define all possible candidate variants, 34 Charolais-sired calves were whole-genome sequenced, including 17 homozygous for a long-range haplotype associated with increased growth (QQ) and 17 homozygous for potential ancestral haplotypes for this region (qq). The Q haplotype was refined to an 814 kb region between chr6:37,199,897-38,014,080 and contained 218 variants not found in qq individuals. These variants include an insertion in an intron of NCAPG, a previously documented mutation in NCAPG (rs109570900), two coding sequence mutations in LCORL (rs109696064 and rs384548488), and 15 variants located within ATAC peaks that were predicted to affect transcription factor binding. Notably, rs384548488 is a frameshift variant likely resulting in loss of function for long isoforms of LCORL. To test the association of the coding sequence variants of LCORL with phenotype, 405 cattle from five populations were genotyped. The two variants were in complete linkage disequilibrium. Statistical analysis of the three populations that contained QQ animals revealed significant (p < 0.05) associations with genotype and birth weight, live weight, carcass weight, hip height, and average daily gain. These findings affirm the link between this locus and growth in beef cattle and describe DNA variants that define the haplotype. However, further studies will be required to define the true causative mutation.

Is CYP2C Haplotype Relevant for Efficacy and Bleeding Risk in Clopidogrel-Treated Patients?

A recently discovered haplotype-CYP2C:TG-determines the ultrarapid metabolism of several CYP2C19 substrates. The platelet inhibitor clopidogrel requires CYP2C19-mediated activation: the risk of ischemic events is increased in patients with a poor (PM) or intermediate (IM) CYP2C19 metabolizer phenotype (vs. normal, NM; rapid, RM; or ultrarapid, UM). We investigated whether the CYP2C:TG haplotype affected efficacy/bleeding risk in clopidogrel-treated patients. Adults (n = 283) treated with clopidogrel over 3-6 months were classified by CYP2C19 phenotype based on the CYP2C19*2*17 genotype, and based on the CYP2C19/CYP2C cluster genotype, and regarding carriage of the CYP2:TG haplotype, and were balanced on a number of covariates across the levels of phenotypes/haplotype carriage. Overall, 45 (15.9%) patients experienced ischemic events, and 49 (17.3%) experienced bleedings. By either classification, the incidence of ischemic events was similarly numerically higher in PM/IM patients (21.6%, 21.8%, respectively) than in mutually similar NM, RM, and UM patients (13.2-14.8%), whereas the incidence of bleeding events was numerically lower (13.1% vs. 16.6-20.5%). The incidence of ischemic events was similar in CYP2C:TG carries and non-carries (14.1% vs. 16.1%), whereas the incidence of bleedings appeared mildly lower in the former (14.9% vs. 20.1%). We observed no signal to suggest a major effect of the CYP2C19/CYP2C cluster genotype or CYP2C:TG haplotype on the clinical efficacy/safety of clopidogrel.

The Landscape of Presence/Absence Variations during the Improvement of Rice.

Rice is one of the most important staple crops in the world; therefore, the improvement of rice holds great significance for enhancing agricultural production and addressing food security challenges. Although there have been numerous studies on the role of single-nucleotide polymorphisms (SNPs) in rice improvement with the development of next-generation sequencing technologies, research on the role of presence/absence variations (PAVs) in the improvement of rice is limited. In particular, there is a scarcity of studies exploring the traits and genes that may be affected by PAVs in rice. Here, we extracted PAVs utilizing resequencing data from 148 improved rice varieties distributed in Asia. We detected a total of 33,220 PAVs and found that the number of variations decreased gradually as the length of the PAVs increased. The number of PAVs was the highest on chromosome 1. Furthermore, we identified a 6 Mb hotspot region on chromosome 11 containing 1091 PAVs in which there were 29 genes related to defense responses. By conducting a genome-wide association study (GWAS) using PAV variation data and phenotypic data for five traits (flowering time, plant height, flag leaf length, flag leaf width, and panicle number) across all materials, we identified 186 significantly associated PAVs involving 20 cloned genes. A haplotype analysis and expression analysis of candidate genes revealed that important genes might be affected by PAVs, such as the flowering time gene OsSFL1 and the flag leaf width gene NAL1. Our work investigated the pattern in PAVs and explored important PAV key functional genes associated with agronomic traits. Consequently, these results provide potential and exploitable genetic resources for rice breeding.

From Caves to the Savannah, the Mitogenome History of Modern Lions (Panthera leo) and Their Ancestors.

Lions (Panthera leo) play a crucial ecological role in shaping and maintaining fragile ecosystems within Africa. Conservation efforts should focus on genetic variability within wild populations when considering reintroduction attempts. We studied two groups of lions from two conservation sites located in Zambia and Zimbabwe to determine their genetic make-up, information that is usually unknown to the sites. In this study, we analysed 17 specimens for cytb and seven microsatellite markers to ascertain family relationships and genetic diversity previously obtained by observational studies. We then produced a standardised haplogroup phylogeny using all available entire mitogenomes, as well as calculating a revised molecular clock. The modern lion lineage diverged ~151 kya and was divided into two subspecies, both containing three distinct haplogroups. We confirm that Panthera leo persica is not a subspecies, but rather a haplogroup of the northern P.l. leo that exited Africa at least ~31 kya. The progenitor to all lions existed ~1.2 Mya, possibly in SE Africa, and later exited Africa and split into the two cave lion lineages ~175 kya. Species demography is correlated to major climactic events. We now have a detailed phylogeny of lion evolution and an idea of their conservation status given the threat of climate change.

No Association of Polymorphisms in the Genes Encoding Interleukin-6 and Interleukin-6 Receptor Subunit Alpha with the Risk of Keloids in Polish Patients.

A keloid is a benign fibroproliferative hypertrophy of scar tissue that extends outside the original wound and invades adjacent healthy skin. Keloid formation is thought to be a complex process including overactivity of the interleukin-6 signaling pathway and genetic susceptibility. The aim of the study was to investigate possible associations between rs1800797, rs1800796, and rs1800795 polymorphisms in the promoter of the IL6 gene encoding interleukin-6 and the rs2228145 polymorphism in the IL6R gene encoding the interleukin-6 receptor subunit alpha with the predisposition to keloids in Polish patients. The genetic polymorphisms were identified either using Polymerase Chain Reaction-Restriction Fragment Length Polymorphism (PCR-RFLP) or sequencing of samples of genomic DNA extracted from blood leukocytes of 86 adult patients with keloids and 100 newborns comprising a control group. No significant differences in the distributions of IL6 or IL6R alleles or genotypes were found between keloid patients and newborn controls. There were also no significant differences between both groups in the distribution of IL6 haplotypes. The IL6 rs1800797, rs1800796 and rs1800795 and IL6R rs2228145 polymorphisms were not found to predispose individuals in the study group to keloids. IL6 promoter haplotypes were not found to be associated with a higher risk of keloids in the studied group.

Genotyping Hepatitis B virus by Next-Generation Sequencing: Detection of Mixed Infections and Analysis of Sequence Conservation.

Our aim was to develop an accurate, highly sensitive method for HBV genotype determination and detection of genotype mixtures. We examined the preS and 5' end of the HBV X gene (5X) regions of the HBV genome using next-generation sequencing (NGS). The 1852 haplotypes obtained were subjected to genotyping via the Distance-Based discrimination method (DB Rule) using two sets of 95 reference sequences of genotypes A-H. In clinical samples from 125 patients, the main genotypes were A, D, F and H in Caucasian, B and C in Asian and A and E in Sub-Saharan patients. Genotype mixtures were identified in 28 (22.40%) cases, and potential intergenotypic recombination was observed in 29 (23.20%) cases. Furthermore, we evaluated sequence conservation among haplotypes classified into genotypes A, C, D, and E by computing the information content. The preS haplotypes exhibited limited shared conserved regions, whereas the 5X haplotypes revealed two groups of conserved regions across the genotypes assessed. In conclusion, we developed an NGS-based HBV genotyping method utilizing the DB Rule for genotype classification. We identified two regions conserved across different genotypes at 5X, offering promising targets for RNA interference-based antiviral therapies.

Pangenome characterization and analysis of the NAC gene family reveals genes for Sclerotinia sclerotiorum resistance in sunflower (Helianthus annuus).

BACKGROUND: Sunflower (Helianthus annuus) is one of the most important economic crops in oilseed production worldwide. The different cultivars exhibit variability in their resistance genes. The NAC transcription factor (TF) family plays diverse roles in plant development and stress responses. With the completion of the H. annuus genome sequence, the entire complement of genes coding for NACs has been identified. However, the reference genome of a single individual cannot cover all the genetic information of the species. RESULTS: Considering only a single reference genome to study gene families will miss many meaningful genes. A pangenome-wide survey and characterization of the NAC genes in sunflower species were conducted. In total, 139 HaNAC genes are identified, of which 114 are core and 25 are variable. Phylogenetic analysis of sunflower NAC proteins categorizes these proteins into 16 subgroups. 138 HaNACs are randomly distributed on 17 chromosomes. SNP-based haplotype analysis shows haplotype diversity of the HaNAC genes in wild accessions is richer than in landraces and modern cultivars. Ten HaNAC genes in the basal stalk rot (BSR) resistance quantitative trait loci (QTL) are found. A total of 26 HaNAC genes are differentially expressed in response to Sclerotinia head rot (SHR). A total of 137 HaNAC genes are annotated in Gene Ontology (GO) and are classified into 24 functional groups. GO functional enrichment analysis reveals that HaNAC genes are involved in various functions of the biological process. CONCLUSIONS: We identified NAC genes in H. annuus (HaNAC) on a pangenome-wide scale and analyzed S. sclerotiorum resistance-related NACs. This study provided a theoretical basis for further genomic improvement targeting resistance-related NAC genes in sunflowers.

The effects of IL-27 and IL-35 gene variation and expression levels on the susceptibility and clinical manifestations of pulmonary tuberculosis.

Inflammatory cytokines have crucial roles in the pathogenesis of tuberculosis (TB), and interleukin (IL)-27 and IL-35 have a pro-inflammatory and anti-inflammatory effect on many diseases, including infectious diseases. Therefore, we evaluated the relationship between IL-27 and IL-35 gene polymorphism, expression levels, and pulmonary TB (PTB) susceptibility. Nine single-nucleotide polymorphisms (SNPs) in the IL-27 gene (rs181206, rs153109, and rs17855750) and the IL-35 gene (rs4740, rs428253, rs9807813, rs2243123, rs2243135, and rs568408) were genotyped by the SNPscan technique in 497 patients with PTB and 501 controls. There was no significant difference regarding the genotype and allele frequencies of the above SNPs in the IL-27 and IL-35 genes between patients with PTB and controls. Haplotype analysis showed that the frequency of the GAC haplotype in the IL-35 gene was significantly decreased in patients with PTB when compared to controls (p = 0.036). Stratified analysis suggested that the frequency of the IL-27 rs17855750 GG genotype was significantly increased in patients with PTB with fever. Moreover, the lower frequency of the IL-35 rs568408 GA genotype was associated with drug-induced liver injury in patients with PTB. The IL-35 rs428253 GC genotype, as well as the rs4740 AA genotype and A allele, showed significant relationships with hypoproteinemia in patients with PTB. When compared with controls, the IL-27 level was significantly increased in patients with PTB. Taken together, IL-35 gene variation might contribute to a protective role on the susceptibility to PTB, and IL-27 and IL-35 gene polymorphisms were associated with several clinical manifestations of patients with PTB.

Association of tuberculosis risk with genetic polymorphisms of the immune checkpoint genes PDCD1, CTLA-4, and TIM3.

The immune checkpoint proteins were reported to involve to host resistance to Mycobacteria tuberculosis (Mtb). Here, we evaluated 11 single nucleotide polymorphisms (SNPs) in PDCD1, CTLA4, and HAVCR2 genes between participants with and without TB infection. Genomic DNA isolated from 285 patients with TB and 270 controls without TB infection were used to perform the genotyping assay. Odds ratios were used to characterize the association of 11 SNPs with TB risk. In this study, the various genotypes of the 11 SNPs did not differ significantly in frequency between the non-TB and TB groups. When patients were stratified by sex, however, men differed significantly from women in genotype frequencies at HAVCR2 rs13170556. Odds ratios indicated that rs2227982, rs13170556, rs231775, and rs231779 were sex-specifically associated with TB risk. In addition, the combinations of rs2227982/rs13170556 GA/TC in men and the A-C-C haplotype of rs231775-rs231777-rs231779 in women were significantly associated with TB risk. Our results indicate that rs2227982 in PDCD1 and rs13170556 in HAVCR2 are associated with increased TB susceptibility in men and that the CTLA4 haplotype appears protective against TB in women.

Multiple polygenic risk scores can improve the prediction of systemic lupus erythematosus in Taiwan.

OBJECTIVE: To identify new genetic variants associated with SLE in Taiwan and establish polygenic risk score (PRS) models to improve the early diagnostic accuracy of SLE. METHODS: The study enrolled 2429 patients with SLE and 48 580 controls from China Medical University Hospital in Taiwan. A genome-wide association study (GWAS) and PRS analyses of SLE and other three SLE markers, namely ANA, anti-double-stranded DNA antibody (dsDNA) and anti-Smith antibody (Sm), were conducted. RESULTS: Genetic variants associated with SLE were identified through GWAS. Some novel genes, which have been previously reported, such as RCC1L and EGLN3, were revealed to be associated with SLE in Taiwan. Multiple PRS models were established, and optimal cut-off points for each PRS were determined using the Youden Index. Combining the PRSs for SLE, ANA, dsDNA and Sm yielded an area under the curve of 0.64 for the optimal cut-off points. An analysis of human leucocyte antigen (HLA) haplotypes in SLE indicated that individuals with HLA-DQA1*01:01 and HLA-DQB1*05:01 were at a higher risk of being classified into the SLE group. CONCLUSIONS: The use of PRSs to predict SLE enables the identification of high-risk patients before abnormal laboratory data were obtained or symptoms were manifested. Our findings underscore the potential of using PRSs and GWAS in identifying SLE markers, offering promise for early diagnosis and prediction of SLE.

Lytic and Latent Genetic Diversity of the Epstein-Barr Virus Reveals Raji-Related Variants from Southeastern Brazil Associated with Recombination Markers.

Epstein-Barr virus (EBV) is a ubiquitous gammaherpesvirus etiologically associated with benign and malignant diseases. Since the pathogenic mechanisms of EBV are not fully understood, understanding EBV genetic diversity is an ongoing goal. Therefore, the present work describes the genetic diversity of the lytic gene BZLF1 in a sampling of 70 EBV-positive cases from southeastern Brazil. Additionally, together with the genetic regions previously characterized, the aim of the present study was to determine the impact of viral genetic factors that may influence EBV genetic diversity. Accordingly, the phylogenetic analysis of the BZLF1 indicated two main clades with high support, BZ-A and BZ-B (PP > 0.85). Thus, the BZ-A clade was the most diverse clade associated with the main polymorphisms investigated, including the haplotype Type 1 + V3 (p < 0.001). Furthermore, the multigene phylogenetic analysis (MLA) between BZLF1 and the oncogene LMP1 showed specific clusters, revealing haplotypic segregation that previous single-gene phylogenies from both genes failed to demonstrate. Surprisingly, the LMP1 Raji-related variant clusters were shown to be more diverse, associated with BZ-A/B and the Type 2/1 + V3 haplotypes. Finally, due to the high haplotypic diversity of the Raji-related variants, the number of DNA recombination-inducing motifs (DRIMs) was evaluated within the different clusters defined by the MLA. Similarly, the haplotype BZ-A + Raji was shown to harbor a greater number of DRIMs (p < 0.001). These results call attention to the high haplotype diversity of EBV in southeast Brazil and strengthen the hypothesis of the recombinant potential of South American Raji-related variants via the LMP1 oncogene.

Differential Activation of TAS2R4 May Recover Ability to Taste Propylthiouracil for Some TAS2R38 AVI Homozygotes.

Bitterness from phenylthiocarbamide and 6-n-propylthiouracil (PROP) varies with polymorphisms in the TAS2R38 gene. Three SNPs form two common (AVI, PAV) and four rare haplotypes (AAI, AAV, PVI, and PAI). AVI homozygotes exhibit higher detection thresholds and lower suprathreshold bitterness for PROP compared to PAV homozygotes and heterozygotes, and these differences may influence alcohol and vegetable intake. Within a diplotype, substantial variation in suprathreshold bitterness persists, and some AVI homozygotes report moderate bitterness at high concentrations. A second receptor encoded by a gene containing a functional polymorphism may explain this. Early work has suggested that PROP might activate TAS2R4 in vitro, but later work did not replicate this. Here, we identify three TAS2R4 SNPs that result in three diplotypes-SLN/SLN, FVS/SLN, and FVS/FVS-which make up 25.1%, 44.9%, and 23.9% of our sample. These TAS2R4 haplotypes show minimal linkage disequilibrium with TAS2R38, so we examined the suprathreshold bitterness as a function of both. The participants (n = 243) rated five PROP concentrations in duplicate, interleaved with other stimuli. As expected, the TAS2R38 haplotypes explained ~29% (p < 0.0001) of the variation in the bitterness ratings, with substantial variation within the haplotypes (AVI/AVI, PAV/AVI, and PAV/PAV). Notably, the TAS2R4 diplotypes (independent of the TAS2R38 haplotypes) explained ~7-8% of the variation in the bitterness ratings (p = 0.0001). Given this, we revisited if PROP could activate heterologously expressed TAS2R4 in HEK293T cells, and calcium imaging indicated 3 mM PROP is a weak TAS2R4 agonist. In sum, our data are consistent with the second receptor hypothesis and may explain the recovery of the PROP tasting phenotype in some AVI homozygotes; further, this finding may potentially help explain the conflicting results on the TAS2R38 diplotype and food intake.

Virus Quasispecies Rarefaction: Subsampling with or without Replacement?

In quasispecies diversity studies, the comparison of two samples of varying sizes is a common necessity. However, the sensitivity of certain diversity indices to sample size variations poses a challenge. To address this issue, rarefaction emerges as a crucial tool, serving to normalize and create fairly comparable samples. This study emphasizes the imperative nature of sample size normalization in quasispecies diversity studies using next-generation sequencing (NGS) data. We present a thorough examination of resampling schemes using various simple hypothetical cases of quasispecies showing different quasispecies structures in the sense of haplotype genomic composition, offering a comprehensive understanding of their implications in general cases. Despite the big numbers implied in this sort of study, often involving coverages exceeding 100,000 reads per sample and amplicon, the rarefaction process for normalization should be performed with repeated resampling without replacement, especially when rare haplotypes constitute a significant fraction of interest. However, it is noteworthy that different diversity indices exhibit distinct sensitivities to sample size. Consequently, some diversity indicators may be compared directly without normalization, or instead may be resampled safely with replacement.

Haplotype-resolved genome assembly and implementation of VitExpress, an open interactive transcriptomic platform for grapevine.

Haplotype-resolved genome assemblies were produced for Chasselas and Ugni Blanc, two heterozygous Vitis vinifera cultivars by combining high-fidelity long-read sequencing and high-throughput chromosome conformation capture (Hi-C). The telomere-to-telomere full coverage of the chromosomes allowed us to assemble separately the two haplo-genomes of both cultivars and revealed structural variations between the two haplotypes of a given cultivar. The deletions/insertions, inversions, translocations, and duplications provide insight into the evolutionary history and parental relationship among grape varieties. Integration of de novo single long-read sequencing of full-length transcript isoforms (Iso-Seq) yielded a highly improved genome annotation. Given its higher contiguity, and the robustness of the IsoSeq-based annotation, the Chasselas assembly meets the standard to become the annotated reference genome for V. vinifera. Building on these resources, we developed VitExpress, an open interactive transcriptomic platform, that provides a genome browser and integrated web tools for expression profiling, and a set of statistical tools (StatTools) for the identification of highly correlated genes. Implementation of the correlation finder tool for MybA1, a major regulator of the anthocyanin pathway, identified candidate genes associated with anthocyanin metabolism, whose expression patterns were experimentally validated as discriminating between black and white grapes. These resources and innovative tools for mining genome-related data are anticipated to foster advances in several areas of grapevine research.

Population genetic structure of Pomacea canaliculata in China based on the COI and ITS1 genes.

Comprehending the phylogeography of invasive organisms enhances our insight into their distribution dynamics, which is instrumental for the development of effective prevention and management strategies. In China, Pomacea canaliculata and Pomacea maculata are the two most widespread and damaging species of the non-native Pomacea spp.. Given this species' rapid spread throughout country, it is urgent to investigate the genetic diversity and structure of its different geographic populations, a task undertaken in the current study using the COI and ITS1 mitochondrial and ribosomal DNA genes, respectively. The result of this study, based on a nationwide systematic survey, a collection of Pomacea spp., and the identification of cryptic species, showed that there is a degree of genetic diversity and differentiation in P. canaliculata, and that all of its variations are mainly due to differences between individuals within different geographical populations. Indeed, this species contains multiple haplotypes, but none of them form a systematic geographical population structure. Furthermore, the COI gene exhibits higher genetic diversity than the ITS1 gene. Our study further clarifies the invasive pathways and dispersal patterns of P. canaliculata in China to provide a theoretical basis.

Genetic diversity and haplotypes of Cysticercus tenuicollis isolates from slaughtered sheep and goats in Elazig and Bingol provinces of Turkey.

BACKGROUND: The cestode Taenia hydatigena uses canids, primarily dogs, as definitive hosts, whereas the metacestode larval stage cysticercus infects a range of intermediate hosts, including domestic animals such as goats, sheep and pigs. OBJECTIVES: In this study, we aimed to determine the genetic differences and haplotypes of Cysticercus tenuicollis isolates obtained from sheep and goats slaughtered in the Bingol and Elazig provinces of Turkey. METHODS: C. tenuicollis isolates were collected from 44 sheep and 26 goats slaughtered in slaughterhouses in Bingol and Elazig. After the isolation of total genomic DNA from C. tenuicollis isolates, the genetic characterization of the partial mitochondrial cytochrome oxidase 1 (CO1) gene region (866 bp) was amplified using specific primers by polymerase chain reaction (PCR), the products were then sequenced, and haplotype and genetic diversity analyses were carried out. RESULTS: As a result of the haplotype network analyses, 34 different haplotypes were detected around the main haplotype (Hap02) arranged in a star-like configuration and separated from other haplotypes by 1-28 mutation steps and covering 22.85% (16/70) of all isolates. Twenty-seven polymorphic fields were detected, 77.77% (21/27) of which were parsimony-informative, and secondary haplotype and nucleotide diversity were observed. Additionally, we detected high intraspecies haplotype diversity (hd: 0.933) and high nucleotide diversity (π: 0.00383), with 27 different nucleotide variation positions among the haplotypes of the isolates. Tajima's D value was negative, indicating population expansion and/or selection purification. The significantly negative Fu's Fs values indicated recent population expansion or the presence of expected rare haplotypes. CONCLUSION: The results of this study confirmed that C. tenuicollis isolates clustered in one lineage and were closely related to the relevant reference sequences in different countries, confirming the circulation of C. tenuicollis in different geographical regions.

Identification and map-based cloning of an EMS-induced mutation in wheat gene TaSP1 related to spike architecture.

KEY MESSAGE: A candidate gene TaSP1 related to spike shape was cloned, and the gene-specific marker was developed to efficiently track the superior haplotype in common wheat. Spike shape, an important factor that affects wheat grain yield, is mainly defined by spike length (SPL), spikelet number (SPN), and compactness. Zhoumai32 mutant 1160 (ZM1160), a mutant obtained from ethyl methane sulfonate (EMS) treatment of hexaploid wheat variety Zhoumai32, was used to identify and clone the candidate gene that conditioned the spike shape. Genetic analysis of an F2 population derived from a cross of ZM1160 and Bainong207 suggested that the compact spike shape in ZM1160 was controlled by a single recessive gene, and therefore, the mutated gene was designated as Tasp1. With polymorphic markers identified through bulked segregant analysis (BSA), the gene was mapped to a 2.65-cM interval flanked by markers YZU0852 and MIS46239 on chromosome 7D, corresponding to a 0.42-Mb physical interval of Chinese spring (CS) reference sequences (RefSeq v1.0). To fine map TaSP1, 15 and seven recombinants were, respectively, screened from 1599 and 1903 F3 plants derived from the heterozygous F2 plants. Finally, TaSP1 was delimited to a 21.9 Kb (4,870,562 to 4,892,493 bp) Xmis48123-Xmis48104 interval. Only one high-confidence gene TraesCS7D02G010200 was annotated in this region, which encodes an unknown protein with a putative vWA domain. Quantitative reverse transcription PCR (qRT-PCR) analysis showed that TraesCS7D02G010200 was mainly expressed in the spike. Haplotype analysis of 655 wheat cultivars using the candidate gene-specific marker Xg010200p2 identified a superior haplotype TaSP1b with longer spike and more spikelet number. TaSP1 is beneficial to the improvement in wheat spike shape.

A fully haplotype-resolved and nearly gap-free genome assembly of wheat stripe rust fungus.

Stripe rust fungus Puccinia striiformis f. sp. tritici (Pst) is a destructive pathogen of wheat worldwide. Pst has a macrocyclic-heteroecious lifecycle, in which one-celled urediniospores are dikaryotic, each nucleus containing one haploid genome. We successfully generated the first fully haplotype-resolved and nearly gap-free chromosome-scale genome assembly of Pst by combining PacBio HiFi sequencing and trio-binning strategy. The genome size of the two haploid assemblies was 75.59 Mb and 75.91 Mb with contig N50 of 4.17 Mb and 4.60 Mb, and both had 18 pseudochromosomes. The high consensus quality values of 55.57 and 59.02 for both haplotypes confirmed the correctness of the assembly. Of the total 18 chromosomes, 15 and 16 were gapless while there were only five and two gaps for the remaining chromosomes of the two haplotypes, respectively. In total, 15,046 and 15,050 protein-coding genes were predicted for the two haplotypes, and the complete BUSCO scores achieved 97.7% and 97.9%, respectively. The genome will lay the foundation for further research on genetic variations and the evolution of rust fungi.