Major quantitative trait locus qLA3.1 is related to tomato leaf angle by regulating cell length at the petiole base.

KEY MESSAGE: qLA3.1, controlling leaf angle in tomato, was fine-mapped to an interval of 4.45 kb on chromosome A03, and one gene encoding auxin response factor was identified as a candidate gene. Leaf angle is a crucial trait in plant architecture that plays an important role in achieving optimal plant structure. However, there are limited reports on gene localization, cloning, and the function of plant architecture in horticultural crops, particularly regarding leaf angle. In this study, we selected 'Z3' with erect leaves and 'Heinz1706' with horizontal leaves as the phenotype and cytological observation. We combined bulked segregant analysis and fine genetic mapping to identify a candidate gene, known as, i.e., qLA3.1, which was related to tomato leaf angle. Through multiple analyses, we found that Solyc03g113410 was the most probably candidate for qLA3.1, which encoded the auxin response factor SlARF11 in tomato and was homologous to OsARF11 related to leaf angle in rice. We discovered that silencing SlARF11 resulted in upright leaves, while plants with over-expressed SlARF11 exhibited horizontal leaves. We also found that cultivars with erect leaves had a mutation from base G to base A. Moreover, quantitative analysis of plants treated with hormones indicated that SlARF11 might participate in cell elongation and the activation of genes related to auxin and brassinosteroid pathways. Transcriptome analysis further validated that SlARF11 may regulate leaf angle through hormone signaling pathways. These data support the idea that the auxin response factor SlARF11 may have an important function in tomato leaf petiole angles.

Validation and identification of promising gene specific markers governing foliar disease resistance in groundnut (Arachis hypogaea L.).

BACKGROUND: Groundnut is vulnerable to the major foliar fungal disease viz., late leaf spot (LLS) and rust in kharif season, which results in severe yield losses. Until now, LLS and rust resistance linked markers were developed based on GPBD 4 as a major donor source and were validated in its derivatives only, which restricted their use in marker assisted selection (MAS) involving other donors. METHODS AND RESULTS: The current study focused to validate LLS and rust resistance linked markers employing advanced breeding lines of F6 generation, derived from nine different crosses involving nine diverse parents, to identify potential markers for marker-assisted breeding of LLS and rust resistance in groundnut. Out of 28-trait linked markers used for validation, 8 were polymorphic (28.57%). Marker-trait association (MTA) and Single Marker Analysis (SMA) revealed that the SSR marker pPGPseq5D05 is significantly associated with both LLS (15.8% PVE) and rust (17.5% PVE) resistance, whereas, the marker IPAHM103 is tightly linked with rust resistance (26.8% PVE) alone. In silico analysis revealed that the marker gene for IPAHM103 is a zinc finger protein and the marker gene for pPGPseq5D05 is an ADP-ribosylation factor GTPase-activating protein. Both these protein products impart resistance or tolerance to biotic stress in crop plants. Two other markers namely, GMLQ975 and pPGPseq13A10 were also found to be associated with LLS resistance explaining MTA up to 60%. CONCLUSION: These gene specific markers will enable us to screen more number of germplasm lines or newly developed lines in MAS schemes for LLS and rust resistance using a wide range of resistant sources.

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.

Genome-wide association study and expression of candidate genes for Fe and Zn concentration in sorghum grains.

Sorghum germplasm showed grain Fe and Zn genetic variability, but a few varieties were biofortified with these minerals. This work contributes to narrowing this gap. Fe and Zn concentrations along with 55,068 high-quality GBS SNP data from 140 sorghum accessions were used in this study. Both micronutrients exhibited good variability with respective ranges of 22.09-52.55 ppm and 17.92-43.16 ppm. Significant marker-trait associations were identified on chromosomes 1, 3, and 5. Two major effect SNPs (S01_72265728 and S05_58213541) explained 35% and 32% of Fe and Zn phenotypic variance, respectively. The SNP S01_72265728 was identified in the cytochrome P450 gene and showed a positive effect on Fe accumulation in the kernel, while S05_58213541 was intergenic near Sobic.005G134800 (zinc-binding ribosomal protein) and showed negative effect on Zn. Tissue-specific in silico expression analysis resulted in higher levels of Sobic.003G350800 gene product in several tissues such as leaf, root, flower, panicle, and stem. Sobic.005G188300 and Sobic.001G463800 were expressed moderately at grain maturity and anthesis in leaf, root, panicle, and seed tissues. The candidate genes expressed in leaves, stems, and grains will be targeted to improve grain and stover quality. The haplotypes identified will be useful in forward genetics breeding.

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.

