Genomic associations with bill length and disease reveal drift and selection across island bird populations

Island species provide excellent models for investigating how selection and drift operate in wild populations, and for determining how these processes act to influence local adaptation and speciation. Here, we examine the role of selection and drift in shaping genomic and phenotypic variation across recently separated populations of Berthelot's pipit (Anthus berthelotii), a passerine bird endemic to three archipelagos in the Atlantic. We first characterized genetic diversity and population structuring that supported previous inferences of a history of recent colonizations and bottlenecks. We then tested for regions of the genome associated with the ecologically important traits of bill length and malaria infection, both of which vary substantially across populations in this species. We identified a SNP associated with variation in bill length among individuals, islands, and archipelagos; patterns of variation at this SNP suggest that both phenotypic and genotypic variation in bill length is largely shaped by founder effects. Malaria was associated with SNPs near/within genes involved in the immune response, but this relationship was not consistent among archipelagos, supporting the view that disease resistance is complex and rapidly evolving. Although we found little evidence for divergent selection at candidate loci for bill length and malaria resistance, genome scan analyses pointed to several genes related to immunity and metabolism as having important roles in divergence and adaptation. Our findings highlight the utility and challenges involved with combining association mapping and population genetic analysis in nonequilibrium populations, to disentangle the effects of drift and selection on shaping genotypes and phenotypes

Modeling QTL-by-environment interactions for multi-parent populations

Multi-parent populations (MPPs) are attractive for genetic and breeding studies because they combine genetic diversity with an easy-to-control population structure. Most methods for mapping QTLs in MPPs focus on the detection of QTLs in single environments. Little attention has been given to mapping QTLs in multienvironment trials (METs) and to detecting and modeling QTL-by-environment interactions (QEIs). We present mixed model approaches for the detection and modeling of consistent versus environment-dependent QTLs, i.e., QTL-by-environment interaction (QEI). QTL effects are assumed to be normally distributed with variances expressing consistency or dependence on environments and families. The entries of the corresponding design matrices are functions of identity-by-descent (IBD) probabilities between parents and offspring and follow from the parental origin of offspring DNA. A polygenic effect is added to the models to account for background genetic variation. We illustrate the wide applicability of our method by analyzing several public MPP datasets with observations from METs. The examples include diallel, nested association mapping (NAM), and multi-parent advanced inter-cross (MAGIC) populations. The results of our approach compare favorably with those of previous studies that used tailored methods

Sex-specific changes in the aphid DNA methylation landscape

Aphids present an ideal system to study epigenetics as they can produce diverse, but genetically identical, morphs in response to environmental stimuli. Here, using whole genome bisulphite sequencing and transcriptome sequencing of the green peach aphid (Myzus persicae), we present the first detailed analysis of cytosine methylation in an aphid and investigate differences in the methylation and transcriptional landscapes of male and asexual female morphs. We find that methylation primarily occurs in a CG dinucleotide (CpG) context and that exons are highly enriched for methylated CpGs, particularly at the 3' end of genes. Methylation is positively associated with gene expression, and methylated genes are more stably expressed than un-methylated genes. Male and asexual female morphs have distinct methylation profiles. Strikingly, these profiles are divergent between the sex chromosome and the autosomes; autosomal genes are hypo-methylated in males compared to asexual females, whereas genes belonging to the sex chromosome, which is haploid in males, are hyper-methylated. Overall, we find correlated changes in methylation and gene expression between males and asexual females, and this correlation is particularly strong for genes located on the sex chromosome. Our results suggest that differential methylation of sex-biased genes plays a role in aphid sexual differentiation

Genomic erosion in the assessment of species extinction risk and recovery potential

Many species are facing unprecedented population size declines and deterioration of their environment. This exposes species to genomic erosion, which we define here as the damage inflicted to a species’ genome or gene pool due to a loss of genetic diversity, an increase in expressed genetic load, maladaptation, and/or genetic introgression. The International Union for Conservation of Nature (IUCN) bases its extinction risk assessments on direct threats to population size and habitat. However, it does not assess the long-term impacts of genomic erosion, and hence, it is likely to underestimate the extinction risk of many species. High-quality whole genome sequence data that is currently being generated could help improve extinction risk assessments. Genomic data contains information about a species’ past demography, its genome-wide genetic diversity, the incidence of genetic introgression, as well as the genetic load of deleterious mutations. Computer modelling of these data enables forecasting of population trajectories under different management scenarios. In this Perspective, we discuss the threats posed by genomic erosion. Using evolutionary genomic simulations, we argue that whole genome sequence data provides critical information for assessing the extinction risk and recovery potential of species. Genomics-informed assessments of the extinction risk complement the IUCN Red List, and such genomics-informed conservation is invaluable in guiding species recovery programs in the UN’s Decade on Ecosystem Restoration and beyond

