Specific Quantitative Trait Loci Research Articles

Molecular mechanisms of resistance to bovine mastitis

Identification of disease-causing genes that underlie complex traits, such as susceptibility to disease, is the primary goal of genetic and biomedical researches. Bovine mastitis is the most prevalent disease with substantial economic losses and animal welfare issues in dairy industry. Quantitative trait loci (QTLs) and single-nucleotide polymorphism or gene association data as biomarkers were used in the selection programs in this industry. In the current microarray-based study, the gene expression profiles of cows selected for high or low mastitis susceptibility affected by a specific QTL on Bos taurus autosome 18 (BTA18) with induced E. coli mastitis in acute phase (time point 24h) were used to perform functional genomics. The results of a functional enrichment analysis showed that the Q genotype (favorable allele) helps to correct suppression of the metabolic process during infection and help to moderate signaling pathways related to immune response and inflammation. The most important pathways suppressed by Q genotype are related to triggering the inflammatory and immune responses which cause moderate responses to infection. Furthermore, the pathways related to growth, cell proliferation and reproductive process enriched by down-regulated expressed genes. The results of network analysis showed that transcription factors TP53, SP1 ligands such as INS, IL1B, IFNG, TGFB1, EGF and protein kinase MAPK1 are the common regulatory factors in the constructed networks.These findings regardless of revealing the new biological process of gene regulatory network analysis of differentially expressed genes (DE)s in bovine mastitis, could also introduce the important key regulatory elements in resistance to mastitis.

Disentangling group specific QTL allele effects from genetic background epistasis using admixed individuals in GWAS: An application to maize flowering.

When handling a structured population in association mapping, group-specific allele effects may be observed at quantitative trait loci (QTLs) for several reasons: (i) a different linkage disequilibrium (LD) between SNPs and QTLs across groups, (ii) group-specific genetic mutations in QTL regions, and/or (iii) epistatic interactions between QTLs and other loci that have differentiated allele frequencies between groups. We present here a new genome-wide association (GWAS) approach to identify QTLs exhibiting such group-specific allele effects. We developed genetic materials including admixed progeny from different genetic groups with known genome-wide ancestries (local admixture). A dedicated statistical methodology was developed to analyze pure and admixed individuals jointly, allowing one to disentangle the factors causing the heterogeneity of allele effects across groups. This approach was applied to maize by developing an inbred “Flint-Dent” panel including admixed individuals that was evaluated for flowering time. Several associations were detected revealing a wide range of configurations of allele effects, both at known flowering QTLs (Vgt1, Vgt2 and Vgt3) and new loci. We found several QTLs whose effect depended on the group ancestry of alleles while others interacted with the genetic background. Our GWAS approach provides useful information on the stability of QTL effects across genetic groups and can be applied to a wide range of species.

Construction of a single nucleotide polymorphism marker based QTL map and validation of resistance loci to bacterial wilt caused by Ralstonia solanacearum species complex in tomato

