Empirical threshold values for quantitative trait mapping

Immunological roads diverged: mapping tuberculosis outcomes in mice

The journey from phenotypic observation to causal genetic mechanism is a long and challenging road. For pathogens like Mycobacterium tuberculosis (Mtb), which causes tuberculosis (TB), host-pathogen coevolution has spanned millennia, costing millions of human lives. Mammalian models can systematically recapitulate host genetic variation, producing a spectrum of disease outcomes. Leveraging genome sequences and deep phenotyping data from infected mouse genetic reference populations (GRPs), quantitative trait locus (QTL) mapping approaches have successfully identified host genomic regions associated with TB phenotypes. Here, we review the ongoing optimization of QTL mapping study design alongside advances in mouse GRPs. These next-generation resources and approaches have enabled identification of novel host-pathogen interactions governing one of the most prevalent infectious diseases in the world today.

A Plasmodium falciparum genetic cross reveals the contributions of pfcrt and plasmepsin II/III to piperaquine drug resistance

Piperaquine (PPQ) is widely used in combination with dihydroartemisinin as a first-line treatment against malaria. Multiple genetic drivers of PPQ resistance have been reported, including mutations in the Plasmodium falciparum chloroquine resistance transporter (pfcrt) and increased copies of plasmepsin II/III (pm2/3). We generated a cross between a Cambodia-derived multidrug-resistant KEL1/PLA1 lineage isolate (KH004) and a drug-susceptible Malawian parasite (Mal31). Mal31 harbors a wild-type (3D7-like) pfcrt allele and a single copy of pm2/3, while KH004 has a chloroquine-resistant (Dd2-like) pfcrt allele with an additional G367C substitution and multiple copies of pm2/3. We recovered 104 unique recombinant parasites and examined a targeted set of progeny representing all possible combinations of variants at pfcrt and pm2/3. We performed a detailed analysis of competitive fitness and a range of PPQ susceptibility phenotypes with these progenies, including PPQ survival assay, area under the dose response curve, and a limited point IC50. We find that inheritance of the KH004 pfcrt allele is required for reduced PPQ sensitivity, whereas copy number variation in pm2/3 further decreases susceptibility but does not confer resistance in the absence of additional mutations in pfcrt. A deep investigation of genotype-phenotype relationships demonstrates that progeny clones from experimental crosses can be used to understand the relative contributions of pfcrt, pm2/3, and parasite genetic background to a range of PPQ-related traits. Additionally, we find that the resistance phenotype associated with parasites inheriting the G367C substitution in pfcrt is consistent with previously validated PPQ resistance mutations in this transporter.IMPORTANCEResistance to piperaquine, used in combination with dihydroartemisinin, has emerged in Cambodia and threatens to spread to other malaria-endemic regions. Understanding the causal mutations of drug resistance and their impact on parasite fitness is critical for surveillance and intervention and can also reveal new avenues to limiting the evolution and spread of drug resistance. An experimental genetic cross is a powerful tool for pinpointing the genetic determinants of key drug resistance and fitness phenotypes and has the distinct advantage of quantifying the effects of naturally evolved genetic variation. Our study was strengthened since the full range of copies of KH004 pm2/3 was inherited among the progeny clones, allowing us to directly test the role of the pm2/3 copy number on resistance-related phenotypes in the context of a unique pfcrt allele. Our multigene model suggests an important role for both loci in the evolution of this multidrug-resistant parasite lineage.

Quantitative trait locus analysis of gray leaf spot resistance in the maize IBM Syn10 DH population

KEY MESSAGE

The exploration and dissection of a set of QTLs and candidate genes for gray leaf spot disease resistance using two fully assembled parental genomes may help expedite maize resistance breeding. The fungal disease of maize known as gray leaf spot (GLS), caused by Cercospora zeae-maydis and Cercospora zeina, is a significant concern in China, Southern Africa, and the USA. Resistance to GLS is governed by multiple genes with an additive effect and is influenced by both genotype and environment. The most effective way to reduce the cost of production is to develop resistant hybrids. In this study, we utilized the IBM Syn 10 Doubled Haploid (IBM Syn10 DH) population to identify quantitative trait loci (QTLs) associated with resistance to gray leaf spot (GLS) in multiple locations. Analysis of seven distinct environments revealed a total of 58 QTLs, 49 of which formed 12 discrete clusters distributed across chromosomes 1, 2, 3, 4, 8 and 10. By comparing these findings with published research, we identified colocalized QTLs or GWAS loci within eleven clustering intervals. By integrating transcriptome data with genomic structural variations between parental individuals, we identified a total of 110 genes that exhibit both robust disparities in gene expression and structural alterations. Further analysis revealed 19 potential candidate genes encoding conserved resistance gene domains, including putative leucine-rich repeat receptors, NLP transcription factors, fucosyltransferases, and putative xyloglucan galactosyltransferases. Our results provide a valuable resource and linked loci for GLS marker resistance selection breeding in maize.

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

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

Evidence for an Association Between a pH-Dependent Potassium Channel, TWIK-1, and the Accuracy of Smooth Pursuit Eye Movements

Purpose

Within the healthy population there is a large variation in the ability to perform smooth pursuit eye movements. Our purpose was to investigate the genetic and physiological bases for this variation.

Methods

We carried out a whole-genome association study, recording smooth pursuit movements for 1040 healthy volunteers by infrared oculography. The primary phenotypic measure was root mean square error (RMSE) of eye position relative to target position. Secondary measures were pursuit gain, frequency of catch-up saccades, and frequency of anticipatory saccades. Ten percent of participants, chosen randomly, were tested twice, giving estimates of test-retest reliability.

Results

No significant association was found with three genes previously identified as candidate genes for variation in smooth pursuit: DRD3, COMT, NRG1. A strong association (P = 3.55 × 10-11) was found between RMSE and chromosomal region 1q42.2. The most strongly associated marker (rs701232) lies in an intron of KCNK1, which encodes a two-pore-domain potassium ion channel TWIK-1 (or K2P1) that affects cell excitability. Each additional copy of the A allele decreased RMSE by 0.29 standard deviation. When a psychophysical test of visually perceived motion was used as a covariate in the regression analysis, the association with rs701232 did not weaken (P = 5.38 × 10-12).

Conclusions

Variation in the sequence or the expression of the pH-dependent ion channel TWIK-1 is a likely source of variance in smooth pursuit. The variance associated with TWIK-1 appears not to arise from sensory mechanisms, because the use of a perceptual covariate left the association intact.

Discovery of genomic loci for liver health and steatosis reveals overlap with glutathione redox genetics

Non-alcoholic fatty liver disease (NAFLD) is the most common chronic liver condition in the United States, encompassing a wide spectrum of liver pathologies including steatosis, steatohepatitis, fibrosis, and cirrhosis. Despite its high prevalence, there are no medications currently approved by the Food and Drug Administration for the treatment of NAFLD. Recent work has suggested that NAFLD has a strong genetic component and identifying causative genes will improve our understanding of the molecular mechanisms contributing to NAFLD and yield targets for future therapeutic investigations. Oxidative stress is known to play an important role in NAFLD pathogenesis, yet the underlying mechanisms accounting for disturbances in redox status are not entirely understood. To better understand the relationship between the glutathione redox system and signs of NAFLD in a genetically-diverse population, we measured liver weight, serum biomarkers aspartate aminotransferase (AST) and alanine aminotransferase (ALT), and graded liver pathology in a large cohort of Diversity Outbred mice. We compared hepatic endpoints to those of the glutathione redox system previously measured in the livers and kidneys of the same mice, and we screened for statistical and genetic associations using the R/qtl2 software. We discovered several novel genetic loci associated with markers of liver health, including loci that were associated with both liver steatosis and glutathione redox status. Candidate genes within each locus point to possible new mechanisms underlying the complex relationship between NAFLD and the glutathione redox system, which could have translational implications for future studies targeting NAFLD pathology.

Combining genome-wide association study and linkage mapping in the genetic dissection of amylose content in maize (Zea mays L.)

KEY MESSAGE

Integrated linkage and association analysis revealed genetic basis across multiple environments. The genes Zm00001d003102 and Zm00001d015905 were further verified to influence amylose content using gene-based association study. Maize kernel amylose is an important source of human food and industrial raw material. However, the genetic basis underlying maize amylose content is still obscure. Herein, we used an intermated B73 × Mo17 (IBM) Syn10 doubled haploid population composed of 222 lines and a germplasm set including 305 inbred lines to uncover the genetic control for amylose content under four environments. Linkage mapping detected 16 unique QTL, among which four were individually repeatedly identified across multiple environments. Genome-wide association study revealed 17 significant (P = 2.24E-06) single-nucleotide polymorphisms, of which two (SYN19568 and PZE-105090500) were located in the intervals of the mapped QTL (qAC2 and qAC5-3), respectively. According to the two population co-localized loci, 20 genes were confirmed as the candidate genes for amylose content. Gene-based association analysis indicated that the variants in Zm00001d003102 (Beta-16-galactosyltransferase GALT29A) and Zm00001d015905 (Sugar transporter 4a) affected amylose content across multi-environment. Tissue expression analysis showed that the two genes were specifically highly expressed in the ear and stem, respectively, suggesting that they might participate in sugar transport from source to sink organs. Our study provides valuable genetic information for breeding maize varieties with high amylose.