RicePilaf: a post-GWAS/QTL dashboard to integrate pangenomic, coexpression, regulatory, epigenomic, ontology, pathway, and text-mining information to provide functional insights into rice QTLs and GWAS loci.

BACKGROUND: As the number of genome-wide association study (GWAS) and quantitative trait locus (QTL) mappings in rice continues to grow, so does the already long list of genomic loci associated with important agronomic traits. Typically, loci implicated by GWAS/QTL analysis contain tens to hundreds to thousands of single-nucleotide polmorphisms (SNPs)/genes, not all of which are causal and many of which are in noncoding regions. Unraveling the biological mechanisms that tie the GWAS regions and QTLs to the trait of interest is challenging, especially since it requires collating functional genomics information about the loci from multiple, disparate data sources. RESULTS: We present RicePilaf, a web app for post-GWAS/QTL analysis, that performs a slew of novel bioinformatics analyses to cross-reference GWAS results and QTL mappings with a host of publicly available rice databases. In particular, it integrates (i) pangenomic information from high-quality genome builds of multiple rice varieties, (ii) coexpression information from genome-scale coexpression networks, (iii) ontology and pathway information, (iv) regulatory information from rice transcription factor databases, (v) epigenomic information from multiple high-throughput epigenetic experiments, and (vi) text-mining information extracted from scientific abstracts linking genes and traits. We demonstrate the utility of RicePilaf by applying it to analyze GWAS peaks of preharvest sprouting and genes underlying yield-under-drought QTLs. CONCLUSIONS: RicePilaf enables rice scientists and breeders to shed functional light on their GWAS regions and QTLs, and it provides them with a means to prioritize SNPs/genes for further experiments. The source code, a Docker image, and a demo version of RicePilaf are publicly available at https://github.com/bioinfodlsu/rice-pilaf.

Genetic control of root morphology in response to nitrogen across rapeseed diversity.

Rapeseed (Brassica napus L.) is an oil-containing crop of great economic value but with considerable nitrogen requirement. Breeding root systems that efficiently absorb nitrogen from the soil could be a driver to ensure genetic gains for more sustainable rapeseed production. The aim of this study is to identify genomic regions that regulate root morphology in response to nitrate availability. The natural variability offered by 300 inbred lines was screened at two experimental locations. Seedlings grew hydroponically with low or elevated nitrate levels. Fifteen traits related to biomass production and root morphology were measured. On average across the panel, a low nitrate level increased the root-to-shoot biomass ratio and the lateral root length. A large phenotypic variation was observed, along with important heritability values and genotypic effects, but low genotype-by-nitrogen interactions. Genome-wide association study and bulk segregant analysis were used to identify loci regulating phenotypic traits. The first approach nominated 319 SNPs that were combined into 80 QTLs. Three QTLs identified on the A07 and C07 chromosomes were stable across nitrate levels and/or experimental locations. The second approach involved genotyping two groups of individuals from an experimental F2 population created by crossing two accessions with contrasting lateral root lengths. These individuals were found in the tails of the phenotypic distribution. Co-localized QTLs found in both mapping approaches covered a chromosomal region on the A06 chromosome. The QTL regions contained some genes putatively involved in root organogenesis and represent selection targets for redesigning the root morphology of rapeseed.

Novel putative causal mutations associated with fat traits in Nellore cattle uncovered by eQTLs located in open chromatin regions.

Intramuscular fat (IMF) and backfat thickness (BFT) are critical economic traits impacting meat quality. However, the genetic variants controlling these traits need to be better understood. To advance knowledge in this area, we integrated RNA-seq and single nucleotide polymorphisms (SNPs) identified in genomic and transcriptomic data to generate a linkage disequilibrium filtered panel of 553,581 variants. Expression quantitative trait loci (eQTL) analysis revealed 36,916 cis-eQTLs and 14,408 trans-eQTLs. Association analysis resulted in three eQTLs associated with BFT and 24 with IMF. Functional enrichment analysis of genes regulated by these 27 eQTLs revealed noteworthy pathways that can play a fundamental role in lipid metabolism and fat deposition, such as immune response, cytoskeleton remodeling, iron transport, and phospholipid metabolism. We next used ATAC-Seq assay to identify and overlap eQTL and open chromatin regions. Six eQTLs were in regulatory regions, four in predicted insulators and possible CCCTC-binding factor DNA binding sites, one in an active enhancer region, and the last in a low signal region. Our results provided novel insights into the transcriptional regulation of IMF and BFT, unraveling putative regulatory variants.

