Genome-wide DNA methylation dynamics following recent polyploidy in the allotetraploid Tragopogon miscellus (Asteraceae).

Polyploidy is an important evolutionary force, yet epigenetic mechanisms, such as DNA methylation, that regulate genome-wide expression of duplicated genes remain largely unknown. Here, we use Tragopogon (Asteraceae) as a model system to discover patterns and temporal dynamics of DNA methylation in recently formed polyploids. The naturally occurring allotetraploid Tragopogon miscellus formed in the last 95-100 yr from parental diploids Tragopogon dubius and T. pratensis. We profiled the DNA methylomes of these three species using whole-genome bisulfite sequencing. Genome-wide methylation levels in T. miscellus were intermediate between its diploid parents. However, nonadditive CG and CHG methylation occurred in transposable elements (TEs), with variation among TE types. Most differentially methylated regions (DMRs) showed parental legacy, but some novel DMRs were detected in the polyploid. Differentially methylated genes (DMGs) were also identified and characterized. This study provides the first assessment of both overall and locus-specific patterns of DNA methylation in a recent natural allopolyploid and shows that novel methylation variants can be generated rapidly after polyploid formation. Together, these results demonstrate that mechanisms to regulate duplicate gene expression may arise soon after allopolyploid formation and that these mechanisms vary among genes.

Identification of pleiotropic loci mediating structural and non-structural carbohydrate accumulation within the sorghum bioenergy association panel using high-throughput markers.

Molecular characterization of diverse germplasm can contribute to breeding programs by increasing genetic gain for sorghum [Sorghum bicolor (L.) Moench] improvement. Identifying novel marker-trait associations and candidate genes enriches the existing genomic resources and can improve bioenergy-related traits using genomic-assisted breeding. In the current scenario, identifying the genetic loci underlying biomass and carbon partitioning is vital for ongoing efforts to maximize each carbon sink's yield for bioenergy production. Here, we have processed a high-density genomic marker (22 466 550) data based on whole-genome sequencing (WGS) using a set of 365 accessions from the bioenergy association panel (BAP), which includes ~19.7 million (19 744 726) single nucleotide polymorphism (SNPs) and 2.7 million (~2 721 824) insertion deletions (indels). A set of high-quality filtered SNP (~5.48 million) derived markers facilitated the assessment of population structure, genetic diversity, and genome-wide association studies (GWAS) for various traits related to biomass and its composition using the BAP. The phenotypic traits for GWAS included seed color (SC), plant height (PH), days to harvest (DTH), fresh weight (FW), dry weight (DW), brix content % (BRX), neutral detergent fiber (NDF), acid detergent fiber (ADF), non-fibrous carbohydrate (NFC), and lignin content. Several novel loci and candidate genes were identified for bioenergy-related traits, and some well-characterized genes for plant height (Dw1 and Dw2) and the YELLOW SEED1 locus (Y1) were validated. We further performed a multi-variate adaptive shrinkage analysis to identify pleiotropic QTL, which resulted in several shared marker-trait associations among bioenergy and compositional traits. Significant marker-trait associations with pleiotropic effects can be used to develop molecular markers for trait improvement using a marker-assisted breeding approach. Significant nucleotide diversity and heterozygosity were observed between photoperiod-sensitive and insensitive individuals of the panel. This diverse bioenergy panel with genomic resources will provide an excellent opportunity for further genetic studies, including selecting parental lines for superior hybrid development to improve biomass-related traits in sorghum.

Investigation of regulatory divergence between homoeologs in the recently formed allopolyploids, Tragopogon mirus and T. miscellus (Asteraceae).

