Mapping QTL associated with resistance to Pseudomonas syringae pv. actinidiae in kiwifruit (Actinidia chinensis var. chinensis)

Pseudomonas syringae pv. actinidiae (Psa) is a bacterial pathogen of kiwifruit. This pathogen causes leaf-spotting, cane dieback, wilting, cankers (lesions), and in severe cases, plant death. Families of diploid A. chinensis seedlings grown in the field show a range of susceptibilities to the disease with up to 100% of seedlings in some families succumbing to Psa. But the effect of selection for field resistance to Psa on the alleles that remain in surviving seedlings has not been assessed. The objective of this work was to analyse, the effect of plant removal from Psa on the allele frequency of an incomplete-factorial-cross population. This population was founded using a range of genotypically distinct diploid A. chinensis var. chinensis parents to make 28 F1 families. However, because of the diversity of these families, low numbers of surviving individuals, and a lack of samples from dead individuals, standard QTL mapping approaches were unlikely to yield good results. Instead, a modified bulk segregant analysis (BSA) overcame these drawbacks while reducing the costs of sampling and sample processing, and the complexity of data analysis. Because the method was modified, part one of this work was used to determine the signal strength required for a QTL to be detected with BSA. Once QTL detection accuracy was known, part two of this work analysed the 28 families from the incomplete-factorial-cross population that had multiple individuals removed due to Psa infection. Each family was assigned to one of eight bulks based on a single parent that contributed to the families. DNA was extracted in bulk by grinding sampled leaf discs together before DNA extraction. Each sample bulk was compared against a bulk made up of WGS data from the parents contributing to the sample bulk. The deviation in allele frequency from the expected allele frequency within surviving populations using the modified BSA method was able to identify 11 QTLs for Psa that were present in at least two analyses. The identification of these Psa resistance QTL will enable marker development to selectively breed for resistance to Psa in future kiwifruit breeding programs.

Reconstructing queen genotypes by pool sequencing colonies in eusocial insects: Statistical Methods and their application to honeybee

Eusocial insects are crucial to many ecosystems, and particularly the honeybee (Apis mellifera). One approach to facilitate their study in molecular genetics, is to consider whole-colony genotyping by combining DNA of multiple individuals in a single pool sequencing experiment. Cheap and fast, this technique comes with the drawback of producing data requiring dedicated methods to be fully exploited. Despite this limitation, pool sequencing data have been shown to be informative and cost-effective when working on random mating populations. Here, we present new statistical methods for exploiting pool sequencing of eusocial colonies in order to reconstruct the genotypes of the queen of such colony. This leverages the possibility to monitor genetic diversity, perform genomic-based studies or implement selective breeding. Using simulations and honeybee real data, we show that the new methods allow for a fast and accurate estimation of the queen's genetic ancestry, with correlations of about 0.9 to that obtained from individual genotyping. Also, it allows for an accurate reconstruction of the queen genotypes, with about 2% genotyping error. We further validate these inferences using experimental data on colonies with both pool sequencing and individual genotyping of drones. In brief, in this study we present statistical models to accurately estimate the genetic ancestry and reconstruct the genotypes of the queen from pool sequencing data from workers of an eusocial colony. Such information allows to exploit pool sequencing for traditional population genetics analyses, association studies and for selective breeding. While validated in Apis mellifera, these methods are applicable to other eusocial hymenopterans.

Contribution of Rare and Low-Frequency Variants to Multiple Sclerosis Susceptibility in the Italian Continental Population

Genome-wide association studies identified over 200 risk loci for multiple sclerosis (MS) focusing on common variants, which account for about 50% of disease heritability. The goal of this study was to investigate whether low-frequency and rare functional variants, located in MS-established associated loci, may contribute to disease risk in a relatively homogeneous population, testing their cumulative effect (burden) with gene-wise tests. We sequenced 98 genes in 588 Italian patients with MS and 408 matched healthy controls (HCs). Variants were selected using different filtering criteria based on allelic frequency and in silico functional impacts. Genes showing a significant burden (n = 17) were sequenced in an independent cohort of 504 MS and 504 HC. The highest signal in both cohorts was observed for the disruptive variants (stop-gain, stop-loss, or splicing variants) located in EFCAB13, a gene coding for a protein of an unknown function (p < 10^-4). Among these variants, the minor allele of a stop-gain variant showed a significantly higher frequency in MS versus HC in both sequenced cohorts (p = 0.0093 and p = 0.025), confirmed by a meta-analysis on a third independent cohort of 1298 MS and 1430 HC (p = 0.001) assayed with an SNP array. Real-time PCR on 14 heterozygous individuals for this variant did not evidence the presence of the stop-gain allele, suggesting a transcript degradation by non-sense mediated decay, supported by the evidence that the carriers of the stop-gain variant had a lower expression of this gene (p = 0.0184). In conclusion, we identified a novel low-frequency functional variant associated with MS susceptibility, suggesting the possible role of rare/low-frequency variants in MS as reported for other complex diseases.

