Genetic population structure and relatedness in the narrow‐striped mongoose (Mungotictis decemlineata), a social Malagasy carnivore with sexual segregation

Information on the genetic structure of animal populations can allow inferences about mechanisms shaping their social organization, dispersal, and mating system. The mongooses (Herpestidae) include some of the best‐studied mammalian systems in this respect, but much less is known about their closest relatives, the Malagasy carnivores (Eupleridae), even though some of them exhibit unusual association patterns. We investigated the genetic structure of the Malagasy narrow‐striped mongoose (Mungotictis decemlineata), a small forest‐dwelling gregarious carnivore exhibiting sexual segregation. Based on mtDNA and microsatellite analyses, we determined population‐wide haplotype structure and sex‐specific and within‐group relatedness. Furthermore, we analyzed parentage and sibship relationships and the level of reproductive skew. We found a matrilinear population structure, with several neighboring female units sharing identical haplotypes. Within‐group female relatedness was significantly higher than expected by chance in the majority of units. Haplotype diversity of males was significantly higher than in females, indicating male‐biased dispersal. Relatedness within the majority of male associations did not differ from random, not proving any kin‐directed benefits of male sociality in this case. We found indications for a mildly promiscuous mating system without monopolization of females by males, and low levels of reproductive skew in both sexes based on parentages of emergent young. Low relatedness within breeding pairs confirmed immigration by males and suggested similarities with patterns in social mongooses, providing a starting point for further investigations of mate choice and female control of reproduction and the connected behavioral mechanisms. Our study contributes to the understanding of the determinants of male sociality in carnivores as well as the mechanisms of female competition in species with small social units.     

Cross‐platform compatibility of de novo‐aligned SNPs in a nonmodel butterfly genus

High‐throughput sequencing methods for genotyping genome‐wide markers are being rapidly adopted for phylogenetics of nonmodel organisms in conservation and biodiversity studies. However, the reproducibility of SNP genotyping and degree of marker overlap or compatibility between datasets from different methodologies have not been tested in nonmodel systems. Using double‐digest restriction site‐associated DNA sequencing, we sequenced a common set of 22 specimens from the butterfly genus Speyeria on two different Illumina platforms, using two variations of library preparation. We then used a de novo approach to bioinformatic locus assembly and SNP discovery for subsequent phylogenetic analyses. We found a high rate of locus recovery despite differences in library preparation and sequencing platforms, as well as overall high levels of data compatibility after data processing and filtering. These results provide the first application of NGS methods for phylogenetic reconstruction in Speyeria and support the use and long‐term viability of SNP genotyping applications in nonmodel systems.     

Polymorphisms in a desaturase 2 ortholog associate with cuticular hydrocarbon and male mating success variation in a natural population of Drosophila serrata

Elucidating the nature of genetic variation underlying both sexually selected traits and the fitness components of sexual selection is essential to understanding the broader consequences of sexual selection as an evolutionary process. To date, there have been relatively few attempts to connect the genetic variance in sexually selected traits with segregating DNA sequence polymorphisms. We set out to address this in a well‐characterized sexual selection system – the cuticular hydrocarbons (CHCs) of Drosophila serrata – using an indirect association study design that allowed simultaneous estimation of the genetic variance in CHCs, sexual fitness and single nucleotide polymorphism (SNP) effects in an outbred population. We cloned and sequenced an ortholog of the D. melanogaster desaturase 2 gene, previously shown to affect CHC biosynthesis in D. melanogaster, and associated 36 SNPs with minor allele frequencies > 0.02 with variance in CHCs and sexual fitness. Three SNPs had significant multivariate associations with CHC phenotype (q‐value < 0.05). At these loci, minor alleles had multivariate effects on CHCs that were weakly associated with the multivariate direction of sexual selection operating on these traits. Two of these SNPs had pleiotropic associations with male mating success, suggesting these variants may underlie responses to sexual selection due to this locus. There were 15 significant male mating success associations (q‐value < 0.1), and interestingly, we detected a nonrandom pattern in the relationship between allele frequency and direction of effect on male mating success. The minor‐frequency allele usually reduced male mating success, suggesting a positive association between male mating success and total fitness at this locus.     