Polyploidy is an important evolutionary process throughout eukaryotes, particularly in flowering plants. Duplicated gene pairs (homoeologs) in allopolyploids provide additional genetic resources for changes in molecular, biochemical, and physiological mechanisms that result in evolutionary novelty. Therefore, understanding how divergent genomes and their regulatory networks reconcile is vital for unraveling the role of polyploidy in plant evolution. Here, we compared the leaf transcriptomes of recently formed natural allotetraploids (Tragopogon mirus and T. miscellus) and their diploid parents (T. porrifolius X T. dubius and T. pratensis X T. dubius, respectively). Analysis of 35 400 expressed loci showed a significantly higher level of transcriptomic additivity compared to old polyploids; only 22% were non-additively expressed in the polyploids, with 5.9% exhibiting transgressive expression (lower or higher expression in the polyploids than in the diploid parents). Among approximately 7400 common orthologous regions (COREs), most loci in both allopolyploids exhibited expression patterns that were vertically inherited from their diploid parents. However, 18% and 20.3% of the loci showed novel expression bias patterns in T. mirus and T. miscellus, respectively. The expression changes of 1500 COREs were explained by cis-regulatory divergence (the condition in which the two parental subgenomes do not interact) between the diploid parents, whereas only about 423 and 461 of the gene expression changes represent trans-effects (the two parental subgenomes interact) in T. mirus and T. miscellus, respectively. The low degree of both non-additivity and trans-effects on gene expression may present the ongoing evolutionary processes of the newly formed Tragopogon polyploids (~80-90 years).

Discovering useful genetic variation in the seed parent gene pool for sorghum improvement.

Multi-parent populations contain valuable genetic material for dissecting complex, quantitative traits and provide a unique opportunity to capture multi-allelic variation compared to the biparental populations. A multi-parent advanced generation inter-cross (MAGIC) B-line (MBL) population composed of 708 F6 recombinant inbred lines (RILs), was recently developed from four diverse founders. These selected founders strategically represented the four most prevalent botanical races (kafir, guinea, durra, and caudatum) to capture a significant source of genetic variation to study the quantitative traits in grain sorghum [Sorghum bicolor (L.) Moench]. MBL was phenotyped at two field locations for seven yield-influencing traits: panicle type (PT), days to anthesis (DTA), plant height (PH), grain yield (GY), 1000-grain weight (TGW), tiller number per meter (TN) and yield per panicle (YPP). High phenotypic variation was observed for all the quantitative traits, with broad-sense heritabilities ranging from 0.34 (TN) to 0.84 (PH). The entire population was genotyped using Diversity Arrays Technology (DArTseq), and 8,800 single nucleotide polymorphisms (SNPs) were generated. A set of polymorphic, quality-filtered markers (3,751 SNPs) and phenotypic data were used for genome-wide association studies (GWAS). We identified 52 marker-trait associations (MTAs) for the seven traits using BLUPs generated from replicated plots in two locations. We also identified desirable allelic combinations based on the plant height loci (Dw1, Dw2, and Dw3), which influences yield related traits. Additionally, two novel MTAs were identified each on Chr1 and Chr7 for yield traits independent of dwarfing genes. We further performed a multi-variate adaptive shrinkage analysis and 15 MTAs with pleiotropic effect were identified. The five best performing MBL progenies were selected carrying desirable allelic combinations. Since the MBL population was designed to capture significant diversity for maintainer line (B-line) accessions, these progenies can serve as valuable resources to develop superior sorghum hybrids after validation of their general combining abilities via crossing with elite pollinators. Further, newly identified desirable allelic combinations can be used to enrich the maintainer germplasm lines through marker-assisted backcross breeding.

Fatty acid composition and genome-wide associations of a chickpea (Cicer arietinum L.) diversity panel for biofortification efforts.

