Segregation distortion: Utilizing simulated genotyping data to evaluate statistical methods

Genetic Insights Into Perinatal Outcomes of Maternal Antihypertensive Therapy During Pregnancy

Importance

Limited information exists regarding the impact of pharmacotherapy in pregnancy due to ethical concerns of unintended fetal harm. Yet, maternal prescriptive drug use for chronic conditions such as hypertension is common.

Objective

To investigate potential causal relationships between perturbing maternal genetic variants influencing antihypertensive drug targets and perinatal outcomes among offspring using mendelian randomization (MR).

Design, Setting, and Participants

This 2-sample MR study used individual-level single-nucleotide variation (SNV) outcome data from mother-father-offspring trios with complete genetic and phenotypic information from the Norwegian Mother, Father and Child Cohort Study (MoBa) and summary-level SNV exposure data from UK Biobank participants sourced from the Integrative Epidemiology Unit OpenGWAS project. Pregnant individuals were recruited across Norway during their routine ultrasonography examination at 18 weeks' gestation between June 1999 and December 2008, and mothers, fathers, and offspring were followed up after birth. Novel genetic instruments for maternal antihypertensive drug targets that act via systolic blood pressure (SBP) were derived from individual-level data analyzed in January 2018. Two-sample multivariable MR analysis of these maternal drug targets and offspring outcomes were performed between January 2023 and April 2024.

Exposures

Maternal genetic variants associated with drug targets for treatments of hypertension, as specified in the National Health Service dictionary of medicines and devices.

Main Outcomes and Measures

Offspring outcomes were Apgar score at 1 minute and 5 minutes, offspring developmental score at 6 months, birth length, birth weight z score, gestational age, head circumference, and congenital malformation. Maternal hypertensive disorders of pregnancy were a positive control.

Results

The MoBa sample contained 29 849 family trios, with a mean (SD) maternal age of 30.2 (18.6) years and a mean (SD) paternal age of 32.8 (13.1) years; 51.1% of offspring were male. Seven independent SNVs were identified as influencing maternal SBP via the antihypertensive drug target instruments. For higher levels of maternal SBP acting through the CACNB2 calcium channel blocker target, the estimated change in gestational age was 3.99 days (95% CI, 0.02-7.96 days) per 10-mm Hg decrease in SBP. There was no evidence of differential risk for measured perinatal outcomes from maternal SBP acting through drug targets for multiple hypertensive subclasses, such as between the ADRB1 β-adrenoceptor-blocking target and risk of congenital malformation (estimated odds ratio, 0.28 [95% CI, 0.02-4.71] per 10-mm Hg decrease in SBP). Maternal and paternal SBP acting through the EDNRA vasodilator antihypertensive target did not have a potential causal effect on birth weight z score, with respective β estimates of 0.71 (95% CI, -0.09 to 1.51) and 0.72 (95% CI, -0.08 to 1.53) per 10-mm Hg decrease in SBP.

Conclusions and Relevance

The findings provided little evidence to indicate that perturbation of maternal genetic variants for SBP that influence antihypertensive drug targets had potential causal relationships with measures of perinatal development and health within this study. These findings may be triangulated with existing literature to guide physicians and mothers in decisions about antihypertensive use during pregnancy.

A High-Resolution Linkage Map Construction and QTL Analysis for Morphological Traits in Anthurium (Anthurium andraeanum Linden)

Anthurium andraeanum Linden is a prominent ornamental plant belonging to the family Araceae and is cultivated worldwide. The morphology characteristics are crucial because they significantly impact ornamental values, commercial properties, and the efficiency of space utilization in production. However, only a few related investigations have been conducted in anthurium to date. In this study, an F1 genetic segregation population containing 160 progenies was generated through hybridization between potted and cut anthurium varieties. Fifteen morphological traits were assessed and revealed substantial levels of genetic variation and widespread positive correlation. Based on specific length amplified fragment (SLAF) sequencing technology, 8171 single nucleotide polymorphism (SNP) markers were developed, and the high-density linkage map of 2202.27 cM in length distributed on 15 linkage groups was constructed successfully, with an average distance of 0.30 cM. Using the inclusive composite interval mapping (ICIM) method, 59 QTLs related to 15 key morphological traits were successfully identified, which explained phenotypic variance (PVE) ranging from 6.21% to 17.74%. Thirty-three of those associated with 13 traits were designated as major QTLs with PVE > 10%. These findings offer valuable insights into the genetic basis of quantitative traits and are beneficial for molecular marker-assisted selection (MAS) in anthurium breeding.

