Quantitative genetic analysis of interactions in the pepper-Phytophthora capsici pathosystem

ABSTRACTThe development of pepper cultivars with durable resistance to the oomycete Phytophthora capsici has been challenging due to differential interactions between the species that allow certain pathogen isolates to cause disease on otherwise resistant host genotypes. Currently, little is known about the pathogen genes that are involved in these interactions. To investigate the genetic basis of P. capsici virulence on individual pepper genotypes, we inoculated sixteen pepper accessions – representing commercial varieties, sources of resistance, and host differentials – with 117 isolates of P. capsici, for a total of 1,864 host-pathogen combinations. Analysis of disease outcomes revealed a significant effect of inter-species genotype-by-genotype interactions, although these interactions were quantitative rather than qualitative in scale. Isolates were classified into five pathogen subpopulations, as determined by their genotypes at over 60,000 single-nucleotide polymorphisms (SNPs). While absolute virulence levels on certain pepper accessions significantly differed between subpopulations, a multivariate phenotype reflecting relative virulence levels on certain pepper genotypes compared to others showed the strongest association with pathogen subpopulation. A genome-wide association study (GWAS) identified four pathogen loci significantly associated with virulence, two of which colocalized with putative RXLR effector genes and another with a polygalacturonase gene cluster. All four loci appeared to represent broad-spectrum virulence genes, as significant SNPs demonstrated consistent effects regardless of the host genotype tested. Host genotype-specific virulence variants in P. capsici may be difficult to map via GWAS, perhaps controlled by many genes of small effect or by multiple alleles that have arisen independently at the same loci.

Modeling the efficacy of CRISPR gene drive for schistosomiasis control

CRISPR gene drives could revolutionize the control of infectious diseases by accelerating the spread of engineered traits that limit parasite transmission in wild populations. While much effort has been spent developing gene drives in mosquitoes, gene drive technology in molluscs has received little attention despite the role of freshwater snails as obligate, intermediate hosts of parasitic flukes causing schistosomiasis -- a disease of poverty affecting more than 200 million people worldwide. A successful drive in snails must overcome self-fertilization, which prevents a drive's spread. Simultaneous hermaphroditism is a feature of snails -- distinct from gene drive model organisms -- and is not yet incorporated in gene drive models of disease control. Here we developed a novel population genetic model accounting for snails' sexual and asexual reproduction, susceptibility to parasite infection regulated by multiple alleles, fitness differences between genotypes, and a range of drive characteristics. We then integrated this model with an epidemiological model of schistosomiasis transmission and snail population dynamics. Simulations showed that gene drive establishment can be hindered by a variety of biological and ecological factors, including selfing. However, our model suggests that, under a range of conditions, gene drive mediated immunity in snails could maintain rapid disease reduction achieved by annual chemotherapy treatment of the human population, leading to long-term elimination. These results indicate that gene drives, in coordination with existing public health measures, may become a useful tool to reduce schistosomiasis burden in selected transmission settings with effective CRISPR construct design and close evaluation of the genetic and ecological landscape.

Remote-Sensing-Combined Haplotype Analysis Using Multi-Parental Advanced Generation Inter-Cross Lines Reveals Phenology QTLs for Canopy Height in Rice

High-throughput phenotyping systems with unmanned aerial vehicles (UAVs) enable observation of crop lines in the field. In this study, we show the ability of time-course monitoring of canopy height (CH) to identify quantitative trait loci (QTLs) and to characterise their pleiotropic effect on various traits. We generated a digital surface model from low-altitude UAV-captured colour digital images and investigated CH data of rice multi-parental advanced generation inter-cross (MAGIC) lines from tillering and heading to maturation. Genome-wide association studies (GWASs) using the CH data and haplotype information of the MAGIC lines revealed 11 QTLs for CH. Each QTL showed haplotype effects on different features of CH such as stage-specificity and constancy. Haplotype analysis revealed relationships at the QTL level between CH and, vegetation fraction and leaf colour [derived from UAV red–green–blue (RGB) data], and CH and yield-related traits. Noticeably, haplotypes with canopy lowering effects at qCH1-4, qCH2, and qCH10-2 increased the ratio of panicle weight to leaf and stem weight, suggesting biomass allocation to grain yield or others through growth regulation of CH. Allele mining using gene information with eight founders of the MAGIC lines revealed the possibility that qCH1-4 contains multiple alleles of semi-dwarf 1 (sd1), the IR-8 allele of which significantly contributed to the “green revolution” in rice. This use of remote-sensing-derived phenotyping data into genetics using the MAGIC lines gives insight into how rice plants grow, develop, and produce grains in phenology and provides information on effective haplotypes for breeding with ideal plant architecture and grain yield.