Chickpea is a nutritionally dense pulse crop with high levels of protein, carbohydrates, micronutrients and low levels of fats. Chickpea fatty acids are associated with a reduced risk of obesity, blood cholesterol, and cardiovascular diseases in humans. We measured four primary chickpea fatty acids; palmitic acid (PA), linoleic acid (LA), alpha-linolenic acid (ALA), and oleic acid (OA), which are crucial for human health and plant stress responses in a chickpea diversity panel with 256 accessions (Kabuli and desi types). A wide concentration range was found for PA (450.7-912.6 mg/100 g), LA (1605.7-3459.9 mg/100 g), ALA (416.4-864.5 mg/100 g), and OA (1035.5-1907.2 mg/100 g). The percent recommended daily allowances also varied for PA (3.3-6.8%), LA (21.4-46.1%), ALA (34.7-72%), and OA (4.3-7.9%). Weak correlations were found among fatty acids. Genome-wide association studies (GWAS) were conducted using genotyping-by-sequencing data. Five significant single nucleotide polymorphisms (SNPs) were identified for PA. Admixture population structure analysis revealed seven subpopulations based on ancestral diversity in this panel. This is the first reported study to characterize fatty acid profiles across a chickpea diversity panel and perform GWAS to detect associations between genetic markers and concentrations of selected fatty acids. These findings demonstrate biofortification of chickpea fatty acids is possible using conventional and genomic breeding techniques, to develop superior cultivars with better fatty acid profiles for improved human health and plant stress responses.

Functional genomic effects of indels using Bayesian genome-phenome wide association studies in sorghum.

High-throughput genomic and phenomic data have enhanced the ability to detect genotype-to-phenotype associations that can resolve broad pleiotropic effects of mutations on plant phenotypes. As the scale of genotyping and phenotyping has advanced, rigorous methodologies have been developed to accommodate larger datasets and maintain statistical precision. However, determining the functional effects of associated genes/loci is expensive and limited due to the complexity associated with cloning and subsequent characterization. Here, we utilized phenomic imputation of a multi-year, multi-environment dataset using PHENIX which imputes missing data using kinship and correlated traits, and we screened insertions and deletions (InDels) from the recently whole-genome sequenced Sorghum Association Panel for putative loss-of-function effects. Candidate loci from genome-wide association results were screened for potential loss of function using a Bayesian Genome-Phenome Wide Association Study (BGPWAS) model across both functionally characterized and uncharacterized loci. Our approach is designed to facilitate in silico validation of associations beyond traditional candidate gene and literature-search approaches and to facilitate the identification of putative variants for functional analysis and reduce the incidence of false-positive candidates in current functional validation methods. Using this Bayesian GPWAS model, we identified associations for previously characterized genes with known loss-of-function alleles, specific genes falling within known quantitative trait loci, and genes without any previous genome-wide associations while additionally detecting putative pleiotropic effects. In particular, we were able to identify the major tannin haplotypes at the Tan1 locus and effects of InDels on the protein folding. Depending on the haplotype present, heterodimer formation with Tan2 was significantly affected. We also identified major effect InDels in Dw2 and Ma1, where proteins were truncated due to frameshift mutations that resulted in early stop codons. These truncated proteins also lost most of their functional domains, suggesting that these indels likely result in loss of function. Here, we show that the Bayesian GPWAS model is able to identify loss-of-function alleles that can have significant effects upon protein structure and folding as well as multimer formation. Our approach to characterize loss-of-function mutations and their functional repercussions will facilitate precision genomics and breeding by identifying key targets for gene editing and trait integration.

Development and characterization of a sorghum multi-parent advanced generation intercross (MAGIC) population for capturing diversity among seed parent gene pool.