Multi-GWAS reveals significant genomic regions for Mungbean yellow mosaic India virus resistance in urdbean (Vigna mungo (L.) across multiple environments

KEY MESSAGE

Unraveling genetic markers for MYMIV resistance in urdbean, with 8 high-confidence marker-trait associations identified across diverse environments, provides crucial insights for combating MYMIV disease, informing future breeding strategies. Globally, yellow mosaic disease (YMD) causes significant yield losses, reaching up to 100% in favorable environments within major urdbean cultivating regions. The introgression of genomic regions conferring resistance into urdbean cultivars is crucial for combating YMD, including resistance against mungbean yellow mosaic India virus (MYMIV). To uncover the genetic basis of MYMIV resistance, we conducted a genome-wide association study (GWAS) using three multi-locus models in 100 diverse urdbean genotypes cultivated across six individual and two combined environments. Leveraging 4538 high-quality single nucleotide polymorphism (SNP) markers, we identified 28 unique significant marker-trait associations (MTAs) for MYMIV resistance, with 8 MTAs considered of high confidence due to detection across multiple GWAS models and/or environments. Notably, 4 out of 28 MTAs were found in proximity to previously reported genomic regions associated with MYMIV resistance in urdbean and mungbean, strengthening our findings and indicating consistent genomic regions for MYMIV resistance. Among the eight highly significant MTAs, one localized on chromosome 6 adjacent to previously identified quantitative trait loci for MYMIV resistance, while the remaining seven were novel. These MTAs contain several genes implicated in disease resistance, including four common ones consistently found across all eight MTAs: receptor-like serine-threonine kinases, E3 ubiquitin-protein ligase, pentatricopeptide repeat, and ankyrin repeats. Previous studies have linked these genes to defense against viral infections across different crops, suggesting their potential for further basic research involving cloning and utilization in breeding programs. This study represents the first GWAS investigation aimed at identifying resistance against MYMIV in urdbean germplasm.

QTL Analysis of β-Glucan Content and Other Grain Traits in a Recombinant Population of Spring Barley

Barley with high grain β-glucan content is valuable for functional foods. The identification of loci for high β-glucan content is, thus, of great importance for barley breeding. Segregation mapping for the content in β-glucan and other barley grain components (starch, protein, lipid, ash, phosphorous, calcium, sodium) was performed using the progeny of the cross between Glacier AC38, a mutant with high amylose, and CDC Fibar, a high β-glucan waxy cultivar. The offspring of this cross showed transgressive segregation for β-glucan content. Linkage analysis based on single-nucleotide polymorphism (SNP) molecular markers was used for the genotyping of the parents and recombinant inbred lines (RILs). Two Quantitative Trait Loci (QTL) for β-glucan content and several QTL for other grain components were found. The former ones, located on chromosomes 1H and 7H, explained 27.9% and 27.4% of the phenotypic variance, respectively. Glacier AC38 provided the allele for high β-glucan content at the QTL on chromosome 1H, whereas CDC Fibar contributed the allele at the QTL on chromosome 7H. Their recombination resulted in a novel haplotype with higher β-glucan content, up to 18.4%. Candidate genes are proposed for these two QTL: HvCslF9, involved in β-glucan biosynthesis, for the QTL on chromosome 1H; Horvu_PLANET_7H01G069300, a gene encoding an ATP-Binding Cassette (ABC) transporter, for the QTL on chromosome 7H.

A systems genomics and genetics approach to identify the genetic regulatory network for lignin content in Brassica napus seeds

Seed quality traits of oilseed rape, Brassica napus (B. napus), exhibit quantitative inheritance determined by its genetic makeup and the environment via the mediation of a complex genetic architecture of hundreds to thousands of genes. Thus, instead of single gene analysis, network-based systems genomics and genetics approaches that combine genotype, phenotype, and molecular phenotypes offer a promising alternative to uncover this complex genetic architecture. In the current study, systems genetics approaches were used to explore the genetic regulation of lignin traits in B. napus seeds. Four QTL (qLignin_A09_1, qLignin_A09_2, qLignin_A09_3, and qLignin_C08) distributed on two chromosomes were identified for lignin content. The qLignin_A09_2 and qLignin_C08 loci were homologous QTL from the A and C subgenomes, respectively. Genome-wide gene regulatory network analysis identified eighty-three subnetworks (or modules); and three modules with 910 genes in total, were associated with lignin content, which was confirmed by network QTL analysis. eQTL (expression quantitative trait loci) analysis revealed four cis-eQTL genes including lignin and flavonoid pathway genes, cinnamoyl-CoA-reductase (CCR1), and TRANSPARENT TESTA genes TT4, TT6, TT8, as causal genes. The findings validated the power of systems genetics to identify causal regulatory networks and genes underlying complex traits. Moreover, this information may enable the research community to explore new breeding strategies, such as network selection or gene engineering, to rewire networks to develop climate resilience crops with better seed quality.

Identification and Validation of Quantitative Trait Loci Associated with Fruit Puffiness in a Processing Tomato Population

Physiological disorders impact the yield and quality of marketable fruit in tomato. Puffy fruit caused by cavities inside the locule can be problematic for processing and fresh market quality. In this paper, we used a recombinant inbred line (RIL) and three derived processing tomato populations to map and validate quantitative trait loci (QTLs) for fruit puffiness across environments. Binary interval mapping was used for mapping the incidence of fruit puffiness, and non-parametric interval mapping and parametric composite interval mapping were used for mapping severity. Marker-trait regressions were carried out to validate putative QTLs in subsequent crosses. QTLs were detected on chromosome (Chr) 1, 2, and 4. Only the QTL on Chr 1 was validated in progeny from subsequent crosses. This QTL explained up to 22.5% of the variance in the percentage of puffy fruit, with a significant interaction between loci on Chr 2 and 4, increasing the percentage of puffy fruit by an additional 15%. The allele responsible for puffy fruit on Chr 1 was inherited from parent FG02-188 and was dominant towards increased incidence and severity. Marker-assisted selection (MAS) for the QTL on Chr 1 was as efficient as genomic selection (GS) in reducing the incidence and severity of puffy fruit, despite the potential contribution of other loci.

Selection for Phytophthora Root Rot Resistance in Chickpea Crosses Affects Yield Potential of Chickpea × Cicer echinospermum Backcross Derivatives

Phytophthora root rot (PRR) of chickpea (Cicer arietinum) caused by Phytophthora medicaginis is an important disease. Partial resistance to PRR is sourced from Cicer echinospermum. In this study, we evaluated if lines with low levels of PRR foliage symptoms in two contrasting recombinant inbred line (RIL) populations parented by chickpea cultivars (Yorker and Rupali) and 04067-81-2-1-1 (C. echinospermum, interspecific breeding line) had a significant drag on yield parameters. For the Yorker × 04067-81-2-1-1 population with the highest level of PRR resistance, in the absence of PRR, low foliage symptom RIL had significantly later flowering and podding, lower grain yields, and lighter seed and shorter plant phenotypes than high foliage symptom RIL. A quantitative trait locus analysis identified significant QTL for flowering, height, 100-seed weight, and yield, and there was a significantly higher frequency of alleles for the negative agronomic traits (i.e., drag) from the 04067-81-2-1-1 parent in low foliage symptom RIL than in high foliage symptom RIL. For the Rupali × 04067-81-2-1-1 population with lower levels of PRR resistance, in the absence of PRR, low foliage symptom RIL had significantly lighter seed and shorter plants than high foliage symptom RIL. Significant QTL were detected, the majority were for the timing of flowering and podding (n = 18), others were for plant height, yield, and 100-seed weight. For this second population, the frequency of alleles for the negative agronomic traits from the 04067-81-2-1-1 parent did not differ between low and high foliage symptom RIL. The 100 seed weight of RIL under moderate PRR disease pressure showed some promise as a yield component trait to identify phenotypes with both high levels of PRR resistance and grain yield potential for further seed number evaluations. We identified that large population sizes are required to enable selection among chickpea × C. echinospermum crosses for high levels of PRR resistance without a significant drag on yield.

Dissection of Common Rust Resistance in Tropical Maize Multiparent Population through GWAS and Linkage Studies

Common rust (CR), caused by Puccina sorghi, is a major foliar disease in maize that leads to quality deterioration and yield losses. To dissect the genetic architecture of CR resistance in maize, this study utilized the susceptible temperate inbred line Ye107 as the male parent crossed with three resistant tropical maize inbred lines (CML312, D39, and Y32) to generate 627 F7 recombinant inbred lines (RILs), with the aim of identifying maize disease-resistant loci and candidate genes for common rust. Phenotypic data showed good segregation between resistance and susceptibility, with varying degrees of resistance observed across different subpopulations. Significant genotype effects and genotype × environment interactions were observed, with heritability ranging from 85.7% to 92.2%. Linkage and genome-wide association analyses across the three environments identified 20 QTLs and 62 significant SNPs. Among these, seven major QTLs explained 66% of the phenotypic variance. Comparison with six SNPs repeatedly identified across different environments revealed overlap between qRUST3-3 and Snp-203,116,453, and Snp-204,202,469. Haplotype analysis indicated two different haplotypes for CR resistance for both the SNPs. Based on LD decay plots, three co-located candidate genes, Zm00001d043536, Zm00001d043566, and Zm00001d043569, were identified within 20 kb upstream and downstream of these two SNPs. Zm00001d043536 regulates hormone regulation, Zm00001d043566 controls stomatal opening and closure, related to trichome, and Zm00001d043569 is associated with plant disease immune responses. Additionally, we performed candidate gene screening for five additional SNPs that were repeatedly detected across different environments, resulting in the identification of five candidate genes. These findings contribute to the development of genetic resources for common rust resistance in maize breeding programs.

Cloning a novel reduced-height (Rht) gene TaOSCA1.4 from a QTL in wheat

Reducing plant height (PH) is one of the core contents of the "Green Revolution", which began in the 1960s in wheat. A number of 27 reduced-height (Rht) genes have been identified and a great number of quantitative trait loci (QTLs) for PH have been mapped on all 21 chromosomes. Nonetheless, only several genes regulated PH have been cloned. In this study, we found the interval of QTL QPh-1B included an EST-SSR marker swes1079. According to the sequence of swes1079, we cloned the TaOSCA1.4 gene. We developed a CAPS marker to analyze the variation across a natural population. The result showed that the PH was significantly different between the two haplotypes of TaOSCA1.4-1B under most of the 12 environments and the average values of irrigation and rainfed conditions. This result further demonstrated that TaOSCA1.4 was associated with PH. Then, we validated the TaOSCA1.4 via RNAi technology. The average PHs of the wild-type (WT), RNAi lines 1 (Ri-1) and 2 (Ri-2) were 94.6, 83.6 and 79.2 cm, respectively, with significant differences between the WT and Ri-1 and Ri-2. This result indicated that the TaOSCA1.4 gene controls PH. TaOSCA1.4 is a constitutively expressed gene and its protein localizes to the cell membrane. TaOSCA1.4 gene is a member of the OSCA gene family, which regulates intracellular Ca2+ concentration. We hypothesized that knock down mutants of TaOSCA1.4 gene reduced regulatory ability of Ca2+, thus reducing the PH. Furthermore, the cell lengths of the knock down mutants are not significantly different than that of WT. We speculate that TaOSCA1.4 gene is not directly associated with gibberellin (GA), which should be a novel mechanism for a wheat Rht gene.

The contribution of an X chromosome QTL to non-Mendelian inheritance and unequal chromosomal segregation in Auanema freiburgense

Auanema freiburgense is a nematode with males, females, and selfing hermaphrodites. When XO males mate with XX females, they typically produce a low proportion of XO offspring because they eliminate nullo-X spermatids. This process ensures that most sperm carry an X chromosome, increasing the likelihood of X chromosome transmission compared to random segregation. This occurs because of an unequal distribution of essential cellular organelles during sperm formation, likely dependent on the X chromosome. Some sperm components are selectively segregated into the X chromosome's daughter cell, while others are discarded with the nullo-X daughter cell. Intriguingly, the interbreeding of 2 A. freiburgense strains results in hybrid males capable of producing viable nullo-X sperm. Consequently, when these hybrid males mate with females, they yield a high percentage of male offspring. To uncover the genetic basis of nullo-spermatid elimination and X chromosome drive, we generated a genome assembly for A. freiburgense and genotyped the intercrossed lines. This analysis identified a quantitative trait locus spanning several X chromosome genes linked to the non-Mendelian inheritance patterns observed in A. freiburgense. This finding provides valuable clues to the underlying factors involved in asymmetric organelle partitioning during male meiotic division and thus non-Mendelian transmission of the X chromosome and sex ratios.

Naturally segregating genetic variants contribute to thermal tolerance in a Drosophila melanogaster model system

Thermal tolerance is a fundamental physiological complex trait for survival in many species. For example, everyday tasks such as foraging, finding a mate, and avoiding predation are highly dependent on how well an organism can tolerate extreme temperatures. Understanding the general architecture of the natural variants within the genes that control this trait is of high importance if we want to better comprehend thermal physiology. Here, we take a multipronged approach to further dissect the genetic architecture that controls thermal tolerance in natural populations using the Drosophila Synthetic Population Resource as a model system. First, we used quantitative genetics and Quantitative Trait Loci mapping to identify major effect regions within the genome that influences thermal tolerance, then integrated RNA-sequencing to identify differences in gene expression, and lastly, we used the RNAi system to (1) alter tissue-specific gene expression and (2) functionally validate our findings. This powerful integration of approaches not only allows for the identification of the genetic basis of thermal tolerance but also the physiology of thermal tolerance in a natural population, which ultimately elucidates thermal tolerance through a fitness-associated lens.

Collaborative Cross mice have diverse phenotypic responses to infection with Methicillin-resistant Staphylococcus aureus USA300

Staphylococcus aureus (S. aureus) is an opportunistic pathogen causing diseases ranging from mild skin infections to life threatening conditions, including endocarditis, pneumonia, and sepsis. To identify host genes modulating this host-pathogen interaction, we infected 25 Collaborative Cross (CC) mouse strains with methicillin-resistant S. aureus (MRSA) and monitored disease progression for seven days using a surgically implanted telemetry system. CC strains varied widely in their response to intravenous MRSA infection. We identified eight 'susceptible' CC strains with high bacterial load, tissue damage, and reduced survival. Among the surviving strains, six with minimal colonization were classified as 'resistant', while the remaining six tolerated higher organ colonization ('tolerant'). The kidney was the most heavily colonized organ, but liver, spleen and lung colonization were better correlated with reduced survival. Resistant strains had higher pre-infection circulating neutrophils and lower post-infection tissue damage compared to susceptible and tolerant strains. We identified four CC strains with sexual dimorphism: all females survived the study period while all males met our euthanasia criteria earlier. In these CC strains, males had more baseline circulating monocytes and red blood cells. We identified several CC strains that may be useful as new models for endocarditis, myocarditis, pneumonia, and resistance to MRSA infection. Quantitative Trait Locus (QTL) analysis identified two significant loci, on Chromosomes 18 and 3, involved in early susceptibility and late survival after infection. We prioritized Npc1 and Ifi44l genes as the strongest candidates influencing survival using variant analysis and mRNA expression data from kidneys within these intervals.

Artificial selection footprints in indigenous and commercial chicken genomes

BACKGROUND

Although many studies have been done to reveal artificial selection signatures in commercial and indigenous chickens, a limited number of genes have been linked to specific traits. To identify more trait-related artificial selection signatures and genes, we re-sequenced a total of 85 individuals of five indigenous chicken breeds with distinct traits from Yunnan Province, China.

RESULTS

We found 30 million non-redundant single nucleotide variants and small indels (< 50 bp) in the indigenous chickens, of which 10 million were not seen in 60 broilers, 56 layers and 35 red jungle fowls (RJFs) that we compared with. The variants in each breed are enriched in non-coding regions, while those in coding regions are largely tolerant, suggesting that most variants might affect cis-regulatory sequences. Based on 27 million bi-allelic single nucleotide polymorphisms identified in the chickens, we found numerous selective sweeps and affected genes in each indigenous chicken breed and substantially larger numbers of selective sweeps and affected genes in the broilers and layers than previously reported using a rigorous statistical model. Consistent with the locations of the variants, the vast majority (~ 98.3%) of the identified selective sweeps overlap known quantitative trait loci (QTLs). Meanwhile, 74.2% known QTLs overlap our identified selective sweeps. We confirmed most of previously identified trait-related genes and identified many novel ones, some of which might be related to body size and high egg production traits. Using RT-qPCR, we validated differential expression of eight genes (GHR, GHRHR, IGF2BP1, OVALX, ELF2, MGARP, NOCT, SLC25A15) that might be related to body size and high egg production traits in relevant tissues of relevant breeds.

CONCLUSION

We identify 30 million single nucleotide variants and small indels in the five indigenous chicken breeds, 10 million of which are novel. We predict substantially more selective sweeps and affected genes than previously reported in both indigenous and commercial breeds. These variants and affected genes are good candidates for further experimental investigations of genotype-phenotype relationships and practical applications in chicken breeding programs.

Mapping QTL controlling count traits with excess zeros and ones using a zero-and-one-inflated generalized Poisson regression model

The research on the quantitative trait locus (QTL) mapping of count data has aroused the wide attention of researchers. There are frequent problems in applied research that limit the application of the conventional Poisson model in the analysis of count phenotypes, which include the overdispersion and excess zeros and ones. In this article, a novel model, that is, the zero-and-one-inflated generalized Poisson (ZOIGP) model, is proposed to deal with these problems. Based on the proposed model, a score test is performed for the inflation parameter, in which the ZOIGP model with a constant proportion of excess zeros and ones is compared with a standard generalized Poisson model. To illustrate the practicability of the ZOIGP model, we extend it to the QTL interval mapping application that underpins count phenotype with excess zeros and excess ones. The genetic effects are estimated utilizing the expectation-maximization algorithm embedded with the Newton-Raphson algorithm, and the genome-wide scan and likelihood ratio test is performed to map and test the potential QTLs. The statistical properties exhibited by the proposed method are investigated through simulation. Finally, a real data analysis example is used to illustrate the utility of the proposed method for QTL mapping.

Gene expression variation underlying tissue-specific responses to copper stress in Drosophila melanogaster

Copper is one of a handful of biologically necessary heavy metals that is also a common environmental pollutant. Under normal conditions, copper ions are required for many key physiological processes. However, in excess, copper results in cell and tissue damage ranging in severity from temporary injury to permanent neurological damage. Because of its biological relevance, and because many conserved copper-responsive genes respond to nonessential heavy metal pollutants, copper resistance in Drosophila melanogaster is a useful model system with which to investigate the genetic control of the heavy metal stress response. Because heavy metal toxicity has the potential to differently impact specific tissues, we genetically characterized the control of the gene expression response to copper stress in a tissue-specific manner in this study. We assessed the copper stress response in head and gut tissue of 96 inbred strains from the Drosophila Synthetic Population Resource using a combination of differential expression analysis and expression quantitative trait locus mapping. Differential expression analysis revealed clear patterns of tissue-specific expression. Tissue and treatment specific responses to copper stress were also detected using expression quantitative trait locus mapping. Expression quantitative trait locus associated with MtnA, Mdr49, Mdr50, and Sod3 exhibited both genotype-by-tissue and genotype-by-treatment effects on gene expression under copper stress, illuminating tissue- and treatment-specific patterns of gene expression control. Together, our data build a nuanced description of the roles and interactions between allelic and expression variation in copper-responsive genes, provide valuable insight into the genomic architecture of susceptibility to metal toxicity, and highlight candidate genes for future functional characterization.

Plant necrotrophic bacterial disease resistance phenotypes, QTL, and metabolites identified through integrated genetic mapping and metabolomics in Solanum species

Most food crops are susceptible to necrotrophic bacteria that cause rotting and wilting diseases in fleshy organs and foods. All varieties of cultivated potato (Solanum tuberosum L.) are susceptible to diseases caused by Pectobacterium species, but resistance has been demonstrated in wild potato relatives including S. chacoense. Previous studies demonstrated that resistance is in part mediated by antivirulence activity of phytochemicals in stems and tubers. Little is known about the genetic basis of antivirulence traits, and the potential for inheritance and introgression into cultivated potato is unclear. Here, the metabolites and genetic loci associated with antivirulence traits in S. chacoense were elucidated by screening a sequenced S. tuberosum x S. chacoense recombinant inbred line (RIL) population for antivirulence traits of its metabolite extracts. Metabolite extracts from the RILs exhibited a quantitative distribution for two antivirulence traits that were positively correlated: quorum sensing inhibition and exo-protease inhibition, with some evidence of transgressive segregation, supporting the role of multiple loci and metabolites regulating these resistance-associated systems. Metabolomics was performed on the highly resistant and susceptible RILs that revealed 30 metabolites associated with resistance, including several alkaloids and terpenes. Specifically, several prenylated metabolites were more abundant in resistant RILs. We constructed a high-density linkage map with 795 SNPs mapped to 12 linkage groups, spanning a length of 1,507 cM and a density of 1 marker per 1.89 cM. Genetic mapping of the antivirulence and metabolite data identified five quantitative trait loci (QTLs) related to quorum sensing inhibition that explained 8-28% of the phenotypic variation and two QTLs for protease activity inhibition that explained 14-19% of the phenotypic variation. Several candidate genes including alkaloid, and secondary metabolite biosynthesis that are related to disease resistance were identified within these QTLs. Taken together, these data support that quorum sensing inhibition and exo-protease inhibition assays may serve as breeding targets to improve resistance to nectrotrophic bacterial pathogens in potato and other plants. The identified candidate genes and metabolites can be utilized in marker assisted selection and genomic selection to improve soft- rot and blackleg disease resistance.

Reproducible Quantitative Trait Loci for Resistance to Soft Rot Caused by Dickeya dianthicola Derived from the Wild Potato Solanum microdontum (PI 458355) Are Located on Chromosomes 1, 3, and 5

The potato wild relative Solanum microdontum is a breeder-friendly source of genetic resistance to soft rot. Our objectives were to (i) identify loci associated with soft rot resistance in S. microdontum germplasm and (ii) develop bi-parental populations in a self-compatible S. tuberosum genetic background to recover segregating F2 progenies, construct a linkage map, and identify quantitative trait loci (QTLs). Under objective (i), tubers from 103 S. microdontum genotypes from the United States Potato Genebank were inoculated with a high virulence strain of Dickeya dianthicola, and lesion size was measured after a 24-h incubation period at 30°C. Association analysis using 3,490 polymorphic Infinium array SNP markers identified soft rot resistance loci on chromosomes 1, 2, 3, 5, 7, 8, 11, and 12. Under objective (ii), a resistant S. microdontum accession PI 458355 was crossed with a highly fertile, self-compatible, diploid S. tuberosum pollen parent (PI 654351) to generate segregating F2 populations. Composite interval mapping was conducted using a genetic linkage map with 970 GBS-based SNP markers. Reproducible QTLs were detected on chromosomes 1, 3, and 5, explaining 11, 13, and 23% of the phenotypic variation, respectively. Homozygous S. microdontum alleles at the QTL on chromosome 5 and heterozygous or homozygous S. microdontum alleles at QTLs on chromosomes 1 and 3 significantly decrease lesion size compared with the homozygous S. tuberosum parent. The germplasm created in these studies provides a resource for studying traits from S. microdontum, and we can use the advanced F2 selections for future potato improvement. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Soybean type-B response regulator GmRR1 mediates phosphorus uptake and yield by modifying root architecture

Phosphorus (P) plays a pivotal role in plant growth and development. Low P stress can greatly hamper plant growth. Here, we identified a QTL (named QPH-9-1), which is associated with P efficiency across multiple environments through linkage analysis and genome-wide association study. Furthermore, we successfully cloned the underlying soybean (Glycine max) gene GmRR1 (a soybean type-B Response Regulator 1) that encodes a type-B response regulator protein. Knockout of GmRR1 resulted in a substantial increase in plant height, biomass, P uptake efficiency, and yield-related traits due to the modification of root structure. In contrast, overexpression of GmRR1 in plants resulted in a decrease in these phenotypes. Further analysis revealed that knockout of GmRR1 substantially increased the levels of auxin and ethylene in roots, thereby promoting root hair formation and growth by promoting the formation of root hair primordium and lengthening the root apical meristem. Yeast two-hybrid, bimolecular fluorescence complementation, and dual-luciferase assays demonstrated an interaction between GmRR1 and Histidine-containing Phosphotransmitter protein 1. Expression analysis suggested that these proteins coparticipated in response to low P stress. Analysis of genomic sequences showed that GmRR1 underwent a selection during soybean domestication. Taken together, this study provides further insights into how plants respond to low P stress by modifying root architecture through phytohormone pathways.

Epistasis between mutator alleles contributes to germline mutation spectrum variability in laboratory mice

Maintaining germline genome integrity is essential and enormously complex. Although many proteins are involved in DNA replication, proofreading, and repair, mutator alleles have largely eluded detection in mammals. DNA replication and repair proteins often recognize sequence motifs or excise lesions at specific nucleotides. Thus, we might expect that the spectrum of de novo mutations - the frequencies of C>T, A>G, etc. - will differ between genomes that harbor either a mutator or wild-type allele. Previously, we used quantitative trait locus mapping to discover candidate mutator alleles in the DNA repair gene Mutyh that increased the C>A germline mutation rate in a family of inbred mice known as the BXDs (Sasani et al., 2022, Ashbrook et al., 2021). In this study we developed a new method to detect alleles associated with mutation spectrum variation and applied it to mutation data from the BXDs. We discovered an additional C>A mutator locus on chromosome 6 that overlaps Ogg1, a DNA glycosylase involved in the same base-excision repair network as Mutyh (David et al., 2007). Its effect depends on the presence of a mutator allele near Mutyh, and BXDs with mutator alleles at both loci have greater numbers of C>A mutations than those with mutator alleles at either locus alone. Our new methods for analyzing mutation spectra reveal evidence of epistasis between germline mutator alleles and may be applicable to mutation data from humans and other model organisms.

A novel QTL conferring Fusarium crown rot resistance on chromosome 2A in a wheat EMS mutant

KEY MESSAGE

A novel QTL on chromosome 2A for Fusarium crown rot resistance was identified and validated in wheat. Fusarium crown rot (FCR) is a fungal disease that causes significant yield losses in many cereal growing regions in the world. In this study, genetic analysis was conducted for a wheat EMS mutant C549 which showed stable resistance to FCR at seedling stage. A total of 10 QTL were detected on chromosomes 1A, 2A, 3B, 4A, 6B, and 7B using a population of 138 F7 recombinant inbred lines (RILs) derived from a cross between C549 and a Chinese germplasm 3642. A novel locus Qfcr.cau-2A, which accounted for up to 24.42% of the phenotypic variation with a LOD value of 12.78, was consistently detected across all six trials conducted. Furthermore, possible effects of heading date (HD) and plant height on FCR severity were also investigated in the mapping population. While plant height had no effects on FCR resistance, a weak and negative association between FCR resistance and HD was observed. A QTL for HD (Qhd.cau-2A.2) was coincident with Qfcr.cau-2A. Conditional QTL mapping indicated that although Qfcr.cau-2A and Qhd.cau-2A.2 had significant interactions, Qfcr.cau-2A remained significant after the effects of HD was removed. It is unlikely that genes underlying these two loci are same. Nevertheless, the stable expression of Qfcr.cau-2A in the validation population of 148 F7 RILs developed between C549 and its wild parent Chuannong 16 demonstrated the potential value of this locus in FCR resistance breeding programs.

Identification of genetic drivers of plasma lipoprotein size in the Diversity Outbred mouse population

Despite great progress in understanding lipoprotein physiology, there is still much to be learned about the genetic drivers of lipoprotein abundance, composition, and function. We used ion mobility spectrometry to survey 16 plasma lipoprotein subfractions in 500 Diversity Outbred mice maintained on a Western-style diet. We identified 21 quantitative trait loci (QTL) affecting lipoprotein abundance. To refine the QTL and link them to disease risk in humans, we asked if the human homologs of genes located at each QTL were associated with lipid traits in human genome-wide association studies. Integration of mouse QTL with human genome-wide association studies yielded candidate gene drivers for 18 of the 21 QTL. This approach enabled us to nominate the gene encoding the neutral ceramidase, Asah2, as a novel candidate driver at a QTL on chromosome 19 for large HDL particles (HDL-2b). To experimentally validate Asah2, we surveyed lipoproteins in Asah2-/- mice. Compared to wild-type mice, female Asah2-/- mice showed an increase in several lipoproteins, including HDL. Our results provide insights into the genetic regulation of circulating lipoproteins, as well as mechanisms by which lipoprotein subfractions may affect cardiovascular disease risk in humans.

Genomic analysis and characterization of new loci associated with seed protein and oil content in soybeans

Breeding for increased protein without a reduction in oil content in soybeans [Glycine max (L.) Merr.] is a challenge for soybean breeders but an expected goal. Many efforts have been made to develop new soybean varieties with high yield in combination with desirable protein and/or oil traits. An elite line, R05-1415, was reported to be high yielding, high protein, and low oil. Several significant quantitative trait loci (QTL) for protein and oil were reported in this line, but many of them were unstable across environments or genetic backgrounds. Thus, a new study under multiple field environments using the Infinium BARCSoySNP6K BeadChips was conducted to detect and confirm stable genomic loci for these traits. Genetic analyses consistently detected a single major genomic locus conveying these two traits with remarkably high phenotypic variation explained (R2 ), varying between 24.2% and 43.5%. This new genomic locus is located between 25.0 and 26.7 Mb, distant from the previously reported QTL and did not overlap with other commonly reported QTL and the recently cloned gene Glyma.20G085100. Homolog analysis indicated that this QTL did not result from the paracentric chromosome inversion with an adjacent genomic fragment that harbors the reported QTL. The pleiotropic effect of this QTL could be a challenge for improving protein and oil simultaneously; however, a further study of four candidate genes with significant expressions in the seed developmental stages coupled with haplotype analysis may be able to pinpoint causative genes. The functionality and roles of these genes can be determined and characterized, which lay a solid foundation for the improvement of protein and oil content in soybeans.

Data-adaptive and pathway-based tests for association studies between somatic mutations and germline variations in human cancers

Cancer is a disease driven by a combination of inherited genetic variants and somatic mutations. Recently available large-scale sequencing data of cancer genomes have provided an unprecedented opportunity to study the interactions between them. However, previous studies on this topic have been limited by simple, low statistical power tests such as Fisher's exact test. In this paper, we design data-adaptive and pathway-based tests based on the score statistic for association studies between somatic mutations and germline variations. Previous research has shown that two single-nucleotide polymorphism (SNP)-set-based association tests, adaptive sum of powered score (aSPU) and data-adaptive pathway-based (aSPUpath) tests, increase the power in genome-wide association studies (GWASs) with a single disease trait in a case-control study. We extend aSPU and aSPUpath to multi-traits, that is, somatic mutations of multiple genes in a cohort study, allowing extensive information aggregation at both SNP and gene levels.p $p$ -values from different parameters assuming varying genetic architecture are combined to yield data-adaptive tests for somatic mutations and germline variations. Extensive simulations show that, in comparison with some commonly used methods, our data-adaptive somatic mutations/germline variations tests can be applied to multiple germline SNPs/genes/pathways, and generally have much higher statistical powers while maintaining the appropriate type I error. The proposed tests are applied to a large-scale real-world International Cancer Genome Consortium whole genome sequencing data set of 2583 subjects, detecting more significant and biologically relevant associations compared with the other existing methods on both gene and pathway levels. Our study has systematically identified the associations between various germline variations and somatic mutations across different cancer types, which potentially provides valuable utility for cancer risk prediction, prognosis, and therapeutics.

Epistasis between mutator alleles contributes to germline mutation spectra variability in laboratory mice

Maintaining germline genome integrity is essential and enormously complex. Although many proteins are involved in DNA replication, proofreading, and repair [1], mutator alleles have largely eluded detection in mammals. DNA replication and repair proteins often recognize sequence motifs or excise lesions at specific nucleotides. Thus, we might expect that the spectrum of de novo mutations - the frequencies of C>T, A>G, etc. - will differ between genomes that harbor either a mutator or wild-type allele. Previously, we used quantitative trait locus mapping to discover candidate mutator alleles in the DNA repair gene Mutyh that increased the C>A germline mutation rate in a family of inbred mice known as the BXDs [2,3]. In this study we developed a new method to detect alleles associated with mutation spectrum variation and applied it to mutation data from the BXDs. We discovered an additional C>A mutator locus on chromosome 6 that overlaps Ogg1, a DNA glycosylase involved in the same base-excision repair network as Mutyh [4]. Its effect depended on the presence of a mutator allele near Mutyh, and BXDs with mutator alleles at both loci had greater numbers of C>A mutations than those with mutator alleles at either locus alone. Our new methods for analyzing mutation spectra reveal evidence of epistasis between germline mutator alleles and may be applicable to mutation data from humans and other model organisms.

Amoeba predation of Cryptococcus: A quantitative and population genomic evaluation of the accidental pathogen hypothesis

The "Amoeboid Predator-Fungal Animal Virulence Hypothesis" posits that interactions with environmental phagocytes shape the evolution of virulence traits in fungal pathogens. In this hypothesis, selection to avoid predation by amoeba inadvertently selects for traits that contribute to fungal escape from phagocytic immune cells. Here, we investigate this hypothesis in the human fungal pathogens Cryptococcus neoformans and Cryptococcus deneoformans. Applying quantitative trait locus (QTL) mapping and comparative genomics, we discovered a cross-species QTL region that is responsible for variation in resistance to amoeba predation. In C. neoformans, this same QTL was found to have pleiotropic effects on melanization, an established virulence factor. Through fine mapping and population genomic comparisons, we identified the gene encoding the transcription factor Bzp4 that underlies this pleiotropic QTL and we show that decreased expression of this gene reduces melanization and increases susceptibility to amoeba predation. Despite the joint effects of BZP4 on amoeba resistance and melanin production, we find no relationship between BZP4 genotype and escape from macrophages or virulence in murine models of disease. Our findings provide new perspectives on how microbial ecology shapes the genetic architecture of fungal virulence, and suggests the need for more nuanced models for the evolution of pathogenesis that account for the complexities of both microbe-microbe and microbe-host interactions.

Genetic analysis of basal stalk rot resistance introgressed from wild Helianthus petiolaris into cultivated sunflower (Helianthus annuus L.) using an advanced backcross population

Introduction

Sclerotinia sclerotiorum is a serious pathogen causing severe basal stalk rot (BSR) disease on cultivated sunflower (Helianthus annuus L.) that leads to significant yield losses due to insufficient resistance. The wild annual sunflower species H. petiolaris, commonly known as prairie sunflower is known for its resistance against this pathogen. Sunflower resistance to BSR is quantitative and determined by many genes with small effects on the resistance phenotype. The objective of this study was to identify loci governing BSR resistance derived from H. petiolaris using a quantitative trait loci (QTL) mapping approach.

Methods

BSR resistance among lines of an advanced backcross population (AB-QTL) with 174 lines developed from a cross of inbred line HA 89 with H. petiolaris PI 435843 was determined in the field during 2017-2019, and in the greenhouse in 2019. AB-QTL lines and the HA 89 parent were genotyped using genotyping-by-sequencing and a genetic linkage map was developed spanning 997.51 cM and using 1,150 SNP markers mapped on 17 sunflower chromosomes.

Results and discussion

Highly significant differences (p<0.001) for BSR response among AB-QTL lines were observed disease incidence (DI) in all field seasons, as well as disease rating (DR) and area under the disease progress curve (AUDPC) in the greenhouse with a moderately high broad-sense heritability (H 2) of 0.61 for the tested resistance parameters. A total of 14 QTL associated with BSR resistance were identified on nine chromosomes, each explaining a proportion of the phenotypic variation ranging from 3.5% to 28.1%. Of the 14 QTL, eight were detected for BSR resistance in the field and six were detected under greenhouse conditions. Alleles conferring increased BSR resistance were contributed by the H. petiolaris parent at 11 of the 14 QTL.

Integrating Association Mapping, Linkage Mapping, Fine Mapping with RNA Seq Conferring Seedling Vigor Improvement for Successful Crop Establishment in Deep Sown Direct-Seeded Rice

BACKGROUND

Ongoing large-scale shift towards direct seeded rice (DSR) necessitates a convergence of breeding and genetic approaches for its sustenance and harnessing natural resources and environmental benefits. Improving seedling vigour remains key objective for breeders working with DSR. The present study aims to understand the genetic control of seedling vigour in deep sown DSR. Combined genome-wide association mapping, linkage mapping, fine mapping, RNA-sequencing to identify candidate genes and validation of putative candidate genes were performed in the present study.

RESULTS

Significant phenotypic variations were observed among genotypes in both F3:4:5 and BC2F2:3 populations. The mesocotyl length showed significant positive correlation with %germination, root and shoot length. The 881 kb region on chromosome 7 reported to be associated with mesocotyl elongation. RNA-seq data and RT-PCR results identified and validated seven potential candidate genes. The four promising introgression lines free from linkage drag and with longer mesocotyl length, longer root length, semi-dwarf plant height have been identified.

CONCLUSION

The study will provide rice breeders (1) the pre breeding material in the form of anticipated DSR adapted introgression lines possessing useful traits and alleles improving germination under deep sown DSR field conditions (2) the base for the studies involving functional characterization of candidate genes. The development and utilization of improved introgression lines and molecular markers may play an important role in genomics-assisted breeding (GAB) during the pyramiding of valuable genes providing adaptation to rice under DSR. Our results offer a robust and reliable package that can contribute towards enhancing genetic gains in direct seeded rice breeding programs.

Identification of quantitative trait loci for drought tolerance in Bukoba/Kijivu Andean mapping population of common bean

KEY MESSAGE

Quantitative Trait Loci "hotspots" for drought tolerance were identified on chromosomes Pv06, Pv07 and Pv10 of common bean. Drought is a major production constraint of common bean (Phaseolus vulgaris L.) worldwide. The objective of this study was to identify the Quantitative Trait Loci (QTL) for drought tolerance in an Andean population of Recombinant Inbred Lines (RILs). A total of 155 F5:7 RILs derived from a cross between Kijivu (drought tolerant) and Bukoba (drought susceptible) were evaluated for drought tolerance in field and pot experiments. Four field experiments were conducted at three locations in Zambia in 2020 and 2021. All field trials were conducted in the dry season under irrigation. The 155 RILs were genotyped with 11,292 SNPs, and composite interval mapping was conducted to identify QTL for drought tolerance. Seed yield for Kijivu under drought stress was consistently higher than for Bukoba across all four field trials. A total of 60 QTL were identified for morphological, agronomic, and physiological traits under drought stress and non-stress conditions. However, the majority of these QTL were specific to drought stress. QTL "hotspots" for drought tolerance were identified on chromosomes Pv06, Pv07, and Pv10. Extensive co-localizations for agronomic and morpho-physiological traits under drought stress were observed at the three drought-tolerance QTL hotspots. Additionally, these three QTL hotspots overlapped with previously identified QTL for drought tolerance, while several others identified QTL are novel. The three identified QTL hotspots could be used in future marker-assisted selection for drought tolerance in common bean.

Machine learning reveals genetic modifiers of the immune microenvironment of cancer

Heritability in the immune tumor microenvironment (iTME) has been widely observed yet remains largely uncharacterized. Here, we developed a machine learning approach to map iTME modifiers within loci from genome-wide association studies (GWASs) for breast cancer (BrCa) incidence. A random forest model was trained on a positive set of immune-oncology (I-O) targets, and then used to assign I-O target probability scores to 1,362 candidate genes in linkage disequilibrium with 155 BrCa GWAS loci. Cluster analysis of the most probable candidates revealed two subfamilies of genes related to effector functions and adaptive immune responses, suggesting that iTME modifiers impact multiple aspects of anticancer immunity. Two of the top ranking BrCa candidates, LSP1 and TLR1, were orthogonally validated as iTME modifiers using BrCa patient biopsies and comparative mapping studies, respectively. Collectively, these data demonstrate a robust and flexible framework for functionally fine-mapping GWAS risk loci to identify translatable therapeutic targets.

Towards carotenoid biofortification in wheat: identification of XAT-7A1, a multicopy tandem gene responsible for carotenoid esterification in durum wheat

Yellow pigment content, mainly due to the accumulation of carotenoids, is a quality trait in durum wheat grain as it confers the bright yellow color to pasta preferred by consumers. Also, carotenoids are essential nutrients exerting important biological functions in human health. Consequently, biofortification strategies have been developed in many crops to increase carotenoid content. In this context, carotenoid esterification is emerging as a new breeding target for wheat biofortification, as carotenoid esters have been found to promote both carotenoid accumulation and stability. Until recently, no carotenoid esters have been identified in significant proportions in durum wheat grains, and interspecific breeding programs have been started to transfer esterification ability from common wheat and Hordeum chilense.In this work, XAT-7A1 is identified as the gene responsible for carotenoid esterification in durum wheat. Sequencing, copy number variation and mapping results show that XAT-7A1 is organized as tandem or proximal GDSL esterase/lipase copies in chromosome 7A. Three XAT-7A1 haplotypes are described: Type 1 copies, associated with high levels of carotenoid esters (diesters and monoesters) production and high expression in grain development; Type 2 copies, present in landraces with low levels of carotenoid esters (monoesters) or no esters; and Type 3 copies, without the signal peptide, resulting in zero-ester phenotypes.The identification of XAT-7A1 is a necessary step to make the carotenoid esterification ability available for durum and bread wheat breeding, which should be focused on the Type 1 XAT-7A1 haplotype, which may be assessed as a single gene since XAT-7A1 copies are inherited together.

Thresholding Gini variable importance with a single-trained random forest: An empirical Bayes approach

Random forests (RFs) are a widely used modelling tool capable of feature selection via a variable importance measure (VIM), however, a threshold is needed to control for false positives. In the absence of a good understanding of the characteristics of VIMs, many current approaches attempt to select features associated to the response by training multiple RFs to generate statistical power via a permutation null, by employing recursive feature elimination, or through a combination of both. However, for high-dimensional datasets these approaches become computationally infeasible. In this paper, we present RFlocalfdr, a statistical approach, built on the empirical Bayes argument of Efron, for thresholding mean decrease in impurity (MDI) importances. It identifies features significantly associated with the response while controlling the false positive rate. Using synthetic data and real-world data in health, we demonstrate that RFlocalfdr has equivalent accuracy to currently published approaches, while being orders of magnitude faster. We show that RFlocalfdr can successfully threshold a dataset of 106 datapoints, establishing its usability for large-scale datasets, like genomics. Furthermore, RFlocalfdr is compatible with any RF implementation that returns a VIM and counts, making it a versatile feature selection tool that reduces false discoveries.

Genetics of evolved load resistance in the skeletons of unusually large mice from Gough Island

A primary function of the skeleton is to resist the loads imparted by body weight. Genetic analyses have identified genomic regions that contribute to differences in skeletal load resistance between laboratory strains of mice, but these studies are usually restricted to 1 or 2 bones and leave open the question of how load resistance evolves in natural populations. To address these challenges, we examined the genetics of bone structure using the largest wild house mice on record, which live on Gough Island (GI). We measured structural traits connected to load resistance in the femur, tibia, scapula, humerus, radius, ulna, and mandible of GI mice, a smaller-bodied reference strain from the mainland, and 760 of their F2s. GI mice have bone geometries indicative of greater load resistance abilities but show no increase in bone mineral density compared to the mainland strain. Across traits and bones, we identified a total of 153 quantitative trait loci (QTL) that span all but one of the autosomes. The breadth of QTL detection ranges from a single bone to all 7 bones. Additive effects of QTL are modest. QTL for bone structure show limited overlap with QTL for bone length and width and QTL for body weight mapped in the same cross, suggesting a distinct genetic architecture for load resistance. Our findings provide a rare genetic portrait of the evolution of load resistance in a natural population with extreme body size.

QTL mapping and stability analysis of trichome density in zucchini (Cucurbita pepo L.)

Trichomes provide an excellent model for studying cell differentiation and proliferation. The aboveground tissues of plants with long dense trichomes (LDTs) can cause skin itching in people working in a zucchini field, in which management, pollination, and fruit harvesting are difficult. In this study, an F2 population was constructed with the LDT inbred line "16" and the sparse micro trichome (SMT) inbred line "63" for QTL analysis of type I and II trichome density. Two QTLs were identified on chromosomes 3 and 15 using the QTL-seq method. Additionally, 191 InDel markers were developed on 20 chromosomes, a genetic map was constructed for QTL mapping, and three QTLs were identified on chromosomes 3, 6, and 15. Two QTLs, CpTD3.1 and CpTD15.1, were identified in both QTL-seq and genetic map-based QTL analyses, and CpTD15.1 was the major-effect QTL. The stability of CpTD3.1 and CpTD15.1 was confirmed using data from F2 plants under different environmental conditions. The major-effect QTL CpTD15.1 was located between markers chr15-4991349 and chr15-5766791, with a physical distance of 775.44 kb, and explained 12.71%-29.37% of the phenotypic variation observed in the three environments. CpTD3.1 was located between markers chr3-218350 and chr3-2891236, in a region with a physical distance of 2,672.89 kb, and explained 5.00%-10.64% of the phenotypic variation observed in the three environments. The functional annotations of the genes within the CpTD15.1 region were predicted, and five genes encoding transcription factors regulating trichome development were selected. Cp4.1LG15g04400 encoded zinc finger protein (ZFP) and harbored nonsynonymous SNPs in the conserved ring finger domain between the two parental lines. There were significant differences in Cp4.1LG15g04400 expression between "16" and "63", and a similar pattern was found between germplasm resources of LDT lines and SMT lines. It was presumed that Cp4.1LG15g04400 might regulate trichome density in zucchini. These results lay a foundation for better understanding the density of multicellular nonglandular trichomes and the regulatory mechanism of trichome density in zucchini.

Linked candidate genes of different functions for white mold resistance in common bean (Phaseolus vulgaris L) are identified by multiple QTL mapping approaches

White mold (WM) is a major disease in common bean (Phaseolus vulgaris L.), and its complex quantitative genetic control limits the development of WM resistant cultivars. WM2.2, one of the nine meta-QTL with a major effect on WM tolerance, explains up to 35% of the phenotypic variation and was previously mapped to a large genomic interval on Pv02. Our objective was to narrow the interval of this QTL using combined approach of classic QTL mapping and QTL-based bulk segregant analysis (BSA), and confirming those results with Khufu de novo QTL-seq. The phenotypic and genotypic data from two RIL populations, 'Raven'/I9365-31 (R31) and 'AN-37'/PS02-029C-20 (Z0726-9), were used to select resistant and susceptible lines to generate subpopulations for bulk DNA sequencing. The QTL physical interval was determined by considering overlapping interval of the identified QTL or peak region in both populations by three independent QTL mapping analyses. Our findings revealed that meta-QTL WM2.2 consists of three regions, WM2.2a (4.27-5.76 Mb; euchromatic), WM 2.2b (12.19 to 17.61 Mb; heterochromatic), and WM2.2c (23.01-25.74 Mb; heterochromatic) found in both populations. Gene models encoding for gibberellin 2-oxidase 8, pentatricopeptide repeat, and heat-shock proteins are the likely candidate genes associated with WM2.2a resistance. A TIR-NBS-LRR class of disease resistance protein (Phvul.002G09200) and LRR domain containing family proteins are potential candidate genes associated with WM2.2b resistance. Nine gene models encoding disease resistance protein [pathogenesis-related thaumatin superfamily protein and disease resistance-responsive (dirigent-like protein) family protein etc] found within the WM2.2c QTL interval are putative candidate genes. WM2.2a region is most likely associated with avoidance mechanisms while WM2.2b and WM2.2c regions trigger physiological resistance based on putative candidate genes.

Application of SVR-Mediated GWAS for Identification of Durable Genetic Regions Associated with Soybean Seed Quality Traits

Soybean (Glycine max L.) is an important food-grade strategic crop worldwide because of its high seed protein and oil contents. Due to the negative correlation between seed protein and oil percentage, there is a dire need to detect reliable quantitative trait loci (QTL) underlying these traits in order to be used in marker-assisted selection (MAS) programs. Genome-wide association study (GWAS) is one of the most common genetic approaches that is regularly used for detecting QTL associated with quantitative traits. However, the current approaches are mainly focused on estimating the main effects of QTL, and, therefore, a substantial statistical improvement in GWAS is required to detect associated QTL considering their interactions with other QTL as well. This study aimed to compare the support vector regression (SVR) algorithm as a common machine learning method to fixed and random model circulating probability unification (FarmCPU), a common conventional GWAS method in detecting relevant QTL associated with soybean seed quality traits such as protein, oil, and 100-seed weight using 227 soybean genotypes. The results showed a significant negative correlation between soybean seed protein and oil concentrations, with heritability values of 0.69 and 0.67, respectively. In addition, SVR-mediated GWAS was able to identify more relevant QTL underlying the target traits than the FarmCPU method. Our findings demonstrate the potential use of machine learning algorithms in GWAS to detect durable QTL associated with soybean seed quality traits suitable for genomic-based breeding approaches. This study provides new insights into improving the accuracy and efficiency of GWAS and highlights the significance of using advanced computational methods in crop breeding research.

Naturally segregating genetic variants contribute to thermal tolerance in a D. melanogaster model system

Thermal tolerance is a fundamental physiological complex trait for survival in many species. For example, everyday tasks such as foraging, finding a mate, and avoiding predation, are highly dependent on how well an organism can tolerate extreme temperatures. Understanding the general architecture of the natural variants of the genes that control this trait is of high importance if we want to better comprehend how this trait evolves in natural populations. Here, we take a multipronged approach to further dissect the genetic architecture that controls thermal tolerance in natural populations using the Drosophila Synthetic Population Resource (DSPR) as a model system. First, we used quantitative genetics and Quantitative Trait Loci (QTL) mapping to identify major effect regions within the genome that influences thermal tolerance, then integrated RNA-sequencing to identify differences in gene expression, and lastly, we used the RNAi system to 1) alter tissue-specific gene expression and 2) functionally validate our findings. This powerful integration of approaches not only allows for the identification of the genetic basis of thermal tolerance but also the physiology of thermal tolerance in a natural population, which ultimately elucidates thermal tolerance through a fitness-associated lens.

