Megabase Length Hypermutation Accompanies Human Structural Variation at 17p11.2.

DNA rearrangements resulting in human genome structural variants (SVs) are caused by diverse mutational mechanisms. We used long- and short-read sequencing technologies to investigate end products of de novo chromosome 17p11.2 rearrangements and query the molecular mechanisms underlying both recurrent and non-recurrent events. Evidence for an increased rate of clustered single-nucleotide variant (SNV) mutation in cis with non-recurrent rearrangements was found. Indel and SNV formation are associated with both copy-number gains and losses of 17p11.2, occur up to ∼1 Mb away from the breakpoint junctions, and favor C > G transversion substitutions; results suggest that single-stranded DNA is formed during the genesis of the SV and provide compelling support for a microhomology-mediated break-induced replication (MMBIR) mechanism for SV formation. Our data show an additional mutational burden of MMBIR consisting of hypermutation confined to the locus and manifesting as SNVs and indels predominantly within genes.

Exome-wide assessment of isolated biliary atresia: A report from the National Birth Defects Prevention Study using child-parent trios and a case-control design to identify novel rare variants.

The etiology of biliary atresia (BA) is unknown, but recent studies suggest a role for rare protein-altering variants (PAVs). Exome sequencing data from the National Birth Defects Prevention Study on 54 child-parent trios, one child-mother duo, and 1513 parents of children with other birth defects were analyzed. Most (91%) cases were isolated BA. We performed (1) a trio-based analysis to identify rare de novo, homozygous, and compound heterozygous PAVs and (2) a case-control analysis using a sequence kernel-based association test to identify genes enriched with rare PAVs. While we replicated previous findings on PKD1L1, our results do not suggest that recurrent de novo PAVs play important roles in BA susceptibility. In fact, our finding in NOTCH2, a disease gene associated with Alagille syndrome, highlights the difficulty in BA diagnosis. Notably, IFRD2 has been implicated in other gastrointestinal conditions and warrants additional study. Overall, our findings strengthen the hypothesis that the etiology of BA is complex.

Long read sequencing and expression studies of AHDC1 deletions in Xia-Gibbs syndrome reveal a novel genetic regulatory mechanism.