Identification of novel genes associated with herbicide tolerance in Lentil (Lens culinaris ssp. culinaris Medik.).

Weeds pose a major constraint in lentil cultivation, leading to decrease farmers' revenues by reducing the yield and increasing the management costs. The development of herbicide tolerant cultivars is essential to increase lentil yield. Even though herbicide tolerant lines have been identified in lentils, breeding efforts are still limited and lack proper validation. Marker assisted selection (MAS) can increase selection accuracy at early generations. Total 292 lentil accessions were evaluated under different dosages of two herbicides, metribuzin and imazethapyr, during two seasons at Marchouch, Morocco and Terbol, Lebanon. Highly significant differences among accessions were observed for days to flowering (DF) and maturity (DM), plant height (PH), biological yield (BY), seed yield (SY), number of pods per plant (NP), as well as the reduction indices (RI) for PH, BY, SY and NP. A total of 10,271 SNPs markers uniformly distributed along the lentil genome were assayed using Multispecies Pulse SNP chip developed at Agriculture Victoria, Melbourne. Meta-GWAS analysis was used to detect marker-trait associations, which detected 125 SNPs markers associated with different traits and clustered in 85 unique quantitative trait loci. These findings provide valuable insights for initiating MAS programs aiming to enhance herbicide tolerance in lentil crop.

Characterizing stay-green in barley across diverse environments: unveiling novel haplotypes.

KEY MESSAGE: There is variation in stay-green within barley breeding germplasm, influenced by multiple haplotypes and environmental conditions. The positive genetic correlation between stay-green and yield across multiple environments highlights the potential as a future breeding target. Barley is considered one of the most naturally resilient crops making it an excellent candidate to dissect the genetics of drought adaptive component traits. Stay-green, is thought to contribute to drought adaptation, in which the photosynthetic machinery is maintained for a longer period post-anthesis increasing the photosynthetic duration of the plant. In other cereal crops, including wheat, stay-green has been linked to increased yield under water-limited conditions. Utilizing a panel of diverse barley breeding lines from a commercial breeding program we aimed to characterize stay-green in four environments across two years. Spatiotemporal modeling was used to accurately model senescence patterns from flowering to maturity characterizing the variation for stay-green in barley for the first time. Environmental effects were identified, and multi-environment trait analysis was performed for stay-green characteristics during grain filling. A consistently positive genetic correlation was found between yield and stay-green. Twenty-two chromosomal regions with large effect haplotypes were identified across and within environment types, with ten being identified in multiple environments. In silico stacking of multiple desirable haplotypes showed an opportunity to improve the stay-green phenotype through targeted breeding. This study is the first of its kind to model barley stay-green in a large breeding panel and has detected novel, stable and environment specific haplotypes. This provides a platform for breeders to develop Australian barley with custom senescence profiles for improved drought adaptation.

Disentangling genetic effects on transcriptional and post-transcriptional gene regulation through integrating exon and intron expression QTLs.

Expression quantitative trait loci (eQTL) studies typically consider exon expression of genes and discard intronic RNA sequencing reads despite their information on RNA metabolism. Here, we quantify genetic effects on exon and intron levels of genes and their ratio in lymphoblastoid cell lines, revealing thousands of cis-QTLs of each type. While genetic effects are often shared between cis-QTL types, 7814 (47%) are not detected as top cis-QTLs at exon levels. We show that exon levels preferentially capture genetic effects on transcriptional regulation, while exon-intron-ratios better detect those on co- and post-transcriptional processes. Considering all cis-QTL types substantially increases (by 71%) the number of colocalizing variants identified by genome-wide association studies (GWAS). It further allows dissecting the potential gene regulatory processes underlying GWAS associations, suggesting comparable contributions by transcriptional (50%) and co- and post-transcriptional regulation (46%) to complex traits. Overall, integrating intronic RNA sequencing reads in eQTL studies expands our understanding of genetic effects on gene regulatory processes.

Genetic influence on within-person longitudinal change in anthropometric traits in the UK Biobank.