Molecular Mapping to Discover Reliable Salinity-Resilient QTLs from the Novel Landrace Akundi in Two Bi-Parental Populations Using SNP-Based Genome-Wide Analysis in Rice

Achieving high-yield potential is always the ultimate objective of any breeding program. However, various abiotic stresses such as salinity, drought, cold, flood, and heat hampered rice productivity tremendously. Salinity is one of the most important abiotic stresses that adversely affect rice grain yield. The present investigation was undertaken to dissect new genetic loci, which are responsible for salt tolerance at the early seedling stage in rice. A bi-parental mapping population (F_2:3) was developed from the cross between BRRI dhan28/Akundi, where BRRI dhan28 (BR28) is a salt-sensitive irrigated (boro) rice mega variety and Akundi is a highly salinity-tolerant Bangladeshi origin indica rice landrace that is utilized as a donor parent. We report reliable and stable QTLs for salt tolerance from a common donor (Akundi) irrespective of two different genetic backgrounds (BRRI dhan49/Akundi and BRRI dhan28/Akundi). A robust 1k-Rice Custom Amplicon (1k-RiCA) SNP marker genotyping platform was used for genome-wide analysis of this bi-parental population. After eliminating markers with high segregation distortion, 886 polymorphic SNPs built a genetic linkage map covering 1526.5 cM of whole rice genome with an average SNP density of 1.72 cM for the 12 genetic linkage groups. A total of 12 QTLs for nine different salt tolerance-related traits were identified using QGene and inclusive composite interval mapping of additive and dominant QTL (ICIM-ADD) under salt stress on seven different chromosomes. All of these 12 new QTLs were found to be unique, as no other map from the previous study has reported these QTLs in the similar chromosomal location and found them different from extensively studied Saltol, SKC1, OsSalT, and salT locus. Twenty-eight significant digenic/epistatic interactions were identified between chromosomal regions linked to or unlinked to QTLs. Akundi acts like a new alternate donor source of salt tolerance except for other usually known donors such as Nona Bokra, Pokkali, Capsule, and Hasawi used in salt tolerance genetic analysis and breeding programs worldwide, including Bangladesh. Integration of the seven novel, reliable, stable, and background independent salinity-resilient QTLs (qSES1, qSL1, qRL1, qSUR1, qSL8, qK8, qK1) reported in this investigation will expedite the cultivar development that is highly tolerant to salt stress.