Identification of quantitative trait loci associated with R1-mediated resistance to powdery mildew and sex determination in hop (Humulus lupulus L.)

KEY MESSAGE

Two QTL were identified using linkage mapping approaches, one on hop linkage group 3 (qHl_Chr3.PMR1) associated with powdery mildew resistance and a second on linkage group 10 (cqHl_ChrX.SDR1) associated with sex determination. Hop (Humulus lupulus L.) is a dioecious species cultivated for use in beer. Hop powdery mildew, caused by Podosphaera macularis, is a constraint in many growing regions. Thus, identifying markers associated with powdery mildew resistance and sex provides the opportunity to pyramid R-genes and select female plants as seedlings, respectively. Our objectives were to characterize the genetic basis of R1-mediated resistance in the cultivar Zenith which provides resistance to pathogen races in the US, identify quantitative trait loci (QTL) associated with R1 and sex, and develop markers for molecular breeding-based approaches. Phenotypic evaluation of the population indicated that R1-based resistance and sex are inherited monogenically. We constructed a genetic map using 1339 single nucleotide polymorphisms (SNPs) based upon genotype-by-sequencing of 128 F₁ progeny derived from a Zenith × USDA 21058M biparental population. SNPs were assigned to 10 linkage groups comprising a map length of 1204.97 cM with an average density of 0.94 cM/marker. Quantitative trait locus mapping identified qHl_Chr3.PMR1, associated with R1 on linkage group 3 (LOD = 23.57, R² = 57.2%), and cqHl_ChrX.SDR1, associated with sex on linkage group 10 (LOD = 5.42, R² = 25.0%). Kompetitive allele-specific PCR (KASP) assays were developed for both QTL and assessed against diverse germplasm. Our results indicate that KASP markers associated with R1 may be limited to materials that are pedigree-related to Zenith, whereas markers associated with sex may be transferable across populations. The high-density map, QTL, and associated KASP markers will enable selecting for sex and R1-mediated resistance in hop.

Identification of Stable and Multiple Environment Interaction QTLs and Candidate Genes for Fiber Productive Traits Under Irrigated and Water Stress Conditions Using Intraspecific RILs of Gossypium hirsutum var. MCU5 X TCH1218

Cotton productivity under water-stressed conditions is controlled by multiple quantitative trait loci (QTL). Enhancement of these productivity traits under water deficit stress is crucial for the genetic improvement of upland cotton, Gossypium hirsutum. In the present study, we constructed a genetic map with 504 single nucleotide polymorphisms (SNPs) covering a total span length of 4,416 cM with an average inter-marker distance of 8.76 cM. A total of 181 intra-specific recombinant inbred lines (RILs) were derived from a cross between G. hirsutum var. MCU5 and TCH1218 were used. Although 2,457 polymorphic SNPs were detected between the parents using the CottonSNP50K assay, only 504 SNPs were found to be useful for the construction of the genetic map. In the SNP genotyping, a large number of SNPs showed either >20% missing data, duplication, or segregation distortion. However, the mapped SNPs of this study showed collinearity with the physical map of the reference genome (G. hirsutum var.TM-1), indicating that there was no chromosomal rearrangement within the studied mapping population. RILs were evaluated under multi-environments and seasons for which the phenotypic data were acquired. A total of 53 QTL controlling plant height (PH), number of sympodial branches, boll number (BN), and boll weight (BW) were dissected by QTL analysis under irrigated and water stress conditions. Additionally, it was found that nine QTL hot spots not only co-localized for more than one investigated trait but were also stable with major QTL, i.e., with > 10% of phenotypic variation. One QTL hotspot on chromosome 22 flanked by AX-182254626-AX-182264770 with a span length of 89.4 cM co-localized with seven major and stable QTL linked to a number of sympodial branches both under irrigated and water stress conditions. In addition, putative candidate genes associated with water stress in the QTL hotspots were identified. Besides, few QTL from the hotspots were previously reported across various genetic architects in cotton validating the potential applications of these identified QTL for cotton breeding and improvement. Thus, the major and stable QTL identified in the present study would improve the cotton productivity under water-limited environments through marker-assisted selection.

