Empirical threshold values for quantitative trait mapping

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.

Major genetic locus with pleiotropism determined seed-related traits in cultivated and wild soybeans

KEY MESSAGE

Here, a novel pleiotropic QTL qSS14 simultaneously regulating four seed size traits and two consistently detected QTLs qSW17 and qSLW02 were identified across multiple years. Seed-related traits were the key agronomic traits that have been artificially selected during the domestication of wild soybean. Identifying the genetic loci and genes that regulate seed size could clarify the genetic variations in seed-related traits and provide novel insights into high-yield soybean breeding. In this study, we used a high-density genetic map constructed by F₁₀ RIL populations from a cross between Glycine max and Glycine soja to detect additive QTLs for seven seed-related traits over the last three years. As a result, we identified one novel pleiotropic QTL, qSS14, that simultaneously controlled four seed size traits (100-seed weight, seed length, seed width, and seed thickness) and two consistently detected QTLs, qSW17, and qSLW02, in multiple years of phenotypic data. Furthermore, we predicted two, two and three candidate genes within these three critical loci based on the parental resequencing data and gene function annotations. And the relative expression of four candidate genes GLYMA_14G155100, GLYMA_17G061000, GLYMA_02G273100, and GLYMA_02G273300 showed significant differences among parents and the extreme materials through qRT-PCR analysis. These findings could facilitate the determination of beneficial genes in wild soybean and contribute to our understanding of the soybean domestication process.

Evaluating the Genetic Background Effect on Dissecting the Genetic Basis of Kernel Traits in Reciprocal Maize Introgression Lines

Maize yield is mostly determined by its grain size. Although numerous quantitative trait loci (QTL) have been identified for kernel-related traits, the application of these QTL in breeding programs has been strongly hindered because the populations used for QTL mapping are often different from breeding populations. However, the effect of genetic background on the efficiency of QTL and the accuracy of trait genomic prediction has not been fully studied. Here, we used a set of reciprocal introgression lines (ILs) derived from 417F × 517F to evaluate how genetic background affects the detection of QTLassociated with kernel shape traits. A total of 51 QTL for kernel size were identified by chromosome segment lines (CSL) and genome-wide association studies (GWAS) methods. These were subsequently clustered into 13 common QTL based on their physical position, including 7 genetic-background-independent and 6 genetic-background-dependent QTL, respectively. Additionally, different digenic epistatic marker pairs were identified in the 417F and 517F ILs. Therefore, our results demonstrated that genetic background strongly affected not only the kernel size QTL mapping via CSL and GWAS but also the genomic prediction accuracy and epistatic detection, thereby enhancing our understanding of how genetic background affects the genetic dissection of grain size-related traits.

Genetic dissection of the pluripotent proteome through multi-omics data integration

Genetic background drives phenotypic variability in pluripotent stem cells (PSCs). Most studies to date have used transcript abundance as the primary molecular readout of cell state in PSCs. We performed a comprehensive proteogenomics analysis of 190 genetically diverse mouse embryonic stem cell (mESC) lines. The quantitative proteome is highly variable across lines, and we identified pluripotency-associated pathways that were differentially activated in the proteomics data that were not evident in transcriptome data from the same lines. Integration of protein abundance to transcript levels and chromatin accessibility revealed broad co-variation across molecular layers as well as shared and unique drivers of quantitative variation in pluripotency-associated pathways. Quantitative trait locus (QTL) mapping localized the drivers of these multi-omic signatures to genomic hotspots. This study reveals post-transcriptional mechanisms and genetic interactions that underlie quantitative variability in the pluripotent proteome and provides a regulatory map for mESCs that can provide a basis for future mechanistic studies.

Which mouse multiparental population is right for your study? The Collaborative Cross inbred strains, their F1 hybrids, or the Diversity Outbred population

Multiparental populations (MPPs) encompass greater genetic diversity than traditional experimental crosses of two inbred strains, enabling broader surveys of genetic variation underlying complex traits. Two such mouse MPPs are the Collaborative Cross (CC) inbred panel and the Diversity Outbred (DO) population, which are descended from the same eight inbred strains. Additionally, the F1 intercrosses of CC strains (CC-RIX) have been used and enable study designs with replicate outbred mice. Genetic analyses commonly used by researchers to investigate complex traits in these populations include characterizing how heritable a trait is, i.e. its heritability, and mapping its underlying genetic loci, i.e. its quantitative trait loci (QTLs). Here we evaluate the relative merits of these populations for these tasks through simulation, as well as provide recommendations for performing the quantitative genetic analyses. We find that sample populations that include replicate animals, as possible with the CC and CC-RIX, provide more efficient and precise estimates of heritability. We report QTL mapping power curves for the CC, CC-RIX, and DO across a range of QTL effect sizes and polygenic backgrounds for samples of 174 and 500 mice. The utility of replicate animals in the CC and CC-RIX for mapping QTLs rapidly decreased as traits became more polygenic. Only large sample populations of 500 DO mice were well-powered to detect smaller effect loci (7.5-10%) for highly complex traits (80% polygenic background). All results were generated with our R package musppr, which we developed to simulate data from these MPPs and evaluate genetic analyses from user-provided genotypes.

Genetic mapping of some key plant architecture traits in Brassica juncea using a doubled haploid population derived from a cross between two distinct lines: vegetable type Tumida and oleiferous Varuna

Genetic mapping of some key plant architectural traits in a vegetable type and an oleiferous B. juncea cross revealed QTL and candidate genes for breeding more productive ideotypes. Brassica juncea (AABB, 2n = 36), commonly called mustard, is an allopolyploid crop of recent origin but contains considerable morphological and genetic variation. An F₁-derived doubled haploid population developed from a cross between an Indian oleiferous line, Varuna, and a Chinese stem type vegetable mustard, Tumida showed significant variability for some key plant architectural traits-four stem strength-related traits, stem diameter (Dia), plant height (Plht), branch initiation height (Bih), number of primary branches (Pbr), and days to flowering (Df). Multi-environment QTL analysis identified twenty Stable QTL for the above-mentioned nine plant architectural traits. Though Tumida is ill-adapted to the Indian growing conditions, it was found to contribute favorable alleles in Stable QTL for five architectural traits-press force, Dia, Plht, Bih, and Pbr; these QTL could be used to breed superior ideotypes in the oleiferous mustard lines. A QTL cluster on LG A10 contained Stable QTL for seven architectural traits that included major QTL (phenotypic variance ≥ 10%) for Df and Pbr, with Tumida contributing the trait-enhancing alleles for both. Since early flowering is critical for the cultivation of mustard in the Indian subcontinent, this QTL cannot be used for the improvement of Pbr in the Indian gene pool lines. Conditional QTL analysis for Pbr, however, identified other QTL which could be used for the improvement of Pbr without affecting Df. The Stable QTL intervals were mapped on the genome assemblies of Tumida and Varuna for the identification of candidate genes.

