Identification of candidate gene associated with maize northern leaf blight resistance in a multi-parent population.

KEY MESSAGE: QTL mapping combined with genome-wide association studies, revealed a potential candidate gene for  resistance to northern leaf blight in the tropical CATETO-related maize line YML226, providing a basis for marker-assisted selection of maize varieties Northern leaf blight (NLB) is a foliar disease that can cause severe yield losses in maize. Identifying and utilizing NLB-resistant genes is the most effective way to prevent and control this disease. In this study, five important inbred lines of maize were used as parental lines to construct a multi-parent population for the identification of NLB-resistant loci. QTL mapping and GWAS analysis revealed that QTL qtl_YML226_1, which had the largest phenotypic variance explanation (PVE) of 9.28%, and SNP 5-49,193,921 were co-located in the CATETO-related line YML226. This locus was associated with the candidate gene Zm00001d014471, which encodes a pentatricopeptide repeat (PPR) protein. In the coding region of Zm00001d014471, YML226 had more specific SNPs than the other parental lines. qRT-PCR showed that the relative expressions of Zm00001d014471 in inoculated and uninoculated leaves of YML226 were significantly higher, indicating that the expression of the candidate gene was correlated with NLB resistance. The analysis showed that the higher expression level in YML226 might be caused by SNP mutations. This study identified NLB resistance candidate loci and genes in the tropical maize inbred line YML226 derived from the CATETO germplasm, thereby providing a theoretical basis for using modern marker-assisted breeding techniques to select genetic resources resistant to NLB.

Identification of candidate genes controlling cold tolerance at the early seedling stage from Dongxiang wild rice by QTL mapping, BSA-Seq and RNA-Seq.

BACKGROUND: The cold tolerance of rice is closely related to its production and geographic distribution. The identification of cold tolerance-related genes is of important significance for developing cold-tolerant rice. Dongxiang wild rice (Oryza rufipogon Griff.) (DXWR) is well-adapted to the cold climate of northernmost-latitude habitats ever found in the world, and is one of the most valuable rice germplasms for cold tolerance improvement. RESULTS: Transcriptome analysis revealed genes differentially expressed between Xieqingzao B (XB; a cold sensitive variety) and 19H19 (derived from an interspecific cross between DXWR and XB) in the room temperature (RT), low temperature (LT), and recovery treatments. The results demonstrated that chloroplast genes might be involved in the regulation of cold tolerance in rice. A high-resolution SNP genetic map was constructed using 120 BC5F2 lines derived from a cross between 19H19 and XB based on the genotyping-by-sequencing (GBS) technique. Two quantitative trait loci (QTLs) for cold tolerance at the early seedling stage (CTS), qCTS12 and qCTS8, were detected. Moreover, a total of 112 candidate genes associated with cold tolerance were identified based on bulked segregant analysis sequencing (BSA-seq). These candidate genes were divided into eight functional categories, and the expression trend of candidate genes related to 'oxidation-reduction process' and 'response to stress' differed between XB and 19H19 in the RT, LT and recovery treatments. Among these candidate genes, the expression level of LOC_Os12g18729 in 19H19 (related to 'response to stress') decreased in the LT treatment but restored and enhanced during the recovery treatment whereas the expression level of LOC_Os12g18729 in XB declined during recovery treatment. Additionally, XB contained a 42-bp deletion in the third exon of LOC_Os12g18729, and the genotype of BC5F2 individuals with a survival percentage (SP) lower than 15% was consistent with that of XB. Weighted gene coexpression network analysis (WGCNA) and modular regulatory network learning with per gene information (MERLIN) algorithm revealed a gene interaction/coexpression network regulating cold tolerance in rice. In the network, differentially expressed genes (DEGs) related to 'oxidation-reduction process', 'response to stress' and 'protein phosphorylation' interacted with LOC_Os12g18729. Moreover, the knockout mutant of LOC_Os12g18729 decreased cold tolerance in early rice seedling stage signifcantly compared with that of wild type. CONCLUSIONS: In general, study of the genetic basis of cold tolerance of rice is important for the development of cold-tolerant rice varieties. In the present study, QTL mapping, BSA-seq and RNA-seq were integrated to identify two CTS QTLs qCTS8 and qCTS12. Furthermore, qRT-PCR, genotype sequencing and knockout analysis indicated that LOC_Os12g18729 could be the candidate gene of qCTS12. These results are expected to further exploration of the genetic mechanism of CTS in rice and improve cold tolerance of cultivated rice by introducing the cold tolerant genes from DXWR through marker-assisted selection.