Multiparent advanced eneration inter-cross (MAGIC) populations improve the precision of quantitative trait loci (QTL) mapping over biparental populations by incorporating increased diversity and opportunities to reduce linkage disequilibrium among variants. Here, we describe the development of a MAGIC B-Line (MBL) population from an inter-cross among 4 diverse founders of grain sorghum [Sorghum bicolor (L.) Moench] across different races (kafir, guinea, durra, and caudatum). These founders were selected based on genetic uniqueness and several distinct qualitative features including panicle architecture, plant color, seed color, endosperm texture, and awns. A whole set of MBL (708 F6) recombinant inbred lines along with their founders were genotyped using Diversity Arrays Technology (DArTseq) and 5,683 single-nucleotide polymorphisms (SNPs) were generated. A genetic linkage map was constructed using a set of polymorphic, quality-filtered markers (2,728 SNPs) for QTL interval-mapping. For population validation, 3 traits (seed color, plant color, and awns) were used for QTL mapping and genome-wide association study (GWAS). QTL mapping and GWAS identified 4 major genomic regions located across 3 chromosomes (Chr1, Chr3, and Chr6) that correspond to known genetic loci for the targeted traits. Founders of this population consist of the fertility maintainer (A/B line) gene pool and derived MBL lines could serve as female/seed parents in the cytoplasmic male sterility breeding system. The MBL population will serve as a unique genetic and genomic resource to better characterize the genetics of complex traits and potentially identify superior alleles for crop improvement efforts to enrich the seed parent gene pool.

A Robust Methodology for Assessing Homoeolog-Specific Expression.

Angiosperm evolution is marked by numerous, recurring polyploidization events. While hybridization and polyploidization have greatly increased the degree of genetic and phenotypic diversity in plants, the mechanisms underlying changes in the genotype-to-phenotype relationships remain unclear. As the field of natural sciences continues to expand during the post-genomic era, large datasets are becoming increasingly common. However, the development of tools and workflows available to robustly assess these changes have lagged behind data production. A robust homoeolog-specific expression analysis strongly depends upon proper homoeolog calling, the ability to account for reference sequence biases, flexible and accurate methods for dealing with residual bias, and a reproducible workflow. To that end, this chapter aims to provide a detailed description of the potential pitfalls encountered while estimating homoeolog-specific expression as well as provide a workflow that allows for robust inferences based on precise estimates of expression changes.

Dissecting the Genetic Architecture of Carbon Partitioning in Sorghum Using Multiscale Phenotypes.

Carbon partitioning in plants may be viewed as a dynamic process composed of the many interactions between sources and sinks. The accumulation and distribution of fixed carbon is not dictated simply by the sink strength and number but is dependent upon the source, pathways, and interactions of the system. As such, the study of carbon partitioning through perturbations to the system or through focus on individual traits may fail to produce actionable developments or a comprehensive understanding of the mechanisms underlying this complex process. Using the recently published sorghum carbon-partitioning panel, we collected both macroscale phenotypic characteristics such as plant height, above-ground biomass, and dry weight along with microscale compositional traits to deconvolute the carbon-partitioning pathways in this multipurpose crop. Multivariate analyses of traits resulted in the identification of numerous loci associated with several distinct carbon-partitioning traits, which putatively regulate sugar content, manganese homeostasis, and nitrate transportation. Using a multivariate adaptive shrinkage approach, we identified several loci associated with multiple traits suggesting that pleiotropic and/or interactive effects may positively influence multiple carbon-partitioning traits, or these overlaps may represent molecular switches mediating basal carbon allocating or partitioning networks. Conversely, we also identify a carbon tradeoff where reduced lignin content is associated with increased sugar content. The results presented here support previous studies demonstrating the convoluted nature of carbon partitioning in sorghum and emphasize the importance of taking a holistic approach to the study of carbon partitioning by utilizing multiscale phenotypes.

Sorghum Association Panel whole-genome sequencing establishes cornerstone resource for dissecting genomic diversity.