Xia-Gibbs syndrome (XGS; MIM# 615829) is a rare mendelian disorder characterized by Development Delay (DD), intellectual disability (ID), and hypotonia. Individuals with XGS typically harbor de novo protein-truncating mutations in the AT-Hook DNA binding motif containing 1 (AHDC1) gene, although some missense mutations can also cause XGS. Large de novo heterozygous deletions that encompass the AHDC1 gene have also been ascribed as diagnostic for the disorder, without substantial evidence to support their pathogenicity. We analyzed 19 individuals with large contiguous deletions involving AHDC1, along with other genes. One individual bore the smallest known contiguous AHDC1 deletion (∼350 Kb), encompassing eight other genes within chr1p36.11 (Feline Gardner-Rasheed, IFI6, FAM76A, STX12, PPP1R8, THEMIS2, RPA2, SMPDL3B) and terminating within the first intron of AHDC1. The breakpoint junctions and phase of the deletion were identified using both short and long read sequencing (Oxford Nanopore). Quantification of RNA expression patterns in whole blood revealed that AHDC1 exhibited a mono-allelic expression pattern with no deficiency in overall AHDC1 expression levels, in contrast to the other deleted genes, which exhibited a 50% reduction in mRNA expression. These results suggest that AHDC1 expression in this individual is compensated by a novel regulatory mechanism and advances understanding of mutational and regulatory mechanisms in neurodevelopmental disorders.

Identification and Fine-Mapping of a Novel QTL, qMrdd2, That Confers Resistance to Maize Rough Dwarf Disease.

Maize rough dwarf disease (MRDD) is caused by a virus and seriously affects maize quality and yield worldwide. MRDD can be most effectively controlled with disease-resistant hybrids of corn. Here, MRDD-resistant (Qi319) and -susceptible (Ye478) parental inbred maize lines and their 314 recombinant inbred lines (RILs) that were derived from a cross between them were evaluated across three environments. A stable resistance QTL, qMrdd2, was identified and mapped using best linear unbiased prediction (BLUP) values to a 0.55-Mb region between the markers MK807 and MK811 on chromosome 2 (B73 RefGen_v3) and was found to explain 8.6 to 11.0% of the total phenotypic variance in MRDD resistance. We validated the effect of qMrdd2 using a chromosome segment substitution line (CSSL) that was derived from a cross between maize inbred Qi319 as the MRDD resistance donor and Ye478 as the recipient. Disease severity index of the CSSL haplotype II harboring qMrdd2 was significantly lower than that of the susceptible parent Ye478. Subsequently, we fine-mapped qMrdd2 to a 315-kb region flanked by the markers RD81 and RD87, thus testing recombinant-derived progeny using selfed backcrossed families. In this study, we identified a novel QTL for MRDD resistance by combining the RIL and CSSL populations, thus providing important genetic information that can be used for breeding MRDD-resistant varieties of maize.

Phenotypic and protein localization heterogeneity associated with AHDC1 pathogenic protein-truncating alleles in Xia-Gibbs syndrome.

Xia-Gibbs syndrome (XGS) is a rare Mendelian disease typically caused by de novo stop-gain or frameshift mutations in the AT-hook DNA binding motif containing 1 (AHDC1) gene. Patients usually present in early infancy with hypotonia and developmental delay and later exhibit intellectual disability (ID). The overall presentation is variable, however, and the emerging clinical picture is still evolving. A detailed phenotypic analysis of 34 XGS individuals revealed five core phenotypes (delayed motor milestones, speech delay, low muscle tone, ID, and hypotonia) in more than 80% of individuals and an additional 12 features that occurred more variably. Seizures and scoliosis were more frequently associated with truncations that arise before the midpoint of the protein although the occurrence of most features could not be predicted by the mutation position. Transient expression of wild type and different patient truncated AHDC1 protein forms in human cell lines revealed abnormal patterns of nuclear localization including a diffuse distribution of a short truncated form and nucleolar aggregation in mid-protein truncated forms. Overall, both the occurrence of variable phenotypes and the different distribution of the expressed protein reflect the heterogeneity of this syndrome.

Transmission event of SARS-CoV-2 delta variant reveals multiple vaccine breakthrough infections.

BACKGROUND: This study aims to identify the causative strain of SARS-CoV-2 in a cluster of vaccine breakthroughs. Vaccine breakthrough by a highly transmissible SARS-CoV-2 strain is a risk to global public health. METHODS: Nasopharyngeal swabs from suspected vaccine breakthrough cases were tested for SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) by qPCR (quantitative polymerase chain reaction) for Wuhan-Hu1 and alpha variant. Positive samples were then sequenced by Swift Normalase Amplicon Panels to determine the causal variant. GATK (genome analysis toolkit) variants were filtered with allele fraction ≥80 and min read depth 30x. RESULTS: Viral sequencing revealed an infection cluster of 6 vaccinated patients infected with the delta (B.1.617.2) SARS-CoV-2 variant. With no history of vaccine breakthrough, this suggests the delta variant may possess immune evasion in patients that received the Pfizer BNT162b2, Moderna mRNA-1273, and Covaxin BBV152. CONCLUSIONS: Delta variant may pose the highest risk out of any currently circulating SARS-CoV-2 variants, with previously described increased transmissibility over alpha variant and now, possible vaccine breakthrough. FUNDING: Parts of this work was supported by the National Institute of Allergy and Infectious Diseases (1U19AI144297) and Baylor College of Medicine internal funding.

Neptune: an environment for the delivery of genomic medicine.

PURPOSE: Genomic medicine holds great promise for improving health care, but integrating searchable and actionable genetic data into electronic health records (EHRs) remains a challenge. Here we describe Neptune, a system for managing the interaction between a clinical laboratory and an EHR system during the clinical reporting process. METHODS: We developed Neptune and applied it to two clinical sequencing projects that required report customization, variant reanalysis, and EHR integration. RESULTS: Neptune has been applied for the generation and delivery of over 15,000 clinical genomic reports. This work spans two clinical tests based on targeted gene panels that contain 68 and 153 genes respectively. These projects demanded customizable clinical reports that contained a variety of genetic data types including single-nucleotide variants (SNVs), copy-number variants (CNVs), pharmacogenomics, and polygenic risk scores. Two variant reanalysis activities were also supported, highlighting this important workflow. CONCLUSION: Methods are needed for delivering structured genetic data to EHRs. This need extends beyond developing data formats to providing infrastructure that manages the reporting process itself. Neptune was successfully applied on two high-throughput clinical sequencing projects to build and deliver clinical reports to EHR systems. The software is open source and available at https://gitlab.com/bcm-hgsc/neptune .