Cruz PC, Borromeo TH, Lalusin AG, Mauleon R, McNally KL, Swamy BPM. Molecular dissection of connected rice populations revealed important genomic regions for agronomic and biofortification traits

Breeding staple crops with increased micronutrient concentration is a sustainable approach to address micronutrient malnutrition. We carried out Multi-Cross QTL analysis and Inclusive Composite Interval Mapping for 11 agronomic, yield and biofortification traits using four connected RILs populations of rice. Overall, MC-156 QTLs were detected for agronomic (115) and biofortification (41) traits, which were higher in number but smaller in effects compared to single population analysis. The MC-QTL analysis was able to detect important QTLs viz: qZn_5.2, qFe_7.1, qGY_10.1, qDF_7.1, qPH_1.1, qNT_4.1, qPT_4.1, qPL_1.2, qTGW_5.1, qGL_3.1 , and qGW_6.1 , which can be used in rice genomics assisted breeding. A major QTL (qZn_5.2 ) for grain Zn concentration has been detected on chromosome 5 that accounted for 13% of R². In all, 26 QTL clusters were identified on different chromosomes. qPH_6.1 epistatically interacted with qZn_5.1 and qGY_6.2 . Most of QTLs were co-located with functionally related candidate genes indicating the accuracy of QTL mapping. The genomic region of qZn_5.2 was co-located with putative genes such as OsZIP5, OsZIP9, and LOC_OS05G40490 that are involved in Zn uptake. These genes included polymorphic functional SNPs, and their promoter regions were enriched with cis-regulatory elements involved in plant growth and development, and biotic and abiotic stress tolerance. Major effect QTL identified for biofortification and agronomic traits can be utilized in breeding for Zn biofortified rice varieties.

Multi-omics analysis identifies drivers of protein phosphorylation

BACKGROUND

Phosphorylation of proteins is a key step in the regulation of many cellular processes including activation of enzymes and signaling cascades. The abundance of a phosphorylated peptide (phosphopeptide) is determined by the abundance of its parent protein and the proportion of target sites that are phosphorylated.

RESULTS

We quantified phosphopeptides, proteins, and transcripts in heart, liver, and kidney tissue samples of mice from 58 strains of the Collaborative Cross strain panel. We mapped ~700 phosphorylation quantitative trait loci (phQTL) across the three tissues and applied genetic mediation analysis to identify causal drivers of phosphorylation. We identified kinases, phosphatases, cytokines, and other factors, including both known and potentially novel interactions between target proteins and genes that regulate site-specific phosphorylation. Our analysis highlights multiple targets of pyruvate dehydrogenase kinase 1 (PDK1), a regulator of mitochondrial function that shows reduced activity in the NZO/HILtJ mouse, a polygenic model of obesity and type 2 diabetes.

CONCLUSIONS

Together, this integrative multi-omics analysis in genetically diverse CC strains provides a powerful tool to identify regulators of protein phosphorylation. The data generated in this study provides a resource for further exploration.

High-Density Linkage Maps from Japanese Rice japonica Recombinant Inbred Lines Using Genotyping by Random Amplicon Sequencing-Direct (GRAS-Di)

The genetic dissection of agronomically important traits in closely related Japanese rice cultivars is still in its infancy mainly because of the narrow genetic diversity within japonica rice cultivars. In an attempt to unveil potential polymorphism between closely related Japanese rice cultivars, we used a next-generation-sequencing-based genotyping method: genotyping by random amplicon sequencing-direct (GRAS-Di) to develop genetic linkage maps. In this study, four recombinant inbred line (RIL) populations and their parents were used. A final RIL number of 190 for RIL71, 96 for RIL98, 95 for RIL16, and 94 for RIL91 derived from crosses between a common leading Japanese rice cultivar Koshihikari and Yamadanishiki, Taichung 65, Fujisaka 5, and Futaba, respectively, and the parent plants were subjected to GRAS-Di library construction and sequencing. Approximately 438.7 Mbp, 440 Mbp, 403.1 Mbp, and 392 Mbp called bases covering 97.5%, 97.3%, 98.3%, and 96.1%, respectively, of the estimated rice genome sequence at average depth of 1× were generated. Analysis of genotypic data identified 1050, 1285, 1708, and 1704 markers for each of the above RIL populations, respectively. Markers generated by GRAS-Di were organized into linkage maps and compared with those generated by GoldenGate SNP assay of the same RIL populations; the average genetic distance between markers showed a clear decrease in the four RIL populations when we integrated markers of both linkage maps. Genetic studies using these markers successfully localized five QTLs associated with heading date on chromosomes 3, 6, and 7 and which previously were identified as Hd1, Hd2, Hd6, Hd16, and Hd17. Therefore, GRAS-Di technology provided a low cost and efficient genotyping to overcome the narrow genetic diversity in closely related Japanese rice cultivars and enabled us to generate a high density linkage map in this germplasm.

Linkage and association mapping of ovule number per ovary (ON) in oilseed rape (Brassica napus L.)

Ovule number (ON) produced during flower development determines the maximum number of seeds per silique and thereby affects crop productivity; however, the genetic basis of ON remains poorly understood in oilseed rape (Brassica napus). In this study, we genetically dissected the ON variations in a double haploid (DH) population and in natural population (NP) by linkage mapping and genome-wide association analysis. Phenotypic analysis showed that ON displayed normal distribution in both populations with the broad-sense heritability of 0.861 (DH population) and 0.930 (natural population). Linkage mapping identified 5 QTLs related to ON, including qON-A03, qON-A07, qON-A07-2, qON-A10, and qON-C06. Genome-wide association studies (GWAS) revealed 214, 48, and 40 significant single-nucleotide polymorphisms (SNPs) by individually using the single-locus model GLM and the multiple-locus model MrMLM and FASTMrMLM. The phenotypic variation explained (PVE) by these QTLs and SNPs ranged from 2.00-17.40% to 5.03-7.33%, respectively. Integration of the results from both strategies identified four consensus genomic regions associated with ON from the chromosomes A03, A07, and A10. Our results preliminarily resolved the genetic basis of ON and provides useful molecular markers for plant yield improvement in B. napus.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01355-7.

Identification of the egusi seed trait locus (eg) and its suppressor gene associated with the thin seed coat trait in watermelon

