Photorespiration is an essential process in oxygenic photosynthetic organisms triggered by the oxygenase activity of Rubisco. In peroxisomes, photorespiratory HYDROXYPYRUVATE REDUCTASE1 (HPR1) catalyzes the conversion of hydroxypyruvate to glycerate together with the oxidation of a pyridine nucleotide cofactor. HPR1 regulation remains poorly understood; however, HPR1 phosphorylation at T335 has been reported. By comparing the kinetic properties of phosphomimetic (T335D), nonphosphorylatable (T335A), and wild-type recombinant Arabidopsis (Arabidopsis thaliana) HPR1, it was found that HPR1-T335D exhibits reduced NADH-dependent hydroxypyruvate reductase activity while showing improved NADPH-dependent activity. Complementation of the Arabidopsis hpr1-1 mutant by either wild-type HPR1 or HPR1-T335A fully complemented the photorespiratory growth phenotype of hpr1-1 in ambient air, whereas HPR1-T335D-containing hpr1-1 plants remained smaller and had lower photosynthetic CO₂ assimilation rates. Metabolite analyses indicated that these phenotypes were associated with subtle perturbations in the photorespiratory cycle of HPR1-T335D-complemented hpr1-1 rosettes compared to all other HPR1-containing lines. Therefore, T335 phosphorylation may play a role in the regulation of HPR1 activity in planta, although it was not required for growth under ambient air controlled conditions. Furthermore, improved NADP-dependent HPR1 activities in peroxisomes could not compensate for the reduced NADH-dependent HPR1 activity.

Members of the light-harvesting complex protein family participate in multiple processes connected with light sensing, light absorption, and pigment binding within the thylakoid membrane. Amino acid residues of the light-harvesting chlorophyll a/b-binding proteins involved in pigment binding have been precisely identified through x-ray crystallography experiments. In vitro pigment-binding studies have been performed with LIGHT-HARVESTING-LIKE3 proteins, and the pigment-binding ability of cyanobacterial high-light-inducible proteins has been studied in detail. However, analysis of pigment binding by plant high-light-inducible protein homologs, called ONE-HELIX PROTEINS (OHPs), is lacking. Here, we report on successful in vitro reconstitution of Arabidopsis (Arabidopsis thaliana) OHPs with chlorophylls and carotenoids and show that pigment binding depends on the formation of OHP1/OHP2 heterodimers. Pigment-binding capacity was completely lost in each of the OHPs when residues of the light-harvesting complex chlorophyll-binding motif required for chlorophyll binding were mutated. Moreover, the mutated OHP variants failed to rescue the respective knockout (T-DNA insertion) mutants, indicating that pigment-binding ability is essential for OHP function in vivo. The scaffold protein HIGH CHLOROPHYLL FLUORESCENCE244 (HCF244) is tethered to the thylakoid membrane by the OHP heterodimer. We show that HCF244 stability depends on OHP heterodimer formation and introduce the concept of a functional unit consisting of OHP1, OHP2, and HCF244, in which each protein requires the others. Because of their pigment-binding capacity, we suggest that OHPs function in the delivery of pigments to the D1 subunit of PSII.

The corrinoid B₁₂ is synthesized only by prokaryotes yet is widely required by eukaryotes as an enzyme cofactor. Microalgae have evolved B₁₂ dependence on multiple occasions, and we previously demonstrated that experimental evolution of the non–B₁₂-requiring alga Chlamydomonas reinhardtii in media supplemented with B₁₂ generated a B₁₂-dependent mutant (hereafter metE7). This clone provides a unique opportunity to study the physiology of a nascent B₁₂ auxotroph. Our analyses demonstrate that B₁₂ deprivation of metE7 disrupts C1 metabolism, causes an accumulation of starch and triacylglycerides, and leads to a decrease in photosynthetic pigments, proteins, and free amino acids. B₁₂ deprivation also caused a substantial increase in reactive oxygen species, which preceded rapid cell death. Survival could be improved without compromising growth by simultaneously depriving the cells of nitrogen, suggesting a type of cross protection. Significantly, we found further improvements in survival under B₁₂ limitation and an increase in B₁₂ use efficiency after metE7 underwent a further period of experimental evolution, this time in coculture with a B₁₂-producing bacterium. Therefore, although an early B₁₂-dependent alga would likely be poorly adapted to coping with B₁₂ deprivation, association with B₁₂-producers can ensure long-term survival whilst also providing a suitable environment for evolving mechanisms to tolerate B₁₂ limitation better.