Improved Characterization of Balancing Selection Genome Wide and a Detailed Look at HLA Genes

Balancing selection occurs when different evolutionary pressures impact the fitness of multiple alleles, resulting in increased allelic diversity in the population. A new statistical method was developed to test for selection, improving inference by using efficient Bayesian techniques to test for density and strength of linkage disequilibrium. Evolutionary simulation studies showed that the method consistently outperformed existing methods. Using this methodology, we tested for novel signals of balancing selection genome wide in 500 samples from phased trios. Several novel signals of selection appeared in CYP2A7, GPC6, and CNR2 across multiple ancestries. Additionally, tests in SIRPA demonstrate dramatically strong selection signal, significantly higher than previously observed. Well-known signals around olfactory genes and the MHC, containing HLA genes associated with the immune response, also demonstrated strong signatures of selection. So, utilizing data from the 17th IHIW, a follow up analysis was then performed by leveraging over seven thousand HLA typed samples by NGS; in contrast, the genome wide scan did not include a detailed characterization of the HLA genes. The strongest signals observed in the IHIW samples were in DQA1 and DQB1 in or around exon 2, the portion of the gene responsible for antigen presentation and most likely to be under environmental and evolutionary pressure. Our new statistical approach and analysis suggest novel evolutionary pressure in new regions and additionally highlight the importance of improved sequencing and characterization of variation across the extended MHC and other critical regions.

Tyr is Responsible for the Cctq1a QTL and Links Developmental Environment to Central Corneal Thickness Determination

Central corneal thickness is a quantitative trait with important associations to human health. In a phenotype-driven approach studying corneal thickness of congenic derivatives of C57BLKS/J and SJL/J mice, the critical region for a quantitative trait locus influencing corneal thickness, Cctq1a, was delimited to a 10-gene interval. Exome sequencing, RNAseq, and studying independent mutations eliminated multiple candidate genes and confirmed one. Though the causative gene, Tyr, has no obvious direct function in the transparent cornea, studies with multiple alleles on matched genetic backgrounds, both in isolation and genetic complementation crosses, confirmed allelism of Tyr-Cctq1a; albino mice lacking Tyr function had thin corneas. Albino mice also had increased axial length. Because albinism exposes eyes to increased light, the effect of dark-rearing was tested and found to rescue central corneal thickness. In sum, the results point to an epiphenomenon; developmental light exposure interacts with genotype as an important determinate of adult corneal thickness.

Dog colour patterns explained by modular promoters of ancient canid origin

Distinctive colour patterns in dogs are an integral component of canine diversity. Colour pattern differences are thought to have arisen from mutation and artificial selection during and after domestication from wolves but important gaps remain in understanding how these patterns evolved and are genetically controlled. In other mammals, variation at the ASIP gene controls both the temporal and spatial distribution of yellow and black pigments. Here, we identify independent regulatory modules for ventral and hair cycle ASIP expression, and we characterize their action and evolutionary origin. Structural variants define multiple alleles for each regulatory module and are combined in different ways to explain five distinctive dog colour patterns. Phylogenetic analysis reveals that the haplotype combination for one of these patterns is shared with Arctic white wolves and that its hair cycle-specific module probably originated from an extinct canid that diverged from grey wolves more than 2 million years ago. Natural selection for a lighter coat during the Pleistocene provided the genetic framework for widespread colour variation in dogs and wolves.

The genetics of gene expression in a C. elegans multi parental recombinant inbred line population

