The research team, led by Ruikang Liu, Chiyun Sun and Yonghao Li of Beijing University of Chinese Medicine and the China Academy of Chinese Medical Sciences, together with colleagues at several Chinese hospitals, assembled genome-wide association study (GWAS) summary statistics for heart failure and osteoporosis and subjected them to a battery of complementary statistical tools. Rather than asking whether the two diseases share one or two risk genes, the investigators probed the architecture of the entire genetic landscape: how many variants influence both traits, whether those variants push risk in the same or opposite directions, which genes the signals point to, and where in the body those genes are active.

The first and perhaps most striking finding concerns polygenic overlap. Using MiXeR, a bivariate causal mixture model that estimates the number of genetic variants with non-null effects on each trait, the team calculated that approximately 447 ± 114 variants influence both heart failure and osteoporosis. That shared set corresponds to roughly 52 percent of the estimated non-null variant count for osteoporosis and about 21 percent of that for heart failure. In other words, a substantial fraction of the genetic determinants of bone fragility are also determinants of cardiac dysfunction — a degree of overlap that far exceeds what most single-gene approaches would ever reveal.

Yet here the story takes an unexpected turn. Despite this extensive overlap, the genome-wide genetic correlation between the two diseases was weakly negative, at rg = −0.091 (standard error 0.024, P = 1.07 × 10⁻⁴). Moreover, when the researchers examined the correlation of effect sizes within the shared non-null component, they again found a negative value, rho_beta = −0.295. The interpretation is counterintuitive but biologically plausible: many of the same variants are involved in both conditions, but they frequently act in opposite directions. A variant that raises the risk of heart failure may, on average, lower the risk of osteoporosis, and vice versa. The diseases share their genetic machinery even as that machinery is tuned in opposing ways.

To move from global statistics to specific genomic locations, the team applied conditional and conjunctional false discovery rate (cond/conjFDR) analysis, a framework that borrows statistical strength across related traits to detect associations that neither dataset could confirm alone. This approach identified five pleiotropic loci — genomic regions harboring variants that influence both diseases — and, notably, three novel loci associated with heart failure that had not previously been flagged. The method’s power lies in its ability to elevate weak, sub-threshold signals to significance when they are corroborated across traits, effectively squeezing more discovery out of existing GWAS resources.

The next layer of analysis asked which genes, rather than which markers, the associations implicate. Through transcriptome-wide association studies (TWAS) conducted with the MetaXcan framework, the researchers imputed the genetic regulation of gene expression across tissues and tested whether genetically predicted expression levels associate with disease risk. They then cross-referenced these results with Bayesian colocalization methods — eCAVIAR and fastenloc — which test whether a GWAS signal and an expression quantitative trait locus (eQTL) signal in the same region are driven by the same causal variant rather than by two distinct variants in linkage disequilibrium. The consistent winner of this triangulation was HOXC6, a homeobox transcription factor gene prioritized as a shared transcriptomic candidate across multiple tissues. HOX genes are master regulators of developmental patterning, and their involvement in both a cardiovascular and a skeletal disorder hints at developmental programs that echo into adult disease susceptibility.

Perhaps the most novel component of the study is its spatial dimension. Using gsMap — genetically informed spatial mapping of cells for complex traits — the team projected human disease-associated genetic signals onto a reference atlas of single-cell spatial transcriptomics from E16.5 mouse embryos. This technique asks not just which genes are implicated, but in which anatomical niches their expression is concentrated. The results were strikingly compartmentalized: heart failure-associated genetic signals were enriched in embryonic muscle regions, consistent with the cardiac muscle as the affected tissue, while osteoporosis-associated signals concentrated in cartilage and cartilage-primordium regions, the developmental precursors of the skeletal system. Within this spatial framework, GRB10 emerged as a shared candidate gene. GRB10 is an adaptor protein in insulin-like growth factor signaling, a pathway with well-documented roles in growth, metabolism and tissue size regulation — making its dual appearance in heart and bone contexts a compelling lead for future work.

The authors are careful to frame these results as hypothesis-generating rather than definitive. The analyses rely on summary statistics and reference atlases, and the spatial mapping uses an embryonic mouse reference to interpret human adult disease, an approach that assumes developmental expression patterns illuminate disease-relevant biology. The negative genetic correlation also complicates any simple clinical translation: shared variants acting in opposite directions mean that a therapeutic strategy targeting a shared gene could, in principle, benefit one organ while harming the other. Precisely for that reason, the identification of candidates such as HOXC6 and GRB10 matters — they define concrete molecular entry points whose directionality and tissue context can now be tested experimentally in cellular and animal models.

What the study delivers is a new conceptual map of the heart–bone axis. It quantifies, for the first time at this resolution, just how extensively the genetic architectures of heart failure and osteoporosis interpenetrate, while simultaneously demonstrating that overlap does not equal concordance. The five pleiotropic loci, the three novel heart failure loci, the multi-tissue prioritization of HOXC6 and the spatially resolved nomination of GRB10 together generate a set of testable predictions about regulatory links between cardiovascular and skeletal traits. As the population ages and the clinical burden of both diseases grows, understanding the shared — and antagonistic — genetic wiring beneath them may prove essential for designing interventions that strengthen the heart without weakening the skeleton, and vice versa.

Subject of Research: Shared genetic architecture and spatial transcriptomic correlations between heart failure and osteoporosis

Article Title: Integrating multi-omics analysis reveals shared genetic architecture and spatial correlations between heart failure and osteoporosis

Article References: Integrating multi-omics analysis reveals shared genetic architecture and spatial correlations between heart failure and osteoporosis. (n.d.). https://doi.org/10.1038/s41598-026-74146-7

Image Credits: AI Generated

DOI: 10.1038/s41598-026-74146-7

Keywords: heart failure, osteoporosis, GWAS, polygenic overlap, MiXeR, TWAS, colocalization, HOXC6, GRB10, spatial transcriptomics, genetic correlation, pleiotropy

The team approached the problem with two complementary technologies drawn from statistical genetics. The first was local genetic correlation analysis, carried out with a method called SUPERGNOVA, which estimates how strongly the genetic influences on two traits align or diverge within specific segments of the genome rather than averaging across all roughly three billion base pairs. The researchers applied this tool to sixteen specific cancers plus a pan-cancer phenotype, computing pairwise local correlations for every combination. The second approach relied on polygenic risk scores, or PRSs, which condense the small effects of hundreds of thousands of common genetic variants into a single number estimating an individual’s inherited liability to a disease. By testing whether a PRS built for one cancer predicted risk of another, the researchers could ask whether shared susceptibility extends from broad genomic architecture down to the level of actual prediction.

The genome-wide results set the stage for everything that followed. When the authors computed standard genome-wide genetic correlations across all cancer pairs, they found twenty significantly correlated pairs, a number that confirms earlier reports of modest pan-cancer sharing but also underscores how much of the genome’s influence on each cancer remains disease-specific. Genetic correlation at this scale is a blunt instrument: it asks whether the constellation of variants affecting one cancer tends, on average, to affect another in the same direction. A weak overall correlation does not mean the two cancers share nothing; it may simply mean that shared effects concentrated in a handful of genomic neighborhoods are diluted by thousands of variants acting independently.