Phosphatidylcholine and phosphatidylethanolamine are two major phospholipid classes in eukaryotes. Each biosynthesis pathway starts with the phosphorylation of choline (Cho) or ethanolamine (Etn) catalyzed by either choline or ethanolamine kinase (CEK). Arabidopsis contains four CEK isoforms, but their isozyme-specific roles in metabolism and development are poorly described. Here, we showed that these four CEKs have distinct substrate specificities in vitro. While CEK1 and CEK2 showed substrate preference for Cho over Etn, CEK3 and CEK4 had clear substrate specificity for Cho and Etn, respectively. In vivo, CEK1, CEK2, and CEK3 exhibited kinase activity for Cho but not Etn, although the latter two isoforms showed rather minor contributions to total Cho kinase activity in both shoots and roots. The knockout mutants of CEK2 and CEK3 both affected root growth, and these isoforms had nonoverlapping cell-type-specific expression patterns in the root meristematic zone. In-depth phenotype analysis, as well as chemical and genetic complementation, revealed that CEK3, a Cho-specific kinase, is involved in cell elongation during root development. Phylogenetic analysis of CEK orthologs in Brassicaceae species showed evolutionary divergence between Etn kinases and Cho kinases. Collectively, our results demonstrate the distinct roles of the four CEK isoforms in Cho/Etn metabolism and plant development.

Salicinoids form a specific class of phenolic glycosides characteristic of the Salicaceae. Although salicinoids accumulate in large amounts and have been shown to be involved in plant defense, their biosynthesis is unclear. We identified two sulfated salicinoids, salicin-7-sulfate and salirepin-7-sulfate, in black cottonwood (Populus trichocarpa). Both compounds accumulated in high amounts in above-ground tissues including leaves, petioles, and stems, but were also found at lower concentrations in roots. A survey of salicin-7-sulfate and salirepin-7-sulfate in a subset of poplar (Populus sp.) and willow (Salix sp.) species revealed a broader distribution within the Salicaceae. To elucidate the formation of these compounds, we studied the sulfotransferase (SOT) gene family in P. trichocarpa (PtSOT). One of the identified genes, PtSOT1, was shown to encode an enzyme able to convert salicin and salirepin into salicin-7-sulfate and salirepin-7-sulfate, respectively. The expression of PtSOT1 in different organs of P. trichocarpa matched the accumulation of sulfated salicinoids in planta. Moreover, RNA interference-mediated knockdown of SOT1 in gray poplar (Populus x canescens) resulted in decreased levels of sulfated salicinoids in comparison to wild-type plants, indicating that SOT1 is responsible for their formation in planta. The presence of a nonfunctional SOT1 allele in black poplar (Populus nigra) was shown to correlate with the absence of salicin-7-sulfate and salirepin-7-sulfate in this species. Food choice experiments with leaves from wild-type and SOT1 knockdown trees suggest that sulfated salicinoids do not affect the feeding preference of the generalist caterpillar Lymantria dispar. A potential role of the sulfated salicinoids in sulfur storage and homeostasis is discussed.

The lignin biosynthetic pathway is highly conserved in angiosperms, yet pathway manipulations give rise to a variety of taxon-specific outcomes. Knockout of lignin-associated 4-coumarate:CoA ligases (4CLs) in herbaceous species mainly reduces guaiacyl (G) lignin and enhances cell wall saccharification. Here we show that CRISPR-knockout of 4CL1 in poplar (Populus tremula x alba) preferentially reduced syringyl (S) lignin, with negligible effects on biomass recalcitrance. Concordant with reduced S-lignin was downregulation of ferulate 5-hydroxylases (F5Hs). Lignification was largely sustained by 4CL5, a low-affinity paralog of 4CL1 typically with only minor xylem expression or activity. Levels of caffeate, the preferred substrate of 4CL5, increased in line with significant upregulation of caffeoyl shikimate esterase1. Upregulation of caffeoyl-CoA O-methyltransferase1 and downregulation of F5Hs are consistent with preferential funneling of 4CL5 products toward G-lignin biosynthesis at the expense of S-lignin. Thus, transcriptional and metabolic adaptations to 4CL1-knockout appear to have enabled 4CL5 catalysis at a level sufficient to sustain lignification. Finally, genes involved in sulfur assimilation, the glutathione-ascorbate cycle, and various antioxidant systems were upregulated in the mutants, suggesting cascading responses to perturbed thioesterification in lignin biosynthesis.

Plants require a high concentration of ascorbate as a redox buffer for survival under stress conditions, such as high light. Dehydroascorbate reductases (DHARs) are enzymes that catalyze the reduction of DHA to ascorbate using reduced glutathione (GSH) as an electron donor, allowing rapid ascorbate recycling. However, a recent study using an Arabidopsis (Arabidopsis thaliana) triple mutant lacking all three DHAR genes (herein called dhar) did not find evidence for their role in ascorbate recycling under oxidative stress. To further study the function of DHARs, we generated dhar Arabidopsis plants as well as a quadruple mutant line combining dhar with an additional vtc2 mutation that causes ascorbate deficiency. Measurements of ascorbate in these mutants under low- or high-light conditions indicated that DHARs have a nonnegligible impact on full ascorbate accumulation under high light, but that they are dispensable when ascorbate concentrations are low to moderate. Because GSH itself can reduce DHA nonenzymatically, we used the pad2 mutant that contains ~30% of the wild-type GSH level. The pad2 mutant accumulated ascorbate at a wild-type level under high light; however, when the pad2 mutation was combined with dhar, there was near-complete inhibition of high-light–dependent ascorbate accumulation. The lack of ascorbate accumulation was consistent with a marked increase in the ascorbate degradation product threonate. These findings indicate that ascorbate recycling capacity is limited in dhar pad2 plants, and that both DHAR activity and GSH content set a threshold for high-light–induced ascorbate accumulation.