Bacterial wilt (BW), caused by Ralstonia solanacearum species complex is one of the major biotic factors limiting tomato production in the humid tropics. Pyramiding of resistance genes through marker-assisted selection is an efficient way to develop durable BW resistant cultivars. Tomato line ‘Hawaii 7996’ (H7996) is a stable and robust resistance source against various strains of the species complex. Major BW resistance quantitative trait loci (QTLs) Bwr-12 and Bwr-6, and several minor or strain specific QTLs have been coarse-mapped in this line, but none has been fine-mapped and validated. The objective of the current study was to construct a high density genetic map using single-nucleotide polymorphism (SNP) markers derived from genotyping-by-sequencing, fine-map Bwr-12 and Bwr-6 and determine the effects of these QTLs using a near isogenic line (NIL) population. A high density genetic map using 1604 SNP markers with an average distance of 0.82 cM was developed for 188 F9 recombinant inbred lines derived from the cross H7996 × WVa700. A total of seven QTLs associated with BW resistance to race 1-phylotype I (R. pseudosolanacearum) or/and race 3-phylotype II (R. solanacearum) strains were located on chromosomes 6 (Bwr-6.1, 6.2, 6.3 and 6.4) and 12 (Bwr-12.1, Bwr-12.2 and Bwr-12.3) with logarithm of odds (LOD) scores of 6.2–15.6 and 6.2–31.1, explaining 14.2–33.4% and 15.9–53.9% of the total phenotypic variation contributed from H7996, respectively. To validate the genetic effects of the two QTL regions, a set of 80 BC3F3 NILs containing different sections of Bwr-6 with or without Bwr-12 was phenotyped for disease severity after challenge with either race 1-phylotype I Pss4 or race 3-phylotype II Pss1632 BW strains over two seasons. Bwr-6.1 specific to Pss4 and Bwr-6.3 specific to Pss1632 were mapped to an interval of 5.0 cM (P < 0.05) between 6_33,444,000_SLM6-47 and 6_33,868,000_SLM6-124 SNP marker, and to 2.7 cM (P < 0.01) between positions 6_35,949,000 _SLM6-107 to 6_36,750,000_SLM6-82 marker, respectively. In addition, the specific effect of Bwr-12 for resistance to Pss4 (LOD score of 5.8-16.1, P < 0.01) was confirmed.

Genotyping-by-sequencing and multilocation evaluation of two interspecific backcross populations identify QTLs for yield-related traits in pigeonpea.

This study has identified single-nucleotide polymorphism (SNP) markers associated with nine yield-related traits in pigeonpea by using two backcross populations (BP) developed through interspecific crosses and evaluating them at two locations and 3years. In both the populations, markers have shown strong segregation distortion; therefore, a quantitative trait locus (QTL) mapping mixed model was used. A total of 86 QTLs explaining 12-21% phenotypic variation were detected in BP-1. On the other hand, 107 QTLs explaining 11-29% phenotypic variation were detected in BP-2. Although most QTLs were environment and trait specific, few stable and consistent QTLs were also detected. Interestingly, 11 QTLs in BP-2 were associated with more than one trait. Among these QTLs, eight QTLs associated with days to 50% flowering and days to 75% maturity were located on CcLG07. One SNP "S7_14185076" marker in BP-2 population has been found associated with four traits, namely days to 50% flowering, days to 75% maturity, primary branches per plant and secondary branches per plant with positive additive effect. Hence, the present study has not only identified QTLs for yield-related traits, but also discovered novel alleles from wild species, which can be used for improvement of traits through genomics-assisted breeding.

Quantitative Trait Locus Analysis of Leaf Morphology Indicates Conserved Shape Loci in Grapevine.

Leaf shape in plants plays important roles in water use, canopy structure, and physiological tolerances to abiotic stresses; all important traits for the future development and sustainability of grapevine cultivation. Historically, researchers have used ampelography, the study of leaf shape in grapevines, to differentiate Vitis species and cultivars based on finite leaf attributes. However, ampelographic measurements have limitations and new methods for quantifying shape are now available. We paired an analysis of finite trait attributes with a 17-point landmark survey and generalized Procrustes analysis (GPA) to reconstruct grapevine leaves digitally from five interspecific hybrid mapping families. Using the reconstructed leaves, we performed three types of quantitative trait loci (QTL) analyses to determine the genetic architecture that defines leaf shape. In the first analysis, we compared several important ampelographic measurements as finite trait QTL. In the second and third analyses, we identified significant shape variation via principal components analysis (PCA) and using a multivariate least squares interval mapping (MLSIM) approach. In total, we identified 271 significant QTL across the three measures of leaf shape and identified specific QTL hotspots in the grape genome which appear to drive major aspects of grapevine leaf shape.

Selection signatures in candidate genes and QTL for reproductive traits in Nellore heifers