That is precisely what the local analysis revealed. Zooming into genomic regions rather than whole chromosomes, Hu and colleagues identified eighty-two distinct regions where genetic signals were shared across sixty-six different cancer pairs. In other words, the number of cancer pairs showing evidence of shared genetic influence more than tripled once the analysis moved from the genome-wide to the regional scale. These regions are not random stretches of DNA. The team annotated them using genetic correlations with non-cancer traits and known associations catalogued in the GWAS Catalog, building a functional picture of where and why cancer susceptibilities converge.

Among the eighty-two shared regions, five stood out as hubs of what the authors call mutually correlated cancer clusters. In these locations, the local genetic effects on multiple cancers were correlated with one another in a structured web, rather than in isolated pairs. When the researchers examined what these hub regions do biologically, enrichment appeared in functional domains connected to inflammatory processes. The finding fits a broader theme in cancer genetics: chronic inflammation is a recognized driver of tumor initiation and progression in many tissues, and inherited variation in immune and inflammatory regulation could plausibly tilt risk for several cancers at once. A single regulatory variant that dampens or amplifies inflammatory signaling, for example, might leave detectable fingerprints in the genetic architecture of lung, colorectal, and other inflammation-associated malignancies simultaneously.

The polygenic risk score analysis added an independent line of evidence. Testing PRSs across cancer pairs, the researchers identified five pairs of cancers with statistically significant cross-cancer associations, meaning that a genetic risk score optimized for one cancer carried predictive information for the other. This is a meaningful step beyond correlation of summary statistics, because PRS association tests whether shared genetic liability is strong enough to translate into measurable prediction in independent data. Still, five significant pairs out of the many tested reinforces the study’s headline conclusion: cross-cancer prediction from genetics exists, but it is the exception rather than the rule, and it demands careful dissection before it can be interpreted biologically.

That dissection came from the study’s most technically inventive component, the decomposition of genome-wide PRSs into functionally informed components. Rather than treating a PRS as an indivisible number, the authors split it into pleiotropy-informed parts, which emphasize variants already known to influence multiple traits, and pathway-specific parts, which weight variants according to membership in defined biological pathways. When the significant cross-cancer PRS associations were decomposed this way, a clear pattern emerged: the shared predictive signal was concentrated in specific pleiotropy groups and in immune-related pathways, not spread diffusely across the genome. The cross-cancer associations, in other words, are not a generic reflection of genetic overlap between complex traits; they trace back to identifiable biological machinery.

The immune pathway enrichment ties the two analytical strands of the study together elegantly. The regional analysis pointed to inflammatory functional domains in the hub regions of shared local genetic correlation, and the PRS decomposition pointed to immune-related pathways as carriers of cross-cancer predictive signal. Two independent methods, operating at different scales of genomic resolution, converged on the same biological story. This kind of methodological triangulation matters in a field where pleiotropy, the phenomenon of one variant influencing many traits, can arise from genuine shared biology or from statistical artifacts such as sample overlap and population stratification. Agreement between local correlation and PRS-based prediction makes a genuine biological basis for the shared risk considerably more convincing.

The implications reach into both research practice and, eventually, clinical risk assessment. For researchers, the study supplies a framework: rather than asking only whether two cancers are genetically correlated overall, investigators can map where in the genome the correlation lives, which functional categories those regions implicate, and whether the shared signal is strong enough to survive PRS-based prediction testing. For clinicians and genetic counselors, the work tempers expectations while pointing to future opportunities. Current polygenic risk scores are largely cancer-specific, and the modest number of significant cross-cancer pairs suggests that a single pan-cancer score is not on the immediate horizon. But for the specific pairs identified, particularly those anchored in immune pathways, multi-cancer risk models that borrow strength across related malignancies could eventually sharpen risk stratification, especially as PRS methods improve and datasets grow more diverse.

The study also carries a cautionary note about scale. Findings that vanish at the genome-wide level but appear at the regional level illustrate a general principle in complex-trait genetics: aggregation can obscure as much as it reveals. Eighty-two shared regions hiding inside twenty significant genome-wide correlations is a quantitative demonstration that the architecture of shared disease risk is heterogeneous, concentrated in pockets of the genome where pleiotropic variants cluster. As biobank-scale datasets expand and methods like SUPERGNOVA and PRS decomposition mature, the expectation is that more of these hidden threads connecting cancers will come into view, each one a potential clue to mechanisms that tumor types as different as melanoma and leukemia might, unexpectedly, share. For now, Hu and colleagues have shown that the genetic conversation between cancers is quieter than a pan-cancer model would hope, but far more organized than a genome-wide average would suggest.

Subject of Research: Shared polygenic genetic risk across multiple cancer types

Article Title: Identifying shared polygenic risk across cancers

Article References: Hu, J., Muheyati, M., Xu, L., DeWan, A., & Zhao, H. (2026). Identifying shared polygenic risk across cancers. PLOS Genetics, 22(9), e1012308. https://doi.org/10.1371/journal.pgen.1012308

Image Credits: AI Generated

DOI: 10.1371/journal.pgen.1012308

Keywords: polygenic risk scores, genetic correlation, cancer genetics, pleiotropy, SUPERGNOVA, immune pathways, inflammation, GWAS, cross-cancer susceptibility, genomic regions, PLOS Genetics, risk prediction

The research team, led by Katrine E. Dailey and Nicholas C. Wu, set out to map the evolutionary landscape that shaped CR9114. This antibody is encoded by the IGHV1-69 germline gene, a segment of the human antibody repertoire that is particularly prone to producing stem-directed influenza antibodies. CR9114 acquired its remarkable breadth not at birth but through somatic hypermutation, the process by which B cells randomly mutate their antibody genes after encountering an antigen, followed by selection for variants that bind more tightly. The investigators asked a deceptively simple question: if breadth is so valuable, why has the immune system so rarely stumbled upon it?

To answer this, the researchers performed four separate deep mutational scanning experiments, a high-throughput technique that allows scientists to test tens of thousands of antibody variants simultaneously. They compared the binding affinity landscapes of both the germline precursor and the mature, somatically mutated CR9114 against three distinct hemagglutinin targets: the H1 subtype of influenza A, the H3 subtype of influenza A, and the hemagglutinin of influenza B virus. By systematically mutating every position in the antibody’s binding regions and measuring how each change affected binding to each of the three targets, they generated an unprecedentedly detailed picture of the trade-offs that govern antibody evolution.

The central finding is stark: many mutations that have little effect on binding to H1 hemagglutinin, or that even improve it, are actively detrimental to binding the H3 hemagglutinin or the influenza B hemagglutinin. In other words, the three binding tasks are not compatible at the level of individual amino acid positions. An antibody mutation that helps the immune system track one viral subtype can silently sabotage its ability to recognize another. This creates a rugged fitness landscape in which the path toward broader recognition is littered with valleys of reduced function, valleys that natural selection within a single infected or vaccinated individual has little incentive to cross.

This incompatibility between binding targets is compounded by two well-known but often underappreciated genetic phenomena: epistasis and pleiotropy. Epistasis refers to situations in which the effect of one mutation depends on the presence of other mutations elsewhere in the protein. In the context of IGHV1-69-derived stem antibodies, the study revealed that epistasis is pervasive. A mutation that appears neutral in the germline background may become beneficial or harmful once other somatic mutations have accumulated, meaning that the evolutionary consequences of any single change cannot be predicted in isolation. Pleiotropy, meanwhile, describes the situation in which a single mutation affects multiple traits simultaneously, here the binding affinities to different hemagglutinin types, forcing unavoidable trade-offs.