Egusi watermelon has a unique egusi seed type, which could be useful for breeding both edible seeds and edible flesh in watermelon. However, the genetic basis of the unique egusi seed type is not clear. In the present study, we first reported that at least two genes with inhibitory epistasis were responsible for the thin seed coat (unique egusi seed type) in watermelon. Inheritance analysis of five populations, including F₂, BC, and BCF₂, suggested that the thin seed coat trait was controlled by a suppressor gene together with the egusi seed locus (eg) in egusi watermelon. Based on high-throughput sequencing technology, two quantitative trait loci located on chromosome 1 and chromosome 6 were identified for the thin seed coat trait in watermelon. One of the loci, the eg locus on chromosome 6, was finely mapped to a genomic region of 15.7 kb, which contained only one candidate gene. Comparative transcriptome analysis highlighted differentially expressed genes involved in cellulose and lignin synthesis between watermelon genotypes varying in the thickness of the seed coat and provided several potential candidate genes for the thin seed coat trait. Taken together, our data suggest that at least two genes are complementarily involved in the thin seed coat trait and will be useful for cloning novel genes. The results presented here provide a new reference for uncovering egusi seed genetic mechanisms and valuable information for marker-assisted selection in seed coat breeding.

Identification of major quantitative trait loci and candidate genes for seed weight in soybean

Four major quantitative trait loci for 100-seed weight were identified in a soybean RIL population under five environments, and the most likely candidate genes underlying these loci were identified. Seed weight is an important target of soybean breeding. However, the genes underlying the major quantitative trait loci (QTL) controlling seed weight remain largely unknown. In this study, a soybean population of 300 recombinant inbred lines (RILs) derived from a cross between PI595843 (PI) and WH was used to map the QTL and identify candidate genes for seed weight. The RIL population was genotyped through whole genome resequencing, and phenotyped for 100-seed weight under five environments. A total of 38 QTL were detected, and four major QTL, each explained at least 10% of the variation in 100-seed weight, were identified. Six candidate genes within these four major QTL regions were identified by analyses of their tissue expression patterns, gene annotations, and differential gene expression levels in soybean seeds during four developmental stages between two parental lines. Further sequence variation analyses revealed a C to T substitution in the first exon of the Glyma.19G143300, resulting in an amino acid change between PI and WH, and thus leading to a different predicted kinase domain, which might affect its protein function. Glyma.19G143300 is highly expressed in soybean seeds and encodes a leucine-rich repeat receptor-like protein kinase (LRR-RLK). Its predicted protein has typical domains of LRR-RLK family, and phylogenetic analyses reveled its similarity with the known LRR-RLK protein XIAO (LOC_Os04g48760), which is involved in controlling seed size. The major QTL and candidate genes identified in this study provide useful information for molecular breeding of new soybean cultivars with desirable seed weight.

Candidate loci for breeding compact plant-type soybean varieties

Plant height and node number are important agronomic traits that influence yield in soybean (Glycine max L.). Here, to better understand the genetic basis of the traits, we used two recombinant inbred line (RIL) populations to detect quantitative trait loci (QTLs) associated with plant height and node number in different environments. This analysis detected 9 and 21 QTLs that control plant height and node number, respectively. Among them, we identified two genomic regions that overlap with Determinate stem 1 (Dt1) and Dt2, which are known to influence both plant height and node number. Furthermore, different combinations of Dt1 and Dt2 alleles were enriched in distinct latitudes. In addition, we determined that the QTLs qPH-13-SE and qPH-13-DW in the two RIL populations overlap with genomic intervals associated with plant height and the QTL qNN-04-DW overlaps with an interval associated with node number. Combining the dwarf allele of qPH-13-SE/qPH-13-DW and the multiple-node allele of qNN-04-DW produced plants with ideal plant architecture, i.e., shorter main stems with more nodes. This plant type may help increase yield at high planting density. This study thus provides candidate loci for breeding elite soybean cultivars for plant height and node number.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01352-2.

Genetic dissection of QTLs for oil content in four maize DH populations

Oil is one of the main components in maize kernels. Increasing the total oil content (TOC) is favorable to optimize feeding requirement by improving maize quality. To better understand the genetic basis of TOC, quantitative trait loci (QTL) in four double haploid (DH) populations were explored. TOC exhibited continuously and approximately normal distribution in the four populations. The moderate to high broad-sense heritability (67.00-86.60%) indicated that the majority of TOC variations are controlled by genetic factors. A total of 16 QTLs were identified across all chromosomes in a range of 3.49-30.84% in term of phenotypic variation explained. Among them, six QTLs were identified as the major QTLs that explained phenotypic variation larger than 10%. Especially, qOC-1-3 and qOC-2-3 on chromosome 9 were recognized as the largest effect QTLs with 30.84% and 21.74% of phenotypic variance, respectively. Seventeen well-known genes involved in fatty acid metabolic pathway located within QTL intervals. These QTLs will enhance our understanding of the genetic basis of TOC in maize and offer prospective routes to clone candidate genes regulating TOC for breeding program to cultivate maize varieties with the better grain quality.

Large-scale genomic and transcriptomic profiles of rice hybrids reveal a core mechanism underlying heterosis

BACKGROUND

Heterosis is widely used in agriculture. However, its molecular mechanisms are still unclear in plants. Here, we develop, sequence, and record the phenotypes of 418 hybrids from crosses between two testers and 265 rice varieties from a mini-core collection.

RESULTS

Phenotypic analysis shows that heterosis is dependent on genetic backgrounds and environments. By genome-wide association study of 418 hybrids and their parents, we find that nonadditive QTLs are the main genetic contributors to heterosis. We show that nonadditive QTLs are more sensitive to the genetic background and environment than additive ones. Further simulations and experimental analysis support a novel mechanism, homo-insufficiency under insufficient background (HoIIB), underlying heterosis. We propose heterosis in most cases is not due to heterozygote advantage but homozygote disadvantage under the insufficient genetic background.

CONCLUSION

The HoIIB model elucidates that genetic background insufficiency is the intrinsic mechanism of background dependence, and also the core mechanism of nonadditive effects and heterosis. This model can explain most known hypotheses and phenomena about heterosis, and thus provides a novel theory for hybrid rice breeding in future.

QTLs Related to Rice Callus Regeneration Ability: Localization and Effect Verification of qPRR3

Mature and efficient tissue culture systems are already available for most japonica rice varieties (Oryza sativa ssp. geng). However, it remains challenging to regenerate the majority of indica rice varieties (Oryza sativa ssp. xian). In this study, quantitative trait loci (QTLs) associated with rice callus regeneration ability were identified based on the plant regeneration rate (PRR) and total green plant rate (TGPR) of the 93-11 × Nip recombinant inbred line population. Significant positive correlations were found between PRR and TGPR. A total of three QTLs (one for PRR and two for TGPR) were identified. qPRR3 (located on chromosome 3) was detected for both traits, which could explain 13.40% and 17.07% of the phenotypic variations of PRR and TGPR, respectively. Subsequently, the effect of qPRR3 on callus regeneration ability was validated by cryptographically tagged near-isogenic lines (NILs), and the QTL was narrowed to an interval of approximately 160 kb. The anatomical structure observation of the regenerated callus of the NILs revealed that qPRR3 can improve the callus regeneration ability by promoting the regeneration of shoots.