The causes of temporal fluctuations in adult traits are poorly understood. Here, we investigate the genetic determinants of within-person trait variability of 8 repeatedly measured anthropometric traits in 50,117 individuals from the UK Biobank. We found that within-person (non-directional) variability had a SNP-based heritability of 2-5% for height, sitting height, body mass index (BMI) and weight (P ≤ 2.4 × 10-3). We also analysed longitudinal trait change and show a loss of both average height and weight beyond about 70 years of age. A variant tracking the Alzheimer's risk APOE- E 4 allele (rs429358) was significantly associated with weight loss ( ß = -0.047 kg per yr, s.e. 0.007, P = 2.2 × 10-11), and using 2-sample Mendelian Randomisation we detected a relationship consistent with causality between decreased lumbar spine bone mineral density and height loss (bxy = 0.011, s.e. 0.003, P = 3.5 × 10-4). Finally, population-level variance quantitative trait loci (vQTL) were consistent with within-person variability for several traits, indicating an overlap between trait variability assessed at the population or individual level. Our findings help elucidate the genetic influence on trait-change within an individual and highlight disease risks associated with these changes.

Identification and characterization of whole blood gene expression and splicing quantitative trait loci during early to mid-lactation of dairy cattle.

BACKGROUND: Characterization of regulatory variants (e.g., gene expression quantitative trait loci, eQTL; gene splicing QTL, sQTL) is crucial for biologically interpreting molecular mechanisms underlying loci associated with complex traits. However, regulatory variants in dairy cattle, particularly in specific biological contexts (e.g., distinct lactation stages), remain largely unknown. In this study, we explored regulatory variants in whole blood samples collected during early to mid-lactation (22-150 days after calving) of 101 Holstein cows and analyzed them to decipher the regulatory mechanisms underlying complex traits in dairy cattle. RESULTS: We identified 14,303 genes and 227,705 intron clusters expressed in the white blood cells of 101 cattle. The average heritability of gene expression and intron excision ratio explained by cis-SNPs is 0.28 ± 0.13 and 0.25 ± 0.13, respectively. We identified 23,485 SNP-gene expression pairs and 18,166 SNP-intron cluster pairs in dairy cattle during early to mid-lactation. Compared with the 2,380,457 cis-eQTLs reported to be present in blood in the Cattle Genotype-Tissue Expression atlas (CattleGTEx), only 6,114 cis-eQTLs (P < 0.05) were detected in the present study. By conducting colocalization analysis between cis-e/sQTL and the results of genome-wide association studies (GWAS) from four traits, we identified a cis-e/sQTL (rs109421300) of the DGAT1 gene that might be a key marker in early to mid-lactation for milk yield, fat yield, protein yield, and somatic cell score (PP4 > 0.6). Finally, transcriptome-wide association studies (TWAS) revealed certain genes (e.g., FAM83H and TBC1D17) whose expression in white blood cells was significantly (P < 0.05) associated with complex traits. CONCLUSIONS: This study investigated the genetic regulation of gene expression and alternative splicing in dairy cows during early to mid-lactation and provided new insights into the regulatory mechanisms underlying complex traits of economic importance.

Removing unwanted variation between samples in Hi-C experiments.

Hi-C data are commonly normalized using single sample processing methods, with focus on comparisons between regions within a given contact map. Here, we aim to compare contact maps across different samples. We demonstrate that unwanted variation, of likely technical origin, is present in Hi-C data with replicates from different individuals, and that properties of this unwanted variation change across the contact map. We present band-wise normalization and batch correction, a method for normalization and batch correction of Hi-C data and show that it substantially improves comparisons across samples, including in a quantitative trait loci analysis as well as differential enrichment across cell types.

A comprehensive genetic map of cytokine responses in Lyme borreliosis.

The incidence of Lyme borreliosis has risen, accompanied by persistent symptoms. The innate immune system and related cytokines are crucial in the host response and symptom development. We characterized cytokine production capacity before and after antibiotic treatment in 1,060 Lyme borreliosis patients. We observed a negative correlation between antibody production and IL-10 responses, as well as increased IL-1Ra responses in patients with disseminated disease. Genome-wide mapping the cytokine production allowed us to identify 34 cytokine quantitative trait loci (cQTLs), with 31 novel ones. We pinpointed the causal variant at the TLR1-6-10 locus and validated the regulation of IL-1Ra responses at transcritpome level using an independent cohort. We found that cQTLs contribute to Lyme borreliosis susceptibility and are relevant to other immune-mediated diseases. Our findings improve the understanding of cytokine responses in Lyme borreliosis and provide a genetic map of immune function as an expanded resource.