The practical consequence of these constraints is that the evolutionary trajectory available to a developing B cell is far narrower than the raw diversity of mutations might suggest. Somatic hypermutation generates plenty of variants, but selection can only act on those that maintain or improve binding to the antigen currently driving the immune response. If improving binding to H1 hemagglutinin simultaneously degrades binding to H3 and influenza B hemagglutinins, then a B cell expanding during an H1 infection will be steered away from breadth even if breadth would be valuable in the long run. The immune system, in effect, optimizes locally rather than globally, and the local optima it reaches are frequently narrow.

CR9114, then, is best understood not as a representative example of what stem antibodies can achieve but as an evolutionary outlier, a rare combination of mutations that happened to satisfy the conflicting requirements of binding across both influenza A and B viruses. The deep mutational scanning data show that the mature antibody sits at a point in sequence space that is difficult to reach step by step, because many of the intermediate forms along plausible paths would have lost binding to one or more targets and been eliminated by selection. The rarity of CR9114-like antibodies is therefore not a failure of the immune repertoire to generate diversity but a predictable outcome of the geometry of the binding landscape.

These findings carry significant implications for vaccine design. Many current efforts to build a universal influenza vaccine aim to drive the immune system toward stem-directed antibodies by using engineered immunogens such as chimeric hemagglutinins, in which the head domain is repeatedly swapped to redirect the response toward the conserved stem. The new results suggest that such strategies must contend with the epistatic and pleiotropic constraints documented here. Sequential vaccination with different hemagglutinin subtypes may help, because it can apply selection pressure from multiple targets at once, but the study indicates that the order and combination of immunogens matter enormously, since mutations favored by one antigen may undermine recognition of the next.

The work also provides a mechanistic framework for evaluating antibody repertoires more broadly. If epistasis is prevalent in IGHV1-69-derived stem antibodies, then computational models that treat each mutation independently will systematically mispredict which antibody variants are viable. Deep mutational scanning of this kind, applied across multiple antibody families and viral targets, could help vaccine designers identify which evolutionary paths are accessible and which immunization regimens are most likely to shepherd B cells toward genuinely broad neutralization rather than toward narrow local optima.

Ultimately, the study reframes the quest for a universal influenza vaccine as a problem of navigating a constrained evolutionary landscape rather than simply eliciting more antibodies. The scarcity of antibodies like CR9114 is not an accident but a consequence of the fundamental genetics of antibody maturation, in which the same mutations that confer one binding capability often compromise another, and in which the value of any mutation depends on the company it keeps. Understanding these constraints in molecular detail brings vaccine developers one step closer to designing immunization strategies that can work with, rather than against, the evolutionary rules that govern the human immune response to influenza.

Subject of Research: Evolutionary constraints on broadly neutralizing influenza stem antibodies

Article Title: Epistasis and pleiotropy constrain the evolution of cross-reactive breadth in a broadly neutralizing influenza antibody

Article References: Dailey, K. E., Wang, Y., Teo, Q. W., Wang, C., Lv, H., Lei, R., Tong, M., Rodriguez, L. A., Wang, L., & Wu, N. C. (2026). Epistasis and pleiotropy constrain the evolution of cross-reactive breadth in a broadly neutralizing influenza antibody. PLOS Biology, 24(9), e3004021. https://doi.org/10.1371/journal.pbio.3004021

Image Credits: AI Generated

DOI: 10.1371/journal.pbio.3004021

Keywords: influenza, broadly neutralizing antibodies, CR9114, hemagglutinin stem, epistasis, pleiotropy, deep mutational scanning, IGHV1-69, somatic hypermutation, universal vaccine, PLOS Biology, antibody evolution

The story began to unfold in recent years when large-scale genetic studies pointed an unexpected finger at ITSN1. One analysis combining the deCODE, UK Biobank, and Accelerating Medicines Partnership Parkinson’s Disease datasets found that heterozygous loss-of-function variants in the gene, meaning changes that cripple one of the two inherited copies, confer a substantially elevated risk of Parkinson’s disease. A second study described five patients from three families carrying distinct frameshift variants in ITSN1, documenting their clinical pictures in detail. A third effort, drawing on the UK Biobank, the AMP-PD cohort, and the All of Us program, reported a similar increase in risk and even supplied functional evidence from a fruit fly model, showing that the protein encoded by ITSN1, called Intersectin-1, interacts with alpha-synuclein, the sticky protein that clumps in the brains of Parkinson’s patients. What remained unclear was whether patients carrying these variants show a recognizable clinical signature.

To answer that question, the ROPAD team turned to their own trove of data. They analyzed 8,660 reportedly unrelated Parkinson’s disease patients recruited between April 2019 and August 2024, drawn from Europe (68 percent), North America (16 percent), Israel (10 percent), Turkey (4 percent), and South America (2 percent). Crucially, all of these patients had undergone exome or genome sequencing and had no relevant variants in known Parkinson’s genes. That exclusion criterion matters: the researchers removed anyone carrying a pathogenic or likely pathogenic variant in dominant Parkinson’s genes such as LRRK2, SNCA, VPS35, RAB32, or CHCHD2, anyone with two heterozygous or one homozygous pathogenic variant in recessive genes including PARK7, PINK1, and PRKN, and anyone with a pathogenic or risk-factor variant in GBA1, the most common genetic contributor to the disease. By stripping out known genetic explanations, the team could ask whether ITSN1 variants stand on their own as a risk factor.

As a comparison group, the researchers examined 210,657 exome- or genome-sequenced individuals from the CENTOGENE Biodatabank, people referred for routine diagnostic testing, healthy parents in trio analyses, or participants in other observational studies, none of whom had been referred for or diagnosed with Parkinson’s disease. Loss-of-function variants were defined rigorously as start-loss, stop-gain, stop-loss, splice-donor, splice-acceptor, or frameshift changes, identified through a validated in-house bioinformatic pipeline. The team then mined electronic case report forms for clinical details, including year of birth, sex, age at disease onset, age at diagnosis, age at enrollment, family history, and motor scores from a partial Unified Parkinson’s Disease Rating Scale examination.

The result was striking. Thirteen of the 8,660 patients, or 0.15 percent, carried heterozygous loss-of-function variants in ITSN1, compared with just 0.025 percent of the 210,657 comparison individuals. That difference translates into an odds ratio of 6.09, with a 95 percent confidence interval of 3.31 to 11.19 and a p-value of 1.34 times ten to the minus six, a level of statistical significance that is difficult to dismiss. The allele-level comparison tells the same story: the proportion of ITSN1 loss-of-function alleles in the Parkinson’s patients, 0.075 percent, is roughly six times higher than the proportion found in non-Finnish European individuals in the Genome Aggregation Database, the reference catalog that population geneticists use as a baseline for how common variants should be. The enrichment held against every other population in that database as well.

Perhaps most convincing is how neatly these numbers overlap with the earlier studies. The confidence interval from the ROPAD analysis largely overlaps those reported in the two prior association studies, which found odds ratios of 7.3 and 10.5 respectively. The carrier frequency of 0.15 percent in ROPAD matches the 0.15 percent reported in the deCODE-led meta-analysis and the 0.14 percent in the family study. Within the ROPAD cohort itself, ITSN1 loss-of-function variants rank as the sixth most common disease-relevant genetic finding, trailing only variants in GBA1, LRRK2, PRKN, SNCA, and GCH1. For a gene that was essentially invisible to Parkinson’s researchers a few years ago, that is a rapid ascent.