Identification and validation of a key genomic region on chromosome 6 for resistance to Fusarium stalk rot in tropical maize

A key genomic region was identified for resistance to FSR at 168 Mb on chromosome 6 in GWAS and haplotype regression analysis, which was validated by QTL mapping in two populations. Fusarium stalk rot (FSR) of maize is an economically important post-flowering stalk rot (PFSR) disease caused by Fusarium verticillioides. The pathogen invades the plant individually, or in combination with other stalk rot pathogens or secondary colonizers, thereby making it difficult to make accurate selection for resistance. For identification and validation of genomic regions associated with FSR resistance, a genome-wide association study (GWAS) was conducted with 342 maize lines. The panel was screened for FSR in three environments using standard artificial inoculation methodology. GWAS using the mixed linear model corrected for population structure and kinship was done, in which 290,626 SNPs from genotyping-by-sequencing were used. A total of 7 SNPs, five on chromosome 6 showing strong LD at 168 Mb, were identified to be associated with FSR. Haplotype regression analysis identified 32 haplotypes with a significant effect on the trait. In a QTL mapping experiment in two populations for validating the identified variants, QTLs were identified with confidence intervals having overlapped physical coordinates in both the populations on chromosome 6, which was closely located to the GWAS-identified variants on chromosome 6. It makes this genomic region a crucial one to further investigate the possibility of developing trait markers for deployment in breeding pipelines. It was noted that previously reported QTLs for other stalk rots in maize mapped within the same physical intervals of several haplotypes identified for FSR resistance in this study. The possibility of QTLs controlling broad-spectrum resistance for PFSR in general requires further investigation.

An Integration of Linkage Mapping and GWAS Reveals the Key Genes for Ear Shank Length in Maize

Ear shank length (ESL) has significant effects on grain yield and kernel dehydration rate in maize. Herein, linkage mapping and genome-wide association study were combined to reveal the genetic architecture of maize ESL. Sixteen quantitative trait loci (QTL) were identified in the segregation population, among which five were repeatedly detected across multiple environments. Meanwhile, 23 single nucleotide polymorphisms were associated with the ESL in the association panel, of which four were located in the QTL identified by linkage mapping and were designated as the population-common loci. A total of 42 genes residing in the linkage disequilibrium regions of these common variants and 12 of them were responsive to ear shank elongation. Of the 12 genes, five encode leucine-rich repeat receptor-like protein kinases, proline-rich proteins, and cyclin11, respectively, which were previously shown to regulate cell division, expansion, and elongation. Gene-based association analyses revealed that the variant located in Cyclin11 promoter affected the ESL among different lines. Cyclin11 showed the highest expression in the ear shank 15 days after silking among diverse tissues of maize, suggesting its role in modulating ESL. Our study contributes to the understanding of the genetic mechanism underlying maize ESL and genetic modification of maize dehydration rate and kernel yield.

A large-effect QTL introgressed from ricebean imparts resistance to Mungbean yellow mosaic India virus in blackgram (Vigna mungo (L.) Hepper)

Here, we report identification of a large effect QTL conferring Mungbean yellow mosaic India virus resistance introgressed from ricebean in blackgram variety Mash114. The tightly linked KASP markers would assist in marker-assisted-transfer of this region into Vigna species infected by MYMIV. Until recently, precise location of genes and marker-assisted selection was long thought in legumes such as blackgram due to lack of dense molecular maps. However, advances in next-generation sequencing based on high-throughput genotyping technologies such as QTL-seq have revolutionized trait mapping in marker-orphan crops. Using QTL-seq approach, we have identified a large-effect QTL for resistance to Mungbean yellow mosaic India virus (MYMIV) in blackgram variety Mash114. MYMIV is devastating disease responsible for huge yield losses in blackgram, greengram and other legumes. Mash114 showed consistent and high level of resistance to MYMIV since last nine years. Whole genome re-sequencing of MYMIV-resistant and susceptible bulks derived from RILs of cross KUG253 X Mash114 identified a large-effect QTL (qMYMIV6.1.1) spanning 3.4 Mb on chromosome 6 explaining 70% of total phenotypic variation. This region was further identified as an inter-specific introgression from ricebean. Linkage mapping using KASP markers developed from potent candidate genes involved in virus resistance identified the 500 kb genomic region equaling 1.9 cM on genetic map linked with MYMIV. The three KASP markers closely associated with MYMIV originated from serine threonine kinase, UBE2D2 and BAK1/BRI1-ASSOCIATED RECEPTOR KINASE genes. These KASPs can be used for marker-assisted transfer of introgressed segment into suitable backgrounds of Vigna species.

Systematic trait dissection in oilseed rape provides a comprehensive view, further insight, and exact roadmap for yield determination

Background

Yield is the most important and complex trait that is influenced by numerous relevant traits with very complicated interrelations. While there are a large number of studies on the phenotypic relationship and genetic basis of yield traits, systematic studies with further dissection focusing on yield are limited. Therefore, there is still lack of a comprehensive and in-depth understanding of the determination of yield.

Results

In this study, yield was systematically dissected at the phenotypic, genetic to molecular levels in oilseed rape (Brassica napus L.). The analysis of correlation, network, and principal component for 21 traits in BnaZN-RIL population showed that yield was determined by a complex trait network with key contributors. The analysis of the constructed high-density single nucleotide polymorphism (SNP) linkage map revealed the concentrated distribution of distorted and heterozygous markers, likely due to selection on genes controlling the growth period and yield heterosis. A total of 134 consensus quantitative trait loci (QTL) were identified for 21 traits, of which all were incorporated into an interconnecting QTL network with dozens of hub-QTL. Four representative hub-QTL were further dissected to the target or candidate genes that governed the causal relationships between the relevant traits.

Conclusions

The highly consistent results at the phenotypic, genetic, and molecular dissecting demonstrated that yield was determined by a multilayer composite network that involved numerous traits and genes showing complex up/down-stream and positive/negative regulation. This provides a systematic view, further insight, and exact roadmap for yield determination, which represents a significant advance toward the understanding and dissection of complex traits.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13068-022-02134-w.

Identification of two quantitative genes controlling soybean flowering using bulked-segregant analysis and genetic mapping

Photoperiod responsiveness is important to soybean production potential and adaptation to local environments. Varieties from temperate regions generally mature early and exhibit extremely low yield when grown under inductive short-day (SD) conditions. The long-juvenile (LJ) trait is essentially a reduction and has been introduced into soybean cultivars to improve yield in tropical environments. In this study, we used next-generation sequencing (NGS)-based bulked segregant analysis (BSA) to simultaneously map qualitative genes controlling the LJ trait in soybean. We identified two genomic regions on scaffold_32 and chromosome 18 harboring loci LJ32 and LJ18, respectively. Further, we identified LJ32 on the 228.7-kb scaffold_32 as the soybean pseudo-response-regulator gene Tof11 and LJ18 on a 301-kb region of chromosome 18 as a novel PROTEIN FLOWERING LOCUS T-RELATED gene, Glyma.18G298800. Natural variants of both genes contribute to LJ trait regulation in tropical regions. The molecular identification and functional characterization of Tof11 and LJ18 will enhance understanding of the molecular mechanisms underlying the LJ trait and provide useful genetic resources for soybean molecular breeding in tropical regions.