Genome-wide patterns of genetic diversity, population structure and demographic history in mānuka (Leptospermum scoparium) growing on indigenous Māori land

Leptospermum scoparium J. R. Forst et G. Forst, known as mānuka by Māori, the indigenous people of Aotearoa (New Zealand), is a culturally and economically significant shrub species, native to New Zealand and Australia. Chemical, morphological and phylogenetic studies have indicated geographical variation of mānuka across its range in New Zealand, and genetic differentiation between New Zealand and Australia. We used pooled whole genome re-sequencing of 76 L. scoparium and outgroup populations from New Zealand and Australia to compile a dataset totalling ~2.5 million SNPs. We explored the genetic structure and relatedness of L. scoparium across New Zealand, and between populations in New Zealand and Australia, as well as the complex demographic history of this species. Our population genomic investigation suggests there are five geographically distinct mānuka gene pools within New Zealand, with evidence of gene flow occurring between these pools. Demographic modelling suggests three of these gene pools have undergone expansion events, whilst the evolutionary histories of the remaining two have been subjected to contractions. Furthermore, mānuka populations in New Zealand are genetically distinct from populations in Australia, with coalescent modelling suggesting these two clades diverged ~9-12 million years ago. We discuss the evolutionary history of this species and the benefits of using pool-seq for such studies. Our research will support the management and conservation of mānuka by landowners, particularly Māori, and the development of a provenance story for the branding of mānuka based products.

Optimizing sampling design and sequencing strategy for the genomic analysis of quantitative traits in natural populations

Mapping the genes underlying ecologically relevant traits in natural populations is fundamental to develop a molecular understanding of species adaptation. Current sequencing technologies enable the characterization of a species' genetic diversity across the landscape or even over its whole range. The relevant capture of the genetic diversity across the landscape is critical for a successful genetic mapping of traits and there are no clear guidelines on how to achieve an optimal sampling and which sequencing strategy to implement. Here we determine, through simulation, the sampling scheme that maximizes the power to map the genetic basis of a complex trait in an outbreeding species across an idealized landscape and draw genomic predictions for the trait, comparing individual and pool sequencing strategies. Our results show that quantitative trait locus detection power and prediction accuracy are higher when more populations over the landscape are sampled and this is more cost-effectively done with pool sequencing than with individual sequencing. Additionally, we recommend sampling populations from areas of high genetic diversity. As progress in sequencing enables the integration of trait-based functional ecology into landscape genomics studies, these findings will guide study designs allowing direct measures of genetic effects in natural populations across the environment.

Comparison of anadromous and landlocked Atlantic salmon genomes reveals signatures of parallel and relaxed selection across the Northern Hemisphere

Most Atlantic salmon (Salmo salar L.) populations follow an anadromous life cycle, spending early life in freshwater, migrating to the sea for feeding, and returning to rivers to spawn. At the end of the last ice age ~10,000 years ago, several populations of Atlantic salmon became landlocked. Comparing their genomes to their anadromous counterparts can help identify genetic variation related to either freshwater residency or anadromy. The objective of this study was to identify consistently divergent loci between anadromous and landlocked Atlantic salmon strains throughout their geographical distribution, with the long-term aim of identifying traits relevant for salmon aquaculture, including fresh and seawater growth, omega-3 metabolism, smoltification, and disease resistance. We used a Pool-seq approach (n = 10-40 individuals per population) to sequence the genomes of twelve anadromous and six landlocked Atlantic salmon populations covering a large part of the Northern Hemisphere and conducted a genomewide association study to identify genomic regions having been under different selection pressure in landlocked and anadromous strains. A total of 28 genomic regions were identified and included cadm1 on Chr 13 and ppargc1a on Chr 18. Seven of the regions additionally displayed consistently reduced heterozygosity in fish obtained from landlocked populations, including the genes gpr132, cdca4, and sertad2 on Chr 15. We also found 16 regions, including igf1 on Chr 17, which consistently display reduced heterozygosity in the anadromous populations compared to the freshwater populations, indicating relaxed selection on traits associated with anadromy in landlocked salmon. In conclusion, we have identified 37 regions which may harbor genetic variation relevant for improving fish welfare and quality in the salmon farming industry and for understanding life-history traits in fish.