Digging into the variants themselves, the team identified 12 distinct ITSN1 loss-of-function changes among the 13 patients: four nonsense mutations, six frameshift variants, and two splice site changes. Eleven of the twelve appeared only once, and one, a small deletion designated c.2842_2843del, was found in two patients. When the researchers tested whether those two German patients were related, analysis of extended haplotypes found no evidence of close relatedness or a shared founder chromosome, suggesting the variant arose independently in each case. This pattern of largely private, one-off mutations contrasts sharply with the genetics of other dominant Parkinson’s genes, where a handful of recurrent variants, such as the famous LRRK2 G2019S change, account for most cases. It also hints that the mutational spectrum of ITSN1 resembles that of genes where virtually any disruption of the protein is harmful.

There is an intriguing twist. More than 40 distinct ITSN1 loss-of-function variants have already been reported in individuals with autism and neurodevelopmental disorders, and four of the variants found in the ROPAD patients had appeared in that literature. This overlap raises the possibility of pleiotropy, the phenomenon in which a single gene produces different outcomes depending on context. Intersectin-1 is a scaffolding protein involved in clathrin-mediated endocytosis, synaptic vesicle recycling, and actin and CDC42-related signaling, processes essential both for building neural circuits during development and for maintaining them across a lifetime. Losing one functional copy, a state called haploinsufficiency, might disturb neurodevelopmental circuitry in childhood or undermine neuronal maintenance later in life. Which outcome prevails may depend on genetic modifiers, the residual expression produced by splice-altering variants, environmental stressors, and other factors scientists have yet to untangle.

On the clinical side, the ROPAD analysis found that patients carrying ITSN1 variants look, in most respects, like ordinary Parkinson’s patients. The male-to-female ratio, the proportion with a positive family history, the age at diagnosis, the age at enrollment, and the disease duration all failed to differ significantly between carriers and non-carriers. A multiple linear regression of motor scores adjusted for disease duration also found no significant difference, though the p-value of 0.075 leaves open the possibility of a subtle effect. Interestingly, the two earlier large-cohort studies did report significantly younger ages at onset among carriers, with mean or median onsets in the early sixties versus around seventy in non-carriers. The ROPAD cohort as a whole skews young, with a mean and median onset of 55 years, which likely explains why the same signal did not emerge here. Notably, the median onset of 53 years among ROPAD carriers is comparable to the 47-year average onset in the five family-study patients, and one ROPAD participant carried both an ITSN1 deletion and a mild GBA1 variant, a combination the researchers excluded from their main analysis because the GBA1 change may have modified her disease risk, echoing what has been documented in people carrying both LRRK2 and GBA1 variants.

Taken together, the findings cement ITSN1 as a genuine Parkinson’s disease gene, one that produces a form of the illness clinically indistinguishable from idiopathic disease. The authors propose the designation PARK-ITSN1 for this genetic entity and call for further research into its natural history and treatment responses, questions that will matter as genetic testing becomes routine and as pharmaceutical companies design trials stratified by genotype. The study also carries a broader lesson about the power of deeply curated clinical and genetic datasets: when a diagnostic laboratory pairs sequencing with systematic phenotyping across thousands of patients, even vanishingly rare risk factors, present in fewer than two people per thousand, can be detected, validated, and connected to the biology of one of the world’s most common neurodegenerative diseases.

Subject of Research: Association between ITSN1 loss-of-function variants and Parkinson's disease risk in the ROPAD cohort

Article Title: Data from the ROPAD study corroborate an association between ITSN1 loss-of-function variants and Parkinson’s disease

Article References: Ganoza, C. A., Westenberger, A., Paul, J. J., Curado, F., Somerville, E. N., Rennecke, J., Bauer, P., & Beetz, C. (2026). Data from the ROPAD study corroborate an association between ITSN1 loss-of-function variants and Parkinson’s disease. npj Parkinson's Disease, 12(1), Article 237. https://doi.org/10.1038/s41531-026-01562-x

Image Credits: AI Generated

DOI: 10.1038/s41531-026-01562-x

Keywords: Parkinson's disease, ITSN1, genetics, loss-of-function variants, ROPAD study, haploinsufficiency, Intersectin-1, neurodegeneration, exome sequencing, age at onset, GBA1, pleiotropy

The foundation of the work is a large and diverse genetic resource: 353 core germplasm accessions of cultivated peanut, Arachis hypogaea L., sampled from around the world. Each accession was genotyped with high-quality single nucleotide polymorphism data, the tiny DNA letter differences that distinguish one variety from another. From these variants the team constructed 51,669 linkage disequilibrium blocks, stretches of neighboring SNPs that are inherited together as haplotypes. Working with haplotype blocks rather than individual SNPs matters because the functional unit of inheritance is often a combination of variants traveling together, and grouping them reduces noise while preserving the signal that breeders actually manipulate when they cross two lines. These blocks became the raw material for a systematic dissection of the genetic control of kernel and pod weight.

To connect haplotypes with measurable traits, the researchers applied a multi-locus genome-wide association model, a statistical framework that evaluates many chromosome regions simultaneously rather than testing each one in isolation. This approach is better suited to complex quantitative traits, where many loci each contribute a modest effect and single-locus scans can miss real signals or inflate false ones. The analysis identified 29 main-effect quantitative trait loci containing 192 alleles for 100-kernel weight, distributed across 13 of the peanut chromosomes. Together these loci explained 84.55 percent of the phenotypic variation in kernel weight, an unusually high proportion that suggests the panel captured most of the relevant genetic variation. The single largest-effect QTL sat on chromosome Arahy.16, accounting for 36.95 percent of the variation on its own, while chromosome Arahy.12 harbored the largest number of kernel-weight QTLs, seven in total.

The picture for 100-pod weight was similarly rich. The team detected 28 QTLs carrying 208 alleles, again spread over 13 chromosomes and jointly explaining 85.85 percent of the phenotypic variation. Pod weight showed its own hotspots, with clusters of four QTLs each on chromosomes Arahy.12, Arahy.14 and Arahy.19, and the largest-effect locus on Arahy.14 explaining 18.67 percent of the variation. The fact that different chromosomes dominate each trait tells breeders something important: kernel weight and pod weight, though correlated in the field, are not simply two views of the same genetics. Improving one does not automatically improve the other, which is precisely why simultaneous improvement has been so difficult and why the ability to predict crosses that advance both traits at once is valuable.

Perhaps the most striking single finding is a shared pleiotropic QTL on chromosome Arahy.05, a region that influences both traits at the same time. This one locus explained 29.58 percent of the variation in 100-kernel weight and 15.63 percent of the variation in 100-pod weight, making it a genetic fulcrum on which both yield components partially balance. Pleiotropic loci like this are double-edged: the same haplotype that boosts kernel weight may also shift pod weight, sometimes favorably and sometimes not. Identifying such regions explicitly means breeders can track them with DNA markers and choose allele combinations that push both traits in the desired direction rather than discovering trade-offs years later in the field.

