Impact of the inaccessible genome on genotype imputation and genome-wide association studies.

Genotype imputation is widely used in genome-wide association studies (GWAS). However, both the genotyping chips and imputation reference panels are dependent on next-generation sequencing (NGS). Due to the nature of NGS, some regions of the genome are inaccessible to sequencing. To date, there has been no complete evaluation of these regions and their impact on the identification of associations in GWAS remains unclear. In this study, we systematically assess the extent to which variants in inaccessible regions are underrepresented on genotyping chips and imputation reference panels, in GWAS results and in variant databases. We also determine the proportion of genes located in inaccessible regions and compare the results across variant masks defined by the 1000 Genomes Project and the TOPMed program. Overall, fewer variants were observed in inaccessible regions in all categories analyzed. Depending on the mask used and normalized for region size, only 4%-17% of the genotyped variants are located in inaccessible regions and 52 to 581 genes were almost completely inaccessible. From the Cooperative Health Research in South Tyrol (CHRIS) study, we present a case study of an association located in an inaccessible region that is driven by genotyped variants and cannot be reproduced by imputation in GRCh37. We conclude that genotyping, NGS, genotype imputation and downstream analyses such as GWAS and fine mapping are systematically biased in inaccessible regions, due to missed variants and spurious associations. To help researchers assess gene and variant accessibility, we provide an online application (https://gab.gm.eurac.edu).

Deqformer: high-definition and scalable deep learning probe design method.

Target enrichment sequencing techniques are gaining widespread use in the field of genomics, prized for their economic efficiency and swift processing times. However, their success depends on the performance of probes and the evenness of sequencing depth among each probe. To accurately predict probe coverage depth, a model called Deqformer is proposed in this study. Deqformer utilizes the oligonucleotides sequence of each probe, drawing inspiration from Watson-Crick base pairing and incorporating two BERT encoders to capture the underlying information from the forward and reverse probe strands, respectively. The encoded data are combined with a feed-forward network to make precise predictions of sequencing depth. The performance of Deqformer is evaluated on four different datasets: SNP panel with 38 200 probes, lncRNA panel with 2000 probes, synthetic panel with 5899 probes and HD-Marker panel for Yesso scallop with 11 000 probes. The SNP and synthetic panels achieve impressive factor 3 of accuracy (F3acc) of 96.24% and 99.66% in 5-fold cross-validation. F3acc rates of over 87.33% and 72.56% are obtained when training on the SNP panel and evaluating performance on the lncRNA and HD-Marker datasets, respectively. Our analysis reveals that Deqformer effectively captures hybridization patterns, making it robust for accurate predictions in various scenarios. Deqformer leads to a novel perspective for probe design pipeline, aiming to enhance efficiency and effectiveness in probe design tasks.

A HIV Superinfection Diagnosed in a Patient on Intramuscular Long-acting Combination of Cabotegravir and Rilpivirine.

A case of a male with human immunodeficiency virus with plasma genotyping detecting no resistance and a CRF02_AG subtype had a controlled HIV RNA on antiretroviral therapy since 2010. We introduced intramuscular therapy with cabotegravir and rilpivirine. One month later, his HIV RNA was 1500 copies/mL; genotyping found a subtype B with many mutations.

RHC genotyping in Chinese Han population.

BACKGROUND: The Rh blood group system is characterized by its complexity and polymorphism, encompassing 56 different antigens. Accurately predicting the presence of the C antigen using genotyping methods has been challenging. The objective of this study was to evaluate the accuracy of various genotyping methods for predicting the Rh C and to identify a suitable method for the Chinese Han population. METHODS: In total, 317 donors, consisting 223 D+ (including 20 with the Del phenotype) and 94 D- were randomly selected. For RHC genotyping, 48C and 109bp insertion were detected on the Real-time PCR platform and -292 substitution was analyzed via restriction fragment length polymorphism (RFLP). Moreover, the promoter region of the RHCE gene was sequenced to search for other nucleotide substitutions between RHC and RHc. Agreement between prediction methods was evaluated using the Kappa statistic, and comparisons between methods were conducted via the χ2 test. RESULTS: The analysis revealed that the 48C allele, 109bp insertion, a specific pattern observed in RFLP results, and wild-type alleles of seven single nucleotide polymorphisms (SNPs) were in strong agreement with the Rh C, with Kappa coefficients exceeding 0.8. However, there were instances of false positives or false negatives (0.6% false negative rate for 109bp insertion and 5.4-8.2% false positive rates for other methods). The 109bp insertion method exhibited the highest accuracy in predicting the Rh C, at 99.4%, compared to other methods (P values≤0.001). Although no statistical differences were found among other methods for predicting Rh C (P values>0.05), the accuracies in descending order were 48C (94.6%) > rs586178 (92.7%) > rs4649082, rs2375313, rs2281179, rs2072933, rs2072932, and RFLP (92.4%) > rs2072931 (91.8%). CONCLUSIONS: None of the methods examined can independently and accurately predict the Rh C. However, the 109bp insertion test demonstrated the highest accuracy for predicting the Rh C in the Chinese Han population. Utilizing the 109bp insertion test in combination with other methods may enhance the accuracy of Rh C prediction.

High-Resolution Genotyping of Formalin-Fixed Tissue Accurately Estimates Polygenic Risk Scores in Human Diseases.

Formalin-fixed paraffin-embedded (FFPE) tissues stored in biobanks and pathology archives are a vast but underutilized source for molecular studies on different diseases. Beyond being the "gold standard" for preservation of diagnostic human tissues, FFPE samples retain similar genetic information as matching blood samples, which could make FFPE samples an ideal resource for genomic analysis. However, research on this resource has been hindered by the perception that DNA extracted from FFPE samples is of poor quality. Here, we show that germline disease-predisposing variants and polygenic risk scores (PRS) can be identified from FFPE normal tissue (FFPE-NT) DNA with high accuracy. We optimized the performance of FFPE-NT DNA on a genome-wide array containing 657,675 variants. Via a series of testing and validation phases, we established a protocol for FFPE-NT genotyping with results comparable with blood genotyping. The median call rate of FFPE-NT samples in the validation phase was 99.85% (range 98.26%-99.94%) and median concordance with matching blood samples was 99.79% (range 98.85%-99.9%). We also demonstrated that a rare pathogenic PALB2 genetic variant predisposing to cancer can be correctly identified in FFPE-NT samples. We further imputed the FFPE-NT genotype data and calculated the FFPE-NT genome-wide PRS in 3 diseases and 4 disease risk variables. In all cases, FFPE-NT and matching blood PRS were highly concordant (all Pearson's r > 0.95). The ability to precisely genotype FFPE-NT on a genome-wide array enables translational genomics applications of archived FFPE-NT samples with the possibility to link to corresponding phenotypes and longitudinal health data.

Molecular Confirmation of Anopheles stephensi Mosquitoes in the Al Hudaydah Governorate, Yemen, 2021 and 2022.

We detected malaria vector Anopheles stephensi mosquitoes in the Al Hudaydah governorate in Yemen by using DNA sequencing. We report 2 cytochrome c oxidase subunit I haplotypes, 1 previously found in Ethiopia, Somalia, Djibouti, and Yemen. These findings provide insight into invasive An. stephensi mosquitoes in Yemen and their connection to East Africa.

Emergence of Bluetongue Virus Serotype 3, the Netherlands, September 2023.

Since 1998, notifiable bluetongue virus (BTV) serotypes 1-4, 6, 8, 9, 11, and 16 have been reported in Europe. In August 2006, a bluetongue (BT) outbreak caused by BTV serotype 8 began in northwestern Europe. The Netherlands was declared BT-free in February 2012, and annual monitoring continued. On September 3, 2023, typical BT clinical manifestations in sheep were notified to the Netherlands Food and Product Safety Consumer Authority. On September 6, we confirmed BTV infection through laboratory diagnosis; notifications of clinical signs in cattle were also reported. We determined the virus was serotype 3 by whole-genome sequencing. Retrospective analysis did not reveal BTV circulation earlier than September. The virus source and introduction route into the Netherlands remains unknown. Continuous monitoring and molecular diagnostic testing of livestock will be needed to determine virus spread, and new prevention strategies will be required to prevent BTV circulation within the Netherlands and Europe.

Long-term prediction by DNA methylation of high-grade cervical intraepithelial neoplasia: Results of the ARTISTIC cohort.

Methylation markers have shown potential for triaging high-risk HPV-positive (hrHPV+) women to identify those at increased risk of invasive cervical cancer (ICC). Our aim was to assess the performance of the S5 DNA methylation classifier for predicting incident high-grade cervical intraepithelial neoplasia (CIN) and ICC among hrHPV+ women in the ARTISTIC screening trial cohort. The S5 classifier, comprising target regions of tumour suppressor gene EPB41L3 and L1 and L2 regions of HPV16, HPV18, HPV31, and HPV33, was assayed by pyrosequencing in archived hrHPV+ liquid-based samples from 343 women with high-grade disease (139 CIN2, 186 CIN3, and 18 ICC) compared to 800 hrHPV+ controls. S5 DNA methylation correlated directly with increasing severity of disease and inversely with lead time to diagnosis. S5 could discriminate between hrHPV+ women who developed CIN3 or ICC and hrHPV+ controls (p <.0001) using samples taken on average 5 years before diagnosis. This relationship was independent of cytology at baseline. The S5 test showed much higher sensitivity than HPV16/18 genotyping for identifying prevalent CIN3 (93% vs. 61%, p = .01) but lower specificity (50% vs. 66%, p <.0001). The S5 classifier identified most women at high risk of developing precancer and missed very few prevalent advanced lesions thus appearing to be an objective test for triage of hrHPV+ women. The combination of methylation of host and HPV genes enables S5 to combine the predictive power of methylation with HPV genotyping to identify hrHPV-positive women who are at highest risk of developing CIN3 and ICC in the future.

Clinical performance of the novel full-genotyping OncoPredict HPV Quantitative Typing assay using the VALGENT framework.

Clinical validation of human papillomavirus (HPV) assays according to international criteria is prerequisite for their implementation in cervical cancer screening. OncoPredict HPV Quantitative Typing (QT) assay (Hiantis Srl, Milan, Italy) is a novel full-genotyping multiplex real-time PCR quantitative assay targeting E6/E7 genes, allowing individual viral load determination of 12 high-risk (HR) HPV types. Quality controls for sample adequacy, efficiency of nucleic acid extraction and PCR inhibition are included in the assay. Clinical performance of OncoPredict HPV QT test was assessed as part of the "Validation of HPV Genotyping Tests" (VALGENT-2) framework, consisting of 1300 cervical liquid-based cytology (LBC) samples of women aged between 20 and 60 years who had originally attended for routine cervical screening in Scotland. The clinical accuracy of the OncoPredict HPV QT (index test) for the detection of CIN2+ was assessed relative to the GP5+/6+ Enzyme ImmunoAssay (GP5+/6+ EIA) (comparator test), using noninferiority criteria. Intra- and interlaboratory reproducibility of the assay was assessed on a subpopulation, comprising 526 samples. The relative sensitivity and specificity for OncoPredict HPV QT vs GP5+/6+-PCR-EIA were 1.01 (95% CI: 0.99-1.03) and 1.03 (95% CI: 1.0-1.06) respectively. The P-values for noninferiority were ≤0.001. The intra- and inter-laboratory reproducibility demonstrated a high concordance (>98.7%) with kappas for individual types ranging from 0.66 to 1.00. OncoPredict HPV QT fulfills the international validation criteria for the use of HPV tests in cervical cancer screening.

Application of the TaqMan ARMS-PCR Approach for Genotyping Drug-Induced Hearing Loss Using Dried Blood Samples.

A single nucleotide variant in mitochondrial DNA (mtDNA) 1555A>G is associated with drug-induced hearing loss. For the 1555A>G mutation site, 1555A wild-type and 1555G mutant-type plasmids were constructed, respectively. In this study, a PCR method based on the TaqMan amplification refractory mutation system was proposed to detect mtDNA 1555A>G. A common upstream primer, a common TaqMan probe, and two downstream allele-specific primers with mismatched bases were designed. One-step amplification and detection of the wild-type and mutant type at the 1555 site were realized for the deafness-related gene through two reactions. Based on this detection method, the minimum detection limit of the wild-type and mutant type detection systems for plasmids was 50 copies/µL. The minimum sensitivity for the detection of nucleic acids in real dried blood spot (DBS) samples was 0.1 ng/µL. In the normal DBS DNA sample, the detection limit of the mutation abundance reached 0.78%. The specificity of the detection method was 100%, and the coefficient of variation was less than 3.36%. This approach was validated using clinical DNA extracted from 113 DBS samples of newborns. Additionally, it showed 100% agreement with bi-directional Sanger sequencing. It can be used as an optional method for the clinical detection of deafness-related genes.

A simple and effective genotyping workflow for rapid detection of CRISPR genome editing.

Genetically engineered mouse models play a pivotal role in the modeling of diseases, exploration of gene functions, and the development of novel therapies. In recent years, clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPR-associated protein 9 (Cas9)-mediated genome editing technology has revolutionized the process of developing such models by enabling precise genome modifications of the multiple interested genes simultaneously. Following genome editing, an efficient genotyping methodology is crucial for subsequent characterization. However, current genotyping methods are laborious, time-consuming, and costly. Here, using targeting the mouse trypsinogen genes as an example, we introduced common applications of CRISPR-Cas9 editing and a streamlined cost-effective genotyping workflow for CRISPR-edited mouse models, in which Sanger sequencing is required only at the initial steps. In the F0 mice, we focused on identifying the presence of positive editing by PCR followed by Sanger sequencing without the need to know the exact sequences, simplifying the initial screening. In the F1 mice, Sanger sequencing and algorithms decoding were used to identify the precise editing. Once the edited sequence was established, a simple and effective genotyping strategy was established to distinguish homozygous and heterozygous status by PCR from tail DNA. The genotyping workflow applies to deletions as small as one nucleotide, multiple-gene knockout, and knockin studies. This simplified, efficient, and cost-effective genotyping shall be instructive to new investigators who are unfamiliar with characterizing CRISPR-Cas9-edited mouse strains.NEW & NOTEWORTHY This study presents a streamlined, cost-effective genotyping workflow for clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPR-associated protein 9 (Cas9) edited mouse models, focusing on trypsinogen genes. It simplifies initial F0 mouse screening using PCR and Sanger sequencing without needing exact sequences. For F1 mice, precise editing is identified through Sanger sequencing and algorithm decoding. The workflow includes a novel PCR strategy for distinguishing homozygous and heterozygous statuses in subsequent generations, effective for small deletions, multiple-gene knockouts, and knockins.

Blood group genotype matching for transfusion.

The last decade has seen significant growth in the application of DNA-based methods for extended antigen typing, and the use of gene sequencing to consider variation in blood group genes to guide clinical care. The challenge for the field now lies in educating professionals, expanding accessibility and standardizing the use of genotyping for routine patient care. Here we discuss applications of genotyping when transfusion is not straightforward including when compatibility cannot be demonstrated by routine methods, when Rh type is unclear, when allo- and auto-antibodies are encountered in stem cell and organ transplantation, for prenatal testing to determine maternal and foetal risk for complications, and Group A subtyping for kidney and platelet donors. We summarize current commercial testing resources and new approaches to testing including high-density arrays and targeted next-generation sequencing (NGS).

Effect of genotyping errors on linkage map construction based on repeated chip analysis of two recombinant inbred line populations in wheat (Triticum aestivum L.).

Linkage maps are essential for genetic mapping of phenotypic traits, gene map-based cloning, and marker-assisted selection in breeding applications. Construction of a high-quality saturated map requires high-quality genotypic data on a large number of molecular markers. Errors in genotyping cannot be completely avoided, no matter what platform is used. When genotyping error reaches a threshold level, it will seriously affect the accuracy of the constructed map and the reliability of consequent genetic studies. In this study, repeated genotyping of two recombinant inbred line (RIL) populations derived from crosses Yangxiaomai × Zhongyou 9507 and Jingshuang 16 × Bainong 64 was used to investigate the effect of genotyping errors on linkage map construction. Inconsistent data points between the two replications were regarded as genotyping errors, which were classified into three types. Genotyping errors were treated as missing values, and therefore the non-erroneous data set was generated. Firstly, linkage maps were constructed using the two replicates as well as the non-erroneous data set. Secondly, error correction methods implemented in software packages QTL IciMapping (EC) and Genotype-Corrector (GC) were applied to the two replicates. Linkage maps were therefore constructed based on the corrected genotypes and then compared with those from the non-erroneous data set. Simulation study was performed by considering different levels of genotyping errors to investigate the impact of errors and the accuracy of error correction methods. Results indicated that map length and marker order differed among the two replicates and the non-erroneous data sets in both RIL populations. For both actual and simulated populations, map length was expanded as the increase in error rate, and the correlation coefficient between linkage and physical maps became lower. Map quality can be improved by repeated genotyping and error correction algorithm. When it is impossible to genotype the whole mapping population repeatedly, 30% would be recommended in repeated genotyping. The EC method had a much lower false positive rate than did the GC method under different error rates. This study systematically expounded the impact of genotyping errors on linkage analysis, providing potential guidelines for improving the accuracy of linkage maps in the presence of genotyping errors.

Insights into the genetic architecture of Phytophthora capsici root rot resistance in chile pepper (Capsicum spp.) from multi-locus genome-wide association study.

BACKGROUND: Phytophthora root rot, a major constraint in chile pepper production worldwide, is caused by the soil-borne oomycete, Phytophthora capsici. This study aimed to detect significant regions in the Capsicum genome linked to Phytophthora root rot resistance using a panel consisting of 157 Capsicum spp. genotypes. Multi-locus genome wide association study (GWAS) was conducted using single nucleotide polymorphism (SNP) markers derived from genotyping-by-sequencing (GBS). Individual plants were separately inoculated with P. capsici isolates, 'PWB-185', 'PWB-186', and '6347', at the 4-8 leaf stage and were scored for disease symptoms up to 14-days post-inoculation. Disease scores were used to calculate disease parameters including disease severity index percentage, percent of resistant plants, area under disease progress curve, and estimated marginal means for each genotype. RESULTS: Most of the genotypes displayed root rot symptoms, whereas five accessions were completely resistant to all the isolates and displayed no symptoms of infection. A total of 55,117 SNP markers derived from GBS were used to perform multi-locus GWAS which identified 330 significant SNP markers associated with disease resistance. Of these, 56 SNP markers distributed across all the 12 chromosomes were common across the isolates, indicating association with more durable resistance. Candidate genes including nucleotide-binding site leucine-rich repeat (NBS-LRR), systemic acquired resistance (SAR8.2), and receptor-like kinase (RLKs), were identified within 0.5 Mb of the associated markers. CONCLUSIONS: Results will be used to improve resistance to Phytophthora root rot in chile pepper by the development of Kompetitive allele-specific markers (KASP®) for marker validation, genomewide selection, and marker-assisted breeding.

Genome-wide association study, population structure, and genetic diversity of the tea plant in Guizhou Plateau.

BACKGROUND: Guizhou Plateau, as one of the original centers of tea plant, has a profound multi-ethnic cultural heritage and abundant tea germplasm resources. However, the impact of indigenous community factors on the genetic diversity, population structure and geographical distribution of tea plant is still unclear. RESULTS: Using the genotyping-by-sequencing (GBS) approach, we collected 415 tea plant accessions from the study sites, estimated genetic diversity, developed a core collection, and conducted a genome-wide association study (GWAS) based on 99,363 high-quality single-nucleotide polymorphisms (SNPs). A total of 415 tea accessions were clustered into six populations (GP01, GP02, GP03, GP04, GP05 and GP06), and the results showed that GP04 and GP05 had the highest and lowest genetic diversity (Pi = 0.214 and Pi = 0.145, respectively). Moreover, 136 tea accessions (33%) were selected to construct the core set that can represent the genetic diversity of the whole collection. By analyzing seven significant SNP markers associated with the traits such as the germination period of one bud and two leaves (OTL) and the germination period of one bud and three leaves (OtL), four candidate genes possibly related to OTL and OtL were identified. CONCLUSIONS: This study revealed the impact of indigenous communities on the population structure of 415 tea accessions, indicating the importance of cultural practices for protection and utilization of tea plant genetic resources. Four potential candidate genes associated with the OTL and OtL of tea plant were also identified, which will facilitate genetic research, germplasm conservation, and breeding.

Spectrum of germline pathogenic variants in Brazilian hereditary breast/ovarian cancer cases.

PURPOSE: To define the spectrum of germline pathogenic variants (PVs) and copy number variant (CNV) in cancer susceptibility genes to the burden of breast and ovarian cancer (BC, OvC) in high-risk Brazilians in Minas Gerais with health insurance, southeast Brazil, undergoing multigene panel testing (MGPT). METHODS: Genotyping eligible individuals with health insurance in the Brazilian healthcare system for Hereditary Breast and Ovarian Cancer Syndrome to undergo molecular testing for 44 or 141-gene panels, a decision that was insurance driven. RESULTS: Overall, 701 individuals clinically defined as high BC/OvC risk, underwent MGPT from 1/2021 to 10/2022, with ~ 50% genotyped with a 44-gene panel and the rest with a 141-gene panel. Overall, 16.4% and 22.6% of genotyped individuals harbored PVs using 44-gene and the 141 gene panel, respectively. The most frequently mutated genes were: BRCA2 (3.7%); BRCA1 (3.6%) and monoallelic MUTYH (3.1%). CONCLUSION: The rate of PVs detected in high-risk individuals in this study was twice the 10% threshold used in Brazilian health guidelines. MGPT doubled the detection rate of PVs in cancer susceptibility genes in high-risk individuals compared with BRCA1/BRCA2 genotyping alone. The spectrum of PVs in Southern Brazil is diverse, with few recurring variants such as TP53 (0.6%), suggesting regional founder effects. The use of MGPT in hereditary cancer in Minas Gerais significantly increased the detection rate of P/LPVs compared to existing guidelines and should be considered as the primary genotyping modality in assessing hereditary cancer risk in Brazil.

Increased circulation of human adenovirus in 2023: an investigation of the circulating genotypes, upper respiratory viral loads, and hospital admissions in a large academic medical center.

IMPORTANCE: The circulation of human adenoviruses (HAdV) increased in 2023. In this manuscript, we show that HAdV-B3 was predominant in 2023 in a cohort characterized by the Johns Hopkins Hospital System. We also show that HAdV-B3 was associated with an increase in viral loads in respiratory samples and provide a correlation with the clinical presentations and outcomes.

Performance of the Applied Biosystems HIV-1 Genotyping Kit with Integrase.

HIV genotyping is used to assess HIV susceptibility to antiretroviral drugs. The Applied Biosystems HIV-1 Genotyping Kit with Integrase (AB kit, Thermo Fisher Scientific) detects resistance-associated mutations (RAMs) in HIV protease (PR), reverse transcriptase (RT), and integrase (IN). We compared results from the AB kit with results obtained previously with the ViroSeq HIV-1 Genotyping System. DNA amplicons from the AB kit were also analyzed using next-generation sequencing (NGS). HIV RNA was extracted using the MagNA Pure 24 instrument (Roche Diagnostics; 96 plasma samples, HIV subtype B, viral load range: 530-737,741 copies/mL). FASTA files were generated from AB kit data using Exatype (Hyrax Biosciences). DNA amplicons from the AB kit were also analyzed by NGS using the Nextera XT kit (Illumina). Drug resistance was predicted using the Stanford HIV Drug Resistance Database. The mean genetic distance for sequences from ViroSeq and the AB kit was 0.02% for PR/RT and 0.04% for IN; 103 major RAMs were detected by both methods. Four additional major RAMs were detected by the AB kit only. These four major RAMs were also detected by NGS (detected in 18.1%-38.2% of NGS reads). NGS detected 27 major RAMs that were not detected with either of the Sanger sequencing-based kits. All major RAMs detected with ViroSeq were detected with the AB kit; additional RAMs were detected with the AB kit only. DNA amplicons from the AB kit can be used for NGS for more sensitive detection of RAMs.

Consensus genetic linkage map and QTL mapping allow to capture the genomic regions associated with agronomic traits in pearl millet.

MAIN CONCLUSION: A genetic linkage map representing the pearl millet genome was constructed with SNP markers. Major and stable QTL associated with flowering, number of productive tillers, ear head length, and test weight were mapped on chromosomes 1 and 3. Pearl millet (Pennisetum glaucum) is a major cereal and fodder crop in arid and semi-arid regions of Asia and Africa. Agronomic traits are important traits in pearl millet breeding and genetic and environmental factors highly influence them. In the present study, an F9 recombinant inbred line (RIL) population derived from a cross between PT6029 and PT6129 was evaluated for agronomic traits in three environments. Utilizing a genotyping by sequencing approach, a dense genetic map with 993 single nucleotide polymorphism markers covering a total genetic distance of 1035.4 cM was constructed. The average interval between the markers was 1.04 cM, and the seven chromosomes varied from 115.39 to 206.72 cM. Quantitative trait loci (QTL) mapping revealed 35 QTL for seven agronomic traits, and they were distributed on all pearl millet chromosomes. These QTL individually explained 11.35 to 26.71% of the phenotypic variation, with LOD values ranging from 2.74 to 5.80. Notably, four QTL (qDFF1.1, qNPT3.1, qEHL3.1, and qTW1.1) associated with days to fifty percent flowering, the number of productive tillers, ear head length, and test weight were found to be major and stable QTL located on chromosomes 1 and 3. Collectively, our results provide an important base for understanding the genetic architecture of agronomic traits in pearl millet, which is useful for accelerating the genetic gain toward crop improvement.

Self-(in)compatibility in Tunisian apple accessions [Malus domestica. Borkh]: S-genotypes identification and pollen tube growth analysis.

MAIN CONCLUSION: Self-incompatibility studies have revealed a potential use of Tunisian apple resources for crop improvement and modern breeding programs and a likely correlation between the pollen tube growth and flowering period. Apples [Malus domestica. Borkh] exhibit an S-RNase-based gametophytic self-incompatibility (GSI) system. Four primer combinations were used to S-genotype eighteen Tunisian local apple accessions and twelve introduced accessions that served as references. Within the Tunisian local accessions, S2, S3, S7, and S28 S-alleles were the most frequent and were assigned to 14 S-genotypes; among them, S7S28, S3S7, S2S5, and S2S3 were the most abundant. PCA plot showed that population structuring was affected by the S-alleles frequencies and revealed a modern origin of the Tunisian varieties rather than being ancient ones. Nonetheless, the results obtained with 17 SSR markers showed a separate grouping of local Tunisian accessions that calls into question the hypothesis discussed. Pollination experiments showed that the pollen started to germinate within 24 h of pollination but 48 h after pollination in the "El Fessi" accession. The first pollen tubes arrived in the styles within 36 h of pollination in two early flowering accessions known as "Arbi" and "Bokri", and after 72 h of pollination in late flowering "El Fessi" and 48 h after pollination in remaining accessions. The first pollen tube arrests were observed in accessions "Arbi" and "Bokri" within 84 h of pollination, within 108 h of pollination in "El Fessi" and within 108 h of pollination in remaining accessions. In the apple accession called "Boutabgaya," the pollen tubes reached the base of the style within 120 h of pollination without being aborted. Nevertheless, the self-compatible nature of "Boutabgaya" needs more studies to be confirmed. However, our results revealed the malfunction of the female component of the GSI in this accession. To conclude, this work paved the path for further studies to enhance the insight (i) into the relation between the flowering period and the pollen tube growth, (ii) self-compatible nature of "Boutabgaya", and (iii) the origin of the Tunisian apple.