The identification of selection signature genes may help to detect genomic regions that underwent artificial selection and contributed to phenotypic diversity. The aim of this study, therefore, was to detect selection signatures in candidate genes and quantitative trait locus (QTL) for reproductive traits in a Nellore population being selected for sexual precocity. A total of 2035 Nellore heifers, sourced from breeding programs focused on sexual precocity, were used. Candidate genes and some specific QTL related to reproductive traits were chosen based on published literature and Animal QTL databases, respectively, for investigation whether these regions were affected by selection. Selection signature DNA sequences were detected in the selected regions using the extended haplotype homozygosity (EHH) and relative extended haplotype homozygosity (REHH) methods. From 22,241 single nucleotide polymorphisms (SNPs) located in the candidate genes and QTL, 17,312 SNPs generated 2756 haplotype blocks. A total of 7518 EHH tests were analyzed using haplotypes with a frequency of more than 25%, for which there were 39 tests that were significant for REHH (P<0.01). Selection signature DNA sequences were detected that contained several QTLs for important reproductive traits in cattle, suggesting that reproductive traits may have been affected by selection for sexual precocity in this population. Forty-six genes were located in the selection signature regions, whereas 24 genes participated in important biological processes or pathways that may underlie sexual precocity. These results indicate there are possible molecular mechanisms related to sexual precocity in the Nellore breed.

The Effect of Fhb1 and Fhb5 Quantitative Trait Loci in Hard Red Spring Wheat Does Not Depend on Fungicide Use for Managing Fusarium Head Blight in Wheat.

Fusarium head blight (FHB), caused by Fusarium graminearum, is one of the most damaging diseases that affect wheat in Canada. The disease is best managed by integrating host resistance and fungicides, mainly demethylation inhibitors. Research has shown that the effect of fungicides may be dependent on the level of resistance of the cultivar. However, whether the performance of genotypes carrying specific Sumai 3-derived major FHB quantitative trait loci is dependent on fungicide application has not been explored. In our study, the performance of near-isogenic lines (NILs; <1.0% genome/alleles from the resistance donor), carrying Fhb1 and Fhb5 in a hard red spring wheat cultivar CDC Go background compared with a moderately susceptible (MS) genotype, was evaluated with and without one application of metconazole during full flowering. Field experiments were conducted at five site-years in Saskatchewan, Canada, between 2016 and 2017. In both the individual and combined analysis (all trials), we found that the effect of NILs and metconazole in suppressing FHB symptoms and deoxynivalenol (DON) accumulation in the grain was additive. FHB severity was generally low and fungicide efficacy levels, relative to the untreated control, were increased in the MS cultivar than in the NILs carrying Fhb1 and Fhb5, which were least affected by the disease. The results confirm the importance of integrating fungicides with cultivar resistance to reduce FHB and DON, regardless of the presence of those well-characterized resistant genes.

An integrated whole genome analysis of Mycobacterium tuberculosis reveals insights into relationship between its genome, transcriptome and methylome

Human tuberculosis disease (TB), caused by Mycobacterium tuberculosis (Mtb), is a complex disease, with a spectrum of outcomes. Genomic, transcriptomic and methylation studies have revealed differences between Mtb lineages, likely to impact on transmission, virulence and drug resistance. However, so far no studies have integrated sequence-based genomic, transcriptomic and methylation characterisation across a common set of samples, which is critical to understand how DNA sequence and methylation affect RNA expression and, ultimately, Mtb pathogenesis. Here we perform such an integrated analysis across 22 M. tuberculosis clinical isolates, representing ancient (lineage 1) and modern (lineages 2 and 4) strains. The results confirm the presence of lineage-specific differential gene expression, linked to specific SNP-based expression quantitative trait loci: with 10 eQTLs involving SNPs in promoter regions or transcriptional start sites; and 12 involving potential functional impairment of transcriptional regulators. Methylation status was also found to have a role in transcription, with evidence of differential expression in 50 genes across lineage 4 samples. Lack of methylation was associated with three novel variants in mamA, likely to cause loss of function of this enzyme. Overall, our work shows the relationship of DNA sequence and methylation to RNA expression, and differences between ancient and modern lineages. Further studies are needed to verify the functional consequences of the identified mechanisms of gene expression regulation.