Partitioning and aggregating cross-tissue and tissue-specific genetic effects to identify gene-trait associations.

TWAS have shown great promise in extending GWAS loci to a functional understanding of disease mechanisms. In an effort to fully unleash the TWAS and GWAS information, we propose MTWAS, a statistical framework that partitions and aggregates cross-tissue and tissue-specific genetic effects in identifying gene-trait associations. We introduce a non-parametric imputation strategy to augment the inaccessible tissues, accommodating complex interactions and non-linear expression data structures across various tissues. We further classify eQTLs into cross-tissue eQTLs and tissue-specific eQTLs via a stepwise procedure based on the extended Bayesian information criterion, which is consistent under high-dimensional settings. We show that MTWAS significantly improves the prediction accuracy across all 47 tissues of the GTEx dataset, compared with other single-tissue and multi-tissue methods, such as PrediXcan, TIGAR, and UTMOST. Applying MTWAS to the DICE and OneK1K datasets with bulk and single-cell RNA sequencing data on immune cell types showcases consistent improvements in prediction accuracy. MTWAS also identifies more predictable genes, and the improvement can be replicated with independent studies. We apply MTWAS to 84 UK Biobank GWAS studies, which provides insights into disease etiology.

Unraveling the genetic basis of quantitative resistance to diseases in tomato: a meta-QTL analysis and mining of transcript profiles.

KEY MESSAGE: Whole-genome QTL mining and meta-analysis in tomato for resistance to bacterial and fungal diseases identified 73 meta-QTL regions with significantly refined/reduced confidence intervals. Tomato production is affected by a range of biotic stressors, causing yield losses and quality reductions. While sources of genetic resistance to many tomato diseases have been identified and characterized, stability of the resistance genes or quantitative trait loci (QTLs) across the resources has not been determined. Here, we examined 491 QTLs previously reported for resistance to tomato diseases in 40 independent studies and 54 unique mapping populations. We identified 29 meta-QTLs (MQTLs) for resistance to bacterial pathogens and 44 MQTLs for resistance to fungal pathogens, and were able to reduce the average confidence interval (CI) of the QTLs by 4.1-fold and 6.7-fold, respectively, compared to the average CI of the original QTLs. The corresponding physical length of the CIs of MQTLs ranged from 56 kb to 6.37 Mb, with a median of 921 kb, of which 27% had a CI lower than 500 kb and 53% had a CI lower than 1 Mb. Comparison of defense responses between tomato and Arabidopsis highlighted 73 orthologous genes in the MQTL regions, which were putatively determined to be involved in defense against bacterial and fungal diseases. Intriguingly, multiple genes were identified in some MQTL regions that are implicated in plant defense responses, including PR-P2, NDR1, PDF1.2, Pip1, SNI1, PTI5, NSL1, DND1, CAD1, SlACO, DAD1, SlPAL, Ph-3, EDS5/SID1, CHI-B/PR-3, Ph-5, ETR1, WRKY29, and WRKY25. Further, we identified a number of candidate resistance genes in the MQTL regions that can be useful for both marker/gene-assisted breeding as well as cloning and genetic transformation.

Genome-wide association study reveals marker-trait associations for major agronomic traits in proso millet (Panicum miliaceum L.).