Association mapping panels represent foundational resources for understanding the genetic basis of phenotypic diversity and serve to advance plant breeding by exploring genetic variation across diverse accessions. We report the whole-genome sequencing (WGS) of 400 sorghum (Sorghum bicolor (L.) Moench) accessions from the Sorghum Association Panel (SAP) at an average coverage of 38× (25-72×), enabling the development of a high-density genomic marker set of 43 983 694 variants including single-nucleotide polymorphisms (approximately 38 million), insertions/deletions (indels) (approximately 5 million), and copy number variants (CNVs) (approximately 170 000). We observe slightly more deletions among indels and a much higher prevalence of deletions among CNVs compared to insertions. This new marker set enabled the identification of several novel putative genomic associations for plant height and tannin content, which were not identified when using previous lower-density marker sets. WGS identified and scored variants in 5-kb bins where available genotyping-by-sequencing (GBS) data captured no variants, with half of all bins in the genome falling into this category. The predictive ability of genomic best unbiased linear predictor (GBLUP) models was increased by an average of 30% by using WGS markers rather than GBS markers. We identified 18 selection peaks across subpopulations that formed due to evolutionary divergence during domestication, and we found six Fst peaks resulting from comparisons between converted lines and breeding lines within the SAP that were distinct from the peaks associated with historic selection. This population has served and continues to serve as a significant public resource for sorghum research and demonstrates the value of improving upon existing genomic resources.

Ten new high-quality genome assemblies for diverse bioenergy sorghum genotypes.

Introduction: Sorghum (Sorghum bicolor (L.) Moench) is an agriculturally and economically important staple crop that has immense potential as a bioenergy feedstock due to its relatively high productivity on marginal lands. To capitalize on and further improve sorghum as a potential source of sustainable biofuel, it is essential to understand the genomic mechanisms underlying complex traits related to yield, composition, and environmental adaptations. Methods: Expanding on a recently developed mapping population, we generated de novo genome assemblies for 10 parental genotypes from this population and identified a comprehensive set of over 24 thousand large structural variants (SVs) and over 10.5 million single nucleotide polymorphisms (SNPs). Results: We show that SVs and nonsynonymous SNPs are enriched in different gene categories, emphasizing the need for long read sequencing in crop species to identify novel variation. Furthermore, we highlight SVs and SNPs occurring in genes and pathways with known associations to critical bioenergy-related phenotypes and characterize the landscape of genetic differences between sweet and cellulosic genotypes. Discussion: These resources can be integrated into both ongoing and future mapping and trait discovery for sorghum and its myriad uses including food, feed, bioenergy, and increasingly as a carbon dioxide removal mechanism.

Genome-wide association studies of mineral and phytic acid concentrations in pea (Pisum sativum L.) to evaluate biofortification potential.

Pea (Pisum sativum L.) is an important cool season food legume for sustainable food production and human nutrition due to its nitrogen fixation capabilities and nutrient-dense seed. However, minimal breeding research has been conducted to improve the nutritional quality of the seed for biofortification, and most genomic-assisted breeding studies utilize small populations with few single nucleotide polymorphisms (SNPs). Genomic resources for pea have lagged behind those of other grain crops, but the recent release of the Pea Single Plant Plus Collection (PSPPC) and the pea reference genome provide new tools to study nutritional traits for biofortification. Calcium, phosphorus, potassium, iron, zinc, and phytic acid concentrations were measured in a study population of 299 different accessions grown under greenhouse conditions. Broad phenotypic variation was detected for all parameters except phytic acid. Calcium exhibited moderate broad-sense heritability (H2) estimates, at 50%, while all other minerals exhibited low heritability. Of the accessions used, 267 were previously genotyped in the PSPPC release by the USDA, and we mapped the genotyping data to the pea reference genome for the first time. This study generated 54,344 high-quality SNPs used to investigate the population structure of the PSPPC and perform a genome-wide association study to identify genomic loci associated with mineral concentrations in mature pea seed. Overall, we were able to identify multiple significant SNPs and candidate genes for iron, phosphorus, and zinc. These results can be used for genetic improvement in pea for nutritional traits and biofortification, and the candidate genes provide insight into mineral metabolism.

Trajectories of Homoeolog-Specific Expression in Allotetraploid Tragopogon castellanus Populations of Independent Origins.