Genetic testing in ambulatory cardiology clinics reveals high rate of findings with clinical management implications.

PURPOSE: Cardiovascular disease (CVD) is the leading cause of death in adults in the United States, yet the benefits of genetic testing are not universally accepted. METHODS: We developed the "HeartCare" panel of genes associated with CVD, evaluating high-penetrance Mendelian conditions, coronary artery disease (CAD) polygenic risk, LPA gene polymorphisms, and specific pharmacogenetic (PGx) variants. We enrolled 709 individuals from cardiology clinics at Baylor College of Medicine, and samples were analyzed in a CAP/CLIA-certified laboratory. Results were returned to the ordering physician and uploaded to the electronic medical record. RESULTS: Notably, 32% of patients had a genetic finding with clinical management implications, even after excluding PGx results, including 9% who were molecularly diagnosed with a Mendelian condition. Among surveyed physicians, 84% reported medical management changes based on these results, including specialist referrals, cardiac tests, and medication changes. LPA polymorphisms and high polygenic risk of CAD were found in 20% and 9% of patients, respectively, leading to diet, lifestyle, and other changes. Warfarin and simvastatin pharmacogenetic variants were present in roughly half of the cohort. CONCLUSION: Our results support the use of genetic information in routine cardiovascular health management and provide a roadmap for accompanying research.

Exome sequencing in children with clinically suspected maturity-onset diabetes of the young.

OBJECTIVE: Commercial gene panels identify pathogenic variants in as low as 27% of patients suspected to have MODY, suggesting the role of yet unidentified pathogenic variants. We sought to identify novel gene variants associated with MODY. RESEARCH DESIGN AND METHODS: We recruited 10 children with a clinical suspicion of MODY but non-diagnostic commercial MODY gene panels. We performed exome sequencing (ES) in them and their parents. RESULTS: Mean age at diabetes diagnosis was 10 (± 3.8) years. Six were females; 4 were non-Hispanic white, 5 Hispanic, and 1 Asian. Our variant prioritization analysis identified a pathogenic, de novo variant in INS (c.94G > A, p.Gly32Ser), confirmed by Sanger sequencing, in a proband who was previously diagnosed with "autoantibody-negative type 1 diabetes (T1D)" at 3 y/o. This rare variant, absent in the general population (gnomAD database), has been reported previously in neonatal diabetes. We also identified a frameshift deletion (c.2650delC, p.Gln884AsnfsTer57) in RFX6 in a child with a previous diagnosis of "autoantibody-negative T1D" at 12 y/o. The variant was inherited from the mother, who was diagnosed with "thin type 2 diabetes" at 25 y/o. Heterozygous protein-truncating variants in RFX6 gene have been recently reported in individuals with MODY. CONCLUSIONS: We diagnosed two patients with MODY using ES in children initially classified as "T1D". One has a likely pathogenic novel gene variant not previously associated with MODY. We demonstrate the clinical utility of ES in patients with clinical suspicion of MODY.

Phenotypic expansion in KIF1A-related dominant disorders: A description of novel variants and review of published cases.

KIF1A is a molecular motor for membrane-bound cargo important to the development and survival of sensory neurons. KIF1A dysfunction has been associated with several Mendelian disorders with a spectrum of overlapping phenotypes, ranging from spastic paraplegia to intellectual disability. We present a novel pathogenic in-frame deletion in the KIF1A molecular motor domain inherited by two affected siblings from an unaffected mother with apparent germline mosaicism. We identified eight additional cases with heterozygous, pathogenic KIF1A variants ascertained from a local data lake. Our data provide evidence for the expansion of KIF1A-associated phenotypes to include hip subluxation and dystonia as well as phenotypes observed in only a single case: gelastic cataplexy, coxa valga, and double collecting system. We review the literature and suggest that KIF1A dysfunction is better understood as a single neuromuscular disorder with variable involvement of other organ systems than a set of discrete disorders converging at a single locus.

Marker-assisted selection of qMrdd8 to improve maize resistance to rough dwarf disease.