MAIN CONCLUSION: The pilot-scale genome-wide association study in the US proso millet identified twenty marker-trait associations for five morpho-agronomic traits identifying genomic regions for future studies (e.g. molecular breeding and map-based cloning). Proso millet (Panicum miliaceum L.) is an ancient grain recognized for its excellent water-use efficiency and short growing season. It is an indispensable part of the winter wheat-based dryland cropping system in the High Plains of the USA. Its grains are endowed with high nutritional and health-promoting properties, making it increasingly popular in the global market for healthy grains. There is a dearth of genomic resources in proso millet for developing molecular tools to complement conventional breeding for developing high-yielding varieties. Genome-wide association study (GWAS) is a widely used method to dissect the genetics of complex traits. In this pilot study of the first-ever GWAS in the US proso millet, 71 globally diverse genotypes of 109 the US proso millet core collection were evaluated for five major morpho-agronomic traits at two locations in western Nebraska, and GWAS was conducted to identify single nucleotide polymorphisms (SNPs) associated with these traits. Analysis of variance showed that there was a significant difference among the genotypes, and all five traits were also found to be highly correlated with each other. Sequence reads from genotyping-by-sequencing (GBS) were used to identify 11,147 high-quality bi-allelic SNPs. Population structure analysis with those SNPs showed stratification within the core collection. The GWAS identified twenty marker-trait associations (MTAs) for the five traits. Twenty-nine putative candidate genes associated with the five traits were also identified. These genomic regions can be used to develop genetic markers for marker-assisted selection in proso millet breeding.

Mapping and annotating genomic loci to prioritize genes and implicate distinct polygenic adaptations for skin color.

Evidence for adaptation of human skin color to regional ultraviolet radiation suggests shared and distinct genetic variants across populations. However, skin color evolution and genetics in East Asians are understudied. We quantified skin color in 48,433 East Asians using image analysis and identified associated genetic variants and potential causal genes for skin color as well as their polygenic interplay with sun exposure. This genome-wide association study (GWAS) identified 12 known and 11 previously unreported loci and SNP-based heritability was 23-24%. Potential causal genes were determined through the identification of nonsynonymous variants, colocalization with gene expression in skin tissues, and expression levels in melanocytes. Genomic loci associated with pigmentation in East Asians substantially diverged from European populations, and we detected signatures of polygenic adaptation. This large GWAS for objectively quantified skin color in an East Asian population improves understanding of the genetic architecture and polygenic adaptation of skin color and prioritizes potential causal genes.

Identifying therapeutic target genes for migraine by systematic druggable genome-wide Mendelian randomization.

BACKGROUND: Currently, the treatment and prevention of migraine remain highly challenging. Mendelian randomization (MR) has been widely used to explore novel therapeutic targets. Therefore, we performed a systematic druggable genome-wide MR to explore the potential therapeutic targets for migraine. METHODS: We obtained data on druggable genes and screened for genes within brain expression quantitative trait locis (eQTLs) and blood eQTLs, which were then subjected to two-sample MR analysis and colocalization analysis with migraine genome-wide association studies data to identify genes highly associated with migraine. In addition, phenome-wide research, enrichment analysis, protein network construction, drug prediction, and molecular docking were performed to provide valuable guidance for the development of more effective and targeted therapeutic drugs. RESULTS: We identified 21 druggable genes significantly associated with migraine (BRPF3, CBFB, CDK4, CHD4, DDIT4, EP300, EPHA5, FGFRL1, FXN, HMGCR, HVCN1, KCNK5, MRGPRE, NLGN2, NR1D1, PLXNB1, TGFB1, TGFB3, THRA, TLN1 and TP53), two of which were significant in both blood and brain (HMGCR and TGFB3). The results of phenome-wide research showed that HMGCR was highly correlated with low-density lipoprotein, and TGFB3 was primarily associated with insulin-like growth factor 1 levels. CONCLUSIONS: This study utilized MR and colocalization analysis to identify 21 potential drug targets for migraine, two of which were significant in both blood and brain. These findings provide promising leads for more effective migraine treatments, potentially reducing drug development costs.

A genome-wide association study reveals novel loci and candidate genes associated with plant height variation in Medicago sativa.