Polyploidization can have a significant ecological and evolutionary impact by providing substantially more genetic material that may result in novel phenotypes upon which selection may act. While the effects of polyploidization are broadly reviewed across the plant tree of life, the reproducibility of these effects within naturally occurring, independently formed polyploids is poorly characterized. The flowering plant genus Tragopogon (Asteraceae) offers a rare glimpse into the intricacies of repeated allopolyploid formation with both nascent (< 90 years old) and more ancient (mesopolyploids) formations. Neo- and mesopolyploids in Tragopogon have formed repeatedly and have extant diploid progenitors that facilitate the comparison of genome evolution after polyploidization across a broad span of evolutionary time. Here, we examine four independently formed lineages of the mesopolyploid Tragopogon castellanus for homoeolog expression changes and fractionation after polyploidization. We show that expression changes are remarkably similar among these independently formed polyploid populations with large convergence among expressed loci, moderate convergence among loci lost, and stochastic silencing. We further compare and contrast these results for T. castellanus with two nascent Tragopogon allopolyploids. While homoeolog expression bias was balanced in both nascent polyploids and T. castellanus, the degree of additive expression was significantly different, with the mesopolyploid populations demonstrating more non-additive expression. We suggest that gene dosage and expression noise minimization may play a prominent role in regulating gene expression patterns immediately after allopolyploidization as well as deeper into time, and these patterns are conserved across independent polyploid lineages.

Genome-wide association mapping of lentil (Lens culinaris Medikus) prebiotic carbohydrates toward improved human health and crop stress tolerance.

Lentil, a cool-season food legume, is rich in protein and micronutrients with a range of prebiotic carbohydrates, such as raffinose-family oligosaccharides (RFOs), fructooligosaccharides (FOSs), sugar alcohols (SAs), and resistant starch (RS), which contribute to lentil's health benefits. Beneficial microorganisms ferment prebiotic carbohydrates in the colon, which impart health benefits to the consumer. In addition, these carbohydrates are vital to lentil plant health associated with carbon transport, storage, and abiotic stress tolerance. Thus, lentil prebiotic carbohydrates are a potential nutritional breeding target for increasing crop resilience to climate change with increased global nutritional security. This study phenotyped a total of 143 accessions for prebiotic carbohydrates. A genome-wide association study (GWAS) was then performed to identify associated variants and neighboring candidate genes. All carbohydrates analyzed had broad-sense heritability estimates (H2) ranging from 0.22 to 0.44, comparable to those reported in the literature. Concentration ranges corresponded to percent recommended daily allowances of 2-9% SAs, 7-31% RFOs, 51-111% RS, and 57-116% total prebiotic carbohydrates. Significant SNPs and associated genes were identified for numerous traits, including a galactosyltransferase (Lcu.2RBY.1g019390) known to aid in RFO synthesis. Further studies in multiple field locations are necessary. Yet, these findings suggest the potential for molecular-assisted breeding for prebiotic carbohydrates in lentil to support human health and crop resilience to increase global food security.

Islet sympathetic innervation and islet neuropathology in patients with type 1 diabetes.

Dysregulation of glucagon secretion in type 1 diabetes (T1D) involves hypersecretion during postprandial states, but insufficient secretion during hypoglycemia. The sympathetic nervous system regulates glucagon secretion. To investigate islet sympathetic innervation in T1D, sympathetic tyrosine hydroxylase (TH) axons were analyzed in control non-diabetic organ donors, non-diabetic islet autoantibody-positive individuals (AAb), and age-matched persons with T1D. Islet TH axon numbers and density were significantly decreased in AAb compared to T1D with no significant differences observed in exocrine TH axon volume or lengths between groups. TH axons were in close approximation to islet α-cells in T1D individuals with long-standing diabetes. Islet RNA-sequencing and qRT-PCR analyses identified significant alterations in noradrenalin degradation, α-adrenergic signaling, cardiac ß-adrenergic signaling, catecholamine biosynthesis, and additional neuropathology pathways. The close approximation of TH axons at islet α-cells supports a model for sympathetic efferent neurons directly regulating glucagon secretion. Sympathetic islet innervation and intrinsic adrenergic signaling pathways could be novel targets for improving glucagon secretion in T1D.