Assessment of Genetic Biodiversity of Some Traits Using SSR Markers in Rice (Oryza sativa L.)

A set of germplasm of ninety accessions of rice (Oryza sativa L.) were collected in Egypt to assess genetic biodiversity. Using SSR markers some specific quantitative trait loci associated with tolerance to abiotic stress, resistance to blast disease and enhance the grain yield were found. Phylogenic data was collected. Eight primers (RM 315, RM 318, RM 166, RM 302, RM 201, RM 234, RM 526 and RM 144) showed different polymorphism frequencies. The polymorphism was nearly 80 %. The size of detected fragments ranged from 105–325 bp. A total of 186 bands were scored from the amplification products with the ten SSR primers. Genetic diversity analyses were conducted on the basis of the scores with 176 unique bands. Phylogenic trees for the fifteen rice accessions from each group and subgroup were established according to the molecular data and based on ten SSRs. A marked genetic diversity was observed in the five groups and the origin genotypes from the five groups were (Sakha 101, IR 03N137, IR 83142-12, IR 87856-10-AJY-1-B, HHZ 12-Y4-DT1-Y2 and IR 1552), respectively which showed higher layers of diversity and hence can be used as donors for the effective conservation, utilization and providing favorable genes and valuable germplasm in rice breeding.

Alveolar bone healing in mice genetically selected in the maximum (AIRmax) or minimum (AIRmin) inflammatory reaction

The exact role of inflammatory immune response in bone healing process is still unclear, but the success of the alveolar bone healing process seems to be associated with a moderate and transitory inflammatory response, while insufficient or exacerbated responses seems to have a detrimental influence in the healing outcome. In this context, we performed a comparative analysis of mice strains genetically selected for maximum (AIRmax) or minimum (AIRmin) acute inflammatory response to address the influence of inflammation genes in alveolar bone healing outcome. Experimental groups comprised 8-week-old male or female AIRmax and AIRmin submitted to extraction of upper right incisor, and evaluated at 0, 3, 7, 14 and 21 days after upper incision extraction by micro-computed tomography (μCT), histomorphometry, birefringence, immunohistochemistry and molecular (PCRArray) analysis. Overall, the results demonstrate a similar successful bone healing outcome at the endpoint was evidenced in both AIRmin and AIRmax strains. The histormophometric analysis reveal a slight but significant decrease in blood clot and inflammatory cells density, as well a delay in the bone formation in AIRmax strain in the early times, associated with a decreased expression of BMP2, BMP4, BMP7, TGFb1, RUNX2, and ALP. The evaluation of inflammatory cells nature reveals increased GR1+ cells counts in AIRmax strain at 3d, associated with increased levels of neutrophil chemoattractants such as CXCL1 and CXCL2, and its receptor CXCR1, while F4/80+ cell prevails in AIRmin strain at 7d. Also, our results demonstrate a relative predominance of M2 macrophages in AIRmin strain, associated with an increased expression of ARG1, IL10, TGFb, while M1 macrophages prevail in AIRmax, which parallel with increased IL-1B, IL-6 and TNF expression. At late repair stage, AIRmax presents evidences of increased bone remodeling, characterized by increased density of blood vessels and osteoclasts in parallel with decreased bone matrix density, as well increased levels of MMPs, osteoclastogenic and osteocyte markers. In the view of contrasting inflammatory and healing phenotypes of AIRmin and AIRmax strains in other models, the unpredicted phenotype observed suggests the existence of specific QTLs (Quantitative trait loci) responsible for the regulation ‘sterile’ inflammation and bone healing events. Despite the similar endpoint healing, AIRmax strain delayed repair was associated with increased presence of neutrophils and M1 macrophages, supporting the association of M2 cells with faster bone healing. Further studies are required to clarify the elements responsible for the regulation of inflammatory events at bone healing sites, as well the determinants of bone healing outcome.