Uneven Distribution of Mating‐Type Alleles Among Togninia minima Isolates,One of the Causal Agents of Leaf Stripe Disease on Grapevines in Northwest Iran

Togninia minima is the main fungal species associated with grapevine leaf stripe disease worldwide. This species is mainly known from its asexual state in nature; nevertheless, a biallelic heterothallic mating strategy has been confirmed for this species based on in vitro crossing studies. There are no data available on the incidence of an active sexual cycle within the populations of this species in many grapevine‐producing countries as well as Iran. The possibility of a clandestine sexual cycle within the Iranian isolates of T. minima was evaluated by analysing the distribution and frequency of the mating‐type alleles on a microspatial and a macrogeographical scales. Towards this aim, a total of 90 T. minima isolates were recovered from grapevines with esca disease from the vineyards in north and north‐western Iran. A multiplex PCR method previously designed by authors was applied for simultaneous identification and determination of the mating‐type alleles in T. minima populations. The results on the screening of mating‐type alleles using multiplex PCR method revealed the mating‐type identity of 77 isolates as Mat1‐2 and 23 isolates as Mat1‐1. Our results showed that both Mat1‐1 and Mat1‐2 isolates are present in a single vineyard and even on single vines. The distribution of mating‐type alleles in the sampled area skewed from the 1 : 1 ratio (77 : 23); however, co‐occurrence of both mating types in a single vineyard and even on single vines is suggestive for the presence of an active sexual cycle for T. minima in north‐western Iran.     

Development of SNP‐genotyping arrays in two shellfish species

Use of SNPs has been favoured due to their abundance in plant and animal genomes, accompanied by the falling cost and rising throughput capacity for detection and genotyping. Here, we present in vitro (obtained from targeted sequencing) and in silico discovery of SNPs, and the design of medium‐throughput genotyping arrays for two oyster species, the Pacific oyster, Crassostrea gigas, and European flat oyster, Ostrea edulis. Two sets of 384 SNP markers were designed for two Illumina GoldenGate arrays and genotyped on more than 1000 samples for each species. In each case, oyster samples were obtained from wild and selected populations and from three‐generation families segregating for traits of interest in aquaculture. The rate of successfully genotyped polymorphic SNPs was about 60% for each species. Effects of SNP origin and quality on genotyping success (Illumina functionality Score) were analysed and compared with other model and nonmodel species. Furthermore, a simulation was made based on a subset of the C. gigas SNP array with a minor allele frequency of 0.3 and typical crosses used in shellfish hatcheries. This simulation indicated that at least 150 markers were needed to perform an accurate parental assignment. Such panels might provide valuable tools to improve our understanding of the connectivity between wild (and selected) populations and could contribute to future selective breeding programmes.     

Sperm dispersal distances estimated by parentage analysis in a brooding scleractinian coral

Within populations of brooding sessile corals, sperm dispersal constitutes the mechanism by which gametes interact and mating occurs, and forms the first link in the network of processes that determine specieswide connectivity patterns. However, almost nothing is known about sperm dispersal for any internally fertilizing coral. In this study, we conducted a parentage analysis on coral larvae collected from an area of mapped colonies, to measure the distance sperm disperses for the first time in a reef‐building coral and estimated the mating system characteristics of a recently identified putative cryptic species within the Seriatopora hystrix complex (ShA; Warner et al. 2015). We defined consensus criteria among several replicated methods (colony 2.0, cervus 3.0, mltr v3.2) to maximize accuracy in paternity assignments. Thirteen progeny arrays indicated that this putative species produces exclusively sexually derived, primarily outcrossed larvae (mean t_m = 0.999) in multiple paternity broods (mean r_p = 0.119). Self‐fertilization was directly detected at low frequency for all broods combined (2.8%), but comprised 23% of matings in one brood. Although over 82% of mating occurred between colonies within 10 m of each other (mean sperm dispersal = 5.5 m ± 4.37 SD), we found no evidence of inbreeding in the established population. Restricted dispersal of sperm compared to slightly greater larval dispersal appears to limit inbreeding among close relatives in this cryptic species. Our findings establish a good basis for further work on sperm dispersal in brooding corals and provide the first information about the mating system of a newly identified and abundant cryptic species.     