With the QTL-allele systems established, the team converted their statistical results into a practical decision tool. They assembled QTL-allele matrices, tables recording which favorable and unfavorable alleles each of the 353 accessions carries at every locus for both traits. These matrices function like a parts inventory for breeding: each parent line is characterized by the specific haplotypes it can contribute to offspring, and a cross is evaluated by the range of allele combinations its progeny could plausibly assemble. The researchers then applied three distinct cross-selection strategies. The first prioritized 100-kernel weight, the second prioritized 100-pod weight, and the third sought a balanced improvement of both traits simultaneously, reflecting the different goals a breeding program might pursue depending on market demands and regional preferences.

The predicted outcomes of these simulated crosses are remarkable. Under the kernel-weight-priority strategy, the thirty optimal crosses showed predicted recombination potentials reaching 162.35 to 170.21 grams for 100-kernel weight. Under the pod-weight-priority strategy, predicted values ranged from 422.17 to 432.21 grams for 100-pod weight. The balanced strategy produced crosses predicted to deliver 136.00 to 154.03 grams for kernel weight together with 371.09 to 413.21 grams for pod weight. In every case these predictions substantially exceed the best values actually observed among the 353 accessions in the population. That gap between observed and predicted performance is the whole point: it represents untapped genetic potential that exists in the germplasm collectively but in no single variety, waiting to be assembled through the right combination of parents.

This in silico approach addresses one of the oldest frustrations in plant breeding. Traditional crossing programs are largely a numbers game: breeders make hundreds of crosses, grow out thousands of progeny, and hope that favorable alleles recombine into winning combinations. The process works but is slow, expensive, and heavily dependent on chance. By predicting which crosses have the highest probability of stacking favorable haplotypes across all relevant loci, the QTL-allele framework lets breeders concentrate their field resources on a shortlist of parent combinations with the greatest expected payoff. The method also identifies elite-allele donor accessions, the specific varieties carrying the best haplotype at each locus, giving breeders a direct shopping list of parents to draw from when constructing their crossing blocks.

The broader significance of the study lies in the pipeline itself, which the authors describe as a sequence of QTL-allele system construction, in silico progeny simulation, and multi-strategy cross selection. This workflow is not limited to peanut or to yield traits. Any crop with adequate SNP data, a sufficiently diverse germplasm panel, and a measurable trait of interest could be run through the same machinery, from rice and wheat to legumes with similarly complex yield architectures. As genotyping costs continue to fall and global germplasm collections become better characterized, the bottleneck in breeding increasingly shifts from generating data to interpreting it, and frameworks like this one turn that interpretation into concrete, actionable crossing decisions.

For a crop that feeds hundreds of millions of people and supports the livelihoods of smallholder farmers across Asia and Africa, the implications are tangible. Peanut yield gains over recent decades have come largely from conventional selection, and further improvement of kernel and pod weight through phenotype-based methods faces diminishing returns as the easy variation has already been captured. The Chinese team’s work demonstrates that a large reservoir of favorable alleles still exists across the global gene pool, distributed piecemeal among accessions that individually look unremarkable. By systematically cataloging those alleles, quantifying their effects, and predicting the crosses that combine them most effectively, the study offers a credible technical route toward the next generation of high-yield peanut varieties, designed first in a database and proven later in the field.

Subject of Research: Haplotype-based QTL-allele analysis and cross prediction for simultaneous improvement of kernel and pod weight in peanut

Article Title: Prediction of optimal crosses based on haplotype-based QTL-allele systems for simultaneous improvement of 100-kernel weight and 100-pod weight in peanut

Article References: Zhang, X., Su, Y., Yang, Z., Zhang, A., Li, D., Qin, H., Liu, Y., & Xue, Q. (2026). Prediction of optimal crosses based on haplotype-based QTL-allele systems for simultaneous improvement of 100-kernel weight and 100-pod weight in peanut. BMC Plant Biology. https://doi.org/10.1186/s12870-026-10043-5

Image Credits: AI Generated

DOI: 10.1186/s12870-026-10043-5

Keywords: peanut, Arachis hypogaea, QTL, haplotypes, genome-wide association study, 100-kernel weight, 100-pod weight, cross prediction, breeding by design, yield improvement, pleiotropy, germplasm

The study, published in SSM – Population Health by Yeongmi Jeong, Meghan M. Skira, and Robbee Wedow, focuses on two polygenic scores, which are numerical summaries of the many small genetic differences scattered across a person’s genome. The first, called the smoking initiation score, aggregates variants associated with whether a person ever becomes a regular smoker. The second, the cigarettes-per-day or CPD score, aggregates variants associated with how heavily someone smokes once they start. Both scores were built from genome-wide association studies, or GWAS, that scanned the genomes of more than a million people of European ancestry to find single-letter DNA differences, known as single nucleotide polymorphisms, that statistically track with smoking behavior. Each person’s score is calculated by adding up the number of copies of the risk-increasing variants they carry, with each variant weighted by the strength of its estimated association.

Neither score is a crystal ball. Out-of-sample predictions indicate that the initiation score explains only about 4 percent of the variation in whether people start smoking, and the CPD score explains about 4 percent of the variation in how much they smoke. Within the study sample itself, a one standard deviation increase in the initiation score raised the probability of ever smoking by roughly 9 percentage points, while a one standard deviation increase in the CPD score was associated with about two additional cigarettes smoked per day among smokers. These are modest effects, but they are large enough, in a sample of this size, to reveal patterns that would otherwise remain invisible.

To measure health, the researchers constructed a standardized index combining thirteen distinct outcomes: self-reported general health, hospitalization, body mass index, seven doctor-diagnosed chronic conditions including high blood pressure, diabetes, cancer, lung disease, heart problems, stroke, and arthritis, and limitations in basic and instrumental activities of daily living. The analysis covered 48,444 person-wave observations from 9,515 unique participants aged 50 to 65, drawn from survey waves between 1994 and 2018. The age range was chosen deliberately, to minimize survival bias, since people with the worst health may not survive to be surveyed at older ages, and to avoid conflating genetic associations with the health effects of retirement.

The central finding is stark. A one standard deviation increase in the CPD score was associated with a 0.043 standard deviation increase in the health index, where higher values indicate worse health. Adding detailed controls for smoking behavior, including current and lifetime smoking status, the exact number of cigarettes smoked per day, the maximum ever smoked, the age of starting and quitting, and teenage smoking, attenuated the association only modestly, to 0.035 standard deviations. For context, that residual association is roughly 38 percent of the average health gap between ever smokers and never smokers in the sample. In other words, even after accounting for everything the surveys record about smoking, the genetic score still carries substantial information about health.

The most provocative result emerged when the researchers restricted the analysis to never smokers, a group of more than 20,000 observations with no reported smoking behavior to explain away. Among these individuals, a one standard deviation increase in either the CPD score or the initiation score was associated with roughly a 0.03 standard deviation increase in the health index. The CPD score also tracked with higher body mass index and obesity among never smokers, with a one standard deviation increase associated with 0.645 units of additional BMI. Because these individuals never smoked, the association cannot be attributed to their own tobacco use, although it could still reflect secondhand smoke exposure, unmeasured smoking-related environments, or something else entirely.

One obvious candidate is the family environment. Because children inherit genes from their parents, a person’s smoking-related genetic score may be correlated with parental smoking, parental substance use, parental education, and childhood socioeconomic conditions, all of which shape later health. The researchers tested this by adding controls for whether parents smoked during the respondent’s childhood, whether parents drank or used drugs heavily enough to cause family problems, parental mortality timing, parental years of schooling, and self-reported childhood financial status. These factors accounted for a modest share of the association, but more than 75 percent of the CPD score’s relationship with health remained, and about 60 percent of the initiation score’s relationship among never smokers survived the most comprehensive specifications.