Berry Anthocyanin, Acid, and Volatile Trait Analyses in a Grapevine-Interspecific F2 Population Using an Integrated GBS and rhAmpSeq Genetic Map

Increased map density and transferability of markers are essential for the genetic analysis of fruit quality and stress tolerance in interspecific grapevine populations. We used 1449 GBS and 2000 rhAmpSeq markers to develop a dense map for an interspecific F2 population (VRS-F2) that was derived by selfing a single F1 from a Vitis riparia x 'Seyval blanc' cross. The resultant map contained 2519 markers spanning 1131.3 cM and was highly collinear with the Vitis vinifera 'PN40024' genome. Quantitative trait loci (QTL) for berry skin color and flower type were used to validate the map. Four rhAmpSeq transferable markers were identified that can be used in pairs (one pistillate and one hermaphroditic) to predict pistillate and hermaphrodite flower type with ≥99.7% accuracy. Total and individual anthocyanin diglucoside QTL mapped to chromosome 9 near a 5-O-GLUCOSYLTRANSFERASE candidate gene. Malic acid QTL were observed on chromosome 1 and 6 with two MALATE DEHYRDROGENASE CYTOPLASMIC 1 and ALUMINUM-ACTIVATED MALATE TRANSPORTER 2-LIKE (ALMT) candidate genes, respectively. Modeling malic acid identified a potential QTL on chromosome 8 with peak position in proximity of another ALMT. A first-ever reported QTL for the grassy smelling volatile (E)-2-hexenal was found on chromosome 2 with a PHOSPHOLIPID HYDROPEROXIDE GLUTATHIONE PEROXIDASE candidate gene near peak markers.

Male and female recombination landscapes of diploid Arabidopsis arenosa

The number and placement of meiotic crossover events during meiosis have important implications for the fidelity of chromosome segregation as well as patterns of inheritance. Despite the functional importance of recombination, recombination landscapes vary widely among and within species, and this can have a strong impact on evolutionary processes. A good knowledge of recombination landscapes is important for model systems in evolutionary and ecological genetics, since it can improve interpretation of genomic patterns of differentiation and genome evolution, and provides an important starting point for understanding the causes and consequences of recombination rate variation. Arabidopsis arenosa is a powerful evolutionary genetic model for studying the molecular basis of adaptation and recombination rate evolution. Here, we generate genetic maps for 2 diploid A. arenosa individuals from distinct genetic lineages where we have prior knowledge that meiotic genes show evidence of selection. We complement the genetic maps with cytological approaches to map and quantify recombination rates, and test the idea that these populations might have distinct patterns of recombination. We explore how recombination differs at the level of populations, individuals, sexes and genomic regions. We show that the positioning of crossovers along a chromosome correlates with their number, presumably a consequence of crossover interference, and discuss how this effect can cause differences in recombination landscape among sexes or species. We identify several instances of female segregation distortion. We found that averaged genome-wide recombination rate is lower and sex differences subtler in A. arenosa than in Arabidopsis thaliana.

Loss of chromatin remodeler DDM1 causes segregation distortion in Arabidopsis thaliana

In ddm1 mutants, the DNA methylation is primarily affected in the heterochromatic region of the chromosomes, which is associated with the segregation distortion of SNPs in the F2 progenies. Segregation distortion (SD) is common in most genetic mapping experiments and a valuable resource to determine how gene loci induce deviation. Meiotic DNA crossing over and SD are under the control of several types of epigenetic modifications. DNA methylation is an important regulatory epigenetic modification that is inherited across generations. In the present study, we investigated the relationship between SD and DNA methylation. The ecotypes Col-0/C24 and chromatin remodeler mutants ddm1-10/Col and ddm1-15/C24 were reciprocally crossed to obtain F2 generations. A total of 300 plants for each reciprocally crossed plant in the F2 generations were subjected to next-generation sequencing to detect the single-nucleotide polymorphisms (SNPs) as DNA markers. All SNPs were analyzed using the Chi-square test method to determine their segregation ratio in F2 generations. Through the segregation ratio, whole-genome SNPs were classified into 16 classes. In class 10, the SNPs in the reciprocal crosses of wild type showed the expected Mendelian ratio of 1:2:1, while those in the reciprocal crosses of ddm1 mutants showed distortion. In contrast, all SNPs in class 16 displayed a normal 1:2:1 ratio, and class 1 showed SD, regardless of wild type or mutants, as assessed using CAPS (cleaved amplified polymorphic sequences) marker analysis to confirm the next-generation sequencing. In ddm1 mutants, the DNA methylation is highly reduced throughout the whole genome and more significantly in the heterochromatic regions of chromosomes. Our results showed that the ddm1 mutants exhibit low levels of DNA methylation, which facilitates the SD of SNPs primarily located in the heterochromatic region of chromosomes by reducing the heterozygous ratio. The present study will provide a strong base for future research focusing on the impact of DNA methylation on trait segregation and plant evolution.