Independent Origin of XY and ZW Sex Determination Mechanisms in Mosquitofish Sister Species

Fish are known for the outstanding variety of their sex determination mechanisms and sex chromosome systems. The western (Gambusia affinis) and eastern mosquitofish (G. holbrooki) are sister species for which different sex determination mechanisms have been described: ZZ/ZW for G. affinis and XX/XY for G. holbrooki Here, we carried out restriction-site associated DNA (RAD-) and pool sequencing (Pool-seq) to characterize the sex chromosomes of both species. We found that the ZW chromosomes of G. affinis females and the XY chromosomes of G. holbrooki males correspond to different linkage groups, and thus evolved independently from separate autosomes. In interspecific hybrids, the Y chromosome is dominant over the W chromosome, and X is dominant over Z. In G. holbrooki, we identified a candidate region for the Y-linked melanic pigmentation locus, a rare male phenotype that constitutes a potentially sexually antagonistic trait and is associated with other such characteristics, e.g., large body size and aggressive behavior. We developed a SNP-based marker in the Y-linked allele of GIPC PDZ domain containing family member 1 (gipc1), which was linked to melanism in all tested G. holbrooki populations. This locus represents an example for a color locus that is located in close proximity to a putative sex determiner, and most likely substantially contributed to the evolution of the Y.

Genomic dissection of pod shattering in common bean: mutations at non-orthologous loci at the basis of convergent phenotypic evolution under domestication of leguminous species

The complete or partial loss of shattering ability occurred independently during the domestication of several crops. Therefore, the study of this trait can provide an understanding of the link between phenotypic and molecular convergent evolution. The genetic dissection of 'pod shattering' in Phaseolus vulgaris is achieved here using a population of introgression lines and next-generation sequencing techniques. The 'occurrence' of the indehiscent phenotype (indehiscent versus dehiscent) depends on a major locus on chromosome 5. Furthermore, at least two additional genes are associated with the 'level' of shattering (number of shattering pods per plant: low versus high) and the 'mode' of shattering (non-twisting versus twisting pods), with all of these loci contributing to the phenotype by epistatic interactions. Comparative mapping indicates that the major gene identified on common bean chromosome 5 corresponds to one of the four quantitative trait loci for pod shattering in Vigna unguiculata. None of the loci identified comprised genes that are homologs of the known shattering genes in Glycine max. Therefore, although convergent domestication can be determined by mutations at orthologous loci, this was only partially true for P. vulgaris and V. unguiculata, which are two phylogenetically closely related crop species, and this was not the case for the more distant P. vulgaris and G. max. Conversely, comparative mapping suggests that the convergent evolution of the indehiscent phenotype arose through mutations in different genes from the same underlying gene networks that are involved in secondary cell-wall biosynthesis and lignin deposition patterning at the pod level.

Measuring Genetic Differentiation from Pool-seq Data

The advent of high throughput sequencing and genotyping technologies enables the comparison of patterns of polymorphisms at a very large number of markers. While the characterization of genetic structure from individual sequencing data remains expensive for many nonmodel species, it has been shown that sequencing pools of individual DNAs (Pool-seq) represents an attractive and cost-effective alternative. However, analyzing sequence read counts from a DNA pool instead of individual genotypes raises statistical challenges in deriving correct estimates of genetic differentiation. In this article, we provide a method-of-moments estimator of [Formula: see text] for Pool-seq data, based on an analysis-of-variance framework. We show, by means of simulations, that this new estimator is unbiased and outperforms previously proposed estimators. We evaluate the robustness of our estimator to model misspecification, such as sequencing errors and uneven contributions of individual DNAs to the pools. Finally, by reanalyzing published Pool-seq data of different ecotypes of the prickly sculpin Cottus asper, we show how the use of an unbiased [Formula: see text] estimator may question the interpretation of population structure inferred from previous analyses.