Genetic characterization of a Sorghum bicolor multiparent mapping population emphasizing carbon-partitioning dynamics.

Sorghum bicolor, a photosynthetically efficient C4 grass, represents an important source of grain, forage, fermentable sugars, and cellulosic fibers that can be utilized in myriad applications ranging from bioenergy to bioindustrial feedstocks. Sorghum's efficient fixation of carbon per unit time per unit area per unit input has led to its classification as a preferred biomass crop highlighted by its designation as an advanced biofuel by the U.S. Department of Energy. Due to its extensive genetic diversity and worldwide colonization, sorghum has considerable diversity for a range of phenotypes influencing productivity, composition, and sink/source dynamics. To dissect the genetic basis of these key traits, we present a sorghum carbon-partitioning nested association mapping (NAM) population generated by crossing 11 diverse founder lines with Grassl as the single recurrent female. By exploiting existing variation among cellulosic, forage, sweet, and grain sorghum carbon partitioning regimes, the sorghum carbon-partitioning NAM population will allow the identification of important biomass-associated traits, elucidate the genetic architecture underlying carbon partitioning and improve our understanding of the genetic determinants affecting unique phenotypes within Poaceae. We contrast this NAM population with an existing grain population generated using Tx430 as the recurrent female. Genotypic data are assessed for quality by examining variant density, nucleotide diversity, linkage decay, and are validated using pericarp and testa phenotypes to map known genes affecting these phenotypes. We release the 11-family NAM population along with corresponding genomic data for use in genetic, genomic, and agronomic studies with a focus on carbon-partitioning regimes.

Identification of Novel Genomic Associations and Gene Candidates for Grain Starch Content in Sorghum.

Starch accumulated in the endosperm of cereal grains as reserve energy for germination serves as a staple in human and animal nutrition. Unraveling genetic control for starch metabolism is important for breeding grains with high starch content. In this study, we used a sorghum association panel with 389 individuals and 141,557 single nucleotide polymorphisms (SNPs) to fit linear mixed models (LMM) for identifying genomic regions and potential candidate genes associated with starch content. Three associated genomic regions, one in chromosome (chr) 1 and two novel associations in chr-8, were identified using combination of LMM and Bayesian sparse LMM. All significant SNPs were located within protein coding genes, with SNPs ∼ 52 Mb of chr-8 encoding a Casperian strip membrane protein (CASP)-like protein (Sobic.008G111500) and a heat shock protein (HSP) 90 (Sobic.008G111600) that were highly expressed in reproductive tissues including within the embryo and endosperm. The HSP90 is a potential hub gene with gene network of 75 high-confidence first interactors that is enriched for five biochemical pathways including protein processing. The first interactors of HSP90 also showed high transcript abundance in reproductive tissues. The candidates of this study are likely involved in intricate metabolic pathways and represent candidate gene targets for source-sink activities and drought and heat stress tolerance during grain filling.

STAT3 imparts BRCAness by impairing homologous recombination repair in Epstein-Barr virus-transformed B lymphocytes.