The polysaccharide pectin is a major component of the plant cell wall. The pectic glycan homogalacturonan (HG) is a proportionally small but important component of a specialized seed cell wall called mucilage. HG is synthesized in a highly methylesterified form, and, following secretion, is de-methylesterified by pectin methylesterases (PMEs). The degree of methylesterification of HG determines the structural and functional properties of pectin, but how methylesterification is regulated remains largely unknown. Here, we identified two BEL1-Like homeodomain (BLH) transcription factors, BLH2 and BLH4, as positive regulators of HG de-methylesterification in Arabidopsis (Arabidopsis thaliana) seed coat mucilage. BLH2 and BLH4 were significantly expressed in mucilage secretory cells during seed mucilage production. BLH2 and BLH4 single mutants exhibited no obvious mucilage phenotype, but the blh2 blh4 double mutant displayed significantly reduced mucilage adherence to the seed. Reduced mucilage adherence in blh2 blh4 was caused by decreased PME activity in the seed coat, which increased the degree of methylesterification of HG in mucilage. The expression of several PME metabolism-related genes, including PME58, PECTIN METHYLESTERASE INHIBITOR6, SEEDSTICK, and MYB52 was significantly altered in blh2 blh4 seeds. BLH2 and BLH4 directly activated PME58 expression by binding to its TGACAGGT cis-element. Moreover, pme58 mutants exhibited reduced mucilage adherence similar to that of blh2 blh4, and the blh2 blh4 pme58 triple mutant exhibited no additional mucilage adherence defects. Furthermore, overexpression of PME58 in blh2 blh4 rescued the mucilage adherence defect. Together, these results demonstrate that BLH2 and BLH4 redundantly regulate de-methylesterification of HG in seed mucilage by directly activating PME58.

RIPENING INHIBITOR (RIN) is a transcription factor with transcriptional activator activity that plays a major role in regulating fruit ripening in tomato (Solanum lycopersicum). Recent studies have revealed that (1) RIN is indispensable for full ripening but not for the induction of ripening; and (2) the rin mutation, which produces nonripening fruits that never turn red or soften, is not a null mutation but instead converts the encoded transcriptional activator into a repressor. Here, we have uncovered aspects of RIN function by characterizing a series of allelic mutations within this locus that were produced by CRISPR/Cas9. Fruits of RIN-knockout plants, which are characterized by partial ripening and low levels of lycopene but never turn fully red, showed excess flesh softening compared to the wild type. The knockout mutant fruits also showed accelerated cell wall degradation, suggesting that, contrary to the conventional view, RIN represses over-ripening in addition to facilitating ripening. A C-terminal domain-truncated RIN protein, encoded by another allele of the RIN locus (rinG2), did not activate transcription but formed transcription factor complexes that bound to target genomic regions in a manner similar to that observed for wild-type RIN protein. Fruits expressing this truncated RIN protein exhibited extended shelf life, but unlike rin fruits, they accumulated lycopene and appeared orange. The diverse ripening properties of the RIN allelic mutants suggest that substantial phenotypic variation can be produced by tuning the activity of a transcription factor.

Despite the ecological relevance of diatoms, many aspects of their photosynthetic machinery remain poorly understood. Diatoms differ from the green lineage of oxygenic organisms by their photosynthetic pigments and light-harvesting complex (Lhc) proteins, the latter of which are also called fucoxanthin-chlorophyll proteins (FCP). These are composed of three groups of proteins: Lhcf as the main group, Lhcr that are PSI associated, and Lhcx that are involved in photoprotection. The FCP complexes are assembled in trimers and higher oligomers. Several studies have investigated the biochemical properties of purified FCP complexes, but limited knowledge is available about their interaction with the photosystem cores. In this study, isolation of stable supercomplexes from the centric diatom Thalassiosira pseudonana was achieved. To preserve in vivo structure, the separation of thylakoid complexes was performed by native PAGE and sucrose density centrifugation. Different subpopulations of PSI and PSII supercomplexes were isolated and their subunits identified. Analysis of Lhc antenna composition identified Lhc(s) specific for either PSI (Lhcr 1, 3, 4, 7, 10–14, and Lhcf10) or PSII (Lhcf 1–7, 11, and Lhcr2). Lhcx6_1 was reproducibly found in PSII supercomplexes, whereas its association with PSI was unclear. No evidence was found for the interaction between photosystems and higher oligomeric FCPs, comprising Lhcf8 as the main component. Although the subunit composition of the PSII supercomplexes in comparison with that of the trimeric FCP complexes indicated a close mutual association, the higher oligomeric pool is only weakly associated with the photosystems, albeit its abundance in the thylakoid membrane.