Architecture of the outbred brown fat proteome defines regulators of metabolic physiology

Brown adipose tissue (BAT) regulates metabolic physiology. However, nearly all mechanistic studies of BAT protein function occur in a single inbred mouse strain, which has limited the understanding of generalizable mechanisms of BAT regulation over physiology. Here, we perform deep quantitative proteomics of BAT across a cohort of 163 genetically defined diversity outbred mice, a model that parallels the genetic and phenotypic variation found in humans. We leverage this diversity to define the functional architecture of the outbred BAT proteome, comprising 10,479 proteins. We assign co-operative functions to 2,578 proteins, enabling systematic discovery of regulators of BAT. We also identify 638 proteins that correlate with protection from, or sensitivity to, at least one parameter of metabolic disease. We use these findings to uncover SFXN5, LETMD1, and ATP1A2 as modulators of BAT thermogenesis or adiposity, and provide OPABAT as a resource for understanding the conserved mechanisms of BAT regulation over metabolic physiology.

Genetic mapping of the Andean anthracnose resistance gene present in the common bean cultivar BRSMG Realce

The rajado seeded Andean bean (Phaseolus vulgaris L.) cultivar BRSMG Realce (striped seed coat) developed by Embrapa expressed a high level of anthracnose resistance, caused by Colletotrichum lindemuthianum, in field and greenhouse screenings. The main goal of this study was to evaluate the inheritance of anthracnose resistance in BRSMG Realce, map the resistance locus or major gene cluster previously named as Co-Realce, identify resistance-related positional genes, and analyze potential markers linked to the resistance allele. F₂ plants derived from the cross BRSMG Realce × BRS FC104 (Mesoamerican) and from the cross BRSMG Realce × BRS Notável (Mesoamerican) were inoculated with the C. lindemuthianum races 475 and 81, respectively. The BRSMG Realce × BRS FC104 F₂ population was also genotyped using the DArTseq technology. Crosses between BRSMG Realce and BAT 93 (Mesoamerican) were also conducted and resulting F₂ plants were inoculated with the C. lindemuthianum races 65 and 1609, individually. The results shown that anthracnose resistance in BRSMG Realce is controlled by a single locus with complete dominance. A genetic map including 1,118 SNP markers was built and shown 78% of the markers mapped at a distances less than 5.0 cM, with a total genetic length of 4,473.4 cM. A major locus (Co-Realce) explaining 54.6% of the phenotypic variation of symptoms caused by the race 475 was identified in Pv04, flanked by the markers snp1327 and snp12782 and 4.48 cM apart each other. These SNPs are useful for marker-assisted selection, due to an estimated selection efficiency of 99.2%. The identified resistance allele segregates independently of the resistance allele Co-3³ (Pv04) present in BAT 93. The mapped genomic region with 704,867 bp comprising 63 putative genes, 44 of which were related to the pathogen-host interaction. Based on all these results and evidence, anthracnose resistance in BRSMG Realce should be considered as monogenic, useful for breeding purpose. It is proposed that locus Co-Realce is unique and be provisionally designated as CoPv04^R until be officially nominated in accordance with the rules established by the Bean Improvement Cooperative Genetics Committee.

AirMeasurer: open-source software to quantify static and dynamic traits derived from multiseason aerial phenotyping to empower genetic mapping studies in rice

Low-altitude aerial imaging, an approach that can collect large-scale plant imagery, has grown in popularity recently. Amongst many phenotyping approaches, unmanned aerial vehicles (UAVs) possess unique advantages as a consequence of their mobility, flexibility and affordability. Nevertheless, how to extract biologically relevant information effectively has remained challenging. Here, we present AirMeasurer, an open-source and expandable platform that combines automated image analysis, machine learning and original algorithms to perform trait analysis using 2D/3D aerial imagery acquired by low-cost UAVs in rice (Oryza sativa) trials. We applied the platform to study hundreds of rice landraces and recombinant inbred lines at two sites, from 2019 to 2021. A range of static and dynamic traits were quantified, including crop height, canopy coverage, vegetative indices and their growth rates. After verifying the reliability of AirMeasurer-derived traits, we identified genetic variants associated with selected growth-related traits using genome-wide association study and quantitative trait loci mapping. We found that the AirMeasurer-derived traits had led to reliable loci, some matched with published work, and others helped us to explore new candidate genes. Hence, we believe that our work demonstrates valuable advances in aerial phenotyping and automated 2D/3D trait analysis, providing high-quality phenotypic information to empower genetic mapping for crop improvement.

Natural variation of Dt2 determines branching in soybean

Shoot branching is fundamentally important in determining soybean yield. Here, through genome-wide association study, we identify one predominant association locus on chromosome 18 that confers soybean branch number in the natural population. Further analyses determine that Dt2 is the corresponding gene and the natural variations in Dt2 result in significant differential transcriptional levels between the two major haplotypes. Functional characterization reveals that Dt2 interacts with GmAgl22 and GmSoc1a to physically bind to the promoters of GmAp1a and GmAp1d and to activate their transcription. Population genetic investigation show that the genetic differentiation of Dt2 display significant geographic structure. Our study provides a predominant gene for soybean branch number and may facilitate the breeding of high-yield soybean varieties.

High density linkage maps, genetic architecture, and genomic prediction of growth and wood properties in Pinus radiata

Background

The growing availability of genomic resources in radiata pine paves the way for significant advances in fundamental and applied genomic research. We constructed robust high-density linkage maps based on exome-capture genotyping in two F₁ populations, and used these populations to perform quantitative trait locus (QTL) scans, genomic prediction and quantitative analyses of genetic architecture for key traits targeted by tree improvement programmes.

Results

Our mapping approach used probabilistic error correction of the marker data, followed by an iterative approach based on stringent parameters. This approach proved highly effective in producing high-density maps with robust marker orders and realistic map lengths (1285–4674 markers per map, with sizes ranging from c. 1643–2292 cM, and mean marker intervals of 0.7–2.1 cM). Colinearity was high between parental linkage maps, although there was evidence for a large chromosomal rearrangement (affecting ~ 90 cM) in one of the parental maps. In total, 28 QTL were detected for growth (stem diameter) and wood properties (wood density and fibre properties measured by Silviscan) in the QTL discovery population, with 1–3 QTL of small to moderate effect size detected per trait in each parental map. Four of these QTL were validated in a second, unrelated F₁ population. Results from genomic prediction and analyses of genetic architecture were consistent with those from QTL scans, with wood properties generally having moderate to high genomic heritabilities and predictive abilities, as well as somewhat less complex genetic architectures, compared to growth traits.