The researchers then turned to psychology. Prior work has shown that smokers tend to be more impatient, risk-tolerant, and impulsive than non-smokers, and the team asked whether people with higher smoking-related genetic scores differ along similar dimensions. They found little evidence of a link with measured risk preferences, but higher scores were significantly associated with shorter financial planning horizons, greater pessimism, and lower conscientiousness, traits that economists and psychologists have long connected to lower investment in health. When these behavioral and psychological characteristics were added to the models simultaneously, they explained at most 21 percent of the remaining association between the CPD score and health, and 18 percent of the initiation score’s association among never smokers. Even with the fullest set of controls, at least 69 percent of the CPD score’s association with health persisted.

The authors are careful to emphasize what the study does not show. Polygenic scores are not clean measures of direct genetic effects. The GWAS estimates underlying them can absorb indirect genetic effects, such as parental genotypes shaping the environments children experience, as well as confounding from assortative mating, population stratification, collider bias, and case-control ascertainment. There is also the problem of pleiotropy: the same genetic variants that predict smoking may influence other biological processes and traits that affect health directly. Genetic correlation analyses have indeed documented significant overlap between smoking behaviors and a wide range of health, behavioral, and socioeconomic traits. The associations reported here are descriptive, not causal, and the study’s restriction to individuals of European ancestry, a consequence of the limited predictive power of European-derived scores in other populations, means the findings may not generalize more broadly.

Even with those caveats, the implications are significant. Economists studying smoking have long suspected that people who smoke differ from those who do not in ways that extend beyond the behavior itself, and some modeling work has speculated that genetics or correlated risk-taking might explain why smokers face elevated chronic disease and mortality risks beyond what smoking alone explains. This study provides some of the most direct evidence yet that such unobserved heterogeneity is real and partially visible in the genome. Smoking-related polygenic scores, the results suggest, contain information about health that detailed behavioral histories do not capture, pointing toward shared genetic architecture linking tobacco use, personality, and long-run wellbeing. Untangling those threads, the authors argue, will be essential for interpreting genetic risk scores in health research, and for understanding why two people with identical smoking histories can end up with very different health outcomes decades later.

Subject of Research: Associations between smoking-related polygenic scores and health outcomes in midlife

Article Title: Exploring Links Between Health and the Genetic Risk for Smoking

Article References: Jeong, Y., Skira, M. M., & Wedow, R. (2026). Exploring Links Between Health and the Genetic Risk for Smoking. SSM – Population Health, Article 101982. https://doi.org/10.1016/j.ssmph.2026.101982

Image Credits: AI Generated

DOI: 10.1016/j.ssmph.2026.101982

Keywords: polygenic scores, smoking initiation, smoking intensity, Health and Retirement Study, GWAS, pleiotropy, chronic disease, health behavior, genetic risk, conscientiousness, never smokers, population health

The team assembled a natural population of 389 wheat varieties and subjected them to drought stress simulated with 20 percent PEG-6000, a water-soluble polymer that lowers the water potential of the germination medium and mimics the osmotic challenge a seed faces in drying soil. The experiment was repeated across three independent annual environments, a design choice that matters enormously in drought research, because the expression of stress tolerance traits is notoriously sensitive to year-to-year variation in conditions. Rather than relying on a single season, the researchers could ask which genetic signals held up consistently, and which were artifacts of a particular trial.

Six germination-related traits were measured for every variety in every environment: germination rate, shoot length, root length, root number, coleoptile length, and germination potential. From each trait the team calculated a drought resistance coefficient, which expresses how well a variety maintains performance under stress relative to its unstressed baseline. This normalization is important because varieties differ greatly in intrinsic vigor; a coefficient-based approach separates a variety’s inherent capacity from its ability to withstand water deficit specifically.

But the researchers went a step further, addressing a long-standing weakness in stress-tolerance phenotyping. Judging a variety by any single trait can be misleading, since drought tolerance is a composite property that emerges from root growth, shoot elongation, and the speed and completeness of germination acting together. To capture this, the team combined principal component analysis with a membership function method to distill the six traits into a single comprehensive drought tolerance evaluation value, abbreviated as the D-value. This statistical framework, widely used in multi-trait stress evaluations of crops such as maize and cotton, weights each trait according to how much independent information it contributes, producing a ranking that reflects overall tolerance rather than performance on one metric.

The consistency of these rankings across environments was striking. Pairwise correlation coefficients between the D-value rankings from different years all exceeded 0.82, with statistical significance at P below 0.001. In practical terms, a variety that ranked as drought tolerant in one year almost always ranked as tolerant in the others. That level of repeatability is exactly what breeders need, because it suggests the underlying genetic differences are stable rather than being swamped by environmental noise, and it validates the multi-index evaluation system the authors propose as a robust screening tool for germination-stage drought tolerance.

With reliable phenotypes in hand, the team turned to the genome. Using a 660K SNP array, they performed a genome-wide association study, or GWAS, a technique that scans the genome for genetic variants whose frequencies correlate with trait values across the population. The scan identified 270 significantly associated SNPs at a threshold of P less than 1 times 10 to the power of minus 4. Raw association hits, however, often include signals that appear in only one environment or affect only one trait, so the researchers applied a stringent filter: a locus had to be detected in at least two environments or be associated with multiple traits to qualify as an important quantitative trait locus. This filtering narrowed the field to ten important QTLs, of which eight were stably expressed across multiple environments and four showed pleiotropy, meaning a single genomic region influences several different traits simultaneously.

The standout discovery was a stable pleiotropic hotspot on chromosome 5D, comprising two adjacent loci designated qDT-G5D.1 and qDT-G5D.2. This region was associated with root number, germination potential, and germination rate, three traits that together capture much of what determines whether a seedling can establish itself in dry soil. A pleiotropic hotspot of this kind is particularly valuable for breeding, because selecting for the favorable allele at one locus can simultaneously improve several components of tolerance. Coordinated improvement of multiple drought tolerance traits through a single genomic region is far more efficient than pyramiding separate genes for each trait, and the 5D hotspot therefore stands out as a prime candidate for marker-assisted selection.

To move from statistical associations to biological mechanism, the authors integrated previously published transcriptome data, looking for genes within the important QTL intervals that behave differently between drought-tolerant and drought-sensitive varieties. This screen yielded twelve putative candidate genes whose differential expression links them plausibly to the tolerance phenotype. The team then carried out haplotype analysis, examining the natural sequence variants of these candidate genes and testing whether different haplotypes, or combinations of variants inherited together, correspond to different trait performance. The answer was affirmative for ten of the twelve genes: significant phenotypic differences were observed among haplotypes, meaning the sequence variation within these genes tracks with real differences in drought tolerance. That convergence of positional, expression, and haplotype evidence substantially strengthens the case that these genes are functionally involved in germination-stage drought response rather than merely riding along with nearby causal variants.