BACKGROUND: Plant height (PH) is an important agronomic trait influenced by a complex genetic network. However, the genetic basis for the variation in PH in Medicago sativa remains largely unknown. In this study, a comprehensive genome-wide association analysis was performed to identify genomic regions associated with PH using a diverse panel of 220 accessions of M. sativa worldwide. RESULTS: Our study identified eight novel single nucleotide polymorphisms (SNPs) significantly associated with PH evaluated in five environments, explaining 8.59-12.27% of the phenotypic variance. Among these SNPs, the favorable genotype of chr6__31716285 had a low frequency of 16.4%. Msa0882400, located proximal to this SNP, was annotated as phosphate transporter 3;1, and its role in regulating alfalfa PH was supported by transcriptome and candidate gene association analysis. In addition, 21 candidate genes were annotated within the associated regions that are involved in various biological processes related to plant growth and development. CONCLUSIONS: Our findings provide new molecular markers for marker-assisted selection in M. sativa breeding programs. Furthermore, this study enhances our understanding of the underlying genetic and molecular mechanisms governing PH variations in M. sativa.

Systems genetics of metabolic health in the BXD mouse genetic reference population.

Susceptibility to metabolic syndrome (MetS) is dependent on genetics, environment, and gene-by-environment interactions, rendering the study of underlying mechanisms challenging. The majority of experiments in model organisms do not incorporate genetic variation and lack specific evaluation criteria for MetS. Here, we derived a continuous metric, the metabolic health score (MHS), based on standard clinical parameters and defined its molecular signatures in the liver and circulation. In human UK Biobank, the MHS associated with MetS status and was predictive of future disease incidence, even in individuals without MetS. Using quantitative trait locus analyses in mice, we found two MHS-associated genetic loci and replicated them in unrelated mouse populations. Through a prioritization scheme in mice and human genetic data, we identified TNKS and MCPH1 as candidates mediating differences in the MHS. Our findings provide insights into the molecular mechanisms sustaining metabolic health across species and uncover likely regulators. A record of this paper's transparent peer review process is included in the supplemental information.

Single-cell and genome-wide Mendelian randomization identifies causative genes for gout.

BACKGROUND: Gout is a prevalent manifestation of metabolic osteoarthritis induced by elevated blood uric acid levels. The purpose of this study was to investigate the mechanisms of gene expression regulation in gout disease and elucidate its pathogenesis. METHODS: The study integrated gout genome-wide association study (GWAS) data, single-cell transcriptomics (scRNA-seq), expression quantitative trait loci (eQTL), and methylation quantitative trait loci (mQTL) data for analysis, and utilized two-sample Mendelian randomization study to comprehend the causal relationship between proteins and gout. RESULTS: We identified 17 association signals for gout at unique genetic loci, including four genes related by protein-protein interaction network (PPI) analysis: TRIM46, THBS3, MTX1, and KRTCAP2. Additionally, we discerned 22 methylation sites in relation to gout. The study also found that genes such as TRIM46, MAP3K11, KRTCAP2, and TM7SF2 could potentially elevate the risk of gout. Through a Mendelian randomization (MR) analysis, we identified three proteins causally associated with gout: ADH1B, BMP1, and HIST1H3A. CONCLUSION: According to our findings, gout is linked with the expression and function of particular genes and proteins. These genes and proteins have the potential to function as novel diagnostic and therapeutic targets for gout. These discoveries shed new light on the pathological mechanisms of gout and clear the way for future research on this condition.

Shared genetic effect of kidney function on bipolar and major depressive disorders: a large-scale genome-wide cross-trait analysis.