Semi‐quantitative characterisation of mixed pollen samples using MinION sequencing and Reverse Metagenomics (RevMet)

1. The ability to identify and quantify the constituent plant species that make up a mixed‐species sample of pollen has important applications in ecology, conservation, and agriculture. Recently, metabarcoding protocols have been developed for pollen that can identify constituent plant species, but there are strong reasons to doubt that metabarcoding can accurately quantify their relative abundances. A PCR‐free, shotgun metagenomics approach has greater potential for accurately quantifying species relative abundances, but applying metagenomics to eukaryotes is challenging due to low numbers of reference genomes. 2. We have developed a pipeline, RevMet (Reverse Metagenomics) that allows reliable and semi‐quantitative characterization of the species composition of mixed‐species eukaryote samples, such as bee‐collected pollen, without requiring reference genomes. Instead, reference species are represented only by ‘genome skims’: low‐cost, low‐coverage, short‐read sequence datasets. The skims are mapped to individual long reads sequenced from mixed‐species samples using the MinION, a portable nanopore sequencing device, and each long read is uniquely assigned to a plant species. 3. We genome‐skimmed 49 wild UK plant species, validated our pipeline with mock DNA mixtures of known composition, and then applied RevMet to pollen loads collected from wild bees. We demonstrate that RevMet can identify plant species present in mixed‐species samples at proportions of DNA ≥ 1%, with few false positives and false negatives, and reliably differentiate species represented by high versus low amounts of DNA in a sample. 4. RevMet could readily be adapted to generate semi‐quantitative datasets for a wide range of mixed eukaryote samples. Our per‐sample costs were £90 per genome skim and £60 per pollen sample, and new versions of sequencers available now will further reduce these costs

Genetic analysis of wheat sensitivity to the ToxB fungal effector from Pyrenophora tritici-repentis, the causal agent of tan spot.

Genetic mapping of sensitivity to the Pyrenophora tritici-repentis effector ToxB allowed development of a diagnostic genetic marker, and investigation of wheat pedigrees allowed transmission of sensitive alleles to be tracked. Tan spot, caused by the necrotrophic fungal pathogen Pyrenophora tritici-repentis, is a major disease of wheat (Triticum aestivum). Secretion of the P. tritici-repentis effector ToxB is thought to play a part in mediating infection, causing chlorosis of plant tissue. Here, genetic analysis using an association mapping panel (n = 480) and a multiparent advanced generation intercross (MAGIC) population (n founders = 8, n progeny = 643) genotyped with a 90,000 feature single nucleotide polymorphism (SNP) array found ToxB sensitivity to be highly heritable (h2 ≥ 0.9), controlled predominantly by the Tsc2 locus on chromosome 2B. Genetic mapping of Tsc2 delineated a 1921-kb interval containing 104 genes in the reference genome of ToxB-insensitive variety 'Chinese Spring'. This allowed development of a co-dominant genetic marker for Tsc2 allelic state, diagnostic for ToxB sensitivity in the association mapping panel. Phenotypic and genotypic analysis in a panel of wheat varieties post-dated the association mapping panel further supported the diagnostic nature of the marker. Combining ToxB phenotype and genotypic data with wheat pedigree datasets allowed historic sources of ToxB sensitivity to be tracked, finding the variety 'Maris Dove' to likely be the historic source of sensitive Tsc2 alleles in the wheat germplasm surveyed. Exploration of the Tsc2 region gene space in the ToxB-sensitive line 'Synthetic W7984' identified candidate genes for future investigation. Additionally, a minor ToxB sensitivity QTL was identified on chromosome 2A. The resources presented here will be of immediate use for marker-assisted selection for ToxB insensitivity and the development of germplasm with additional genetic recombination within the Tsc2 region

Genetic analysis of wheat sensitivity to the ToxB fungal effector from Pyrenophora tritici-repentis, the causal agent of tan spot