The broader significance of the work lies in what it hands to wheat breeders and to the research community. The multi-index D-value system offers a standardized, repeatable way to score germination-stage drought tolerance that can be adopted in other germplasm screens. The ten important QTLs, and especially the eight stably expressed ones, provide chromosomal landmarks that can be converted into molecular markers for marker-assisted breeding, allowing breeders to track tolerance alleles in seedling nurseries without waiting for full field evaluations. The 5D pleiotropic hotspot offers a single target with the potential to improve several traits at once, and the ten haplotype-validated candidate genes supply a starting point for functional studies aimed at understanding the molecular pathways that let a germinating seed cope with water deficit.

As climate change intensifies drought pressure on wheat-growing regions worldwide, the earliest days of a plant’s life are becoming an increasingly important focus of crop improvement. By combining multi-environment phenotyping, a statistically rigorous composite evaluation, high-density GWAS, and transcriptome-informed candidate gene analysis, this study demonstrates how modern genomic tools can dissect a complex stress response at a life stage that has long been overlooked. The genetic resources it delivers, from stable QTLs to a pleiotropic hotspot on chromosome 5D and a set of validated candidate genes, are the raw material from which more resilient wheat varieties can be built, ensuring that the critical transition from seed to seedling no longer depends quite so heavily on the vagaries of the weather.

Subject of Research: Genetic loci controlling drought tolerance during wheat germination identified through multi-environment GWAS

Article Title: Multi-environment genome-wide association study reveals genetic loci associated with drought tolerance during wheat germination

Article References: Liu, T., Gao, L., Li, J., Qiao, N., Jiang, D., Yang, J., Jia, J., Wang, J., Sun, D., & Li, N. (2026). Multi-environment genome-wide association study reveals genetic loci associated with drought tolerance during wheat germination. Plant Cell Reports, 45(10), Article 318. https://doi.org/10.1007/s00299-026-04014-1

Image Credits: AI Generated

DOI: 10.1007/s00299-026-04014-1

Keywords: wheat, drought tolerance, germination, GWAS, QTL, chromosome 5D, pleiotropy, SNP, candidate genes, haplotype analysis, marker-assisted breeding, PEG-6000

To understand why this matters, it helps to revisit how Mendelian randomization works. The technique uses genetic variants, typically single nucleotide polymorphisms, as instrumental variables to estimate whether an exposure truly causes an outcome. Its validity rests on three assumptions: the instrument must be strongly associated with the exposure, it must be independent of confounders, and it must influence the outcome only through the exposure. In practice, the second and third assumptions are fragile. Horizontal pleiotropy, in which a variant affects the outcome through pathways unrelated to the exposure, violates the exclusion restriction and has spawned an entire industry of robust correction methods. But dynastic effects present a different and arguably more insidious problem, one that cannot be fixed by looking harder at population-level summary statistics alone.

Dynastic effects arise specifically within families. In a trio consisting of two parents and an offspring, parental genotypes can influence offspring traits through unobserved parental phenotypes or other intergenerational mechanisms, such as the family environment parents create or lifestyle behaviors passed down across generations. These pathways create a persistent confounding structure in the association between a SNP and an outcome, violating the independence assumption that Mendelian randomization depends on. Crucially, the bias does not average out with large sample sizes; it simply persists. The only reliable way to block these non-causal pathways is to bring family-level data into the analysis, which is precisely what the authors set out to do in a more comprehensive way than previous designs allowed.

Existing family-based approaches each come with trade-offs. A standard trio-based method fits a linear model regressing the offspring’s phenotype on both offspring and parental genotypes, which effectively controls for dynastic bias but discards parental phenotypic information entirely, limiting statistical power. Mother-offspring pair designs use structural equation models to separate maternal and fetal effects, but they ignore paternal contributions and discard offspring genotype data, confining their use largely to birth-weight research. Sibling-pair designs, increasingly popular as biobanks accumulate sibling data, exploit differences between siblings to cancel out dynastic effects, but they cannot be applied to trios with an only child. Perhaps most importantly, none of these frameworks had been extended to multivariate analysis, leaving correlated traits to be examined one at a time.

FT-SEM fills that gap by modeling multiple correlated phenotypes through a shared latent factor. In the framework, the genetic effect at each SNP is partitioned into paternal, maternal, and offspring components, and the correlations among several traits are assumed to be driven by a common unobserved construct. For the obesity application, that latent factor was built from body mass index, body fat percentage, and waist-to-hip ratio, three moderately correlated measures that each capture a different dimension of adiposity. By estimating SNP effects on the shared factor rather than on each trait separately, the method reduces phenotype-specific measurement error and biological heterogeneity. Latent grandparental genotypes, which are never actually observed, are incorporated as latent variables with variances constrained to match observed genotypes under Hardy-Weinberg equilibrium, and Mendelian transmission paths are fixed at 0.5 to identify the model.

One technical obstacle the authors had to overcome is that the genetic effects of interest are not directly identifiable when the latent factor is estimated jointly with the SNP effect. Their solution is a two-stage procedure: first, a null model excluding all genetic variant components is fitted to obtain reference factor loadings; then the full model is fitted, and a correction factor computed from the ratio of loadings between the two fits is applied to recover accurate estimates of the offspring, paternal, and maternal effects. The offspring effect is the primary target for downstream interpretation and causal inference, because it represents the direct genetic effect on the offspring’s phenotype, while the parental effect estimates serve as diagnostic quantities that characterize and quantify dynastic effects themselves.

The simulation evidence is extensive. Across 10,000 replicates per setting, with sample sizes ranging from 1,000 to 3,000 trios and dynastic effect magnitudes of zero, 0.2 percent, and 0.4 percent of phenotypic variance, FT-SEM produced unbiased effect estimates with well-controlled type I error rates and confidence interval coverage near the nominal 95 percent level. A population-based comparator, IndSEM, was more efficient when no dynastic effects were present, achieving the lowest root mean squared error and the narrowest intervals, but its performance collapsed as dynastic effects increased, producing biased estimates and inflated false positive rates. FT-SEM consistently outperformed PSEM, a multivariate extension of the conventional trio-based linear model, by exhibiting lower root mean squared error and higher statistical power, a direct benefit of incorporating parental phenotypic data. Additional simulations showed the framework remained robust even when the assumed latent-factor structure or cross-generational measurement invariance was deliberately misspecified.

The real-data application turned on 778 complete family trios from the UK Biobank, identified through stringent kinship-based filtering and restricted to individuals of European ancestry to minimize population stratification. After quality control on genotypes and phenotypes, the authors ran genome-wide scans under FT-SEM, PSEM, and IndSEM. No locus reached genome-wide significance, which is unsurprising given the modest trio sample size, but that was not the point. The family-corrected summary statistics for the latent obesity factor served as an unbiased outcome base for downstream two-sample Mendelian randomization, paired with sibling-based exposure data covering 99,998 individuals for the within-family analyses and independent-individual data of 315,133 and 757,601 participants for the comparison arm. Because complete trios are typically excluded from conventional population-based GWAS, the two-sample design involved no sample overlap between exposure and outcome datasets.

The results were revealing. For high-density lipoprotein cholesterol, the population-based IndSEM framework detected significant causal associations under weighted median, inverse-variance weighting, and MRcML methods, with p values around 0.041 to 0.042, and the family-based methods, though underpowered, pointed in the same direction. For systolic blood pressure, the picture was different: IndSEM flagged a significant association under MR-Egger regression, but the family-based FT-SEM and PSEM estimates were substantially attenuated toward zero and nowhere near significance. A sensitivity analysis using a harmonized instrument set drawn from sibling-based GWAS found no significant blood pressure association either, though the authors caution that weaker sibling-based instruments and reduced statistical power likely contributed to that attenuation. Taken together, the pattern suggests the apparent blood pressure to obesity link reported in conventional studies may be spurious, driven largely by familial confounding rather than genuine causation.