BACKGROUND: Epidemiological studies have revealed a significant association between impaired kidney function and certain mental disorders, particularly bipolar disorder (BIP) and major depressive disorder (MDD). However, the evidence regarding shared genetics and causality is limited due to residual confounding and reverse causation. METHODS: In this study, we conducted a large-scale genome-wide cross-trait association study to investigate the genetic overlap between 5 kidney function biomarkers (eGFRcrea, eGFRcys, blood urea nitrogen (BUN), serum urate, and UACR) and 2 mental disorders (MDD, BIP). Summary-level data of European ancestry were extracted from UK Biobank, Chronic Kidney Disease Genetics Consortium, and Psychiatric Genomics Consortium. RESULTS: Using LD score regression, we found moderate but significant genetic correlations between kidney function biomarker traits on BIP and MDD. Cross-trait meta-analysis identified 1 to 19 independent significant loci that were found shared among 10 pairs of 5 kidney function biomarkers traits and 2 mental disorders. Among them, 3 novel genes: SUFU, IBSP, and PTPRJ, were also identified in transcriptome-wide association study analysis (TWAS), most of which were observed in the nervous and digestive systems (FDR < 0.05). Pathway analysis showed the immune system could play a role between kidney function biomarkers and mental disorders. Bidirectional mendelian randomization analysis suggested a potential causal relationship of kidney function biomarkers on BIP and MDD. CONCLUSIONS: In conclusion, the study demonstrated that both BIP and MDD shared genetic architecture with kidney function biomarkers, providing new insights into their genetic architectures and suggesting that larger GWASs are warranted.

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

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

Dual-Model GWAS Analysis and Genomic Selection of Maize Flowering Time-Related Traits.

An appropriate flowering period is an important selection criterion in maize breeding. It plays a crucial role in the ecological adaptability of maize varieties. To explore the genetic basis of flowering time, GWAS and GS analyses were conducted using an associating panel consisting of 379 multi-parent DH lines. The DH population was phenotyped for days to tasseling (DTT), days to pollen-shedding (DTP), and days to silking (DTS) in different environments. The heritability was 82.75%, 86.09%, and 85.26% for DTT, DTP, and DTS, respectively. The GWAS analysis with the FarmCPU model identified 10 single-nucleotide polymorphisms (SNPs) distributed on chromosomes 3, 8, 9, and 10 that were significantly associated with flowering time-related traits. The GWAS analysis with the BLINK model identified seven SNPs distributed on chromosomes 1, 3, 8, 9, and 10 that were significantly associated with flowering time-related traits. Three SNPs 3_198946071, 9_146646966, and 9_152140631 showed a pleiotropic effect, indicating a significant genetic correlation between DTT, DTP, and DTS. A total of 24 candidate genes were detected. A relatively high prediction accuracy was achieved with 100 significantly associated SNPs detected from GWAS, and the optimal training population size was 70%. This study provides a better understanding of the genetic architecture of flowering time-related traits and provides an optimal strategy for GS.

Assessing Falling Number Stability Increases the Genomic Prediction Ability of Pre-Harvest Sprouting Resistance in Common Winter Wheat.

Pre-harvest sprouting (PHS) resistance is a complex trait, and many genes influencing the germination process of winter wheat have already been described. In the light of interannual climate variation, breeding for PHS resistance will remain mandatory for wheat breeders. Several tests and traits are used to assess PHS resistance, i.e., sprouting scores, germination index, and falling number (FN), but the variation of these traits is highly dependent on the weather conditions during field trials. Here, we present a method to assess falling number stability (FNS) employing an after-ripening period and the wetting of the kernels to improve trait variation and thus trait heritability. Different genome-based prediction scenarios within and across two subsequent seasons based on overall 400 breeding lines were applied to assess the predictive abilities of the different traits. Based on FNS, the genome-based prediction of the breeding values of wheat breeding material showed higher correlations across seasons (r=0.505-0.548) compared to those obtained for other traits for PHS assessment (r=0.216-0.501). By weighting PHS-associated quantitative trait loci (QTL) in the prediction model, the average predictive abilities for FNS increased from 0.585 to 0.648 within the season 2014/2015 and from 0.649 to 0.714 within the season 2015/2016. We found that markers in the Phs-A1 region on chromosome 4A had the highest effect on the predictive abilities for FNS, confirming the influence of this QTL in wheat breeding material, whereas the dwarfing genes Rht-B1 and Rht-D1 and the wheat-rye translocated chromosome T1RS.1BL exhibited effects, which are well-known, on FN per se exclusively.