Comparing allele specific expression and local expression quantitative trait loci and the influence of gene expression on complex trait variation in cattle

BackgroundThe mutations changing the expression level of a gene, or expression quantitative trait loci (eQTL), can be identified by testing the association between genetic variants and gene expression in multiple individuals (eQTL mapping), or by comparing the expression of the alleles in a heterozygous individual (allele specific expression or ASE analysis). The aims of the study were to find and compare ASE and local eQTL in 4 bovine RNA-sequencing (RNA-Seq) datasets, validate them in an independent ASE study and investigate if they are associated with complex trait variation.ResultsWe present a novel method for distinguishing between ASE driven by polymorphisms in cis and parent of origin effects. We found that single nucleotide polymorphisms (SNPs) driving ASE are also often local eQTL and therefore presumably cis eQTL. These SNPs often, but not always, affect gene expression in multiple tissues and, when they do, the allele increasing expression is usually the same. However, there were systematic differences between ASE and local eQTL and between tissues and breeds. We also found that SNPs significantly associated with gene expression (p < 0.001) were likely to influence some complex traits (p < 0.001), which means that some mutations influence variation in complex traits by changing the expression level of genes.ConclusionWe conclude that ASE detects phenomenon that overlap with local eQTL, but there are also systematic differences between the SNPs discovered by the two methods. Some mutations influencing complex traits are actually eQTL and can be discovered using RNA-Seq including eQTL in the genes CAST, CAPN1, LCORL and LEPROTL1.

Stable, Environmental Specific and Novel QTL Identification as Well as Genetic Dissection of Fatty Acid Metabolism in Brassica napus.

Fatty acid (FA) composition is the typical quantitative trait in oil seed crops, of which study is not only closely related to oil content, but is also more critical for the quality improvement of seed oil. The double haploid (DH) population named KN with a high density SNP linkage map was applied for quantitative trait loci (QTL) analysis of FA composition in this study. A total of 406 identified QTL were detected for eight FA components with an average confidence interval (CI) of 2.92 cM, the explained phenotypic variation (PV) value ranged from 1.49 to 45.05%. Totally, 204 consensus and 91 unique QTL were further obtained via meta-analysis method for the purpose of detecting multiple environment expressed and pleiotropic QTL, respectively. Of which, 74 stable expressed and 22 environmental specific QTL were also revealed, respectively. In order to make clear the genetic mechanism of FA metabolism at individual QTL level, conditional QTL analysis was also conducted and more than two thousand conditional QTL which could not be detected under the unconditional mapping were detected, which indicated the complex interrelationship of the QTL controlling FA content in rapeseed. Through comparative genomic analysis and homologous gene annotation, 61 candidates related to acyl lipid metabolism were identified underlying the CI of FA QTL. To further visualize the genetic mechanism of FA metabolism, an intuitive and meticulous network about acyl lipid metabolism was constructed and some closely related candidates were positioned. This study provided a more accurate localization for stable and pleiotropic QTL, and a deeper dissection of the molecular regulatory mechanism of FA metabolism in rapeseed.

Genome Reduction in Tetraploid Potato Reveals Genetic Load, Haplotype Variation, and Loci Associated With Agronomic Traits.

The cultivated potato (Solanum tuberosum) has a complex genetic structure due to its autotetraploidy and vegetative propagation which leads to accumulation of mutations and a highly heterozygous genome. A high degree of heterozygosity has been considered to be the main driver of fitness and agronomic trait performance in potato improvement efforts, which is negatively impacted by genetic load. To understand the genetic landscape of cultivated potato, we constructed a gynogenic dihaploid (2n = 2x = 24) population from cv. Superior, prior to development of a high-density genetic map containing 12,753 single nucleotide polymorphisms (SNPs). Common quantitative trait loci (QTL) were identified for tuber traits, vigor and height on chromosomes 2, 4, 7, and 10, while specific QTL for number of inflorescences per plant, and tuber shape were present on chromosomes 4, 6, 10, and 11. Simplex rather than duplex loci were mainly associated with traits. In general, the Q allele (main effect) detected in one or two homologous chromosomes was associated with lower mean trait values suggesting the importance of dosage allelic effects, and the presence of up to two undesired alleles in the QTL region. Loss of heterozygosity has been associated with a lower rate of fitness, yet no correlation between the percent heterozygosity and increased fitness or agronomic performance was observed. Based upon linkage phase, we reconstructed the four homologous chromosome haplotypes of cv. Superior. revealing heterogeneity throughout the genome yet nearly duplicate haplotypes occurring among the homologs of particular chromosomes. These results suggest that the potentially deleterious mutations associated with genetic load in tetraploid potato could be mitigated by multiple loci which is consistent with the theory that epistasis complicates the identification of associations between markers and phenotypic performance.