Studying genetic variation of gene expression provides a powerful way to unravel the molecular components underlying complex traits. Expression QTL studies have been performed in several different model species, yet most of these linkage studies have been based on genetic segregation of two parental alleles. Recently we developed a multi-parental segregating population of 200 recombinant inbred lines (mpRILs) derived from four wild isolates (JU1511, JU1926, JU1931 and JU1941) in the nematode Caenorhabditis elegans. We used RNA-seq to investigate how multiple alleles affect gene expression in these mpRILs. We found 1,789 genes differentially expressed between the parental lines. Transgression, expression beyond any of the parental lines in the mpRILs, was found for 7,896 genes. For expression QTL mapping almost 9,000 SNPs were available. By combining these SNPs and the RNA-seq profiles of the mpRILs, we detected almost 6,800 eQTLs. Most trans-eQTLs (63%) co-locate in six newly identified trans-bands. The trans-eQTLs found in previous 2-parental allele eQTL experiments and this study showed some overlap (17.5%-46.8%), highlighting on the one hand that a large group of genes is affected by polymorphic regulators across populations and conditions, on the other hand it shows that the mpRIL population allows identification of novel gene expression regulatory loci. Taken together, the analysis of our mpRIL population provides a more refined insight into C. elegans complex trait genetics and eQTLs in general, as well as a starting point to further test and develop advanced statistical models for detection of multi-allelic eQTLs and systems genetics studying the genotype-phenotype relationship.

Pseudo-Random Mating with Multiple Alleles

The conditions on the mating matrix associated with a stable equilibrium are specified for an autosomal locus with four alleles. An example illustrates how Hardy–Weinberg proportions are maintained with nonrandom mating. The ABO blood group provides an illustration.

Examining the distribution and impact of single nucleotide polymorphisms in the capsular locus of Streptococcus pneumoniae serotype 19A

Streptococcus pneumoniae serotype 19A prevalence has increased after implementation of PCV7 and PCV10 vaccines. In this study, we have provided, with high accuracy, the genetic diversity of the 19A serotype in a cohort of Dutch invasive pneumococcal disease patients and asymptomatic carriers obtained in the period 2004-2016. Whole genomes of the 338 pneumococcal isolates in this cohort were sequenced and their capsule ( cps ) loci compared to examine the diversity and determine the impact on the production of CPS sugar precursors and CPS shedding. We discovered 79 types with a unique CPS locus sequence. Most variation was observed in the rmlB and rmlD genes of the TDP-Rha synthesis pathway, and in the wzg gene, of unknown function. Interestingly, gene variation in the cps locus was conserved in multiple alleles. Using RmlB and RmlD protein models, we predict that enzymatic function is not affected by the single nucleotide polymorphisms as identified. To determine if RmlB and RmlD function was affected, we analyzed nucleotide sugar levels using UHPLC-MS. CPS precursors differed between 19A cps locus subtypes, including TDP-Rha, but no clear correlation was observed. Also, a significant difference in multiple nucleotide sugar levels was observed between phylogenetically branched groups. Because of indications of a role for Wzg in capsule shedding, we analyzed if this was affected. No clear indication of a direct role in shedding was found. We thus describe genotypic variety in rmlB , rmlD and wzg in serotype 19A the Netherlands, for which we have not discovered an associated phenotype.

CRISPR/Cas9-Mediated Multi-Allelic Gene Targeting in Sugarcane Confers Herbicide Tolerance

Sugarcane is the source of 80% of the sugar and 26% of the bioethanol produced globally. However, its complex, highly polyploid genome (2n = 100 – 120) impedes crop improvement. Here, we report efficient and reproducible gene targeting (GT) in sugarcane, enabling precise co-editing of multiple alleles via template-mediated and homology-directed repair (HDR) of DNA double strand breaks induced by the programmable nuclease CRISPR/Cas9. The evaluation of 146 independently transformed plants from five independent experiments revealed a targeted nucleotide replacement that resulted in both targeted amino acid substitutions W574L and S653I in the acetolactate synthase (ALS) in 11 lines in addition to single, targeted amino acid substitutions W574L or S653I in 25 or 18 lines, respectively. Co-editing of up to three ALS copies/alleles that confer herbicide tolerance was confirmed by Sanger sequencing of cloned long polymerase chain reaction (PCR) amplicons. This work will enable crop improvement by conversion of inferior alleles to superior alleles through targeted nucleotide substitutions.