Maize rough dwarf disease (MRDD) is caused by viruses in the Fijivirus genus in the family Reoviridae. MRDD resistance can be improved by a combination of marker-assisted selection (MAS) and conventional breeding strategies. In our previous study, we fine-mapped a major QTL qMrdd8 and developed the functional Indel marker IDP25K. In the present study, qMrdd8 from the donor parent X178 was introgressed into elite inbred lines derived from the three corn heterotic groups using multi-generation backcrossing and MAS. Recipient lines included Huangzao4, Chang7-2, Ye478, Zheng58, Zhonghuang68, B73, and Ji846. Markers used for foreground selection included IDRQ4, IDRQ47, IDP25K, and IDP27K. Background selection was carried out in the BC3 or BC4 using 107 SSR markers to select lines with the highest rate of recovery of the particular recurrent parent genome. Plants from BC4F2 and BC3F2 that carried the shortest qMrdd8 interval from X178 and those with the highest rate of recovery of the recurrent parent genome were then selected to create converted homozygous inbred lines. In 2017, seven converted inbred lines and five hybrids exhibited enhanced resistance to MRDD, while other agronomic traits were not affected under nonpathogenic stress conditions. Thus, the MRDD resistance allele at the qMrdd8 locus, or IDP25K, should be valuable for maize breeding programs in China.

Mutations in PI3K110δ cause impaired natural killer cell function partially rescued by rapamycin treatment.

BACKGROUND: Heterozygous gain-of-function mutations in PI3K110δ lead to lymphadenopathy, lymphoid hyperplasia, EBV and cytomegalovirus viremia, and sinopulmonary infections. OBJECTIVE: The known role of natural killer (NK) cell function in the control of EBV and cytomegalovirus prompted us to investigate the functional and phenotypic effects of PI3K110δ mutations on NK cell subsets and cytotoxic function. METHODS: Mutations in patients were identified by using whole-exome or targeted sequencing. We performed NK cell phenotyping and functional analysis of patients' cells using flow cytometry, standard Cr51 cytotoxicity assays, and quantitative confocal microscopy. RESULTS: PI3K110δ mutations led to an altered NK cell developmental phenotype and cytotoxic dysfunction. Impaired NK cell cytotoxicity was due to decreased conjugate formation with susceptible target cells and abrogated activation of cell machinery required for target cell killing. These defects were restored partially after initiation of treatment with rapamycin in 3 patients. CONCLUSION: We describe novel NK cell functional deficiency caused by PI3K110δ mutation, which is a likely contributor to the severe viremia observed in these patients. Rapamycin treatment partially restores NK cell function, providing a further rationale for its use in patients with this disease.

The phenotypic spectrum of Xia-Gibbs syndrome.