Conclusions

Despite the economic importance of radiata pine as a plantation forest tree, robust high-density linkage maps constructed from reproducible, sequence-anchored markers have not been published to date. The maps produced in this study will be a valuable resource for several applications, including the selection of marker panels for genomic prediction and anchoring a recently completed de novo whole genome assembly. We also provide the first map-based evidence for a large genomic rearrangement in radiata pine. Finally, results from our QTL scans, genomic prediction, and genetic architecture analyses are informative about the genomic basis of variation in important phenotypic traits.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08950-6.

QTL analysis for sodium removal ability in rice leaf sheaths under salinity using an IR-44595/318 F2 population

Over-accumulation of salt in rice plants is an effect of salt stress which decreases growth and grain yield. Salt removal ability in leaf sheaths is a tolerance mechanism to decrease salt entry and accumulation in leaf blades and maintain photosynthesis under salinity. In this study, a QTL analysis of removal ability of sodium ions (Na⁺) in leaf sheaths and Na⁺ accumulation-related traits, was conducted using F₂ population between two rice varieties, IR-44595 with superior Na⁺ removal ability, and 318 with contrasting Na⁺ removal ability in leaf sheaths under salinity. Suggestive QTLs for Na⁺ removal ability in leaf sheaths were found on chromosomes 4 and 11. The suggestive QTL on chromosome 11 overlapped with other significant QTLs for Na⁺ concentration in shoots, leaf blades and leaf sheaths, and Na⁺/K⁺ ratio in leaf blades. Correlation analysis indicated that Na⁺ removal ability in leaf sheaths is important in reducing Na⁺ accumulation in leaf blades. The varietal difference of Na⁺ removal ability in leaf sheaths at the whole plant level was greater at lower NaCl concentrations and became smaller as the treatment NaCl concentration increased. Although the Na⁺ removal ability in leaf sheath was comparable between IR-44595 and 318 under high salinity at the whole plant level, the younger leaves of IR-44595 still showed a higher Na⁺ sheath-blade ratio than 318, which implied the Na⁺ removal ability functions in the younger leaves in IR-44595 to reduce Na⁺ entry in young leaf blades even under high salinity.

Intergenomic and epistatic interactions control free radical mediated pancreatic β-cell damage

Alloxan (AL)-generated Reactive Oxygen Species (ROS) selectively destroy insulin-producing pancreatic β-cells. A previous genome-wide scan (GWS) using a cohort of 296 F2 hybrids between NOD (AL-sensitive) and ALR (AL-resistant) mice identified linkages contributing to β-cell susceptibility or resistance to AL-induced diabetes on Chromosomes (Chr) 2, 3, 8, and a single nucleotide polymorphism in mt-Nd2 of the mitochondrial genome (mtDNA). AL treatment of congenic and consomic NOD mouse stocks confirmed resistance linked to both the mtDNA and the Chr 8 locus from ALR [NOD.mt^ALR.ALR-(D8Mit293-D8Mit137)]. To identify possible epistatic interactions, the GWS analysis was expanded to 678 F2 mice. ALR-derived diabetes-resistance linkages on Chr 8 as well as the mt-Nd2^a allele were confirmed and novel additional linkages on Chr 4, 5, 6, 7, and 13 were identified. Epistasis was observed between the linkages on Chr 8 and 2 and Chr 8 and 6. Furthermore, the mt-Nd2 genotype affected the epistatic interactions between Chr 8 and 2. These results demonstrate that a combination of nuclear-cytoplasmic genome interactions regulates β-cell sensitivity to ROS-mediated ALD.

Identification and characterization of pleiotropic and epistatic QDRL conferring partial resistance to Pythium irregulare and P. sylvaticum in soybean

Pleiotropic and epistatic quantitative disease resistance loci (QDRL) were identified for soybean partial resistance to different isolates of Pythium irregulare and Pythium sylvaticum. Pythium root rot is an important seedling disease of soybean [Glycine max (L.) Merr.], a crop grown worldwide for protein and oil content. Pythium irregulare and P. sylvaticum are two of the most prevalent and aggressive Pythium species in soybean producing regions in the North Central U.S. Few studies have been conducted to identify soybean resistance for management against these two pathogens. In this study, a mapping population (derived from E13390 x E13901) with 228 F4:5 recombinant inbred lines were screened against P. irregulare isolate MISO 11-6 and P. sylvaticum isolate C-MISO2-2-30 for QDRL mapping. Correlation analysis indicated significant positive correlations between soybean responses to the two pathogens, and a pleiotropic QDRL (qPirr16.1) was identified. Further investigation found that the qPirr16.1 imparts dominant resistance against P. irregulare, but recessive resistance against P. sylvaticum. In addition, two QDRL, qPsyl15.1, and qPsyl18.1 were identified for partial resistance to P. sylvaticum. Further analysis revealed epistatic interactions between qPirr16.1 and qPsyl15.1 for RRW and DRX, whereas qPsyl18.1 contributed resistance to RSE. Marker-assisted resistance spectrum analysis using F6:7 progeny lines verified the resistance of qPirr16.1 against four additional P. irregulare isolates. Intriguingly, although the epistatic interaction of qPirr16.1 and qPsyl15.1 can be confirmed using two additional isolates of P. sylvaticum, the interaction appears to be suppressed for the other two P. sylvaticum isolates. An 'epistatic gene-for-gene' model was proposed to explain the isolate-specific epistatic interactions. The integration of the QDRL into elite soybean lines containing all the desirable alleles has been initiated.

Improved pea reference genome and pan-genome highlight genomic features and evolutionary characteristics

Complete and accurate reference genomes and annotations provide fundamental resources for functional genomics and crop breeding. Here we report a de novo assembly and annotation of a pea cultivar ZW6 with contig N50 of 8.98 Mb, which features a 243-fold increase in contig length and evident improvements in the continuity and quality of sequence in complex repeat regions compared with the existing one. Genome diversity of 118 cultivated and wild pea demonstrated that Pisum abyssinicum is a separate species different from P. fulvum and P. sativum within Pisum. Quantitative trait locus analyses uncovered two known Mendel's genes related to stem length (Le/le) and seed shape (R/r) as well as some candidate genes for pod form studied by Mendel. A pan-genome of 116 pea accessions was constructed, and pan-genes preferred in P. abyssinicum and P. fulvum showed distinct functional enrichment, indicating the potential value of them as pea breeding resources in the future.

Mapping QTLs for Reproductive Stage Salinity Tolerance in Rice Using a Cross between Hasawi and BRRI dhan28