The authors are candid about limitations. The trio sample is small and underpowered for discovery, the framework assumes multivariate normality, additive genetic effects, and a constant genetic architecture across generations, and it is currently restricted to trios with a single offspring and does not model gene-environment or gene-gene interactions. The structural equation model is also computationally heavier than standard linear models, since it must be re-estimated for every SNP, though the authors argue the burden is manageable on modern parallel computing platforms. Still, the broader message is hard to ignore: as Mendelian randomization grows in popularity, susceptibility to dynastic effects has become a critical concern, and discrepancies between population-based and family-based estimates should generally be resolved in favor of the family-based result. FT-SEM offers researchers a way to have both rigor and efficiency, blocking non-causal intergenerational pathways while squeezing maximal information from every member of the family trio, and providing a template for future pedigree-wide analyses that could dissect direct genetic effects, assortative mating, and other intergenerational mechanisms across entire family trees.

Subject of Research: A multivariate GWAS framework using family trio data and parental phenotypes to control dynastic effects in genetic association and Mendelian randomization analyses

Article Title: Multivariate GWAS framework for family trios with parental phenotypes to control dynastic effects

Article References: Zhang, S., Mai, J.-H., Zhong, Q., Zhu, Q.-W., Li, Y.-S., Wu, X.-B., & Zhou, J.-Y. (2026). Multivariate GWAS framework for family trios with parental phenotypes to control dynastic effects. BMC Bioinformatics, 27(1), Article 222. https://doi.org/10.1186/s12859-026-06597-8

Image Credits: AI Generated

DOI: 10.1186/s12859-026-06597-8

Keywords: GWAS, Mendelian randomization, family trios, dynastic effects, structural equation modeling, UK Biobank, obesity, direct genetic effects, parental phenotypes, latent factor, pleiotropy, BMC Bioinformatics

Harvest index sits at the heart of crop productivity. It is the ratio of grain yield to total above-ground biomass, and it captures how efficiently a plant partitions its resources between the parts people eat and the parts they do not. While decades of research have documented how maize domestication restructured the plant’s overall architecture and the architecture of its male and female inflorescences, the question of whether and how domestication genetically shaped harvest index had remained open. Answering that question requires comparing wild, traditional, and modern material in a controlled framework, and that is precisely what a research team led by Li Guo of China Agricultural University set out to do.

The team assembled a diversity panel that spans the full arc of maize evolution: thirty lines of teosinte, the wild ancestor; twenty maize landraces, representing the semi-domesticated forms grown by early farmers; and fifty modern inbred lines, the elite genetic material that underpins contemporary maize breeding. By evaluating harvest index across these three groups, the researchers quantified the cumulative effect of domestication on this trait. The result was striking. On average, domestication increased harvest index by 81 percent, meaning that modern maize devotes a dramatically larger share of its total plant mass to grain than its wild progenitor does.

An 81 percent increase is not a subtle adjustment. It reflects a deep reorganization of how the maize plant allocates carbon, nutrients, and developmental resources. Wild teosinte invests heavily in stems, branches, tassels, and protective structures that help it survive in the wild, but produce little food. Domesticated maize, through the combined action of many genes selected over generations, was redirected toward producing large, exposed ears packed with kernels. Understanding which genetic elements underlie this shift is essential both for reconstructing the history of the crop and for guiding future breeding efforts, particularly as breeders seek to push yield gains further in the face of climate change and growing global demand for grain.

To dissect the genetics behind the shift, Guo and colleagues turned to quantitative trait locus mapping, a classical but powerful technique that links regions of the genome to measurable trait differences in segregating populations. The researchers created two maize-teosinte crossing populations, each mixing chromosomes from the wild ancestor and domesticated maize in different combinations, and then mapped which genomic intervals were associated with variation in harvest index. The scan identified twenty quantitative trait loci, confirming that harvest index is a genetically complex trait shaped by many loci spread across the maize genome, as expected for a trait so tightly connected to overall plant growth and partitioning.

Among those twenty loci, one stood out. A major locus, which the researchers named Harvest Index 1, or HI1, showed consistent and large effects in both crossing populations. Consistency across independent populations is an important signal in genetics: it indicates that the effect is real and robust rather than an artifact of a particular cross. The magnitude of the HI1 effect, together with its reproducibility, marks it as a principal genetic contributor to the harvest index difference between teosinte and maize.

To pin down the function of HI1 more precisely, the team used near-isogenic lines, pairs of plants that are genetically almost identical except for a small chromosome segment carrying either the teosinte or the maize version of the locus. This design isolates the effect of the target gene from the noise of the rest of the genome. When the researchers compared near-isogenic lines carrying the maize allele of HI1 with those carrying the teosinte allele, they confirmed that the maize allele increases harvest index. But the locus also influenced total biomass and grain weight, revealing a pleiotropic effect: a single genetic change simultaneously affecting multiple traits of the plant.

Pleiotropy is a recurring theme in crop domestication genetics. The best-known domestication genes in maize, such as those controlling branching and glume architecture, also produce coordinated suites of changes, because the regulatory variants involved often act on shared developmental pathways rather than on a single trait. The apparent role of HI1 in modulating biomass and grain weight alongside harvest index suggests it may tune the growth and partitioning machinery of the plant more broadly, rather than acting narrowly on grain fill alone. Notably, as the research highlight’s framing emphasizes, HI1 joins a growing list of domestication-selected enhancers, regulatory DNA elements that altered gene expression under human selection rather than changes in the protein-coding sequence itself. This regulatory mode of domestication has emerged as a hallmark of how early farmers, deliberately or not, reshaped crops: subtle changes in when and where genes are turned on can reorganize a plant’s architecture without breaking the genes outright.

The identification of HI1 carries implications well beyond evolutionary history. Harvest index has been a cornerstone of modern yield improvement, and the gains achieved during the Green Revolution of the twentieth century came largely from breeding varieties that invested more heavily in grain and less in straw. Knowing the specific loci that raise harvest index gives breeders molecular targets for the next round of improvement, and the teosinte alleles preserved at loci like HI1 represent a reservoir of functional variation that could be deployed either to boost or, in some contexts, to rebalance trait architecture. As genomic tools make it increasingly feasible to edit or introgress specific regulatory regions, the boundary between the wild ancestor and the modern crop becomes not just a historical record but a practical resource. The study of harvest index during maize domestication, with HI1 as its newly identified centerpiece, illustrates how a single enhancer selected by ancient farmers continues to shape the productivity of one of humanity’s most important food crops, and how decoding the genetic logic of domestication can inform the future of crop improvement.

Subject of Research: Identification of a domestication-selected enhancer, Harvest Index 1, that increased harvest index during maize domestication

Article Title: Yet another domestication-selected enhancer

Article References: Lyu, J. (2026). Yet another domestication-selected enhancer. Nature Plants. https://doi.org/10.1038/s41477-026-02416-3

Image Credits: AI Generated

DOI: 10.1038/s41477-026-02416-3

Keywords: maize domestication, teosinte, harvest index, HI1, enhancer, QTL mapping, near-isogenic lines, plant breeding, biomass allocation, grain yield, pleiotropy, Nature Plants