Genetic Architecture of Nitrogen-Deficiency Tolerance in Wheat Seedlings Based on a Nested Association Mapping (NAM) Population.

Genetic divergence for nitrogen utilization in germplasms is important in wheat breeding programs, especially for low nitrogen input management. In this study, a nested association mapping (NAM) population, derived from “Yanzhan 1” (a Chinese domesticated cultivar) crossed with “Hussar” (a British domesticated cultivar) and another three semi-wild wheat varieties, namely, “Cayazheda 29” (Triticum aestivum ssp. tibetanum Shao), “Yunnan” (T. aestivum ssp. yunnanense King), and “Yutian” (T. aestivum petropavloski Udats et Migusch), was used to detect quantitative trait loci (QTLs) for nitrogen utilization at the seedling stage. An integrated genetic map was constructed using 2,059 single nucleotide polymorphism (SNP) markers from a 90 K SNP chip, with a total coverage of 2,355.75 cM and an average marker spacing of 1.13 cM. A total of 67 QTLs for RDW (root dry weight), SDW (shoot dry weight), TDW (total dry weight), and RSDW (root to shoot ratio) were identified under normal nitrogen conditions (N+) and nitrogen deficient conditions (N−). Twenty-three of these QTLs were only detected under N− conditions. Moreover, 23 favorable QTLs were identified in the domesticated cultivar Yanzhan 1, 15 of which were detected under N+ conditions, while only four were detected under N− conditions. In contrast, the semi-wild cultivars contributed more favorable N−−specific QTLs (eight from Cayazheda 29; nine from Yunnan), which could be further explored for breeding cultivars adapted to nitrogen-deficient conditions. In particular, QRSDW-5A.1 from YN should be further evaluated using high-resolution mapping.

Statistical Approach for Gene Set Analysis with Trait Specific Quantitative Trait Loci