Terpene volatiles are found in many important fruit crops, but their relationship to flavor is poorly understood. Here, we demonstrate using sensory descriptive and discriminant analysis that 1,8-cineole contributes a key floral/eucalyptus note to the aroma of ripe 'Hort16A’ kiwifruit (Actinidia chinensis). Two quantitative trait loci (QTLs) for 1,8-cineole production were identified on linkage groups 27 and 29a in a segregating A. chinensis population, with the QTL on LG29a colocating with a complex cluster of putative terpene synthase (TPS)-encoding genes. Transient expression in Nicotiana benthamiana and analysis of recombinant proteins expressed in Escherichia coli showed four genes in the cluster (AcTPS1a–AcTPS1d) encoded functional TPS enzymes, which produced predominantly sabinene, 1,8-cineole, geraniol, and springene, respectively. The terpene profile produced by AcTPS1b closely resembled the terpenes detected in red-fleshed A. chinensis. AcTPS1b expression correlated with 1,8-cineole content in developing/ripening fruit and also showed a positive correlation with 1,8-cineole content in the mapping population, indicating the basis for segregation is an expression QTL. Transient overexpression of AcTPS1b in Actinidia eriantha fruit confirmed this gene produced 1,8-cineole in Actinidia. Structure-function analysis showed AcTPS1a and AcTPS1b are natural variants at key TPS catalytic site residues previously shown to change enzyme specificity in vitro. Together, our results indicate that AcTPS1b is a key gene for production of the signature flavor terpene 1,8-cineole in ripe kiwifruit. Using a sensory-directed strategy for compound identification provides a rational approach for applying marker-aided selection to improving flavor in kiwifruit as well as other fruits.

Heat stress (HS) has serious effects on plant development, resulting in heavy agricultural losses. A critical transcription factor network is involved in plant adaptation to high temperature. DEHYDRATION RESPONSIVE ELEMENT-BINDING PROTEIN2A (DREB2A) is a key transcription factor that functions in plant thermotolerance. The DREB2A protein is unstable under normal temperature and is degraded by the 26S proteasome; however, the mechanism by which DREB2A protein stability dramatically increases in response to HS remains poorly understood. In this study, we found that the DREB2A protein of Arabidopsis (Arabidopsis thaliana) is stabilized under high temperature by the posttranslational modification SUMOylation. Biochemical data indicated that DREB2A is SUMOylated at K163, a conserved residue adjacent to the negative regulatory domain during HS. SUMOylation of DREB2A suppresses its interaction with BPM2, a ubiquitin ligase component, consequently increasing DREB2A protein stability under high temperature. In addition, analysis of plant heat tolerance and marker gene expression indicated that DREB2A SUMOylation is essential for its function in the HS response. Collectively, our data reveal a role for SUMOylation in the maintenance of DREB2A stability under high temperature, thus improving our understanding of the regulatory mechanisms underlying HS response in plant cells.

Multienvironment trials (METs) are widely used to assess the performance of promising crop germplasm. Though seldom designed to elucidate genetic mechanisms, MET data sets are often much larger than could be duplicated for genetic research and, given proper interpretation, may offer valuable insights into the genetics of adaptation across time and space. The Cooperative Dry Bean Nursery (CDBN) is a MET for common bean (Phaseolus vulgaris) grown for > 70 years in the United States and Canada, consisting of 20–50 entries each year at 10–20 locations. The CDBN provides a rich source of phenotypic data across entries, years, and locations that is amenable to genetic analysis. To study stable genetic effects segregating in this MET, we conducted genome-wide association studies (GWAS) using best linear unbiased predictions derived across years and locations for 21 CDBN phenotypes and genotypic data (1.2 million SNPs) for 327 CDBN genotypes. The value of this approach was confirmed by the discovery of three candidate genes and genomic regions previously identified in balanced GWAS. Multivariate adaptive shrinkage (mash) analysis, which increased our power to detect significant correlated effects, found significant effects for all phenotypes. Mash found two large genomic regions with effects on multiple phenotypes, supporting a hypothesis of pleiotropic or linked effects that were likely selected on in pursuit of a crop ideotype. Overall, our results demonstrate that statistical genomics approaches can be used on MET phenotypic data to discover significant genetic effects and to define genomic regions associated with crop improvement.