Salinity stress is a major constraint to rice production in many coastal regions due to saline groundwater and river sources, especially during the dry season in coastal areas when seawater intrudes further inland due to reduced river flows. Since salinity tolerance is a complex trait, breeding efforts can be assisted by mapping quantitative trait loci (QTLs) for complementary salt tolerance mechanisms, which can then be combined to provide higher levels of tolerance. While an abundance of seedling stage salinity tolerance QTLs have been mapped, few studies have investigated reproductive stage tolerance in rice due to the difficulty of achieving reliable stage-specific phenotyping techniques. In the current study, a BC₁F₂ mapping population consisting of 435 individuals derived from a cross between a salt-tolerant Saudi Arabian variety, Hasawi, and a salt-sensitive Bangladeshi variety, BRRI dhan28, was evaluated for yield components after exposure to EC 10 dS/m salinity stress during the reproductive stage. After selecting tolerant and sensitive progeny, 190 individuals were genotyped by skim sequencing, resulting in 6209 high quality single nucleotide polymorphic (SNP) markers. Subsequently, a total of 40 QTLs were identified, of which 24 were for key traits, including productive tillers, number and percent filled spikelets, and grain yield under stress. Importantly, three yield-related QTLs, one each for productive tillers (qPT3.1), number of filled spikelets (qNFS3.1) and grain yield (qGY3.1) under salinity stress, were mapped at the same position (6.7 Mb or 26.1 cM) on chromosome 3, which had not previously been associated with grain yield under salinity stress. These QTLs can be investigated further to dissect the molecular mechanisms underlying reproductive stage salinity tolerance in rice.

Complex genetics cause and constrain fungal persistence in different parts of the mammalian body

Determining how genetic polymorphisms enable certain fungi to persist in mammalian hosts can improve understanding of opportunistic fungal pathogenesis, a source of substantial human morbidity and mortality. We examined the genetic basis of fungal persistence in mice using a cross between a clinical isolate and the lab reference strain of the budding yeast Saccharomyces cerevisiae. Employing chromosomally encoded DNA barcodes, we tracked the relative abundances of 822 genotyped, haploid segregants in multiple organs over time and performed linkage mapping of their persistence in hosts. Detected loci showed a mix of general and antagonistically pleiotropic effects across organs. General loci showed similar effects across all organs, while antagonistically pleiotropic loci showed contrasting effects in the brain vs the kidneys, liver, and spleen. Persistence in an organ required both generally beneficial alleles and organ-appropriate pleiotropic alleles. This genetic architecture resulted in many segregants persisting in the brain or in nonbrain organs, but few segregants persisting in all organs. These results show complex combinations of genetic polymorphisms collectively cause and constrain fungal persistence in different parts of the mammalian body.

Quantitative trait loci mapping reveals important genomic regions controlling root architecture and shoot biomass under nitrogen, phosphorus, and potassium stress in rapeseed (Brassica napus L.)

Plants rely on root systems for nutrient uptake from soils. Marker-assisted selection helps breeders to select desirable root traits for effective nutrient uptake. Here, 12 root and biomass traits were investigated at the seedling stage under low nitrogen (LN), low phosphorus (LP), and low potassium (LK) conditions, respectively, in a recombinant inbred line (RIL) population, which was generated from Brassica napus L. Zhongshuang11 and 4D122 with significant differences in root traits and nutrient efficiency. Significant differences for all the investigated traits were observed among RILs, with high heritabilities (0.43-0.74) and high correlations between the different treatments. Quantitative trait loci (QTL) mapping identified 57, 27, and 36 loci, explaining 4.1-10.9, 4.6-10.8, and 4.9-17.4% phenotypic variances under LN, LP, and LK, respectively. Through QTL-meta analysis, these loci were integrated into 18 significant QTL clusters. Four major QTL clusters involved 25 QTLs that could be repeatedly detected and explained more than 10% phenotypic variances, including two NPK-common and two specific QTL clusters (K and NK-specific), indicating their critical role in cooperative nutrients uptake of N, P, and K. Moreover, 264 genes within the four major QTL clusters having high expressions in roots and SNP/InDel variations between two parents were identified as potential candidate genes. Thirty-eight of them have been reported to be associated with root growth and development and/or nutrient stress tolerance. These key loci and candidate genes lay the foundation for deeper dissection of the NPK starvation response mechanisms in B. napus.

Compatibility between snails and schistosomes: insights from new genetic resources, comparative genomics, and genetic mapping

The freshwater snail Biomphalaria glabrata is an important intermediate host of the parasite Schistosoma mansoni that causes human intestinal schistosomiasis. To better understand vector snail biology and help advance innovative snail control strategies, we have developed a new snail model consisting of two homozygous B. glabrata lines (iM line and iBS90) with sharply contrasting schistosome-resistance phenotypes. We produced and compared high-quality genome sequences for iM line and iBS90 which were assembled from 255 (N50 = 22.7 Mb) and 346 (N50 = 19.4 Mb) scaffolds, respectively. Using F2 offspring bred from the two lines and the newly generated iM line genome, we constructed 18 linkage groups (representing the 18 haploid chromosomes) covering 96% of the genome and identified three new QTLs (quantitative trait loci), two involved in snail resistance/susceptibility and one relating to body pigmentation. This study provides excellent genomic resources for unveiling complex vector snail biology, reveals genomic difference between resistant and susceptible lines, and offers novel insights into genetic mechanism of the compatibility between snail and schistosome.

Linkage mapping and association analysis to identify a reliable QTL for stigma exsertion rate in rice

The commercialization of hybrid rice has greatly contributed to the increase in rice yield, with the improvement of its seed production capacity having played an important role. The stigma exsertion rate (SER) is a key factor for improving the outcrossing of the sterile line and the hybrid rice seed production. We used the Zhenshan 97B × IRAT109 recombinant inbred population comprising 163 lines and a natural population of 138 accessions to decipher the genetic foundation of SER over 2 years in three environments. Additionally, we detected eight QTLs for SER on chromosomes 1, 2, and 8 via linkage mapping. We also identified seven and 19 significant associations for SER using genome-wide association study in 2016 and 2017, respectively. Interestingly, we located two lead SNPs (sf0803343504 and sf083344610) on chromosome 8 in the qTSE8 QTL region that were significantly associated with total SER. After transcriptomic analysis, quantitative real-time PCR, and haplotype analysis, we found 13 genes within this reliable region as important candidate genes. Our study results will be beneficial to molecular marker-assisted selection of rice lines with high outcrossing rate, thereby improving the efficiency of hybrid seed production.

A Genomic BSAseq Approach for the Characterization of QTLs Underlying Resistance to Fusarium oxysporum in Eggplant