Abstract: Key message: Genetic mapping of sensitivity to the Pyrenophora tritici-repentis effector ToxB allowed development of a diagnostic genetic marker, and investigation of wheat pedigrees allowed transmission of sensitive alleles to be tracked. Abstract: Tan spot, caused by the necrotrophic fungal pathogen Pyrenophora tritici-repentis, is a major disease of wheat (Triticum aestivum). Secretion of the P. tritici-repentis effector ToxB is thought to play a part in mediating infection, causing chlorosis of plant tissue. Here, genetic analysis using an association mapping panel (n = 480) and a multiparent advanced generation intercross (MAGIC) population (n founders = 8, n progeny = 643) genotyped with a 90,000 feature single nucleotide polymorphism (SNP) array found ToxB sensitivity to be highly heritable (h2 ≥ 0.9), controlled predominantly by the Tsc2 locus on chromosome 2B. Genetic mapping of Tsc2 delineated a 1921-kb interval containing 104 genes in the reference genome of ToxB-insensitive variety ‘Chinese Spring’. This allowed development of a co-dominant genetic marker for Tsc2 allelic state, diagnostic for ToxB sensitivity in the association mapping panel. Phenotypic and genotypic analysis in a panel of wheat varieties post-dated the association mapping panel further supported the diagnostic nature of the marker. Combining ToxB phenotype and genotypic data with wheat pedigree datasets allowed historic sources of ToxB sensitivity to be tracked, finding the variety ‘Maris Dove’ to likely be the historic source of sensitive Tsc2 alleles in the wheat germplasm surveyed. Exploration of the Tsc2 region gene space in the ToxB-sensitive line ‘Synthetic W7984’ identified candidate genes for future investigation. Additionally, a minor ToxB sensitivity QTL was identified on chromosome 2A. The resources presented here will be of immediate use for marker-assisted selection for ToxB insensitivity and the development of germplasm with additional genetic recombination within the Tsc2 region

An improved assembly and annotation of the allohexaploid wheat genome identifies complete families of agronomic genes and provides genomic evidence for chromosomal translocations

Advances in genome sequencing and assembly technologies are generating many high-quality genome sequences, but assemblies of large, repeat-rich polyploid genomes, such as that of bread wheat, remain fragmented and incomplete. We have generated a new wheat whole-genome shotgun sequence assembly using a combination of optimized data types and an assembly algorithm designed to deal with large and complex genomes. The new assembly represents >78% of the genome with a scaffold N50 of 88.8 kb that has a high fidelity to the input data. Our new annotation combines strand-specific Illumina RNA-seq and Pacific Biosciences (PacBio) full-length cDNAs to identify 104,091 high-confidence protein-coding genes and 10,156 noncoding RNA genes. We confirmed three known and identified one novel genome rearrangements. Our approach enables the rapid and scalable assembly of wheat genomes, the identification of structural variants, and the definition of complete gene models, all powerful resources for trait analysis and breeding of this key global crop

Heterarchy of Transcription Factors Driving Basal and Luminal Cell Phenotypes in Human Urothelium

Cell differentiation is effected by complex networks of transcription factors that co-ordinate re-organisation of the chromatin landscape. The hierarchies of these relationships can be difficult to dissect. During in vitro differentiation of normal human uro-epithelial cells, formaldehyde-assisted isolation of regulatory elements (FAIRE-seq) and RNA-seq were used to identify alterations in chromatin accessibility and gene expression changes following activation of the nuclear receptor PPARG as a differentiation-initiating event. Regions of chromatin identified by FAIRE-seq as having altered accessibility during differentiation were found to be enriched with sequence-specific binding motifs for transcription factors predicted to be involved in driving basal and differentiated urothelial cell phenotypes, including FOXA1, P63, GRHL2, CTCF and GATA3. In addition, co-occurrence of GATA3 motifs was observed within sub-sets of differentiation-specific peaks containing P63 or FOXA1 after induction of differentiation. Changes in abundance of GRHL2, GATA3, and P63 were observed in immunoblots of chromatin-enriched extracts. Transient siRNA knockdown of P63 revealed that P63 favoured a basal-like phenotype by inhibiting differentiation and promoting expression of basal marker genes. GATA3 siRNA prevented differentiation-associated downregulation of P63 protein and transcript, and demonstrated positive feedback of GATA3 on PPARG transcript, but showed no effect on FOXA1 transcript or protein expression. This approach indicates that as a transcriptionally-regulated programme, urothelial differentiation operates as a heterarchy wherein GATA3 is able to co-operate with FOXA1 to drive expression of luminal marker genes, but that P63 has potential to transrepress expression of the same genes