In plant–pathogen relations, disease symptoms arise from the interaction of the host and pathogen genomes. Host–pathogen functional gene interactions are well described, whereas little is known about how the pathogen genetic variation modulates both organisms’ transcriptomes. To model and generate hypotheses on a generalist pathogen control of gene expression regulation, we used the Arabidopsis thaliana–Botrytis cinerea pathosystem and the genetic diversity of a collection of 96 B. cinerea isolates. We performed expression-based genome-wide association (eGWA) for each of 23,947 measurable transcripts in Arabidopsis (host), and 9267 measurable transcripts in B. cinerea (pathogen). Unlike other eGWA studies, we detected a relative absence of locally acting expression quantitative trait loci (cis-eQTL), partly caused by structural variants and allelic heterogeneity hindering their identification. This study identified several distantly acting trans-eQTL linked to eQTL hotspots dispersed across Botrytis genome that altered only Botrytis transcripts, only Arabidopsis transcripts, or transcripts from both species. Gene membership in the trans-eQTL hotspots suggests links between gene expression regulation and both known and novel virulence mechanisms in this pathosystem. Genes annotated to these hotspots provide potential targets for blocking manipulation of the host response by this ubiquitous generalist necrotrophic pathogen.

Secondary branch number per panicle plays a crucial role in regulating grain number and yield in rice. Here, we report the positional cloning and functional characterization for SECONDARY BRANCH NUMBER7 (qSBN7), a quantitative trait locus affecting secondary branch per panicle and grain number. Our research revealed that the causative variants of qSBN7 are located in the distal promoter region of FRIZZLE PANICLE (FZP), a gene previously associated with the repression of axillary meristem development in rice spikelets. qSBN7 is a novel allele of FZP that causes an ~56% decrease in its transcriptional level, leading to increased secondary branch and grain number, and reduced grain length. Field evaluations showed that qSBN7 increased grain yield by 10.9% in a temperate japonica variety, TN13, likely due to its positive effect on sink capacity. Our findings suggest that incorporation of qSBN7 can increase yield potential and improve the breeding of elite rice varieties.

Many complex human traits exhibit differences between sexes. While numerous factors likely contribute to this phenomenon, growing evidence from genome-wide studies suggest a partial explanation: that males and females from the same population possess differing genetic architectures. Despite this, mapping gene-by-sex (GxS) interactions remains a challenge likely because the magnitude of such an interaction is typically and exceedingly small; traditional genome-wide association techniques may be underpowered to detect such events, due partly to the burden of multiple test correction. Here, we developed a local Bayesian regression (LBR) method to estimate sex-specific SNP marker effects after fully accounting for local linkage-disequilibrium (LD) patterns. This enabled us to infer sex-specific effects and GxS interactions either at the single SNP level, or by aggregating the effects of multiple SNPs to make inferences at the level of small LD-based regions. Using simulations in which there was imperfect LD between SNPs and causal variants, we showed that aggregating sex-specific marker effects with LBR provides improved power and resolution to detect GxS interactions over traditional single-SNP-based tests. When using LBR to analyze traits from the UK Biobank, we detected a relatively large GxS interaction impacting bone mineral density within ABO, and replicated many previously detected large-magnitude GxS interactions impacting waist-to-hip ratio. We also discovered many new GxS interactions impacting such traits as height and body mass index (BMI) within regions of the genome where both male- and female-specific effects explain a small proportion of phenotypic variance (R² < 1 x 10^–4), but are enriched in known expression quantitative trait loci.

Single-cross hybrids have been critical to the improvement of maize (Zea mays L.), but the characterization of their genetic architectures remains challenging. Previous studies of hybrid maize have shown the contribution of within-locus complementation effects (dominance) and their differential importance across functional classes of loci. However, they have generally considered panels of limited genetic diversity, and have shown little benefit from genomic prediction based on dominance or functional enrichments. This study investigates the relevance of dominance and functional classes of variants in genomic models for agronomic traits in diverse populations of hybrid maize. We based our analyses on a diverse panel of inbred lines crossed with two testers representative of the major heterotic groups in the U.S. (1106 hybrids), as well as a collection of 24 biparental populations crossed with a single tester (1640 hybrids). We investigated three agronomic traits: days to silking (DTS), plant height (PH), and grain yield (GY). Our results point to the presence of dominance for all traits, but also among-locus complementation (epistasis) for DTS and genotype-by-environment interactions for GY. Consistently, dominance improved genomic prediction for PH only. In addition, we assessed enrichment of genetic effects in classes defined by genic regions (gene annotation), structural features (recombination rate and chromatin openness), and evolutionary features (minor allele frequency and evolutionary constraint). We found support for enrichment in genic regions and subsequent improvement of genomic prediction for all traits. Our results suggest that dominance and gene annotations improve genomic prediction across diverse populations in hybrid maize.