QTL mapping and identification of candidate genes for cold tolerance at the germination stage in wild rice

BACKGROUND

Cold damage stress significantly affects rice growth (germination and seedling) and causes serious losses in yield in temperate and high-altitude areas around the globe.

OBJECTIVE

This study aimed to explore the cold tolerance (CT) locus of rice and create new cold-tolerant germplasm. We constructed a chromosome segment substitution line (CSSL) with strong CT and fine mapped quantitative trait loci (QTLs) associated with CT by performing the whole-genome resequencing of CSSL with phenotypes under cold treatment.

METHODS

A chromosome CSSL, including 271 lines from a cross between the cold-tolerant wild rice Y11 (Oryza rufipogon Griff.) and the cold-sensitive rice variety GH998, was developed to map QTLs conferring CT at the germination stage. The whole-genome resequencing was performed on CSSL for mapping QTLs of associated with CT at the germination stage.

RESULTS

A high-density linkage map of the CSSLs was developed using the whole-genome resequencing of 1484 bins. The QTL analysis using 615,466 single-nucleotide polymorphisms (SNPs) led to the identification of 2 QTLs related to germination rate at low-temperature on chromosome 8 (qCTG-8) and chromosome 11 (qCTG-11). The qCTG-8 and qCTG-11 explained 14.55% and 14.31% of the total phenotypic variation, respectively. We narrowed down qCTG-8 and qCTG-11 to 195.5 and 78.83-kb regions, respectively. The expression patterns of important candidate genes in different tissues, and of RNA-sequencing (RNA-seq) in CSSLs, were identified based on gene sequences in qCTG-8 and qCTG-11 cold-induced expression analysis. LOC_Os08g01120 and LOC_Os08g01390 were identified as candidate genes in qCTG-8, and LOC_Os11g32880 was identified as a candidate gene in qCTG-11.