Clubroot resistance derived from the European Brassica napus cv. 'Tosca' is not effective against virulent Plasmodiophora brassicae isolates from Alberta, Canada

In this study, clubroot resistance in the resynthesized European winter Brassica napus cv. 'Tosca' was introgressed into a Canadian spring canola line '11SR0099', which was then crossed with the clubroot susceptible spring line '12DH0001' to produce F₁ seeds. The F₁ plants were used to develop a doubled haploid (DH) mapping population. The parents and the DH lines were screened against 'old' pathotypes 2F, 3H, 5I, 6M and 8N of the clubroot pathogen, Plasmodiophora brassicae, as well as against the 'new' pathotypes 5X, 5L, 2B, 3A, 3D, 5G, 8E, 5C, 8J, 5K, 3O and 8P. Genotyping was conducted using a Brassica 15K SNP array. The clubroot screening showed that 'Tosca, '11SR0099' and the resistant DH lines were resistant to three (2F, 3H and 5I) of the five 'old' pathotypes and four (2B, 3O, 8E and 8P) of the 12 'new' pathotypes, while being moderately resistant to the 'old' pathotype 8N and the 'new' pathotypes 3D and 5G. 'Tosca' was susceptible to isolates representing pathotype 3A (the most common among the 'new' pathotypes) as well as pathotypes 6M, 5X, 5L, 5K and 8J. Linkage analysis and QTL mapping identified a ca. 0.88-0.95 Mb genomic region on the A03 chromosome of 'Tosca' as conferring resistance to pathotypes 2F, 3H, 5I, 2B, 3D, 5G, 8E, 3O and 8P. The identified QTL genomic region housed the CRk, Crr3 and CRd gene(s). However, the susceptibility of 'Tosca' to most of the common virulent pathotypes makes it unattractive as a sole CR donor in the breeding of commercial canola varieties in western Canada.

QTL mapping of agronomic traits in wheat using the UK Avalon × Cadenza reference mapping population grown in Kazakhstan

Background

The success of wheat production is largely dependent on local breeding projects that focus on the development of high-yielding cultivars with the use of novel molecular tools. One strategy for improving wheat productivity involves the deployment of diverse germplasms with a high potential yield. An important factor for achieving success involves the dissection of quantitative trait loci (QTLs) for complex agronomic traits, such as grain yield components, in targeted environments for wheat growth.

Methods

In this study, we tested the United Kingdom (UK) spring set of the doubled haploid (DH) reference population derived from the cross between two British cultivars, Avalon (winter wheat) and Cadenza (spring wheat), in the Northern, Central, and Southern regions (Karabalyk, Karaganda, Kyzylorda) of Kazakhstan over three years (2013–2015). The DH population has previously been genotyped by UK scientists using 3647 polymorphic DNA markers. The list of tested traits includes the heading time, seed maturation time, plant height, spike length, productive tillering, number of kernels per spike, number of kernels per meter, thousand kernel weight, and yield per square meter. Windows QTL Cartographer was applied for QTL mapping using the composite interval mapping method.

Results

In total, 83 out of 232 QTLs were identified as stable QTLs from at least two environments. A literature survey suggests that 40 QTLs had previously been reported elsewhere, indicating that this study identified 43 QTLs that are presumably novel marker-trait associations (MTA) for these environments. Hence, the phenotyping of the DH population in new environments led to the discovery of novel MTAs. The identified SNP markers associated with agronomic traits in the DH population could be successfully used in local Kazakh breeding projects for the improvement of wheat productivity.