Bayesian genome-wide TWAS with reference transcriptomic data of brain and blood tissues identified 141 risk genes for Alzheimer's disease dementia.

BACKGROUND: Transcriptome-wide association study (TWAS) is an influential tool for identifying genes associated with complex diseases whose genetic effects are likely mediated through transcriptome. TWAS utilizes reference genetic and transcriptomic data to estimate effect sizes of genetic variants on gene expression (i.e., effect sizes of a broad sense of expression quantitative trait loci, eQTL). These estimated effect sizes are employed as variant weights in gene-based association tests, facilitating the mapping of risk genes with genome-wide association study (GWAS) data. However, most existing TWAS of Alzheimer's disease (AD) dementia are limited to studying only cis-eQTL proximal to the test gene. To overcome this limitation, we applied the Bayesian Genome-wide TWAS (BGW-TWAS) method to leveraging both cis- and trans- eQTL of brain and blood tissues, in order to enhance mapping risk genes for AD dementia. METHODS: We first applied BGW-TWAS to the Genotype-Tissue Expression (GTEx) V8 dataset to estimate cis- and trans- eQTL effect sizes of the prefrontal cortex, cortex, and whole blood tissues. Estimated eQTL effect sizes were integrated with the summary data of the most recent GWAS of AD dementia to obtain BGW-TWAS (i.e., gene-based association test) p-values of AD dementia per gene per tissue type. Then we used the aggregated Cauchy association test to combine TWAS p-values across three tissues to obtain omnibus TWAS p-values per gene. RESULTS: We identified 85 significant genes in prefrontal cortex, 82 in cortex, and 76 in whole blood that were significantly associated with AD dementia. By combining BGW-TWAS p-values across these three tissues, we obtained 141 significant risk genes including 34 genes primarily due to trans-eQTL and 35 mapped risk genes in GWAS Catalog. With these 141 significant risk genes, we detected functional clusters comprised of both known mapped GWAS risk genes of AD in GWAS Catalog and our identified TWAS risk genes by protein-protein interaction network analysis, as well as several enriched phenotypes related to AD. CONCLUSION: We applied BGW-TWAS and aggregated Cauchy test methods to integrate both cis- and trans- eQTL data of brain and blood tissues with GWAS summary data, identifying 141 TWAS risk genes of AD dementia. These identified risk genes provide novel insights into the underlying biological mechanisms of AD dementia and potential gene targets for therapeutics development.

Genetic mapping of regions associated with root system architecture in rice using MutMap QTL-seq.

The root system architecture is an important complex trait in rice. With changing climatic conditions and soil nutrient deficiencies, there is an immediate need to breed nutrient-use-efficient rice varieties with robust root system architectural (RSA) traits. To map the genomic regions associated with crucial component traits of RSA viz. root length and root volume, a biparental F2 mapping population was developed using TI-128, an Ethyl Methane Sulphonate (EMS) mutant of a mega variety BPT-5204 having high root length (RL) and root volume (RV) with wild type BPT-5204. Extreme bulks having high RL and RV and low RL and RV were the whole genome re-sequenced along with parents. Genetic mapping using the MutMap QTL-Seq approach elucidated two genomic intervals on Chr.12 (3.14-3.74 Mb, 18.11-20.85 Mb), and on Chr.2 (23.18-23.68 Mb) as potential regions associated with both RL and RV. The Kompetitive Allele Specific PCR (KASP) assays for SNPs with delta SNP index near 1 were associated with higher RL and RV in the panel of sixty-two genotypes varying in root length and volume. The KASP_SNPs viz. Chr12_S4 (C→T; Chr12:3243938), located in the 3' UTR region of LOC_Os12g06670 encoding a protein kinase domain-containing protein and Chr2_S6 (C→T; Chr2:23181622) present upstream in the regulator of chromosomal condensation protein LOC_Os2g38350. Validation of these genes using qRT-PCR and in-silico studies using various online tools and databases revealed higher expression in TI-128 as compared to BPT- 5204 at the seedling and panicle initiation stages implying the functional role in enhancing RL and RV.