CONCLUSIONS

This study demonstrated a general method that could be used to identify useful loci and genes in wild rice and aid in the future cloning of candidate genes of qCTG-8 and qCTG-11. The CSSLs with strong CT were supported for breeding cold-tolerant rice varieties.

Lipidomic QTL in Diversity Outbred mice identifies a novel function for α/β hydrolase domain 2 (Abhd2) as an enzyme that metabolizes phosphatidylcholine and cardiolipin

We and others have previously shown that genetic association can be used to make causal connections between gene loci and small molecules measured by mass spectrometry in the bloodstream and in tissues. We identified a locus on mouse chromosome 7 where several phospholipids in liver showed strong genetic association to distinct gene loci. In this study, we integrated gene expression data with genetic association data to identify a single gene at the chromosome 7 locus as the driver of the phospholipid phenotypes. The gene encodes α/β-hydrolase domain 2 (Abhd2), one of 23 members of the ABHD gene family. We validated this observation by measuring lipids in a mouse with a whole-body deletion of Abhd2. The Abhd2KO mice had a significant increase in liver levels of phosphatidylcholine and phosphatidylethanolamine. Unexpectedly, we also found a decrease in two key mitochondrial lipids, cardiolipin and phosphatidylglycerol, in male Abhd2KO mice. These data suggest that Abhd2 plays a role in the synthesis, turnover, or remodeling of liver phospholipids.