Eggplant (Solanum melongena L.), similar to many other crops, suffers from soil-borne diseases, including Fusarium oxysporum f. sp. melongenae (Fom), causing wilting and heavy yield loss. To date, the genetic factors underlying plant responses to Fom are not well known. We previously developed a Recombinant Inbred Lines (RILs) population using as a female parent the fully resistant line '305E40' and as a male parent the partially resistant line '67/3'. The fully resistant trait to Fom was introgressed from the allied species S. aethiopicum. In this work, the RIL population was assessed for the responses to Fom and by using a genomic mapping approach, two major QTLs on chromosomes CH02 and CH11 were identified, associated with the full and partial resistance trait to Fom, respectively. A targeted BSAseq procedure in which Illumina reads bulks of RILs grouped according to their resistance score was aligned to the appropriate reference genomes highlighted differentially enriched regions between resistant/susceptible progeny in the genomic regions underlying both QTLs. The characterization of such regions allowed us to identify the most reliable candidate genes for the two resistance traits. With the aim of revealing exclusive species-specific contigs and scaffolds inherited from the allied species and thus associated with the full resistance trait, a draft de-novo assembly of available Illumina sequences of the '305E40' parent was developed to better resolve the non-recombining genomic region on its CH02 carrying the introgressed Fom resistance locus from S. aethiopicum.

Fine mapping of candidate quantitative trait loci for plant and ear height in a maize nested-association mapping population

Plant height (PH) and ear height (EH) are two important traits in maize (Zea mays L.), as they are closely related to lodging resistance and planting density. Our objectives were to (1) investigate single-nucleotide polymorphisms (SNPs) that are associated with PH and EH for detecting quantitative trait loci (QTL) and new gene that determines PH and EH, (2) explore the value of the QTL in maize breeding, and (3) investigate whether the "triangle heterotic group" theory is applicable for lowering PH and EH to increase yield. Seven inbred female parents were crossed with a common founder male parent Ye 107 to create a nested association mapping (NAM) population. The analysis of phenotypic data on PH and EH revealed wide variation among the parents of the NAM population. Genome-wide association study (GWAS) and high-resolution linkage mapping were conducted using the NAM population, which generated 264,694 SNPs by genotyping-by-sequencing. A total of 105 SNPs and 22 QTL were identified by GWAS and found to be significantly associated with PH and EH. A high-confidence QTL for PH, Qtl-chr1-EP, was identified on chromosome 1 via GWAS and confirmed by linkage analysis in two recombinant inbred line (RIL) populations. Results revealed that the SNP variation in the promoter region of the candidate gene Zm00001d031938, located at Qtl-chr1-EP, which encoded UDP-N-acetylglucosamine-peptide N-acetyl-glucosaminyl-transferase, might decrease PH and EH. Furthermore, the triangle heterotic pattern adopted in maize breeding programs by our team is practicable in selecting high-yield crosses based on the low ratio of EH/PH (EP).

Identification of genetic loci conferring seed coat color based on a high-density map in soybean

Seed coat color is a typical evolutionary trait. Identification of the genetic loci that control seed coat color during the domestication of wild soybean could clarify the genetic variations between cultivated and wild soybean. We used 276 F10 recombinant inbred lines (RILs) from the cross between a cultivated soybean (JY47) and a wild soybean (ZYD00321) as the materials to identify the quantitative trait loci (QTLs) for seed coat color. We constructed a high-density genetic map using re-sequencing technology. The average distance between adjacent markers was 0.31 cM on this map, comprising 9,083 bin markers. We identified two stable QTLs (qSC08 and qSC11) for seed coat color using this map, which, respectively, explained 21.933 and 26.934% of the phenotypic variation. Two candidate genes (CHS3C and CHS4A) in qSC08 were identified according to the parental re-sequencing data and gene function annotations. Five genes (LOC100786658, LOC100801691, LOC100806824, LOC100795475, and LOC100787559) were predicted in the novel QTL qSC11, which, according to gene function annotations, might control seed coat color. This result could facilitate the identification of beneficial genes from wild soybean and provide useful information to clarify the genetic variations for seed coat color in cultivated and wild soybean.

The genetic basis of neurocranial size and shape across varied lab mouse populations

Brain and skull tissues interact through molecular signalling and mechanical forces during head development, leading to a strong correlation between the neurocranium and the external brain surface. Therefore, when brain tissue is unavailable, neurocranial endocasts are often used to approximate brain size and shape. Evolutionary changes in brain morphology may have resulted in secondary changes to neurocranial morphology, but the developmental and genetic processes underlying this relationship are not well understood. Using automated phenotyping methods, we quantified the genetic basis of endocast variation across large genetically varied populations of laboratory mice in two ways: (1) to determine the contributions of various genetic factors to neurocranial form and (2) to help clarify whether a neurocranial variation is based on genetic variation that primarily impacts bone development or on genetic variation that primarily impacts brain development, leading to secondary changes in bone morphology. Our results indicate that endocast size is highly heritable and is primarily determined by additive genetic factors. In addition, a non-additive inbreeding effect led to founder strains with lower neurocranial size, but relatively large brains compared to skull size; suggesting stronger canalization of brain size and/or a general allometric effect. Within an outbred sample of mice, we identified a locus on mouse chromosome 1 that is significantly associated with variation in several positively correlated endocast size measures. Because the protein-coding genes at this locus have been previously associated with brain development and not with bone development, we propose that genetic variation at this locus leads primarily to variation in brain volume that secondarily leads to changes in neurocranial globularity. We identify a strain-specific missense mutation within Akt3 that is a strong causal candidate for this genetic effect. Whilst it is not appropriate to generalize our hypothesis for this single locus to all other loci that also contribute to the complex trait of neurocranial skull morphology, our results further reveal the genetic basis of neurocranial variation and highlight the importance of the mechanical influence of brain growth in determining skull morphology.

Uncontrolled transposition following RNAi loss causes hypermutation and antifungal drug resistance in clinical isolates of Cryptococcus neoformans

Cryptococcus neoformans infections cause approximately 15% of AIDS-related deaths owing to a combination of limited antifungal therapies and drug resistance. A collection of clinical and environmental C. neoformans isolates were assayed for increased mutation rates via fluctuation analysis, and we identified two hypermutator C. neoformans clinical isolates with increased mutation rates when exposed to the combination of rapamycin and FK506. Sequencing of drug target genes found that Cnl1 transposon insertions conferred the majority of resistance to rapamycin and FK506 and could also independently cause resistance to 5-fluoroorotic acid and the clinically relevant antifungal 5-flucytosine. Whole-genome sequencing revealed both hypermutator genomes harbour a nonsense mutation in the RNA-interference component ZNF3 and hundreds of Cnl1 elements organized into massive subtelomeric arrays on each of the fourteen chromosomes. Quantitative trait locus mapping in 28 progeny derived from a cross between a hypermutator and wild-type identified a locus associated with hypermutation that included znf3. CRISPR editing of the znf3 nonsense mutation abolished hypermutation and restored small-interfering-RNA production. We conclude that hypermutation and drug resistance in these clinical isolates result from RNA-interference loss and accumulation of Cnl1 elements.