Xia-Gibbs syndrome (XGS: OMIM # 615829) results from de novo truncating mutations within the AT-Hook DNA Binding Motif Containing 1 gene (AHDC1). To further define the phenotypic and molecular spectrum of this disorder, we established an XGS Registry and recruited patients from a worldwide pool of approximately 60 probands. Additional de novo truncating mutations were observed among 25 individuals, extending both the known number of mutation sites and the range of positions within the coding region that were sensitive to alteration. Detailed phenotypic examination of 20 of these patients via clinical records review and data collection from additional surveys showed a wider age range than previously described. Data from developmental milestones showed evidence for delayed speech and that males were more severely affected. Neuroimaging from six available patients showed an associated thinning of the corpus callosum and posterior fossa cysts. An increased risk of both scoliosis and seizures relative to the population burden was also observed. Data from a modified autism screening tool revealed that XGS shares significant overlap with autism spectrum disorders. These details of the phenotypic heterogeneity of XGS implicate specific genotype/phenotype correlations and suggest potential clinical management guidelines.

Pharmacogenetic characterization of naturally occurring germline NT5C1A variants to chemotherapeutic nucleoside analogs.

BACKGROUND: Mutations or alterations in expression of the 5' nucleotidase gene family can lead to altered responses to treatment with nucleoside analogs. While investigating leukemia susceptibility genes, we discovered a very rare p.L254P NT5C1A missense variant in the substrate recognition motif. Given the paucity of cellular drug response data from the NT5C1A germline variation, we characterized p.L254P and eight rare variants of NT5C1A from genomic databases. MATERIALS AND METHODS: Through lentiviral infection, we created HEK293 cell lines that stably overexpress wild-type NT5C1A, p.L254P, or eight NT5C1A variants reported in the National Heart Lung and Blood Institute Exome Variant Server (one truncating and seven missense). IC50 values were determined by cytotoxicity assays after exposure to chemotherapeutic nucleoside analogs (cladribine, gemcitabine, 5-fluorouracil). In addition, we used structure-based homology modeling to generate a three-dimensional model for the C-terminal region of NT5C1A. RESULTS: The p.R180X (truncating), p.A214T, and p.L254P missense changes were the only variants that significantly impaired protein function across all nucleotide analogs tested (>5-fold difference vs. wild-type; P<0.05). Several of the remaining variants individually showed differential effects (both more and less resistant) across the analogs tested. The homology model provided a structural framework to understand the impact of NT5C1A mutants on catalysis and drug processing. The model predicted active site residues within NT5C1A motif III and we experimentally confirmed that p.K314 (not p.K320) is required for NT5C1A activity. CONCLUSION: We characterized germline variation and predicted protein structures of NT5C1A. Individual missense changes showed considerable variation in response to the different nucleoside analogs tested, which may impact patients' responses to treatment.

Dual transcriptome analysis reveals insights into the response to Rice black-streaked dwarf virus in maize.

Maize rough dwarf disease (MRDD) is a viral infection that results in heavy yield losses in maize worldwide, particularly in the summer maize-growing regions of China. MRDD is caused by the Rice black-streaked dwarf virus (RBSDV). In the present study, analyses of microRNAs (miRNAs), the degradome, and transcriptome sequences were used to elucidate the RBSDV-responsive pathway(s) in maize. Genomic analysis indicated that the expression of three non-conserved and 28 conserved miRNAs, representing 17 known miRNA families and 14 novel miRNAs, were significantly altered in response to RBSDV when maize was inoculated at the V3 (third leaf) stage. A total of 99 target transcripts from 48 genes of 10 known miRNAs were found to be responsive to RBSDV infection. The annotations of these target genes include a SQUAMOSA promoter binding (SPB) protein, a P450 reductase, an oxidoreductase, and a ubiquitin-related gene, among others. Characterization of the entire transcriptome suggested that a total of 28 and 1085 differentially expressed genes (DEGs) were detected at 1.5 and 3.0 d, respectively, after artificial inoculation with RBSDV. The expression patterns of cell wall- and chloroplast-related genes, and disease resistance- and stress-related genes changed significantly in response to RBSDV infection. The negatively regulated genes GRMZM2G069316 and GRMZM2G031169, which are the target genes for miR169i-p5 and miR8155, were identified as a nucleolin and a NAD(P)-binding Rossmann-fold superfamily protein in maize, respectively. The gene ontology term GO:0003824, including GRMZM2G031169 and other 51 DEGs, was designated as responsive to RBSDV.

Molecular variation and expansion of a rice black-streaked dwarf virus population based on analysis of segment 1 in Jining, China.

To analyze the variation in rice black-streaked dwarf virus (RBSDV) in an area with high incidence of maize rough dwarf disease (MRDD), the RBSDV S1 segment in a collection of 100 maize isolates (sample population A100) from Jining, Shandong Province, was sequenced. An additional 21 maize and rice isolates (subpopulation B21) that were sampled from nine other geographic locations in China in 2012 and 2013 were used as a control. A total of 914 nucleotide mutations, including 239 singleton variable and 675 parsimony-informative sites were detected among the segment 1 (S1) sequences from A100. A total of 614 nucleotide mutation sites including 164 singleton variable and 450 parsimony-informative sites were detected among the S1 sequences from B21, while 97.55 % of the parsimony-informative sites from B21 were also detected in A100. The nucleotide sequence diversities of A100 (π = 0.0479) and B21 (π = 0.0396) were significantly different (P = 0.0002) but showed similar trends. Phylogenetic analysis showed that the 121 RBSDV isolates could be classified into two groups based on their S1 sequences, independent of subpopulation, with a combination of host species and locations. A100 and B21 were under the same level of negative and purifying selection, with Ka/Ks ratios of 0.0337 and 0.0369, respectively. The combined RBSDV population, including 121 isolates, was expanding, with negative values for Tajima's D, Fu and Li's D, and Fu and Li's F in both A100 and B21, except Tajima's D in A100. Based on S1, the RBSDV population in China has long-term phytogeographic stability, and there do not appear to be any newly-emerging strains.

PacBio-LITS: a large-insert targeted sequencing method for characterization of human disease-associated chromosomal structural variations.

BACKGROUND: Generation of long (>5 Kb) DNA sequencing reads provides an approach for interrogation of complex regions in the human genome. Currently, large-insert whole genome sequencing (WGS) technologies from Pacific Biosciences (PacBio) enable analysis of chromosomal structural variations (SVs), but the cost to achieve the required sequence coverage across the entire human genome is high. RESULTS: We developed a method (termed PacBio-LITS) that combines oligonucleotide-based DNA target-capture enrichment technologies with PacBio large-insert library preparation to facilitate SV studies at specific chromosomal regions. PacBio-LITS provides deep sequence coverage at the specified sites at substantially reduced cost compared with PacBio WGS. The efficacy of PacBio-LITS is illustrated by delineating the breakpoint junctions of low copy repeat (LCR)-associated complex structural rearrangements on chr17p11.2 in patients diagnosed with Potocki-Lupski syndrome (PTLS; MIM#610883). We successfully identified previously determined breakpoint junctions in three PTLS cases, and also were able to discover novel junctions in repetitive sequences, including LCR-mediated breakpoints. The new information has enabled us to propose mechanisms for formation of these structural variants. CONCLUSIONS: The new method leverages the cost efficiency of targeted capture-sequencing as well as the mappability and scaffolding capabilities of long sequencing reads generated by the PacBio platform. It is therefore suitable for studying complex SVs, especially those involving LCRs, inversions, and the generation of chimeric Alu elements at the breakpoints. Other genomic research applications, such as haplotype phasing and small insertion and deletion validation could also benefit from this technology.

Closing the gap: Solving complex medically relevant genes at scale.

Comprehending the mechanism behind human diseases with an established heritable component represents the forefront of personalized medicine. Nevertheless, numerous medically important genes are inaccurately represented in short-read sequencing data analysis due to their complexity and repetitiveness or the so-called 'dark regions' of the human genome. The advent of PacBio as a long-read platform has provided new insights, yet HiFi whole-genome sequencing (WGS) cost remains frequently prohibitive. We introduce a targeted sequencing and analysis framework, Twist Alliance Dark Genes Panel (TADGP), designed to offer phased variants across 389 medically important yet complex autosomal genes. We highlight TADGP accuracy across eleven control samples and compare it to WGS. This demonstrates that TADGP achieves variant calling accuracy comparable to HiFi-WGS data, but at a fraction of the cost. Thus, enabling scalability and broad applicability for studying rare diseases or complementing previously sequenced samples to gain insights into these complex genes. TADGP revealed several candidate variants across all cases and provided insight into LPA diversity when tested on samples from rare disease and cardiovascular disease cohorts. In both cohorts, we identified novel variants affecting individual disease-associated genes (e.g., IKZF1, KCNE1). Nevertheless, the annotation of the variants across these 389 medically important genes remains challenging due to their underrepresentation in ClinVar and gnomAD. Consequently, we also offer an annotation resource to enhance the evaluation and prioritization of these variants. Overall, we can demonstrate that TADGP offers a cost-efficient and scalable approach to routinely assess the dark regions of the human genome with clinical relevance.

GmNAC5, a NAC transcription factor, is a transient response regulator induced by abiotic stress in soybean.

GmNAC5 is a member of NAM subfamily belonging to NAC transcription factors in soybean (Glycine max (L.) Merr.). Studies on NAC transcription factors have shown that this family functioned in the regulation of shoot apical meristem (SAM), hormone signalling, and stress responses. In this study, we examined the expression levels of GmNAC5. GmNAC5 was highly expressed in the roots and immature seeds, especially strongly in immature seeds of 40 days after flowering. In addition, we found that GmNAC5 was induced by mechanical wounding, high salinity, and cold treatments but was not induced by abscisic acid (ABA). The subcellular localization assay suggested that GmNAC5 was targeted at nucleus. Together, it was suggested that GmNAC5 might be involved in seed development and abiotic stress responses in soybean.

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.