CropSNPdb: a database of SNP array data for Brassica crops and hexaploid bread wheat

Advances in sequencing technology have led to a rapid rise in the genomic data available for plants, driving new insights into the evolution, domestication and improvement of crops. Single nucleotide polymorphisms (SNPs) are a major component of crop genomic diversity, and are invaluable as genetic markers in research and breeding programs. High‐throughput SNP arrays, or ‘SNP chips’, can generate reproducible sets of informative SNP markers and have been broadly adopted. Although there are many public repositories for sequencing data, which are routinely uploaded, there are no formal repositories for crop SNP array data. To make SNP array data more easily accessible, we have developed CropSNPdb ( http://snpdb.appliedbioinformatics.com.au ), a database for SNP array data produced by the Illumina Infinium? hexaploid bread wheat (Triticum aestivum) 90K and Brassica 60K arrays. We currently host SNPs from datasets covering 526 Brassica lines and 309 bread wheat lines, and provide search, download and upload utilities for users. CropSNPdb provides a useful repository for these data, which can be applied for a range of genomics and molecular crop‐breeding activities.     

Mapping DNA damage‐dependent genetic interactions in yeast via party mating and barcode fusion genetics

Condition‐dependent genetic interactions can reveal functional relationships between genes that are not evident under standard culture conditions. State‐of‐the‐art yeast genetic interaction mapping, which relies on robotic manipulation of arrays of double‐mutant strains, does not scale readily to multi‐condition studies. Here, we describe barcode fusion genetics to map genetic interactions (BFG‐GI), by which double‐mutant strains generated via en masse “party” mating can also be monitored en masse for growth to detect genetic interactions. By using site‐specific recombination to fuse two DNA barcodes, each representing a specific gene deletion, BFG‐GI enables multiplexed quantitative tracking of double mutants via next‐generation sequencing. We applied BFG‐GI to a matrix of DNA repair genes under nine different conditions, including methyl methanesulfonate (MMS), 4‐nitroquinoline 1‐oxide (4NQO), bleomycin, zeocin, and three other DNA‐damaging environments. BFG‐GI recapitulated known genetic interactions and yielded new condition‐dependent genetic interactions. We validated and further explored a subnetwork of condition‐dependent genetic interactions involving MAG1, SLX4, and genes encoding the Shu complex, and inferred that loss of the Shu complex leads to an increase in the activation of the checkpoint protein kinase Rad53.     

Development of high‐density SNP genotyping arrays for white spruce (Picea glauca) and transferability to subtropical and nordic congeners

High‐density SNP genotyping arrays can be designed for any species given sufficient sequence information of high quality. Two high‐density SNP arrays relying on the Infinium iSelect technology (Illumina) were designed for use in the conifer white spruce (Picea glauca). One array contained 7338 segregating SNPs representative of 2814 genes of various molecular functional classes for main uses in genetic association and population genetics studies. The other one contained 9559 segregating SNPs representative of 9543 genes for main uses in population genetics, linkage mapping of the genome and genomic prediction. The SNPs assayed were discovered from various sources of gene resequencing data. SNPs predicted from high‐quality sequences derived from genomic DNA reached a genotyping success rate of 64.7%. Nonsingleton in silico SNPs (i.e. a sequence polymorphism present in at least two reads) predicted from expressed sequenced tags obtained with the Roche 454 technology and Illumina GAII analyser resulted in a similar genotyping success rate of 71.6% when the deepest alignment was used and the most favourable SNP probe per gene was selected. A variable proportion of these SNPs was shared by other nordic and subtropical spruce species from North America and Europe. The number of shared SNPs was inversely proportional to phylogenetic divergence and standing genetic variation in the recipient species, but positively related to allele frequency in P. glauca natural populations. These validated SNP resources should open up new avenues for population genetics and comparative genetic mapping at a genomic scale in spruce species.     

Ongoing hybridization obscures phylogenetic relationships in the Drosophila subquinaria species complex

Inferring evolutionary relationships among recently diverged lineages is necessary to understand how isolating barriers produce independent lineages. Here, we investigate the phylogenetic relationships between three incompletely isolated and closely related mushroom‐feeding Drosophila species. These species form the Drosophila subquinaria species complex and consist of one Eurasian species (D. transversa) and two widespread North American species (D. subquinaria and D. recens) that are sympatric in central Canada. Although patterns of pre‐ and post‐mating isolation among these species are well characterized, previous work on their phylogenetic relationships is limited and conflicting. In this study, we generated a multi‐locus data set of 29 loci from across the genome sequenced in a population sample from each species, and then, we inferred species relationships and patterns of introgression. We find strong statistical support that D. subquinaria is paraphyletic, showing that samples from the geographic region sympatric with D. recens are most closely related to D. recens, whereas samples from the geographic region allopatric with D. recens are most closely related to D. transversa. We present several lines of evidence that both incomplete lineage sorting and gene flow are causing phylogenetic discordance. We suggest that ongoing gene flow primarily from D. recens into D. subquinaria in the sympatric part of their ranges causes phylogenetic uncertainty in the evolutionary history of these species. Our results highlight how population genetic data can be used to disentangle the sources of phylogenetic discordance among closely related species.     

Facultative pupal mating in Heliconius erato: Implications for mate choice,female preference,and speciation

Mating systems have broad impacts on how sexual selection and mate choice operate within a species, but studies of mating behavior in the laboratory may not reflect how these processes occur in the wild. Here, we examined the mating behavior of the neotropical butterfly Heliconius erato in the field by releasing larvae and virgin females and observing how they mated. H. erato is considered a pupal‐mating species (i.e., males mate with females as they emerge from the pupal case). However, we observed only two teneral mating events, and experimentally released virgins were almost all mated upon recapture. Our study confirms the presence of some pupal‐mating behavior in H. erato, but suggests that adult mating is likely the prevalent mating strategy in this species. These findings have important implications for the role of color pattern and female mate choice in the generation of reproductive isolation in this diverse genus.     

Genomewide single nucleotide polymorphism discovery in Atlantic salmon (Salmo salar): validation in wild and farmed American and European populations

A considerable number of single nucleotide polymorphisms (SNPs) are required to elucidate genotype–phenotype associations and determine the molecular basis of important traits. In this work, we carried out de novo SNP discovery accounting for both genome duplication and genetic variation from American and European salmon populations. A total of 9 736 473 nonredundant SNPs were identified across a set of 20 fish by whole‐genome sequencing. After applying six bioinformatic filtering steps, 200 K SNPs were selected to develop an Affymetrix Axiom^® myDesign Custom Array. This array was used to genotype 480 fish representing wild and farmed salmon from Europe, North America and Chile. A total of 159 099 (79.6%) SNPs were validated as high quality based on clustering properties. A total of 151 509 validated SNPs showed a unique position in the genome. When comparing these SNPs against 238 572 markers currently available in two other Atlantic salmon arrays, only 4.6% of the SNP overlapped with the panel developed in this study. This novel high‐density SNP panel will be very useful for the dissection of economically and ecologically relevant traits, enhancing breeding programmes through genomic selection as well as supporting genetic studies in both wild and farmed populations of Atlantic salmon using high‐resolution genomewide information.     

Genotype‐free estimation of allele frequencies reduces bias and improves demographic inference from RADSeq data

Restriction‐site associated DNA sequencing (RADSeq) facilitates rapid generation of thousands of genetic markers at relatively low cost; however, several sources of error specific to RADSeq methods often lead to biased estimates of allele frequencies and thereby to erroneous population genetic inference. Estimating the distribution of sample allele frequencies without calling genotypes was shown to improve population inference from whole genome sequencing data, but the ability of this approach to account for RADSeq‐specific biases remains unexplored. Here we assess in how far genotype‐free methods of allele frequency estimation affect demographic inference from empirical RADSeq data. Using the well‐studied pied flycatcher (Ficedula hypoleuca) as a study system, we compare allele frequency estimation and demographic inference from whole genome sequencing data with that from RADSeq data matched for samples using both genotype‐based and genotype free methods. The demographic history of pied flycatchers as inferred from RADSeq data was highly congruent with that inferred from whole genome resequencing (WGS) data when allele frequencies were estimated directly from the read data. In contrast, when allele frequencies were derived from called genotypes, RADSeq‐based estimates of most model parameters fell outside the 95% confidence interval of estimates derived from WGS data. Notably, more stringent filtering of the genotype calls tended to increase the discrepancy between parameter estimates from WGS and RADSeq data, respectively. The results from this study demonstrate the ability of genotype‐free methods to improve allele frequency spectrum‐ (AFS‐) based demographic inference from empirical RADSeq data and highlight the need to account for uncertainty in NGS data regardless of sequencing method.     

An empirical comparison of character‐based and coalescent‐based approaches to species delimitation in a young avian complex

The process of discovering species is a fundamental responsibility of systematics. Recently, there has been a growing interest in coalescent‐based methods of species delimitation aimed at objectively identifying species early in the divergence process. However, few empirical studies have compared these new methods with character‐based approaches for discovering species. In this study, we applied both a character‐based and a coalescent‐based approaches to delimit species in a closely related avian complex, the light‐vented/Taiwan bulbul (Pycnonotus sinensis/Pycnonotus taivanus). Population aggregation analyses of plumage, mitochondrial and 13 nuclear intron character data sets produced conflicting species hypotheses with plumage data suggesting three species, mitochondrial data suggesting two species, and nuclear intron data suggesting one species. Such conflict is expected among recently diverged species, and by integrating all sources of data, we delimited three species verified with independently congruent character evidence as well as a more weakly supported fourth species identified by a single character. Attempts to validate species hypothesis using Bayesian Phylogenetics and Phylogeography (BPP), a coalescent‐based method of species delimitation, revealed several issues that can seemingly affect statistical support for species recognition. We found that θ priors had a dramatic impact on speciation probabilities, with lower values consistently favouring splitting and higher values consistently favouring lumping. More resolved guide trees also resulted in overall higher speciation probabilities. Finally, we found suggestive evidence that BPP is sensitive to the divergent effects of nonrandom mating caused by intraspecific processes such as isolation‐with‐distance, and therefore, BPP may not be a conservative method for delimiting independently evolving population lineages. Based on these concerns, we questioned the reliability of BPP results and based our conclusions about species limits exclusively on character data.     

Introgression browser: high‐throughput whole‐genome SNP visualization

Breeding by introgressive hybridization is a pivotal strategy to broaden the genetic basis of crops. Usually, the desired traits are monitored in consecutive crossing generations by marker‐assisted selection, but their analyses fail in chromosome regions where crossover recombinants are rare or not viable. Here, we present the Introgression Browser (iBrowser ), a bioinformatics tool aimed at visualizing introgressions at nucleotide or SNP (Single Nucleotide Polymorphisms) accuracy. The software selects homozygous SNPs from Variant Call Format (VCF) information and filters out heterozygous SNPs, multi‐nucleotide polymorphisms (MNPs) and insertion–deletions (InDels). For data analysis iBrowser makes use of sliding windows, but if needed it can generate any desired fragmentation pattern through General Feature Format (GFF) information. In an example of tomato (Solanum lycopersicum) accessions we visualize SNP patterns and elucidate both position and boundaries of the introgressions. We also show that our tool is capable of identifying alien DNA in a panel of the closely related S. pimpinellifolium by examining phylogenetic relationships of the introgressed segments in tomato. In a third example, we demonstrate the power of the iBrowser in a panel of 597 Arabidopsis accessions, detecting the boundaries of a SNP‐free region around a polymorphic 1.17 Mbp inverted segment on the short arm of chromosome 4. The architecture and functionality of iBrowser makes the software appropriate for a broad set of analyses including SNP mining, genome structure analysis, and pedigree analysis. Its functionality, together with the capability to process large data sets and efficient visualization of sequence variation, makes iBrowser a valuable breeding tool.     

Assessment of genetic diversity and population structure in cultured Australian Pacific oysters

The recent development of Pacific oyster (Crassostrea gigas) SNP genotyping arrays has allowed detailed characterisation of genetic diversity and population structure within and between oyster populations. It also raises the potential of harnessing genomic selection for genetic improvement in oyster breeding programmes. The aim of this study was to characterise a breeding population of Australian oysters through genotyping and analysis of 18 027 SNPs, followed by comparison with genotypes of oyster sampled from Europe and Asia. This revealed that the Australian populations had similar population diversity (H_E) to oysters from New Zealand, the British Isles, France and Japan. Population divergence was assessed using PCA of genetic distance and revealed that Australian oysters were distinct from all other populations tested. Australian Pacific oysters originate from planned introductions sourced from three Japanese populations. Approximately 95% of these introductions were from geographically, and potentially genetically, distinct populations from the Nagasaki oysters assessed in this study. Finally, in preparation for the application of genomic selection in oyster breeding programmes, the strength of LD was evaluated and subsets of loci were tested for their ability to accurately infer relationships. Weak LD was observed on average; however, SNP subsets were shown to accurately reconstitute a genomic relationship matrix constructed using all loci. This suggests that low‐density SNP panels may have utility in the Australian population tested, and the findings represent an important first step towards the design and implementation of genomic approaches for applied breeding in Pacific oysters.     

A genome‐wide significant association on chromosome 2 for footrot resistance/susceptibility in Swiss White Alpine sheep

Footrot is one of the most important causes of lameness in global sheep populations and is characterized by a bacterial infection of the interdigital skin. As a multifactorial disease, its clinical representation depends not only on pathogen factors and environmental components but also on the individual resistance/susceptibility of the host. A genetic component has been shown in previous studies; however, so far no causative genetic variant influencing the risk of developing footrot has been identified. In this study, we genotyped 373 Swiss White Alpine sheep, using the ovine high‐density 600k SNP chip, in order to run a DNA‐based comparison of individuals with known clinical footrot status. We performed a case–control genome‐wide association study, which revealed a genome‐wide significant association for SNP rs418747104 on ovine chromosome 2 at 81.2 Mb. The three best associated SNP markers were located at the MPDZ gene, which codes for the multiple PDZ domain crumbs cell polarity complex component protein, also known as multi‐PDZ domain protein 1 (MUPP1). This protein is possibly involved in maintaining the barrier function and integrity of tight junctions. Therefore, we speculate that individuals carrying MPDZ variants may differ in their footrot resistance/susceptibility due to modified horn and interdigital skin integrity. In conclusion, our study reveals that MPDZ might represent a functional candidate gene, and further research is needed to explore its role in footrot affected sheep.     

Comparison of environmental inference approaches for ecometric analyses: Using hypsodonty to estimate precipitation

1. Ecometrics is the study of community‐level functional trait–environment relationships. We use ecometric analyses to estimate paleoenvironment and to investigate community‐level functional changes through time.
2. We evaluate four methods that have been used or have the potential to be used in ecometric analyses for estimating paleoenvironment to determine whether there have been systematic differences in paleoenvironmental estimation due to choice of the estimation method. Specifically, we evaluated linear regression, polynomial regression, nearest neighbor, and maximum‐likelihood methods to explore the predictive ability of the relationship for a well‐known ecometric dataset of mammalian herbivore hypsodonty metrics (molar tooth crown to root height ratio) and annual precipitation. Each method was applied to 43 Pleistocene fossil sites and compared to annual precipitation from global climate models. Sites were categorized as glacial or interglacial, and paleoprecipitation estimates were compared to the appropriate model.
3. Estimation methods produce results that are highly correlated with log precipitation and estimates from the other methods (p < 0.001). Differences between estimated precipitation and observed precipitation are not significantly different across the four methods, but maximum likelihood produces the most accurate estimates of precipitation. When applied to paleontological sites, paleoprecipitation estimates align more closely with glacial global climate models than with interglacial models regardless of the age of the site.
4. Each method has constraints that are important to consider when designing ecometric analyses to avoid misinterpretations when ecometric relationships are applied to the paleontological record. We show interglacial fauna estimates of paleoprecipitation more closely match glacial global climate models. This is likely because of the anthropogenic effects on community reassembly in the Holocene.