Epstein-Barr virus (EBV) causes lymphomas and epithelial cell cancers. Though generally silent in B lymphocytes, this widely prevalent virus can cause endemic Burkitt lymphoma and post-transplant lymphoproliferative disorders/lymphomas in immunocompromised hosts. By learning how EBV breaches barriers to cell proliferation, we hope to undermine those strategies to treat EBV lymphomas and potentially other cancers. We had previously found that EBV, through activation of cellular STAT3 prevents phosphorylation of Chk1, and thereby, suppresses activation of the intra-S phase cell-cycle checkpoint, a potent barrier to oncogene-driven proliferation. This observation prompted us to examine the consequences on DNA repair since homologous recombination repair, the most error-free form, requires phosphoChk1. We now report that the defect in Chk1 phosphorylation also curtails RAD51 nucleation, and thereby, homologous recombination repair of DNA double strand breaks. The resulting reliance on error-prone microhomology-mediated end-joining (MMEJ) repair makes EBV-transformed cells susceptible to PARP inhibition and simultaneous accrual of genome-wide deletions and insertions resulting from synthesis-dependent MMEJ. Analysis of transcriptomic and drug susceptibility data from hundreds of cancer lines reveals a STAT3-dependent gene-set predictive of susceptibility of cancers to synthetic lethal PARP inhibition. These findings i) demonstrate how the tumor virus EBV re-shapes cellular DNA repair, ii) provide the first genome-wide evidence for insertions resulting from MMEJ in human cells, and iii) expand the range of cancers (EBV-related and -unrelated) that are likely to respond to synthetic lethal inhibitors given the high prevalence of cancers with constitutively active STAT3.

Jasmonate induced alternative splicing responses in Arabidopsis.

Jasmonate is an essential phytohormone regulating plant growth, development, and defense. Alternative splicing (AS) in jasmonate ZIM-domain (JAZ) repressors is well-characterized and plays an important role in jasmonate signaling regulation. However, it is unknown whether other genes in the jasmonate signaling pathway are regulated by AS. We explore the potential for AS regulation in three Arabidopsis genotypes (WT, jaz2, jaz7) in response to methyl jasmonate (MeJA) treatment with respect to: (a) differential AS, (b) differential miRNA targeted AS, and (c) AS isoforms with novel functions. AS events identified from transcriptomic data were validated with proteomic data. Protein interaction networks identified two genes, SKIP and ALY4 whose products have both DNA- and RNA-binding affinities, as potential key regulators mediating jasmonate signaling and AS regulation. We observed cases where AS alone, or AS and transcriptional regulation together, can influence gene expression in response to MeJA. Twenty-one genes contain predicted miRNA target sites subjected to AS, which implies that AS is coupled to miRNA regulation. We identified 30 cases where alternatively spliced isoforms may have novel functions. For example, AS of bHLH160 generates an isoform without a basic domain, which may convert it from an activator to a repressor. Our study identified potential key regulators in AS regulation of jasmonate signaling pathway. These findings highlight the importance of AS regulation in the jasmonate signaling pathway, both alone and in collaboration with other regulators. SIGNIFICANCE STATEMENT: By exploring alternative splicing, we demonstrate its regulation in the jasmonate signaling pathway alone or in collaboration with other posttranscriptional regulations such as nonsense and microRNA-mediated decay. A signal transduction network model for alternative splicing in jasmonate signaling pathway was generated, contributing to our understanding for this important, prevalent, but relatively unexplored regulatory mechanism in plants.

Employing Molecular Phylodynamic Methods to Identify and Forecast HIV Transmission Clusters in Public Health Settings: A Qualitative Study.

Molecular HIV surveillance is a promising public health strategy for curbing the HIV epidemic. Clustering technologies used by health departments to date are limited in their ability to infer/forecast cluster growth trajectories. Resolution of the spatiotemporal dynamics of clusters, through phylodynamic and phylogeographic modelling, is one potential strategy to develop a forecasting tool; however, the projected utility of this approach needs assessment. Prior to incorporating novel phylodynamic-based molecular surveillance tools, we sought to identify possible issues related to their feasibility, acceptability, interpretation, and utility. Qualitative data were collected via focus groups among field experts (n = 17, 52.9% female) using semi-structured, open-ended questions. Data were coded using an iterative process, first through the development of provisional themes and subthemes, followed by independent line-by-line coding by two coders. Most participants routinely used molecular methods for HIV surveillance. All agreed that linking molecular sequences to epidemiological data is important for improving HIV surveillance. We found that, in addition to methodological challenges, a variety of implementation barriers are expected in relation to the uptake of phylodynamic methods for HIV surveillance. The participants identified several opportunities to enhance current methods, as well as increase the usability and utility of promising works-in-progress.