Real geography is continuous, but standard models in population genetics are based on discrete, well-mixed populations. As a result, many methods of analyzing genetic data assume that samples are a random draw from a well-mixed population, but are applied to clustered samples from populations that are structured clinally over space. Here, we use simulations of populations living in continuous geography to study the impacts of dispersal and sampling strategy on population genetic summary statistics, demographic inference, and genome-wide association studies (GWAS). We find that most common summary statistics have distributions that differ substantially from those seen in well-mixed populations, especially when Wright’s neighborhood size is < 100 and sampling is spatially clustered. "Stepping-stone" models reproduce some of these effects, but discretizing the landscape introduces artifacts that in some cases are exacerbated at higher resolutions. The combination of low dispersal and clustered sampling causes demographic inference from the site frequency spectrum to infer more turbulent demographic histories, but averaged results across multiple simulations revealed surprisingly little systematic bias. We also show that the combination of spatially autocorrelated environments and limited dispersal causes GWAS to identify spurious signals of genetic association with purely environmentally determined phenotypes, and that this bias is only partially corrected by regressing out principal components of ancestry. Last, we discuss the relevance of our simulation results for inference from genetic variation in real organisms.

The question of the relative evolutionary roles of adaptive and nonadaptive processes has been a central debate in population genetics for nearly a century. While advances have been made in the theoretical development of the underlying models, and statistical methods for estimating their parameters from large-scale genomic data, a framework for an appropriate null model remains elusive. A model incorporating evolutionary processes known to be in constant operation, genetic drift (as modulated by the demographic history of the population) and purifying selection, is lacking. Without such a null model, the role of adaptive processes in shaping within- and between-population variation may not be accurately assessed. Here, we investigate how population size changes and the strength of purifying selection affect patterns of variation at "neutral" sites near functional genomic components. We propose a novel statistical framework for jointly inferring the contribution of the relevant selective and demographic parameters. By means of extensive performance analyses, we quantify the utility of the approach, identify the most important statistics for parameter estimation, and compare the results with existing methods. Finally, we reanalyze genome-wide population-level data from a Zambian population of Drosophila melanogaster, and find that it has experienced a much slower rate of population growth than was inferred when the effects of purifying selection were neglected. Our approach represents an appropriate null model, against which the effects of positive selection can be assessed.

Positive selection causes beneficial alleles to rise to high frequency, resulting in a selective sweep of the diversity surrounding the selected sites. Accordingly, the signature of a selective sweep in an ancestral population may still remain in its descendants. Identifying signatures of selection in the ancestor that are shared among its descendants is important to contextualize the timing of a sweep, but few methods exist for this purpose. We introduce the statistic SS-H12, which can identify genomic regions under shared positive selection across populations and is based on the theory of the expected haplotype homozygosity statistic H12, which detects recent hard and soft sweeps from the presence of high-frequency haplotypes. SS-H12 is distinct from comparable statistics because it requires a minimum of only two populations, and properly identifies and differentiates between independent convergent sweeps and true ancestral sweeps, with high power and robustness to a variety of demographic models. Furthermore, we can apply SS-H12 in conjunction with the ratio of statistics we term and to further classify identified shared sweeps as hard or soft. Finally, we identified both previously reported and novel shared sweep candidates from human whole-genome sequences. Previously reported candidates include the well-characterized ancestral sweeps at LCT and SLC24A5 in Indo-Europeans, as well as GPHN worldwide. Novel candidates include an ancestral sweep at RGS18 in sub-Saharan Africans involved in regulating the platelet response and implicated in sudden cardiac death, and a convergent sweep at C2CD5 between European and East Asian populations that may explain their different insulin responses.

The vertebrate limb serves as an experimental paradigm to study mechanisms that regulate development of the stereotypical skeletal elements. In this study, we simultaneously inactivated Sall4 using Hoxb6Cre and Plzf in mouse embryos, and found that their combined function regulates development of the proximal-anterior skeletal elements in hindlimbs. The Sall4; Plzf double knockout exhibits severe defects in the femur, tibia, and anterior digits, distinct defects compared to other allelic series of Sall4; Plzf. We found that Sall4 regulates Plzf expression prior to hindlimb outgrowth. Further expression analysis indicated that Hox10 genes and GLI3 are severely downregulated in the Sall4; Plzf double knockout hindlimb bud. In contrast, PLZF expression is reduced but detectable in Sall4; Gli3 double knockout limb buds, and SALL4 is expressed in the Plzf; Gli3 double knockout limb buds. These results indicate that Plzf, Gli3, and Hox10 genes downstream of Sall4, regulate femur and tibia development. In the autopod, we show that Sall4 negatively regulates Hedgehog signaling, which allows for development of the most anterior digit. Collectively, our study illustrates genetic systems that regulate development of the proximal-anterior skeletal elements in hindlimbs.