Identification of a Major-Effect Quantitative Trait Loci Associated with Begomovirus Resistance in Cucurbita moschata

Begomoviruses, viz. squash leaf curl China virus and tomato leaf curl New Delhi virus causative diseases are major concerns of quantitative and qualitative losses in pumpkin (Cucurbita moschata) worldwide. Punjab Agricultural University (PAU) in India has identified a resistant source (PVR-1343) against mixed infection (MI-Sq/To) of these begomoviruses. Introgression of resistance in diverse genetic backgrounds requires the identification of quantitative trait loci (QTLs) associated with MI-Sq/To resistance. Phenotyping of 229 F_2:3 progenies derived from the PVR-1343 × P-135 cross revealed digenic recessive inheritance against MI-Sq/To resistance in PVR-1343. To identify the genomic region, resistant and susceptible bulks were subjected to whole-genome resequencing along with their parents. The whole-genome resequence analysis of parents and bulks using QTLseq/QTLseqr approaches identified an overlapping 1.52 Mb region on chromosome 7 (qMI-Sq/To7.1), while chromosomal region spanning 0.87 Mb on chromosome17 (qMI-Sq/To17.1) was additionally identified by QTLseqr. However, the highest peak value on chromosome 7 with three algorithms {G', ∆(SNP-index) and -log₁₀ (P value)} highlighted the major contribution of qMI-Sq/To7.1 in MI-Sq/To resistance. Nine polymorphic SNPs identified within the highly significant qMI-Sq/To7.1 region were converted into KASP markers. KASP genotyping of F₂ individuals narrowed down the qMI-Sq/To7.1 interval to 103 kb region flanked by two markers, Cmo3914729 and Cmo4018182, which contained 16 annotated genes and accounted for 59.84% of phenotypic variation. The Cmo4018182 KASP marker accurately predicted disease reaction in 91% of diverse Cucurbita genotypes and showed nonsynonym substitutions in the coding region of putative candidate SYNTAXIN-121 gene. These findings pave the way for marker-assisted breeding and elucidating the underlying mechanism of begomovirus resistance in C. moschata.

Identification of Candidate QTLs and Genes for Ear Diameter by Multi-Parent Population in Maize

Ear diameter (ED) is a critical component of grain yield (GY) in maize (Zea mays L.). Studying the genetic basis of ED in maize is of great significance in enhancing maize GY. Against this backdrop, this study was framed to (1) map the ED-related quantitative trait locus (QTL) and SNPs associated with ED; and (2) identify putative functional genes that may affect ED in maize. To accomplish this, an elite maize inbred line, Ye107, which belongs to the Reid heterotic group, was used as a common parent and crossed with seven elite inbred lines from three different heterotic groups (Suwan1, Reid, and nonReid) that exhibited abundant genetic variation in ED. This led to the construction of a multi-parent population consisting of 1215 F₇ recombinant inbred lines (F₇RILs). A genome-wide association study (GWAS) and linkage analysis were then conducted for the multi-parent population using 264,694 high-quality SNPs generated via the genotyping-by-sequencing method. Our study identified a total of 11 SNPs that were significantly associated with ED through the GWAS, and three QTLs were revealed by the linkage analysis for ED. The major QTL on chromosome 1 was co-identified in the region by the GWAS at SNP_143985532. SNP_143985532, located upstream of the Zm00001d030559 gene, encodes a callose synthase that is expressed in various tissues, with the highest expression level in the maize ear primordium. Haplotype analysis indicated that the haplotype B (allele AA) of Zm00001d030559 was positively correlated with ED. The candidate genes and SNPs identified in this study provide crucial insights for future studies on the genetic mechanism of maize ED formation, cloning of ED-related genes, and genetic improvement of ED. These results may help develop important genetic resources for enhancing maize yield through marker-assisted breeding.

Multi-parental fungal mapping population study to detect genomic regions associated with Pyrenophora teres f. teres virulence

In recent years multi-parental mapping populations (MPPs) have been widely adopted in many crops to detect quantitative trait loci (QTLs) as this method can compensate for the limitations of QTL analyses using bi-parental mapping populations. Here we report the first multi-parental nested association mapping (MP-NAM) population study used to detect genomic regions associated with host-pathogenic interactions. MP-NAM QTL analyses were conducted on 399 Pyrenophora teres f. teres individuals using biallelic, cross-specific and parental QTL effect models. A bi-parental QTL mapping study was also conducted to compare the power of QTL detection between bi-parental and MP-NAM populations. Using MP-NAM with 399 individuals detected a maximum of eight QTLs with a single QTL effect model whilst only a maximum of five QTLs were detected with an individual bi-parental mapping population of 100 individuals. When reducing the number of isolates in the MP-NAM to 200 individuals the number of QTLs detected remained the same for the MP-NAM population. This study confirms that MPPs such as MP-NAM populations can be successfully used in detecting QTLs in haploid fungal pathogens and that the power of QTL detection with MPPs is greater than with bi-parental mapping populations.

Genetic Dissection of Extreme Seed-Flooding Tolerance in a Wild Soybean PI342618B by Linkage Mapping and Candidate Gene Analysis

Seed-flooding stress is one of major abiotic constraints that adversely affects soybean production worldwide. Identifying tolerant germplasms and revealing the genetic basis of seed-flooding tolerance are imperative goals for soybean breeding. In the present study, high-density linkage maps of two inter-specific recombinant inbred line (RIL) populations, named NJIRNP and NJIR4P, were utilized to identify major quantitative trait loci (QTLs) for seed-flooding tolerance using three parameters viz., germination rate (GR), normal seedling rate (NSR), and electrical conductivity (EC). A total of 25 and 18 QTLs were detected by composite interval mapping (CIM) and mixed-model-based composite interval mapping (MCIM), respectively, and 12 common QTLs were identified through both methods. All favorable alleles for the tolerance are notably from the wild soybean parent. Moreover, four digenic epistatic QTL pairs were identified, and three of them showed no main effects. In addition, the pigmented soybean genotypes exhibited high seed-flooding tolerance compared with yellow seed coat genotypes in both populations. Moreover, out of five identified QTLs, one major region containing multiple QTLs associated with all three traits was identified on Chromosome 8, and most of the QTLs within this hotspot were major loci (R² > 10) and detectable in both populations and multiple environments. Based on the gene expression and functional annotation information, 10 candidate genes from QTL "hotspot 8-2" were screened for further analysis. Furthermore, the results of qRT-PCR and sequence analysis revealed that only one gene, GmDREB2 (Glyma.08G137600), was significantly induced under flooding stress and displayed a TTC tribasic insertion mutation of the nucleotide sequence in the tolerant wild parent (PI342618B). GmDREB2 encodes an ERF transcription factor, and the subcellular localization analysis using green fluorescent protein (GFP) revealed that GmDREB2 protein was localized in the nucleus and plasma membrane. Furthermore, overexpression of GmDREB2 significantly promoted the growth of soybean hairy roots, which might indicate its critical role in seed-flooding stress. Thus, GmDREB2 was considered as the most possible candidate gene for seed-flooding tolerance.

QTL mapping and improvement of pre-harvest sprouting resistance using japonica weedy rice

The stability of cultivation and production in terms of crop yield has been threatened by climate change due to global warming. Pre-harvest sprouting (PHS) is a threat to crops, especially staple foods, including rice, because of reductions in yield and quality. To address the problem of precocious germination before harvest, we performed quantitative trait loci (QTL) analysis for PHS using F₈ RILs populations derived from japonica weedy rice in Korea. QTL analysis revealed that two stable QTLs, qPH7 and qPH2, associated with PHS resistance were identified on chromosomes 7 and 2, respectively, explaining approximately 38% of the phenotypic variation. The QTL effect in the tested lines significantly decreased the degree of PHS, based on the number of QTLs included. Through fine mapping for main QTL qPH7, the region for the PHS was found to be anchored within 23.575-23.785 Mbp on chromosome 7 using 13 cleaved amplified sequence (CAPS) markers. Among 15 open reading frames (ORFs) within the detected region, one ORF, Os07g0584366, exhibited upregulated expression in the resistant donor, which was approximately nine times higher than that of susceptible japonica cultivars under PHS-inducing conditions. Japonica lines with QTLs related to PHS resistance were developed to improve the characteristics of PHS and design practical PCR-based DNA markers for marker-assisted backcrosses of many other PHS-susceptible japonica cultivars.

Mutant Tof11 alleles are highly accumulated in early planting-adaptable Japanese summer type soybeans

To avoid crop failure because of climate change, soybean (Glycine max (L.) Merrill) cultivars adaptable to early planting are required in western Japan. Because current Japanese cultivars may not be adaptable, genetic resources with high early-planting adaptability, and their genetic information must be developed. In the present study, summer type (ST) soybeans developed for early planting were used as plant materials. We examined their phenological characteristics and short reproductive period as an indicator of early planting adaptability and performed genetic studies. Biparental quantitative trait loci (QTL) analysis of a representative ST cultivar revealed a principal QTL for the reproductive period duration on chromosome 11. The results of resequencing analysis suggested that circadian clock-related Tof11 (soybean orthologue of PRR3) is a candidate QTL. Additionally, all 25 early planting-adaptable germplasms evaluated in this study possessed mutant alleles in Tof11, whereas 15 conventional cultivars only had wild-type alleles. These results suggest that mutant alleles in Tof11 are important genetic factors in the high adaptability to early planting of these soybeans, and thus, these alleles were acquired and accumulated in the ST soybean population.