The analysis of gene sets is usually carried out based on gene ontology terms and known biological pathways. These approaches may not establish any formal relation between genotype and trait specific phenotype. In plant biology and breeding, analysis of gene sets with trait specific Quantitative Trait Loci (QTL) data are considered as great source for biological knowledge discovery. Therefore, we proposed an innovative statistical approach called Gene Set Analysis with QTLs (GSAQ) for interpreting gene expression data in context of gene sets with traits. The utility of GSAQ was studied on five different complex abiotic and biotic stress scenarios in rice, which yields specific trait/stress enriched gene sets. Further, the GSAQ approach was more innovative and effective in performing gene set analysis with underlying QTLs and identifying QTL candidate genes than the existing approach. The GSAQ approach also provided two potential biological relevant criteria for performance analysis of gene selection methods. Based on this proposed approach, an R package, i.e., GSAQ (https://cran.r-project.org/web/packages/GSAQ) has been developed. The GSAQ approach provides a valuable platform for integrating the gene expression data with genetically rich QTL data.

Chromosome conformation signatures define predictive markers of inadequate response to methotrexate in early rheumatoid arthritis

BackgroundThere is a pressing need in rheumatoid arthritis (RA) to identify patients who will not respond to first-line disease-modifying anti-rheumatic drugs (DMARD). We explored whether differences in genomic architecture represented by a chromosome conformation signature (CCS) in blood taken from early RA patients before methotrexate (MTX) treatment could assist in identifying non-response to DMARD and, whether there is an association between such a signature and RA specific expression quantitative trait loci (eQTL).MethodsWe looked for the presence of a CCS in blood from early RA patients commencing MTX using chromosome conformation capture by EpiSwitch™. Using blood samples from MTX responders, non-responders and healthy controls, a custom designed biomarker discovery array was refined to a 5-marker CCS that could discriminate between responders and non-responders to MTX. We cross-validated the predictive power of the CCS by generating 150 randomized groups of 59 early RA patients (30 responders and 29 non-responders) before MTX treatment. The CCS was validated using a blinded, independent cohort of 19 early RA patients (9 responders and 10 non-responders). Last, the loci of the CCS markers were mapped to RA-specific eQTL.ResultsWe identified a 5-marker CCS that could identify, at baseline, responders and non-responders to MTX. The CCS consisted of binary chromosome conformations in the genomic regions of IFNAR1, IL-21R, IL-23, CXCL13 and IL-17A. When tested on a cohort of 59 RA patients, the CCS provided a negative predictive value of 90.0% for MTX response. When tested on a blinded independent validation cohort of 19 early RA patients, the signature demonstrated a true negative response rate of 86 and a 90% sensitivity for detection of non-responders to MTX. Only conformations in responders mapped to RA-specific eQTL.ConclusionsHere we demonstrate that detection of a CCS in blood in early RA is able to predict inadequate response to MTX with a high degree of accuracy. Our results provide a proof of principle that a priori stratification of response to MTX is possible, offering a mechanism to provide alternative treatments for non-responders to MTX earlier in the course of the disease.

The phosphatidylinositide 3-kinase (PI3K) signaling pathway is a determinant of zileuton response in adults with asthma

Variable responsiveness to zileuton, a leukotriene antagonist used to treat asthma, may be due in part to genetic variation. While individual SNPs were previously associated with zileuton-related lung function changes, specific quantitative trait loci (QTLs) and biological pathways that may contribute have not been identified. In this study, we investigated the hypothesis that genetic variation within biological pathways is associated with zileuton response. We performed an integrative QTL mapping and pathway enrichment study to investigate data from a GWAS of zileuton response, in addition to mRNA expression profiles and leukotriene production data from lymphoblastoid cell lines (LCLs) (derived from asthmatics) that were treated with zileuton or ethanol (control). We identified 1060 QTLs jointly associated with zileuton-related differential LTB4 production in LCLs and lung function change in patients taking zileuton, of which eight QTLs were also significantly associated with persistent LTB4 production in LCLs following zileuton treatment (i.e., ‘poor’ responders). Four nominally significant trans-eQTLs were predicted to regulate three candidate genes (SELL, MTF2, and GAL), the expression of which was significantly reduced in LCLs following zileuton treatment. Gene and pathway enrichment analyses of QTL associations identified multiple genes and pathways, predominantly related to phosphatidyl inositol signaling via PI3K. We validated the PI3K pathway activation status in a subset of LCLs demonstrating variable zileuton-related LTB4 production, and show that in contrast to LCLs that responded to zileuton, the PI3K pathway was activated in poor responder LCLs. Collectively, these findings demonstrate a role for the PIK3 pathway and its targets as important determinants of differential responsiveness to zileuton.

Impact of regulatory variation across human iPSCs and differentiated cells.

Induced pluripotent stem cells (iPSCs) are an essential tool for studying cellular differentiation and cell types that are otherwise difficult to access. We investigated the use of iPSCs and iPSC-derived cells to study the impact of genetic variation on gene regulation across different cell types and as models for studies of complex disease. To do so, we established a panel of iPSCs from 58 well-studied Yoruba lymphoblastoid cell lines (LCLs); 14 of these lines were further differentiated into cardiomyocytes. We characterized regulatory variation across individuals and cell types by measuring gene expression levels, chromatin accessibility, and DNA methylation. Our analysis focused on a comparison of inter-individual regulatory variation across cell types. While most cell-type–specific regulatory quantitative trait loci (QTLs) lie in chromatin that is open only in the affected cell types, we found that 20% of cell-type–specific regulatory QTLs are in shared open chromatin. This observation motivated us to develop a deep neural network to predict open chromatin regions from DNA sequence alone. Using this approach, we were able to use the sequences of segregating haplotypes to predict the effects of common SNPs on cell-type–specific chromatin accessibility.

Haploinsufficiency of the Transcription Factor Ets-1 Is Renoprotective in Dahl Salt-Sensitive Rats.

Studies using Dahl salt-sensitive (SS) rats identified specific quantitative trait loci that predispose animals to hypertension-associated albuminuria and kidney injury. We explored the hypothesis that kidney-specific expression of the transcription factor Ets-1, located within one of these loci on chromosome 8, mediates glomerular injury in SS hypertension. During the first week on a high-salt diet, SS rats and SS rats with only one functioning Ets-1 gene (ES rats) demonstrated similar increases in BP. However, serum creatinine concentration, albuminuria, and glomerular expression of ETS-1 and two ETS-1 targets, MCP-1 and MMP2, did not increase as substantially in ES rats as in SS rats. Mean BP subsequently increased further in SS rats and remained higher than that of ES rats for the rest of the study. After 4 weeks of high-salt intake, ES rats still showed a lower mean serum creatinine concentration and less albuminuria, as well as less histologic evidence of glomerular injury and kidney fibrosis, than SS rats did. To investigate the specific contribution of renal Ets-1, we transplanted kidneys from ES or SS rats into salt-resistant SS-Chr 13BN/McwiCrl (SS-13BN) rats. Within 10 days on a high-salt diet, BP increased similarly in ES and SS allograft recipients, becoming significantly higher than the BP of control isograft recipients. However, mean serum creatinine concentration and albuminuria remained lower in ES allograft recipients than in SS allograft recipients at 2 weeks, and ES allografts showed less glomerular injury and interstitial fibrosis. In conclusion, reduced renal expression of ETS-1 prevented hypertension-associated kidney injury in SS rats.

Registration of Two Soybean Germplasm Lines Containing Leaf‐Chewing Insect Resistance QTLs from PI 229358 and PI 227687 Introgressed into ‘Benning’

Two soybean [Glycine max (L.) Merr.] germplasm lines, Benning‐ME (Reg. No. GP‐405, PI 679958) and Benning‐MGHE (Reg. No. GP‐406, PI 679959), were developed by the University of Georgia Agricultural Experiment Stations. Control of insect pests is crucial in soybean production; host‐plant resistance reduces the need for insecticide applications, thus diminishing production costs and pesticide concerns. In soybean, resistance to a broad range of leaf‐chewing insects is found in the Japanese accessions PI 229358 and PI 227687. In PI 229358, resistance is conferred by quantitative trait loci (QTLs) M, H, and G. In PI 227687, resistance is conferred by QTL‐E. To enhance soybean resistance to leaf‐chewing insects, QTLs of PI 229358 and PI 227687 were pyramided in Benning‐ME and Benning‐MGHE, which are near‐isogenic lines of ‘Benning’ obtained through marker‐assisted backcrossing. Under field conditions, Benning‐ME and Benning‐MGHE sustain 67 and 57% less defoliation than Benning, respectively. To determine the QTL introgressions in each line, high‐density single‐nucleotide polymorphism (SNP) genotypes were obtained using the SoySNP50K iSelect BeadChip (Illumina). To facilitate selection of lines carrying a specific QTL pyramid, Kompetitive Allele Specific Polymerase Chain Reaction markers were developed for high‐throughput genotyping. These lines are valuable genetic resources for breeding of host‐plant resistance to insects in soybean. The combination of QTL‐M and QTL‐E provides agriculturally relevant levels of resistance, and with only two loci, the use of this pyramid is feasible in a breeding program.