The Hippo pathway is an evolutionarily conserved signaling network that regulates organ size, cell fate, and tumorigenesis. In the context of organ size control, the pathway incorporates a large variety of cellular cues, such as cell polarity and adhesion, into an integrated transcriptional response. The central Hippo signaling effector is the transcriptional coactivator Yorkie, which controls gene expression in partnership with different transcription factors, most notably Scalloped. When it is not activated by Yorkie, Scalloped can act as a repressor of transcription, at least in part due to its interaction with the corepressor protein Tgi. The mechanism by which Tgi represses transcription is incompletely understood, and therefore we sought to identify proteins that potentially operate together with Tgi. Using an affinity purification and mass-spectrometry approach we identified Pits and CtBP as Tgi-interacting proteins, both of which have been linked to transcriptional repression. Both Pits and CtBP were required for Tgi to suppress the growth of the Drosophila melanogaster eye and CtBP loss suppressed the undergrowth of yorkie mutant eye tissue. Furthermore, as reported previously for Tgi, overexpression of Pits repressed transcription of Hippo pathway target genes. These findings suggest that Tgi might operate together with Pits and CtBP to repress transcription of genes that normally promote tissue growth. The human orthologs of Tgi, CtBP, and Pits (VGLL4, CTBP2, and IRF2BP2) have previously been shown to physically and functionally interact to control transcription, implying that the mechanism by which these proteins control transcriptional repression is conserved throughout evolution.

Repeated alcohol experiences can produce long-lasting memories for sensory cues associated with intoxication. These memories can problematically trigger relapse in individuals recovering from alcohol use disorder (AUD). The molecular mechanisms by which ethanol changes memories to become long-lasting and inflexible remain unclear. New methods to analyze gene expression within precise neuronal cell types can provide further insight toward AUD prevention and treatment. Here, we used genetic tools in Drosophila melanogaster to investigate the lasting consequences of ethanol on transcription in memory-encoding neurons. Drosophila rely on mushroom body (MB) neurons to make associative memories, including memories of ethanol-associated sensory cues. Differential expression analyses revealed that distinct transcripts, but not genes, in the MB were associated with experiencing ethanol alone compared to forming a memory of an odor cue associated with ethanol. Adult MB-specific knockdown of spliceosome-associated proteins demonstrated the necessity of RNA-processing in ethanol memory formation. These findings highlight the dynamic, context-specific regulation of transcription in cue-encoding neurons, and the lasting effect of ethanol on transcript usage during memory formation.

BEAF (Boundary Element-Associated Factor) was originally identified as a Drosophila melanogaster chromatin domain insulator-binding protein, suggesting a role in gene regulation through chromatin organization and dynamics. Genome-wide mapping found that BEAF usually binds near transcription start sites, often of housekeeping genes, suggesting a role in promoter function. This would be a nontraditional role for an insulator-binding protein. To gain insight into molecular mechanisms of BEAF function, we identified interacting proteins using yeast two-hybrid assays. Here, we focus on the transcription factor Serendipity (Sry-). Interactions were confirmed in pull-down experiments using bacterially expressed proteins, by bimolecular fluorescence complementation, and in a genetic assay in transgenic flies. Sry- interacted with promoter-proximal BEAF both when bound to DNA adjacent to BEAF or > 2-kb upstream to activate a reporter gene in transient transfection experiments. The interaction between BEAF and Sry- was detected using both a minimal developmental promoter (y) and a housekeeping promoter (RpS12), while BEAF alone strongly activated the housekeeping promoter. These two functions for BEAF implicate it in playing a direct role in gene regulation at hundreds of BEAF-associated promoters.

Replication initiation in eukaryotic cells occurs asynchronously throughout S phase, yielding early- and late-replicating regions of the genome, a process known as replication timing (RT). RT changes during development to ensure accurate genome duplication and maintain genome stability. To understand the relative contributions that cell lineage, cell cycle, and replication initiation regulators have on RT, we utilized the powerful developmental systems available in Drosophila melanogaster. We generated and compared RT profiles from mitotic cells of different tissues and from mitotic and endocycling cells of the same tissue. Our results demonstrate that cell lineage has the largest effect on RT, whereas switching from a mitotic to an endoreplicative cell cycle has little to no effect on RT. Additionally, we demonstrate that the RT differences we observed in all cases are largely independent of transcriptional differences. We also employed a genetic approach in these same cell types to understand the relative contribution the eukaryotic RT control factor, Rif1, has on RT control. Our results demonstrate that Rif1 can function in a tissue-specific manner to control RT. Importantly, the Protein Phosphatase 1 (PP1) binding motif of Rif1 is essential for Rif1 to regulate RT. Together, our data support a model in which the RT program is primarily driven by cell lineage and is further refined by Rif1/PP1 to ultimately generate tissue-specific RT programs.