Transcriptome-wide association mapping provides insights into the genetic basis and candidate genes governing flowering, maturity and seed weight in rice bean (Vigna umbellata).

BACKGROUND: Rice bean (Vigna umbellata), an underrated legume, adapts to diverse climatic conditions with the potential to support food and nutritional security worldwide. It is used as a vegetable, minor food crop and a fodder crop, being a rich source of proteins, minerals, and essential fatty acids. However, little effort has been made to decipher the genetic and molecular basis of various useful traits in this crop. Therefore, we considered three economically important traits i.e., flowering, maturity and seed weight of rice bean and identified the associated candidate genes employing an associative transcriptomics approach on 100 diverse genotypes out of 1800 evaluated rice bean accessions from the Indian National Genebank. RESULTS: The transcriptomics-based genotyping of one-hundred diverse rice bean cultivars followed by pre-processing of genotypic data resulted in 49,271 filtered markers. The STRUCTURE, PCA and Neighbor-Joining clustering of 100 genotypes revealed three putative sub-populations. The marker-trait association analysis involving various genome-wide association study (GWAS) models revealed significant association of 82 markers on 48 transcripts for flowering, 26 markers on 22 transcripts for maturity and 22 markers on 21 transcripts for seed weight. The transcript annotation provided information on the putative candidate genes for the considered traits. The candidate genes identified for flowering include HSC80, P-II PsbX, phospholipid-transporting-ATPase-9, pectin-acetylesterase-8 and E3-ubiquitin-protein-ligase-RHG1A. Further, the WRKY1 and DEAD-box-RH27 were found to be associated with seed weight. Furthermore, the associations of PIF3 and pentatricopeptide-repeat-containing-gene with maturity and seed weight, and aldo-keto-reductase with flowering and maturity were revealed. CONCLUSION: This study offers insights into the genetic basis of key agronomic traits in rice bean, including flowering, maturity, and seed weight. The identified markers and associated candidate genes provide valuable resources for future exploration and targeted breeding, aiming to enhance the agronomic performance of rice bean cultivars. Notably, this research represents the first transcriptome-wide association study in pulse crop, uncovering the candidate genes for agronomically useful traits.

Investigating the genetic control of plant development in spring barley under speed breeding conditions.

KEY MESSAGE: This study found that the genes, PPD-H1 and ELF3, control the acceleration of plant development under speed breeding, with important implications for optimizing the delivery of climate-resilient crops. Speed breeding is a tool to accelerate breeding and research programmes. Despite its success and growing popularity with breeders, the genetic basis of plant development under speed breeding remains unknown. This study explored the developmental advancements of barley genotypes under different photoperiod regimes. A subset of the HEB-25 Nested Association Mapping population was evaluated for days to heading and maturity under two contrasting photoperiod conditions: (1) Speed breeding (SB) consisting of 22 h of light and 2 h of darkness, and (2) normal breeding (NB) consisting of 16 h of light and 8 h of darkness. GWAS revealed that developmental responses under both conditions were largely controlled by two loci: PPDH-1 and ELF3. Allelic variants at these genes determine whether plants display early flowering and maturity under both conditions. At key QTL regions, domesticated alleles were associated with late flowering and maturity in NB and early flowering and maturity in SB, whereas wild alleles were associated with early flowering under both conditions. We hypothesize that this is related to the dark-dependent repression of PPD-H1 by ELF3 which might be more prominent in NB conditions. Furthermore, by comparing development under two photoperiod regimes, we derived an estimate of plasticity for the two traits. Interestingly, plasticity in development was largely attributed to allelic variation at ELF3. Our results have important implications for our understanding and optimization of speed breeding protocols particularly for introgression breeding and the design of breeding programmes to support the delivery of climate-resilient crops.

Unique genetic architecture of prolificacy in 'Sikkim Primitive' maize unraveled through whole-genome resequencing-based DNA polymorphism.