Quantitative trait loci mapping provides insights into the genetic regulation of dendritic cell numbers in mouse tissues.

To explore the influence of genetics on homeostatic regulation of dendritic cell (DC) numbers, we present a screen of DCs and their progenitors in lymphoid and non-lymphoid tissues in Collaborative Cross (CC) and Diversity Outbred (DO) mice. We report 30 and 71 loci with logarithm of the odds (LOD) scores >8.18 and ranging from 6.67 to 8.19, respectively. The analysis reveals the highly polygenic and pleiotropic architecture of this complex trait, including many of the previously identified genetic regulators of DC development and maturation. Two SNPs in genes potentially underlying variation in DC homeostasis, a splice variant in Gramd4 (rs235532740) and a missense variant in Orai3 (rs216659754), are confirmed by gene editing using CRISPR-Cas9. Gramd4 is a central regulator of DC homeostasis that impacts the entire DC lineage, and Orai3 regulates cDC2 numbers in tissues. Overall, the data reveal a large number of candidate genes regulating DC homeostasis in vivo.

Leaf morphology related genes revealed by integrating Pan-transcriptome, GWAS and eQTL analyses in a Liriodendron population.

Leaf plays an indispensable role in plant development and growth. Although many known genes related to leaf morphology development have been identified, elucidating the complex genetic basis of leaf morphological traits remains a challenge. Liriodendron plants are common ornamental trees due to their unique leaf shapes, while the molecular mechanism underlying Liriodendron leaf morphogenesis has remained unknown. Herein, we firstly constructed a population-level pan-transcriptome of Liriodendron from 81 accessions to explore the expression presence or absence variations (ePAVs), global expression differences at the population level, as well as differentially expressed genes (DEGs) between the Liriodendron chinense and Liriodendron tulipifera accessions. Subsequently, we integrated a genome-wide association study (GWAS), expression quantitative trait loci (eQTL), and transcriptome-wide association study (TWAS) to identify candidate genes related to leaf morphology. Through GWAS analysis, we identified 18 and 17 significant allelic loci in the leaf size and leaf shape modules, respectively. In addition, we discerned 16 candidate genes in relation to leaf morphological traits via TWAS. Further, integrating the co-localization results of GWAS and eQTL, we determined two regulatory hotspot regions, hot88 and hot758, related to leaf size and leaf shape, respectively. Finally, co-expression analysis, eQTL, and linkage mapping together demonstrated that Lchi_4g10795 regulate their own expression levels through cis-eQTL to affect the expression of downstream genes and cooperatively participate in the development of Liriodendron leaf morphology. These findings will improve our understanding of the molecular regulatory mechanism of Liriodendron leaf morphogenesis and will also accelerate molecular breeding of Liriodendron.

Genome-wide expression quantitative trait locus analysis reveals silk-preferential gene regulatory network in maize.

Silk of maize (Zea mays L.) contains diverse metabolites with complicated structures and functions, making it a great challenge to explore the mechanisms of metabolic regulation. Genome-wide identification of silk-preferential genes and investigation of their expression regulation provide an opportunity to reveal the regulatory networks of metabolism. Here, we applied the expression quantitative trait locus (eQTL) mapping on a maize natural population to explore the regulation of gene expression in unpollinated silk of maize. We obtained 3,985 silk-preferential genes that were specifically or preferentially expressed in silk using our population. Silk-preferential genes showed more obvious expression variations compared with broadly expressed genes that were ubiquitously expressed in most tissues. We found that trans-eQTL regulation played a more important role for silk-preferential genes compared to the broadly expressed genes. The relationship between 38 transcription factors and 85 target genes, including silk-preferential genes, were detected. Finally, we constructed a transcriptional regulatory network around the silk-preferential gene Bx10, which was proposed to be associated with response to abiotic stress and biotic stress. Taken together, this study deepened our understanding of transcriptome variation in maize silk and the expression regulation of silk-preferential genes, enhancing the investigation of regulatory networks on metabolic pathways.