During meiotic prophase, concurrent transcription, recombination, and chromosome synapsis place substantial topological strain on chromosomal DNA, but the role of topoisomerases in this context remains poorly defined. Here, we analyzed the roles of topoisomerases I and II (Top1 and Top2) during meiotic prophase in Saccharomyces cerevisiae. We show that both topoisomerases accumulate primarily in promoter-containing intergenic regions of actively transcribing genes, including many meiotic double-strand break (DSB) hotspots. Despite the comparable binding patterns, top1 and top2 mutations have different effects on meiotic recombination. TOP1 disruption delays DSB induction and shortens the window of DSB accumulation by an unknown mechanism. By contrast, temperature-sensitive top2-1 mutants exhibit a marked delay in meiotic chromosome remodeling and elevated DSB signals on synapsed chromosomes. The problems in chromosome remodeling were linked to altered Top2 binding patterns rather than a loss of Top2 catalytic activity, and stemmed from a defect in recruiting the chromosome remodeler Pch2/TRIP13 to synapsed chromosomes. No chromosomal defects were observed in the absence of TOP1. Our results imply independent roles for Top1 and Top2 in modulating meiotic chromosome structure and recombination.

Age-at-onset is one of the critical traits in cohort studies of age-related diseases. Large-scale genome-wide association studies (GWAS) of age-at-onset traits can provide more insights into genetic effects on disease progression and transitions between stages. Moreover, proportional hazards (or Cox) regression models can achieve higher statistical power in a cohort study than a case-control trait using logistic regression. Although mixed-effects models are widely used in GWAS to correct for sample dependence, application of Cox mixed-effects models (CMEMs) to large-scale GWAS is so far hindered by intractable computational cost. In this work, we propose COXMEG, an efficient R package for conducting GWAS of age-at-onset traits using CMEMs. COXMEG introduces fast estimation algorithms for general sparse relatedness matrices including, but not limited to, block-diagonal pedigree-based matrices. COXMEG also introduces a fast and powerful score test for dense relatedness matrices, accounting for both population stratification and family structure. In addition, COXMEG generalizes existing algorithms to support positive semidefinite relatedness matrices, which are common in twin and family studies. Our simulation studies suggest that COXMEG, depending on the structure of the relatedness matrix, is orders of magnitude computationally more efficient than coxme and coxph with frailty for GWAS. We found that using sparse approximation of relatedness matrices yielded highly comparable results in controlling false-positive rate and retaining statistical power for an ethnically homogeneous family-based sample. By applying COXMEG to a study of Alzheimer’s disease (AD) with a Late-Onset Alzheimer’s Disease Family Study from the National Institute on Aging sample comprising 3456 non-Hispanic whites and 287 African Americans, we identified the APOE 4 variant with strong statistical power (P = 1e–101), far more significant than that reported in a previous study using a transformed variable and a marginal Cox model. Furthermore, we identified novel SNP rs36051450 (P = 2e–9) near GRAMD1B, the minor allele of which significantly reduced the hazards of AD in both genders. These results demonstrated that COXMEG greatly facilitates the application of CMEMs in GWAS of age-at-onset traits.

There is increasing interest in developing diagnostics that discriminate individual mutagenic mechanisms in a range of applications that include identifying population-specific mutagenesis and resolving distinct mutation signatures in cancer samples. Analyses for these applications assume that mutagenic mechanisms have a distinct relationship with neighboring bases that allows them to be distinguished. Direct support for this assumption is limited to a small number of simple cases, e.g., CpG hypermutability. We have evaluated whether the mechanistic origin of a point mutation can be resolved using only sequence context for a more complicated case. We contrasted single nucleotide variants originating from the multitude of mutagenic processes that normally operate in the mouse germline with those induced by the potent mutagen N-ethyl-N-nitrosourea (ENU). The considerable overlap in the mutation spectra of these two samples make this a challenging problem. Employing a new, robust log-linear modeling method, we demonstrate that neighboring bases contain information regarding point mutation direction that differs between the ENU-induced and spontaneous mutation variant classes. A logistic regression classifier exhibited strong performance at discriminating between the different mutation classes. Concordance between the feature set of the best classifier and information content analyses suggest our results can be generalized to other mutation classification problems. We conclude that machine learning can be used to build a practical classification tool to identify the mutation mechanism for individual genetic variants. Software implementing our approach is freely available under an open-source license.

The way genes contribute to behavior is complicated. Although there are some single genes with large contributions, most behavioral differences are due to small effects from many interacting genes. This makes it hard to identify the genes that cause behavioral differences. Mutagenesis screens in model organisms, selective breeding experiments in animals, comparisons between related populations with different behaviors, and genome-wide association studies in humans are promising and complementary approaches to understanding the heritable aspects of complex behaviors. To connect genes to behaviors requires measuring behavioral differences, locating correlated genetic changes, determining when, where, and how these candidate genes act, and designing causative confirmatory experiments. This area of research has implications from basic discovery science to human mental health.