KEY MESSAGE: 'Sikkim Primitive' maize landrace, unique for prolificacy (7-9 ears per plant) possesses unique genomic architecture in branching and inflorescence-related gene(s), and locus Zm00001eb365210 encoding glycosyltransferases was identified as the putative candidate gene underlying QTL (qProl-SP-8.05) for prolificacy. The genotype possesses immense usage in breeding high-yielding baby-corn genotypes. 'Sikkim Primitive' is a native landrace of North Eastern Himalayas, and is characterized by having 7-9 ears per plant compared to 1-2 ears in normal maize. Though 'Sikkim Primitive' was identified in the 1960s, it has not been characterized at a whole-genome scale. Here, we sequenced the entire genome of an inbred (MGUSP101) derived from 'Sikkim Primitive' along with three non-prolific (HKI1128, UMI1200, and HKI1105) and three prolific (CM150Q, CM151Q and HKI323) inbreds. A total of 942,417 SNPs, 24,160 insertions, and 27,600 deletions were identified in 'Sikkim Primitive'. The gene-specific functional mutations in 'Sikkim Primitive' were classified as 10,847 missense (54.36%), 402 non-sense (2.015%), and 8,705 silent (43.625%) mutations. The number of transitions and transversions specific to 'Sikkim Primitive' were 666,021 and 279,950, respectively. Among all base changes, (G to A) was the most frequent (215,772), while (C to G) was the rarest (22,520). Polygalacturonate 4-α-galacturonosyltransferase enzyme involved in pectin biosynthesis, cell-wall organization, nucleotide sugar, and amino-sugar metabolism was found to have unique alleles in 'Sikkim Primitive'. The analysis further revealed the Zm00001eb365210 gene encoding glycosyltransferases as the putative candidate underlying QTL (qProl-SP-8.05) for prolificacy in 'Sikkim Primitive'. High-impact nucleotide variations were found in ramosa3 (Zm00001eb327910) and zeaxanthin epoxidase1 (Zm00001eb081460) genes having a role in branching and inflorescence development in 'Sikkim Primitive'. The information generated unraveled the genetic architecture and identified key genes/alleles unique to the 'Sikkim Primitive' genome. This is the first report of whole-genome characterization of the 'Sikkim Primitive' landrace unique for its high prolificacy.

African ancestry-derived APOL1 risk genotypes show proximal epigenetic associations.

Apolipoprotein L1 (APOL1) coding variants, termed G1 and G2, are established genetic risk factors for a growing spectrum of diseases, including kidney disease, in individuals of African ancestry. Evidence suggests that the risk variants, which show a recessive mode of inheritance, lead to toxic gain-of-function changes of the APOL1 protein. Disease occurrence and presentation vary, likely due to modifiers or second hits. To understand the role of the epigenetic landscape in relation to APOL1 risk variants, we performed methylation quantitative trait locus (meQTL) analysis to identify differentially methylated CpGs influenced by APOL1 risk variants in 611 African American individuals. We identified five CpGs that were significantly associated with APOL1 risk alleles in discovery and replication studies, and one CpG-APOL1 association was independent of other genomic variants. Our study highlights proximal DNA methylation alterations that may help explain the variable disease risk and clinical manifestation of APOL1 variants.

Fine mapping of major QTL qshgd1 for spontaneous haploid genome doubling in maize (Zea mays L.).

KEY MESSAGE: A large-effect QTL was fine mapped, which revealed 79 gene models, with 10 promising candidate genes, along with a novel inversion. In commercial maize breeding, doubled haploid (DH) technology is arguably the most efficient resource for rapidly developing novel, completely homozygous lines. However, the DH strategy, using in vivo haploid induction, currently requires the use of mutagenic agents which can be not only hazardous, but laborious. This study focuses on an alternative approach to develop DH lines-spontaneous haploid genome duplication (SHGD) via naturally restored haploid male fertility (HMF). Inbred lines A427 and Wf9, the former with high HMF and the latter with low HMF, were selected to fine-map a large-effect QTL associated with SHGD-qshgd1. SHGD alleles were derived from A427, with novel haploid recombinant groups having varying levels of the A427 chromosomal region recovered. The chromosomal region of interest is composed of 45 megabases (Mb) of genetic information on chromosome 5. Significant differences between haploid recombinant groups for HMF were identified, signaling the possibility of mapping the QTL more closely. Due to suppression of recombination from the proximity of the centromere, and a newly discovered inversion region, the associated QTL was only confined to a 25 Mb region, within which only a single recombinant was observed among ca. 9,000 BC1 individuals. Nevertheless, 79 gene models were identified within this 25 Mb region. Additionally, 10 promising candidate genes, based on RNA-seq data, are described for future evaluation, while the narrowed down genome region is accessible for straightforward introgression into elite germplasm by BC methods.