Construction of relatedness matrices using genotyping-by-sequencing data

Rumen metagenome profiles are heritable and rank the New Zealand national sheep flock for enteric methane emissions

BACKGROUND

Global targets to reduce greenhouse gas emissions to meet international climate change commitments have driven the livestock industry to develop solutions to reduce methane emission in ruminants while maintaining production. Research has shown that selective breeding for low methane emitting ruminants using genomic selection is one viable solution to meet methane targets at a national level. However, this requires obtaining sufficient measures of methane on individual animals across the national herd. In sheep, one affordable method for measuring methane on-farm to rank animals on their methane emissions is portable accumulation chambers (PAC), although this method is not without its challenges. An alternative is to use a proxy trait that is genetically correlated with PAC methane measures. One such trait that has shown promise is rumen metagenome community (RMC) profiles. In this study, we investigate the potential of using RMC profiles as a proxy trait for methane emissions from PAC using a large sheep dataset consisting of 4585 mixed-sex lambs from several flocks and years across New Zealand.

RESULTS

RMC profiles were generated from rumen samples collected on the animals immediately after being measured through PAC using restriction enzyme-reduced representation sequencing. We predicted methane (CH4) and carbon dioxide (CO2) emissions (grams per day), as well as the ratio CH4/(CO2 + CH4) (CH4Ratio), from the RMC profiles and SNP-array genotype data. Heritability and microbiability estimates were similar to values found in the literature for all traits. The correlation of PAC methane with predicted methane was 1.9- to 2.3-fold (CH4) and 1.2- to 1.5-fold (CH4Ratio) greater for RMC profiles compared to host genomics only. The genetic correlation between methane predicted from RMC profiles and PAC methane was 0.75 ± 0.12 for CH4 and 0.64 ± 0.11 for CH4Ratio when using a validation set consisting of the animals with the most recent year of birth in the dataset.

CONCLUSIONS

RMC profiles are predictive of, and genetically correlated, with PAC methane measures. Therefore, RMC profiles are a suitable proxy trait for determining the genetic merit of an animal's methane emissions and could be incorporated into existing breeding programs to facilitate selective breeding for low methane emitting sheep.

Milk fatty acids and rumen composition as proxy measures of enteric methane

Genomic selection is a powerful tool to reduce methane emissions in ruminants. However, it requires large-scale on-farm phenotypic measures of methane. Current technologies to measure methane emissions have several limitations and may not be suitable for lactating animals. Because enteric methane is closely linked to the fermentation process in the rumen, which in turn affects milk composition, breeding for low-methane ruminants may change the rumen microbial and milk composition. Consequently, these compositions may provide proxy measures of methane for use in selective breeding of low-methane ruminants. We investigated the effect on rumen and milk composition in sheep bred for divergent methane yield and the potential for generating proxy measures of methane emissions from rumen or milk samples in lactating ewes. Four hundred genotyped lactating ewes from a sheep research flock bred specifically for high and low-methane emissions had methane measured and rumen and milk samples collected approximately 6 wk post-lambing across 4 lactation years. Rumen samples were processed to generate VFA and metagenomic profiles of the rumen microbial community, and fatty acid profiles and mid-infrared spectra were generated for the milk samples. Although no differences in total fat, protein, and lactose percentages in the milk were found, the milk fatty acid profiles differed between methane selection lines, with higher PUFA and branched-chain fatty acids levels, and lower total SFA contents in ewes from the low-methane line. Lower proportions of acetate relative to propionate were found in the rumen samples from the low-methane ewes. Predictions of methane were obtained from the rumen VFA and metagenomic profiles and the fatty acid profiles and mid-infrared spectra from milk. These predictions formed the proxy methane measures and were heritable (between 0.12 to 0.36) and correlated (between 0.29 and 0.42) with the measured methane values. The genetic correlation between proxies and measured methane was between 0.52 and 0.71. The estimated efficiency of indirect selection for methane was higher for the milk sample proxies (49%-75%) than the rumen metagenomic profiles (45%-47%) and rumen VFA profiles (12%-38%). These results suggest that milk fatty acid, MIR spectroscopic, and rumen microbial composition phenotypes have the potential to be used as proxy measures of methane in lactating ruminants, with the milk-based proxies showing greater promise. Results show that the number of animals with methane proxy measures could be increased substantially and will enable access to breeding technology in countries with limited methane measurement infrastructure.

Single-cell eQTL mapping in yeast reveals a tradeoff between growth and reproduction

Expression quantitative trait loci (eQTLs) provide a key bridge between noncoding DNA sequence variants and organismal traits. The effects of eQTLs can differ among tissues, cell types, and cellular states, but these differences are obscured by gene expression measurements in bulk populations. We developed a one-pot approach to map eQTLs in Saccharomyces cerevisiae by single-cell RNA sequencing (scRNA-seq) and applied it to over 100,000 single cells from three crosses. We used scRNA-seq data to genotype each cell, measure gene expression, and classify the cells by cell-cycle stage. We mapped thousands of local and distant eQTLs and identified interactions between eQTL effects and cell-cycle stages. We took advantage of single-cell expression information to identify hundreds of genes with allele-specific effects on expression noise. We used cell-cycle stage classification to map 20 loci that influence cell-cycle progression. One of these loci influenced the expression of genes involved in the mating response. We showed that the effects of this locus arise from a common variant (W82R) in the gene GPA1, which encodes a signaling protein that negatively regulates the mating pathway. The 82R allele increases mating efficiency at the cost of slower cell-cycle progression and is associated with a higher rate of outcrossing in nature. Our results provide a more granular picture of the effects of genetic variants on gene expression and downstream traits.

Genomic selection shows improved expected genetic gain over phenotypic selection of agronomic traits in allotetraploid white clover

KEY MESSAGE

Genomic selection using white clover multi-year-multi-site data showed predicted genetic gains through integrating among-half-sibling-family phenotypic selection and within-family genomic selection were up to 89% greater than half-sibling-family phenotypic selection alone. Genomic selection, an effective breeding tool used widely in plants and animals for improving low-heritability traits, has only recently been applied to forages. We explored the feasibility of implementing genomic selection in white clover (Trifolium repens L.), a key forage legume which has shown limited genetic improvement in dry matter yield (DMY) and persistence traits. We used data from a training population comprising 200 half-sibling (HS) families evaluated in a cattle-grazed field trial across three years and two locations. Combining phenotype and genotyping-by-sequencing (GBS) data, we assessed different two-stage genomic prediction models, including KGD-GBLUP developed for low-depth GBS data, on DMY, growth score, leaf size and stolon traits. Predictive abilities were similar among the models, ranging from -0.17 to 0.44 across traits, and remained stable for most traits when reducing model input to 100-120 HS families and 5500 markers, suggesting genomic selection is viable with fewer resources. Incorporating a correlated trait with a primary trait in multi-trait prediction models increased predictive ability by 28-124%. Deterministic modelling showed integrating among-HS-family phenotypic selection and within-family genomic selection at different selection pressures estimated up to 89% DMY genetic gain compared to phenotypic selection alone, despite a modest predictive ability of 0.3. This study demonstrates the potential benefits of combining genomic and phenotypic selection to boost genetic gains in white clover. Using cost-effective GBS paired with a prediction model optimized for low read-depth data, the approach can achieve prediction accuracies comparable to traditional models, providing a viable path for implementing genomic selection in white clover.

Single-cell eQTL mapping in yeast reveals a tradeoff between growth and reproduction

Expression quantitative trait loci (eQTLs) provide a key bridge between noncoding DNA sequence variants and organismal traits. The effects of eQTLs can differ among tissues, cell types, and cellular states, but these differences are obscured by gene expression measurements in bulk populations. We developed a one-pot approach to map eQTLs in Saccharomyces cerevisiae by single-cell RNA sequencing (scRNA-seq) and applied it to over 100,000 single cells from three crosses. We used scRNA-seq data to genotype each cell, measure gene expression, and classify the cells by cell-cycle stage. We mapped thousands of local and distant eQTLs and identified interactions between eQTL effects and cell-cycle stages. We took advantage of single-cell expression information to identify hundreds of genes with allele-specific effects on expression noise. We used cell-cycle stage classification to map 20 loci that influence cell-cycle progression. One of these loci influenced the expression of genes involved in the mating response. We showed that the effects of this locus arise from a common variant (W82R) in the gene GPA1, which encodes a signaling protein that negatively regulates the mating pathway. The 82R allele increases mating efficiency at the cost of slower cell-cycle progression and is associated with a higher rate of outcrossing in nature. Our results provide a more granular picture of the effects of genetic variants on gene expression and downstream traits.

Genetic parameters and genotype-by-environment interaction estimates for growth and feed efficiency related traits in Chinook salmon, Oncorhynchus tshawytscha, reared under low and moderate flow regimes

BACKGROUND

A genotype-by-environment (G × E) interaction is defined as genotypes responding differently to different environments. In salmonids, G × E interactions can occur in different rearing conditions, including changes in salinity or temperature. However, water flow, an important variable that can influence metabolism, has yet to be considered for potential G × E interactions, although water flows differ across production stages. The salmonid industry is now manipulating flow in tanks to improve welfare and production performance, and expanding sea pen farming offshore, where flow dynamics are substantially greater. Therefore, there is a need to test whether G × E interactions occur under low and higher flow regimes to determine if industry should consider modifying their performance evaluation and selection criteria to account for different flow environments. Here, we used genotype-by-sequencing to create a genomic-relationship matrix of 37 Chinook salmon, Oncorhynchus tshawytscha, families to assess possible G × E interactions for production performance under two flow environments: a low flow regime (0.3 body lengths per second; bl s-1) and a moderate flow regime (0.8 bl s-1).

RESULTS

Genetic correlations for the same production performance trait between flow regimes suggest there is minimal evidence of a G × E interaction between the low and moderate flow regimes tested in this study, for Chinook salmon reared from 82.9 ± 16.8 g ( x ¯  ± s.d.) to 583.2 ± 117.1 g ( x ¯  ± s.d.). Estimates of genetic and phenotypic correlations between traits did not reveal any unfavorable trait correlations for size- (weight and condition factor) and growth-related traits, regardless of the flow regime, but did suggest measuring feed intake would be the preferred approach to improve feed efficiency because of the strong correlations between feed intake and feed efficiency, consistent with previous studies.

CONCLUSION

This new information suggests that Chinook salmon families do not need to be selected separately for performance across different flow regimes. However, further studies are needed to confirm this across a wider range of fish sizes and flows. This information is key for breeding programs to determine if separate evaluation groups are required for different flow regimes that are used for production (e.g., hatchery, post smolt recirculating aquaculture system, or offshore).

Fitness consequences of population bottlenecks in an invasive blowfly

Invasive species often undergo demographic bottlenecks that cause a decrease in genetic diversity and associated reductions in population fitness. Despite this, they manage to thrive in novel environments. Investigating the effects of inbreeding and genetic bottlenecks on population fitness for invasive species is, therefore, key to understanding how they may survive in new environments. We used the blowfly Calliphora vicina (Sciences, Mathématiques et Physique, 1830, 2, 1), which is native to Europe and was introduced to Australia and New Zealand, to examine the effects of genetic diversity on population fitness. We first collected 59 samples from 15 populations across New Zealand and one in Australia, and used 20,501 biallelic SNPs to investigate population genomic diversity, structure and admixture. We then explored the impacts of repeated experimental bottlenecks on population fitness by creating inbred and outbred lines of C. vicina and measuring a variety of fitness traits. In wild-caught samples, we found low overall genetic diversity, signals of genetic admixture and limited (<3%) genetic differentiation between North and South Island populations, with genetic links between the South Island and Australia. Following experimental bottlenecks, we found significant reductions in fitness for inbred lines. However, fitness effects were not felt equally across all phenotypic traits. Moreover, they were not enough to cause population collapse in any experimental line, suggesting that C. vicina (when under relaxed selection, as in laboratory settings) may be able to compensate for population bottlenecks even when highly inbred. Our results demonstrate the value of a tractable experimental system for investigating processes that may facilitate or hamper biological invasion.

Phenotypic diversity and population structure of Pecan (Carya illinoinensis) collections reveals geographic patterns

Pecan (Carya illinoinensis) is an economically important nut crop known for its genetic diversity and adaptability to various climates. Understanding the growth variability, phenological traits, and population structure of pecan populations is crucial for breeding programs and conservation. In this study, plant growth and phenological traits were evaluated over three consecutive seasons (2015-2017) for 550 genotypes from 26 provenances. Significant variations in plant height, stem diameter, and budbreak were observed among provenances, with Southern provenances exhibiting faster growth and earlier budbreak compared to Northern provenances. Population structure analysis using SNP markers revealed eight distinct subpopulations, reflecting genetic differentiation among provenances. Notably, Southern Mexico collections formed two separate clusters, while Western collections, such as 'Allen 3', 'Allen 4', and 'Riverside', were distinguished from others. 'Burkett' and 'Apache' were grouped together due to their shared maternal parentage. Principal component analysis and phylogenetic tree analysis further supported subpopulation differentiation. Genetic differentiation among the 26 populations was evident, with six clusters highly in agreement with the subpopulations identified by STRUCTURE and fastSTRUCTURE. Principal components analysis (PCA) revealed distinct groups, corresponding to subpopulations identified by genetic analysis. Discriminant analysis of PCA (DAPC) based on provenance origin further supported the genetic structure, with clear separation of provenances into distinct clusters. These findings provide valuable insights into the genetic diversity and growth patterns of pecan populations. Understanding the genetic basis of phenological traits and population structure is essential for selecting superior cultivars adapted to diverse environments. The identified subpopulations can guide breeding efforts to develop resilient rootstocks and contribute to the sustainable management of pecan genetic resources. Overall, this study enhances our understanding of pecan genetic diversity and informs conservation and breeding strategies for the long-term viability of pecan cultivation.

Coat colour in marsupials: genetic variants at the ASIP locus determine grey and black fur of the brushtail possum

The possession of fur or hair is a defining characteristic of mammals and can occur in a variety of colours and patterns. While genetic determinants of coat colour are well described in eutherian 'placental' mammals, the other major mammalian infraclass, marsupials, is grossly understudied. The fur of the common brushtail possum (Trichosurus vulpecula), an iconic native mammal found throughout Australia and introduced into Aotearoa New Zealand, possesses two main colour morphs: grey and black. To identify genetic variants associated with coat colour, we performed a genome-wide association study (GWAS) with genotype by sequencing (GBS) data. Single nucleotide variants (SNVs) on chromosome 3, close to the agouti signalling protein (ASIP) gene that controls the temporal and spatial distribution of pigments in eutherian mammals, were identified. Fine-mapping identified a C>T variant at chr3:100483705 that results in a ASIP:p.Arg115Cys missense substitution, and animals homozygous for this variant have black fur. In addition to uncovering the first genetic determinant of coat colour in a natural marsupial population, comparative analysis of ASIP in divergent marsupial species identified the dasyurids as having accelerated evolution, reflecting their well described diversity of coat colour and pattern.

Low-coverage sequencing and Wahlund effect severely bias estimates of inbreeding, heterozygosity and effective population size in North American wolves

vonHoldt et al. ((2024), Molecular Ecology, 33, e17231) (vH24) used low-coverage (average ~ 7X read depth) restriction site-associated DNA sequence data to estimate individual inbreeding and heterozygosity, and recent effective population size (Ne), in Great Lakes (GL) and Northern Rocky Mountain (RM) wolves. They concluded that RM heterozygosity rapidly declined between 1991 and 2020, and that Ne declined substantially in GL and RM over the last 50 generations. Here, we evaluate the effects of low sequence coverage and sampling strategy on vH24's findings and provide general recommendations for using sequence data to evaluate inbreeding, heterozygosity and Ne. Low-coverage sequencing resulted in downwardly biased estimates of individual inbreeding and heterozygosity, and an erroneous large temporal decline in RM heterozygosity due to declining read depth through time. Additionally, vH24's sampling strategy-which combined individuals from several genetically differentiated populations and across at least eight wolf generations-is expected to have resulted in severe downward bias in estimates of recent Ne for RM. We recommend using high-coverage sequence data (≥ $$ \ge $$ 15-20X) to estimate inbreeding and heterozygosity. Carefully filtering individuals, loci and genotypes, and using genotype imputation or likelihoods can help to minimise bias when low-coverage sequence data must be used. For estimation of contemporary Ne, the marginal benefits of using more than 103-104 loci are small, so aggressive filtering of loci with low average read depth potentially can retain most individuals without sacrificing much precision. Individuals are relatively more valuable than loci because analyses of contemporary Ne should focus on roughly single-generation samples from local breeding populations.

Vocal convergence and social proximity shape the calls of the most basal Passeriformes, New Zealand Wrens

Despite extensive research on avian vocal learning, we still lack a general understanding of how and when this ability evolved in birds. As the closest living relatives of the earliest Passeriformes, the New Zealand wrens (Acanthisitti) hold a key phylogenetic position for furthering our understanding of the evolution of vocal learning because they share a common ancestor with two vocal learners: oscines and parrots. However, the vocal learning abilities of New Zealand wrens remain unexplored. Here, we test for the presence of prerequisite behaviors for vocal learning in one of the two extant species of New Zealand wrens, the rifleman (Acanthisitta chloris). We detect the presence of unique individual vocal signatures and show how these signatures are shaped by social proximity, as demonstrated by group vocal signatures and strong acoustic similarities among distantly related individuals in close social proximity. Further, we reveal that rifleman calls share similar phenotypic variance ratios to those previously reported in the learned vocalizations of the zebra finch, Taeniopygia guttata. Together these findings provide strong evidence that riflemen vocally converge, and though the mechanism still remains to be determined, they may also suggest that this vocal convergence is the result of rudimentary vocal learning abilities.

Construction of relatedness matrices in autopolyploid populations using low-depth high-throughput sequencing data

KEY MESSAGE

An improved estimator of genomic relatedness using low-depth high-throughput sequencing data for autopolyploids is developed. Its outputs strongly correlate with SNP array-based estimates and are available in the package GUSrelate. High-throughput sequencing (HTS) methods have reduced sequencing costs and resources compared to array-based tools, facilitating the investigation of many non-model polyploid species. One important quantity that can be computed from HTS data is the genetic relatedness between all individuals in a population. However, HTS data are often messy, with multiple sources of errors (i.e. sequencing errors or missing parental alleles) which, if not accounted for, can lead to bias in genomic relatedness estimates. We derive a new estimator for constructing a genomic relationship matrix (GRM) from HTS data for autopolyploid species that accounts for errors associated with low sequencing depths, implemented in the R package GUSrelate. Simulations revealed that GUSrelate performed similarly to existing GRM methods at high depth but reduced bias in self-relatedness estimates when the sequencing depth was low. Using a panel consisting of 351 tetraploid potato genotypes, we found that GUSrelate produced GRMs from genotyping-by-sequencing (GBS) data that were highly correlated with a GRM computed from SNP array data, and less biased than existing methods when benchmarking against the array-based GRM estimates. GUSrelate provides researchers with a tool to reliably construct GRMs from low-depth HTS data.

Genome-Wide Association Study Revealed Putative SNPs and Candidate Genes Associated with Growth and Meat Traits in Japanese Quail

The search for SNPs and candidate genes that determine the manifestation of major selected traits is one crucial objective for genomic selection aimed at increasing poultry production efficiency. Here, we report a genome-wide association study (GWAS) for traits characterizing meat performance in the domestic quail. A total of 146 males from an F2 reference population resulting from crossing a fast (Japanese) and a slow (Texas White) growing breed were examined. Using the genotyping-by-sequencing technique, genomic data were obtained for 115,743 SNPs (92,618 SNPs after quality control) that were employed in this GWAS. The results identified significant SNPs associated with the following traits at 8 weeks of age: body weight (nine SNPs), daily body weight gain (eight SNPs), dressed weight (33 SNPs), and weights of breast (18 SNPs), thigh (eight SNPs), and drumstick (three SNPs). Also, 12 SNPs and five candidate genes (GNAL, DNAJC6, LEPR, SPAG9, and SLC27A4) shared associations with three or more traits. These findings are consistent with the understanding of the genetic complexity of body weight-related traits in quail. The identified SNPs and genes can be used in effective quail breeding as molecular genetic markers for growth and meat characteristics for the purpose of genetic improvement.

Population structure and interspecific hybridisation of two invasive blowflies (Diptera: Calliphoridae) following replicated incursions into New Zealand

Rates of biological invasion are increasing globally, with associated negative effects on native biodiversity and ecosystem services. Among other genetic processes, hybridisation can facilitate invasion by producing new combinations of genetic variation that increase adaptive potential and associated population fitness. Yet the role of hybridisation (and resulting gene flow) in biological invasion in invertebrate species is under-studied. Calliphora hilli and Calliphora stygia are blowflies proposed to have invaded New Zealand separately from Australia between 1779 and 1841, and are now widespread throughout the country. Here, we analysed genome-wide single nucleotide polymorphisms (SNPs), generating genotyping-by-sequencing data for 154 individuals collected from 24 populations across New Zealand and Australia to assess the extent of gene flow and hybridisation occurring within and between these blowflies and to better understand their overall population structure. We found that New Zealand populations of both species had weak genetic structure, suggesting high gene flow and an absence of dispersal limitations across the country. We also found evidence that interspecific hybridisation is occurring in the wild between C. hilli and C. stygia in both the native and invasive ranges, and that intraspecific admixture is occurring among populations at appreciable rates. Collectively, these findings provide new insights into the population structure of these two invasive invertebrates and highlight the potential importance of hybridisation and gene flow in biological invasion.

Restriction site associated DNA sequencing for tumour mutation burden estimation and mutation signature analysis

BACKGROUND

Genome-wide measures of genetic disruption such as tumour mutation burden (TMB) and mutation signatures are emerging as useful biomarkers to stratify patients for treatment. Clinicians commonly use cancer gene panels for tumour mutation burden estimation, and whole genome sequencing is the gold standard for mutation signature analysis. However, the accuracy and cost associated with these assays limits their utility at scale.

METHODS

WGS data from 560 breast cancer patients was used for in silico library simulations to evaluate the accuracy of an FDA approved cancer gene panel as well as restriction enzyme associated DNA sequencing (RADseq) libraries for TMB estimation and mutation signature analysis. We also transfected a mouse mammary cell line with APOBEC enzymes and sequenced resulting clones to evaluate the efficacy of RADseq in an experimental setting.

RESULTS

RADseq had improved accuracy of TMB estimation and derivation of mutation profiles when compared to the FDA approved cancer panel. Using simulated immune checkpoint blockade (ICB) trials, we show that inaccurate TMB estimation leads to a reduction in power for deriving an optimal TMB cutoff to stratify patients for immune checkpoint blockade treatment. Additionally, prioritisation of APOBEC hypermutated tumours in these trials optimises TMB cutoff determination for breast cancer. The utility of RADseq in an experimental setting was also demonstrated, based on characterisation of an APOBEC mutation signature in an APOBEC3A transfected mouse cell line.

CONCLUSION

In conclusion, our work demonstrates that RADseq has the potential to be used as a cost-effective, accurate solution for TMB estimation and mutation signature analysis by both clinicians and basic researchers.

Genotyping-by-Sequencing Strategy for Integrating Genomic Structure, Diversity and Performance of Various Japanese Quail (Coturnix japonica) Breeds

Traces of long-term artificial selection can be detected in genomes of domesticated birds via whole-genome screening using single-nucleotide polymorphism (SNP) markers. This study thus examined putative genomic regions under selection that are relevant to the development history, divergence and phylogeny among Japanese quails of various breeds and utility types. We sampled 99 birds from eight breeds (11% of the global gene pool) of egg (Japanese, English White, English Black, Tuxedo and Manchurian Golden), meat (Texas White and Pharaoh) and dual-purpose (Estonian) types. The genotyping-by-sequencing analysis was performed for the first time in domestic quails, providing 62,935 SNPs. Using principal component analysis, Neighbor-Net and Admixture algorithms, the studied breeds were characterized according to their genomic architecture, ancestry and direction of selective breeding. Japanese and Pharaoh breeds had the smallest number and length of homozygous segments indicating a lower selective pressure. Tuxedo and Texas White breeds showed the highest values of these indicators and genomic inbreeding suggesting a greater homozygosity. We revealed evidence for the integration of genomic and performance data, and our findings are applicable for elucidating the history of creation and genomic variability in quail breeds that, in turn, will be useful for future breeding improvement strategies.

Species-wide genomics of kākāpō provides tools to accelerate recovery

The kākāpō is a critically endangered, intensively managed, long-lived nocturnal parrot endemic to Aotearoa New Zealand. We generated and analysed whole-genome sequence data for nearly all individuals living in early 2018 (169 individuals) to generate a high-quality species-wide genetic variant callset. We leverage extensive long-term metadata to quantify genome-wide diversity of the species over time and present new approaches using probabilistic programming, combined with a phenotype dataset spanning five decades, to disentangle phenotypic variance into environmental and genetic effects while quantifying uncertainty in small populations. We find associations for growth, disease susceptibility, clutch size and egg fertility within genic regions previously shown to influence these traits in other species. Finally, we generate breeding values to predict phenotype and illustrate that active management over the past 45 years has maintained both genome-wide diversity and diversity in breeding values and, hence, evolutionary potential. We provide new pathways for informing future conservation management decisions for kākāpō, including prioritizing individuals for translocation and monitoring individuals with poor growth or high disease risk. Overall, by explicitly addressing the challenge of the small sample size, we provide a template for the inclusion of genomic data that will be transformational for species recovery efforts around the globe.

Large-scale analysis of sheep rumen metagenome profiles captured by reduced representation sequencing reveals individual profiles are influenced by the environment and genetics of the host

BACKGROUND

Producing animal protein while reducing the animal's impact on the environment, e.g., through improved feed efficiency and lowered methane emissions, has gained interest in recent years. Genetic selection is one possible path to reduce the environmental impact of livestock production, but these traits are difficult and expensive to measure on many animals. The rumen microbiome may serve as a proxy for these traits due to its role in feed digestion. Restriction enzyme-reduced representation sequencing (RE-RRS) is a high-throughput and cost-effective approach to rumen metagenome profiling, but the systematic (e.g., sequencing) and biological factors influencing the resulting reference based (RB) and reference free (RF) profiles need to be explored before widespread industry adoption is possible.

RESULTS

Metagenome profiles were generated by RE-RRS of 4,479 rumen samples collected from 1,708 sheep, and assigned to eight groups based on diet, age, time off feed, and country (New Zealand or Australia) at the time of sample collection. Systematic effects were found to have minimal influence on metagenome profiles. Diet was a major driver of differences between samples, followed by time off feed, then age of the sheep. The RF approach resulted in more reads being assigned per sample and afforded greater resolution when distinguishing between groups than the RB approach. Normalizing relative abundances within the sampling Cohort abolished structures related to age, diet, and time off feed, allowing a clear signal based on methane emissions to be elucidated. Genus-level abundances of rumen microbes showed low-to-moderate heritability and repeatability and were consistent between diets.

CONCLUSIONS

Variation in rumen metagenomic profiles was influenced by diet, age, time off feed and genetics. Not accounting for environmental factors may limit the ability to associate the profile with traits of interest. However, these differences can be accounted for by adjusting for Cohort effects, revealing robust biological signals. The abundances of some genera were consistently heritable and repeatable across different environments, suggesting that metagenomic profiles could be used to predict an individual's future performance, or performance of its offspring, in a range of environments. These results highlight the potential of using rumen metagenomic profiles for selection purposes in a practical, agricultural setting.

Traces of Human-Mediated Selection in the Gene Pool of Red Deer Populations

In this study, we analysed the effect of human-mediated selection on the gene pool of wild and farmed red deer populations based on genotyping-by-sequencing data. The farmed red deer sample covered populations spread across seven countries and two continents (France, Germany, Hungary, Latvia, New Zealand, Poland, and Slovakia). The Slovak and Spain wild red deer populations (the latter one in a large game estate) were used as control outgroups. The gene flow intensity, relationship and admixture among populations were tested by the Bayesian approach and discriminant analysis of principal components (DAPC). The highest gene diversity (He = 0.19) and the lowest genomic inbreeding (FHOM = 0.04) found in Slovak wild population confirmed our hypothesis that artificial selection accompanied by bottlenecks has led to the increase in overall genomic homozygosity. The Bayesian approach and DAPC consistently identified three separate genetic groups. As expected, the farmed populations were clustered together, while the Slovak and Spanish populations formed two separate clusters. Identified traces of genetic admixture in the gene pool of farmed populations reflected a strong contemporary migration rate between them. This study suggests that even if the history of deer farming has been shorter than traditional livestock species, it may leave significant traces in the genome structure.

Restriction site-associated DNA sequencing technologies as an alternative to low-density SNP chips for genomic selection: a simulation study in layer chickens

BACKGROUND

To reduce the cost of genomic selection, a low-density (LD) single nucleotide polymorphism (SNP) chip can be used in combination with imputation for genotyping selection candidates instead of using a high-density (HD) SNP chip. Next-generation sequencing (NGS) techniques have been increasingly used in livestock species but remain expensive for routine use for genomic selection. An alternative and cost-efficient solution is to use restriction site-associated DNA sequencing (RADseq) techniques to sequence only a fraction of the genome using restriction enzymes. From this perspective, use of RADseq techniques followed by an imputation step on HD chip as alternatives to LD chips for genomic selection was studied in a pure layer line.

RESULTS

Genome reduction and sequencing fragments were identified on reference genome using four restriction enzymes (EcoRI, TaqI, AvaII and PstI) and a double-digest RADseq (ddRADseq) method (TaqI-PstI). The SNPs contained in these fragments were detected from the 20X sequence data of the individuals in our population. Imputation accuracy on HD chip with these genotypes was assessed as the mean correlation between true and imputed genotypes. Several production traits were evaluated using single-step GBLUP methodology. The impact of imputation errors on the ranking of the selection candidates was assessed by comparing a genomic evaluation based on ancestry using true HD or imputed HD genotyping. The relative accuracy of genomic estimated breeding values (GEBVs) was investigated by considering the GEBVs estimated on offspring as a reference. With AvaII or PstI and ddRADseq with TaqI and PstI, more than 10 K SNPs were detected in common with the HD SNP chip, resulting in an imputation accuracy greater than 0.97. The impact of imputation errors on genomic evaluation of the breeders was reduced, with a Spearman correlation greater than 0.99. Finally, the relative accuracy of GEBVs was equivalent.

CONCLUSIONS

RADseq approaches can be interesting alternatives to low-density SNP chips for genomic selection. With more than 10 K SNPs in common with the SNPs of the HD SNP chip, good imputation and genomic evaluation results can be obtained. However, with real data, heterogeneity between individuals with missing data must be considered.

Hidden impacts of conservation management on fertility of the critically endangered kākāpō

Background

Animal conservation often requires intensive management actions to improve reproductive output, yet any adverse effects of these may not be immediately apparent, particularly in threatened species with small populations and long lifespans. Hand-rearing is an example of a conservation management strategy which, while boosting populations, can cause long-term demographic and behavioural problems. It is used in the recovery of the critically endangered kākāpō (Strigops habroptilus), a flightless parrot endemic to New Zealand, to improve the slow population growth that is due to infrequent breeding, low fertility and low hatching success.

Methods

We applied Bayesian mixed models to examine whether hand-rearing and other factors were associated with clutch fertility in kākāpō. We used projection predictive variable selection to compare the relative contributions to fertility from the parents' rearing environment, their age and previous copulation experience, the parental kinship, and the number of mates and copulations for each clutch. We also explored how the incidence of repeated copulations and multiple mates varied with kākāpō density.

Results

The rearing status of the clutch father and the number of mates and copulations of the clutch mother were the dominant factors in predicting fertility. Clutches were less likely to be fertile if the father was hand-reared compared to wild-reared, but there was no similar effect for mothers. Clutches produced by females copulating with different males were more likely to be fertile than those from repeated copulations with one male, which in turn had a higher probability of fertility than those from a single copulation. The likelihood of multiple copulations and mates increased with female:male adult sex ratio, perhaps as a result of mate guarding by females. Parental kinship, copulation experience and age all had negligible associations with clutch fertility.

Conclusions

These results provide a rare assessment of factors affecting fertility in a wild threatened bird species, with implications for conservation management. The increased fertility due to multiple mates and copulations, combined with the evidence for mate guarding and previous results of kākāpō sperm morphology, suggests that an evolutionary mechanism exists to optimise fertility through sperm competition in kākāpō. The high frequency of clutches produced from single copulations in the contemporary population may therefore represent an unnatural state, perhaps due to too few females. This suggests that opportunity for sperm competition should be maximised by increasing population densities, optimising sex ratios, and using artificial insemination. The lower fertility of hand-reared males may result from behavioural defects due to lack of exposure to conspecifics at critical development stages, as seen in other taxa. This potential negative impact of hand-rearing must be balanced against the short-term benefits it provides.

Efficiency of genotyping by sequencing in inferring genomic relatedness and molecular insights into fat tail selection in Tunisian sheep

In developing countries, the use of simple and cost-efficient molecular technology is crucial for genetic characterization of local animal resources and better development of conservation strategies. The genotyping by sequencing (GBS) technique, also called restriction enzyme- reduced representational sequencing, is an efficient, cost-effective method for simultaneous discovery and genotyping of many markers. In the present study, we applied a two-enzyme GBS protocol (PstI/MspI) to discover and genotype SNP markers among 197 Tunisian sheep samples. A total of 100 333 bi-allelic SNPs were discovered and genotyped with an SNP call rate of 0.69 and mean sample depth 3.33. The genomic relatedness between 183 samples grouped the samples perfectly to their populations and pointed out a high genetic relatedness of inbred subpopulation reflecting the current adopted reproductive strategies. The genome-wide association study contrasting fat vs. thin-tailed breeds detected 41 significant variants including a peak positioned on OAR20. We identified FOXC1, GMDS, VEGFA, OXCT1, VRTN and BMP2 as the most promising for sheep tail-type trait. The GBS data have been useful to assess the population structure and improve our understanding of the genomic architecture of distinctive characteristics shaped by selection pressure in local sheep breeds. This study successfully investigates a cost-efficient method to discover genotypes, assign populations and understand insights into sheep adaptation to arid area. GBS could be of potential utility in livestock species in developing/emerging countries.

Outlier analyses and genome-wide association study identify glgC and ERD6-like 4 as candidate genes for foliar water-soluble carbohydrate accumulation in Trifolium repens

Increasing water-soluble carbohydrate (WSC) content in white clover is important for improving nutritional quality and reducing environmental impacts from pastoral agriculture. Elucidation of genes responsible for foliar WSC variation would enhance genetic improvement by enabling molecular breeding approaches. The aim of the present study was to identify single nucleotide polymorphisms (SNPs) associated with variation in foliar WSC in white clover. A set of 935 white clover individuals, randomly sampled from five breeding pools selectively bred for divergent (low or high) WSC content, were assessed with 14,743 genotyping-by-sequencing SNPs, using three outlier detection methods: PCAdapt, BayeScan and KGD-F_ST. These analyses identified 33 SNPs as discriminating between high and low WSC populations and putatively under selection. One SNP was located in the intron of ERD6-like 4, a gene coding for a sugar transporter located on the vacuole membrane. A genome-wide association study using a subset of 605 white clover individuals and 5,757 SNPs, identified a further 12 SNPs, one of which was associated with a starch biosynthesis gene, glucose-1-phosphate adenylyltransferase, glgC. Our results provide insight into genomic regions underlying WSC accumulation in white clover, identify candidate genomic regions for further functional validation studies, and reveal valuable information for marker-assisted or genomic selection in white clover.

Genetic diversity and population structure of wild and cultivated Crotalaria species based on genotyping-by-sequencing

Crotalaria is a plant genus that is found all over the world, with over 700 species of herbs and shrubs. The species are potential alternative food and industrial crops due to their adaptability to different environments. Currently, information on the genetic diversity and population structure of these species is scanty. Genotyping-by-sequencing (GBS) is a cost-effective high-throughput technique in diversity evaluation of plant species that have not been fully sequenced. In the current study, de novo GBS was used to characterize 80 Crotalaria accessions from five geographical regions in Kenya. A total of 9820 single nucleotide polymorphism (SNP) markers were obtained after thinning and filtering, which were then used for the analysis of genetic diversity and population structure in Crotalaria. The proportion of SNPs with a minor allele frequency (maf) > = 0.05 was 45.08%, while the Guanine-Cytosine (GC) content was 0.45, from an average sequence depth of 455,909 reads per base. The transition vs transversion ratio was 1.81 and Heterozygosity (He) ranged between 0.01–0.07 in all the sites and 0.04 to 0.52 in the segregating sites. The mean Tajima’s D value for the population was -0.094, suggesting an excess of rare alleles. The fixation index (Fst) between the different populations based on the Wright Fst (1943) ranged from 0.0119 to 0.066 for the Eastern-Western and Nairobi-Western populations. Model based techniques of population structure analysis including structure, k-means and cross-entropy depicted eight clusters in the study accessions. Non-model based techniques especially DAPC depicted poor population stratification. Correspondence Analysis (CA), Principal coordinate analyses (PCoA) and phylogenetic analysis identified a moderate level of population stratification. Results from this study will help conservationists and breeders understand the genetic diversity of Crotalaria. The study also provides valuable information for genetic improvement of domesticated species.

Identification of Genetic Diversity, Pyrrocidine-Producing Strains and Transmission Modes of Endophytic Sarocladium zeae Fungi from Zea Crops

Genotyping by sequencing (GBS) was used to reveal the inherent genetic variation within the haploid fungi Sarocladium zeae isolated from diverse Zea germplasm, including modern Zea mays and its wild progenitors—the teosintes. In accordance with broad host relationship parameters, GBS analysis revealed significant host lineages of S. zeae genetic diversity, indicating that S. zeae genetic variation may associate with different evolutionary histories of host species or varieties. Based on a recently identified PKS-NRPS gene responsible for pyrrocidine biosynthesis in S. zeae fungi, a novel PCR assay was developed to discriminate pyrrocidine-producing S. zeae strains. This molecular method for screening bioactive strains of S. zeae is complementary to other approaches, such as chemical analyses. An eGFP-labelled S. zeae strain was also developed to investigate the endophytic transmission of S. zeae in Z. mays seedlings, which has further improved our understanding of the transmission modes of S. zeae endophytes in maize tissues.

Genetic diversity and background pollen contamination in Norway spruce and Scots pine seed orchard crops

Seed orchards are the key link between tree breeding and production forest for conifer trees. In Sweden, Scots pine and Norway spruce seed orchards currently supply ca. 85% of seedlings used in annual reforestation. The functionality of these seed orchards is thus crucial for supporting long-term production gain and sustainable diversity. We conducted a large-scale genetic investigation of pine and spruce orchards across Sweden using genotyping-by-sequencing. We genotyped 3,300 seedlings/trees from six orchards and 10 natural stands to gain an overview of mating structure and genetic diversity in orchard crops. We found clear differences in observed heterozygosity (H O) and background pollen contamination (BPC) rates between species, with pine orchard crops showing higher H O and BPC than spruce. BPC in pine crops varied from 87% at young orchard age to 12% at mature age, wherease this rate ranged between 27%-4% in spruce crops. Substantial variance in parental contribution was observed in all orchards with 30%-50% parents contibuting to 80% of the progeny. Selfing was low (2%-6%) in all seed crops. Compared to natural stands, orchard crops had slightly lower H O but no strong signal of inbreeding. Our results provide valuable references for orchard management.

Population genomics of New Zealand pouched lamprey (kanakana; piharau; Geotria australis)

Pouched lamprey (Geotria australis) or kanakana/piharau is a culturally and ecologically significant jawless fish that is distributed throughout Aotearoa New Zealand. Despite its importance, much remains unknown about historical relationships and gene flow between populations of this enigmatic species within New Zealand. To help inform management, we assembled a draft G. australis genome and completed the first comprehensive population genomics analysis of pouched lamprey within New Zealand using targeted gene sequencing (Cyt-b and COI) and restriction site-associated DNA sequencing (RADSeq) methods. Employing 16 000 genome-wide single nucleotide polymorphisms (SNPs) derived from RADSeq (n = 186) and sequence data from Cyt-b (766 bp, n = 94) and COI (589 bp, n = 20), we reveal low levels of structure across 10 sampling locations spanning the species range within New Zealand. F-statistics, outlier analyses, and STRUCTURE suggest a single panmictic population, and Mantel and EEMS tests reveal no significant isolation by distance. This implies either ongoing gene flow among populations or recent shared ancestry among New Zealand pouched lamprey. We can now use the information gained from these genetic tools to assist managers with monitoring effective population size, managing potential diseases, and conservation measures such as artificial propagation programs. We further demonstrate the general utility of these genetic tools for acquiring information about elusive species.

Genomic Prediction of Complex Traits in Forage Plants Species: Perennial Grasses Case

The majority of forage grass species are obligate outbreeders. Their breeding classically consists of an initial selection on spaced plants for highly heritable traits such as disease resistances and heading date, followed by familial selection on swards for forage yield and quality traits. The high level of diversity and heterozygosity, and associated decay of linkage disequilibrium (LD) over very short genomic distances, has hampered the implementation of genomic selection (GS) in these species. However, next generation sequencing technologies in combination with the development of genomic resources have recently facilitated implementation of GS in forage grass species such as perennial ryegrass (Lolium perenne L.), switchgrass (Panicum virgatum L.), and timothy (Phleum pratense L.). Experimental work and simulations have shown that GS can increase significantly the genetic gain per unit of time for traits with different levels of heritability. The main reasons are (1) the possibility to select single plants based on their genomic estimated breeding values (GEBV) for traits measured at sward level, (2) a reduction in the duration of selection cycles, and less importantly (3) an increase in the selection intensity associated with an increase in the genetic variance used for selection. Nevertheless, several factors should be taken into account for the successful implementation of GS in forage grasses. For example, it has been shown that the level of relatedness between the training and the selection population is particularly critical when working with highly structured meta-populations consisting of several genetic groups. A sufficient number of markers should be used to estimate properly the kinship between individuals and to reflect the variability of major QTLs. It is also important that the prediction models are trained for relevant environments when dealing with traits with high genotype × environment interaction (G × E). Finally, in these outbreeding species, measures to reduce inbreeding should be used to counterbalance the high selection intensity that can be achieved in GS.

Development of Epigenetic Clocks for Key Ruminant Species

Robust biomarkers of chronological age have been developed in humans and model mammalian species such as rats and mice using DNA methylation data. The concept of these so-called "epigenetic clocks" has emerged from a large body of literature describing the relationship between genome-wide methylation levels and age. Epigenetic clocks exploit this phenomenon and use small panels of differentially methylated cytosine (CpG) sites to make robust predictions of chronological age, independent of tissue type. Here, we present highly accurate livestock epigenetic clocks for which we have used the custom mammalian methylation array "HorvathMammalMethyl40" to construct the first epigenetic clock for domesticated goat (Capra hircus), cattle (Bos taurus), Red (Cervus elaphus) and Wapiti deer (Cervus canadensis) and composite-breed sheep (Ovis aries). Additionally, we have constructed a 'farm animal clock' for all species included in the study, which will allow for robust predictions to be extended to various breeds/strains. The farm animal clock shows similarly high accuracy to the individual species' clocks (r > 0.97), utilizing only 217 CpG sites to estimate age (relative to the maximum lifespan of the species) with a single mathematical model. We hypothesise that the applications of this livestock clock could extend well beyond the scope of chronological age estimates. Many independent studies have demonstrated that a deviation between true age and clock derived molecular age is indicative of past and/or present health (including stress) status. There is, therefore, untapped potential to utilize livestock clocks in breeding programs as a predictor for age-related production traits.

Genetic diversity and selective sweeps in historical and modern Canadian spring wheat cultivars using the 90K SNP array

Previous molecular characterization studies conducted in Canadian wheat cultivars shed some light on the impact of plant breeding on genetic diversity, but the number of varieties and markers used was small. Here, we used 28,798 markers of the wheat 90K single nucleotide polymorphisms to (a) assess the extent of genetic diversity, relationship, population structure, and divergence among 174 historical and modern Canadian spring wheat varieties registered from 1905 to 2018 and 22 unregistered lines (hereinafter referred to as cultivars), and (b) identify genomic regions that had undergone selection. About 91% of the pairs of cultivars differed by 20-40% of the scored alleles, but only 7% of the pairs had kinship coefficients of < 0.250, suggesting the presence of a high proportion of redundancy in allelic composition. Although the 196 cultivars represented eight wheat classes, our results from phylogenetic, principal component, and the model-based population structure analyses revealed three groups, with no clear structure among most wheat classes, breeding programs, and breeding periods. F_ST statistics computed among different categorical variables showed little genetic differentiation (< 0.05) among breeding periods and breeding programs, but a diverse level of genetic differentiation among wheat classes and predicted groups. Diversity indices were the highest and lowest among cultivars registered from 1970 to 1980 and from 2011 to 2018, respectively. Using two outlier detection methods, we identified from 524 to 2314 SNPs and 41 selective sweeps of which some are close to genes with known phenotype, including plant height, photoperiodism, vernalization, gluten strength, and disease resistance.

Construction of a high-density genetic map and detection of a major QTL of resistance to powdery mildew (Erysiphe necator Sch.) in Caucasian grapes (Vitis vinifera L.)

BACKGROUND

Vitis vinifera L. is the most cultivated grapevine species worldwide. Erysiphe necator Sch., the causal agent of grape powdery mildew, is one of the main pathogens affecting viticulture. V. vinifera has little or no genetic resistances against E. necator and the grape industry is highly dependent on agrochemicals. Some Caucasian V. vinifera accessions have been reported to be resistant to E. necator and to have no genetic relationships to known sources of resistance to powdery mildew. The main purpose of this work was the study and mapping of the resistance to E. necator in the Caucasian grapes 'Shavtsitska' and 'Tskhvedianis tetra'.

RESULTS

The Caucasian varieties 'Shavtsitska' and 'Tskhvedianis tetra' showed a strong partial resistance to E. necator which segregated in two cross populations: the resistant genotypes delayed and limited the pathogen mycelium growth, sporulation intensity and number of conidia generated. A total of 184 seedlings of 'Shavtsitska' x 'Glera' population were genotyped through the Genotyping by Sequencing (GBS) technology and two high-density linkage maps were developed for the cross parents. The QTL analysis revealed a major resistance locus, explaining up to 80.15% of the phenotypic variance, on 'Shavtsitska' linkage group 13, which was associated with a reduced pathogen infection as well as an enhanced plant necrotic response. The genotyping of 105 Caucasian accessions with SSR markers flanking the QTL revealed that the resistant haplotype of 'Shavtsitska' was shared by 'Tskhvedianis tetra' and a total of 25 Caucasian grape varieties, suggesting a widespread presence of this resistance in the surveyed germplasm. The uncovered QTL was mapped in the region where the Ren1 locus of resistance to E. necator, identified in the V. vinifera 'Kishmish vatkana' and related grapes of Central Asia, is located. The genetic analysis conducted revealed that the Caucasian grapes in this study exhibit a resistant haplotype different from that of Central Asian grape accessions.

CONCLUSIONS

The QTL isolated in 'Shavtsitska' and present in the Caucasian V. vinifera varieties could be a new candidate gene of resistance to E. necator to use in breeding programmes. It co-localizes with the Ren1 locus but shows a different haplotype from that of grapevines of Central Asia. We therefore consider that the Caucasian resistance locus, named Ren1.2, contains a member of a cluster of R-genes, of which the region is rich, and to be linked with, or possibly allelic, to Ren1.

Genomic signatures of inbreeding in a critically endangered parrot, the kākāpō

Events of inbreeding are inevitable in critically endangered species. Reduced population sizes and unique life-history traits can increase the severity of inbreeding, leading to declines in fitness and increased risk of extinction. Here, we investigate levels of inbreeding in a critically endangered flightless parrot, the kākāpō (Strigops habroptilus), wherein a highly inbred island population and one individual from the mainland of New Zealand founded the entire extant population. Genotyping-by-sequencing (GBS), and a genotype calling approach using a chromosome-level genome assembly, identified a filtered set of 12,241 single-nucleotide polymorphisms (SNPs) among 161 kākāpō, which together encompass the total genetic potential of the extant population. Multiple molecular-based estimates of inbreeding were compared, including genome-wide estimates of heterozygosity (FH), the diagonal elements of a genomic-relatedness matrix (FGRM), and runs of homozygosity (RoH, FRoH). In addition, we compared levels of inbreeding in chicks from a recent breeding season to examine if inbreeding is associated with offspring survival. The density of SNPs generated with GBS was sufficient to identify chromosomes that were largely homozygous with RoH distributed in similar patterns to other inbred species. Measures of inbreeding were largely correlated and differed significantly between descendants of the two founding populations. However, neither inbreeding nor ancestry was found to be associated with reduced survivorship in chicks, owing to unexpected mortality in chicks exhibiting low levels of inbreeding. Our study highlights important considerations for estimating inbreeding in critically endangered species, such as the impacts of small population sizes and admixture between diverse lineages.

Mixed Mating in a Multi-Origin Population Suggests High Potential for Genetic Rescue in North Island Brown Kiwi, Apteryx mantelli

Reinforcement translocations are increasingly utilised in conservation with the goal of achieving genetic rescue. However, concerns regarding undesirable results, such as genetic homogenisation or replacement, are widespread. One factor influencing translocation outcomes is the rate at which the resident and the introduced individuals interbreed. Consequently, post-release mate choice is a key behaviour to consider in conservation planning. Here we studied mating, and its consequences for genomic admixture, in the North Island brown kiwi Apteryx mantelli population on Ponui Island which was founded by two translocation events over 50 years ago. The two source populations used are now recognised as belonging to two separate management units between which birds differ in size and are genetically differentiated. We examined the correlation between male and female morphometrics for 17 known pairs and quantified the relatedness of 20 pairs from this admixed population. In addition, we compared the genetic similarity and makeup of 106 Ponui Island birds, including 23 known pairs, to birds representing the source populations for the original translocations. We found no evidence for size-assortative mating. On the contrary, genomic SNP data suggested that kiwi of one feather did not flock together, meaning that mate choice resulted in pairing between individuals that were less related than expected by random chance. Furthermore, the birds in the current Ponui Island population were found to fall along a gradient of genomic composition consistent with non-clustered representation of the two parental genomes. These findings indicate potential for successful genetic rescue in future Apteryx reinforcement translocations, a potential that is currently under utilised due to restrictive translocation policies. In light of our findings, we suggest that reconsideration of these policies could render great benefits for the future diversity of this iconic genus in New Zealand.

Peridermal fruit skin formation in Actinidia sp. (kiwifruit) is associated with genetic loci controlling russeting and cuticle formation

BACKGROUND

The skin (exocarp) of fleshy fruit is hugely diverse across species. Most fruit types have a live epidermal skin covered by a layer of cuticle made up of cutin while a few create an outermost layer of dead cells (peridermal layer).

RESULTS

In this study we undertook crosses between epidermal and peridermal skinned kiwifruit, and showed that epidermal skin is a semi-dominant trait. Furthermore, backcrossing these epidermal skinned hybrids to a peridermal skinned fruit created a diverse range of phenotypes ranging from epidermal skinned fruit, through fruit with varying degrees of patches of periderm (russeting), to fruit with a complete periderm. Quantitative trait locus (QTL) analysis of this population suggested that periderm formation was associated with four loci. These QTLs were aligned either to ones associated with russet formation on chromosome 19 and 24, or cuticle integrity and coverage located on chromosomes 3, 11 and 24.

CONCLUSION

From the segregation of skin type and QTL analysis, it appears that skin development in kiwifruit is controlled by two competing factors, cuticle strength and propensity to russet. A strong cuticle will inhibit russeting while a strong propensity to russet can create a continuous dead skinned periderm.

Genetic basis and identification of candidate genes for wooden breast and white striping in commercial broiler chickens

Wooden breast (WB) and white striping (WS) are highly prevalent and economically damaging muscle disorders of modern commercial broiler chickens characterized respectively by palpable firmness and fatty white striations running parallel to the muscle fiber. High feed efficiency and rapid growth, especially of the breast muscle, are believed to contribute to development of such muscle defects; however, their etiology remains poorly understood. To gain insight into the genetic basis of these myopathies, a genome-wide association study was conducted using a commercial crossbred broiler population (n = 1193). Heritability was estimated at 0.5 for WB and WS with high genetic correlation between them (0.88). GWAS revealed 28 quantitative trait loci (QTL) on five chromosomes for WB and 6 QTL on one chromosome for WS, with the majority of QTL for both myopathies located in a ~ 8 Mb region of chromosome 5. This region has highly conserved synteny with a portion of human chromosome 11 containing a cluster of imprinted genes associated with growth and metabolic disorders such as type 2 diabetes and Beckwith-Wiedemann syndrome. Candidate genes include potassium voltage-gated channel subfamily Q member 1 (KCNQ1), involved in insulin secretion and cardiac electrical activity, lymphocyte-specific protein 1 (LSP1), involved in inflammation and immune response.

Genotyping Strategies Using ddRAD Sequencing in Farmed Arctic Charr (Salvelinus alpinus)

Incorporation of genomic technologies into fish breeding programs is a modern reality, promising substantial advances regarding the accuracy of selection, monitoring the genetic diversity and pedigree record verification. Single nucleotide polymorphism (SNP) arrays are the most commonly used genomic tool, but the investments required make them unsustainable for emerging species, such as Arctic charr (Salvelinus alpinus), where production volume is low. The requirement to genotype a large number of animals for breeding practices necessitates cost effective genotyping approaches. In the current study, we used double digest restriction site-associated DNA (ddRAD) sequencing of either high or low coverage to genotype Arctic charr from the Swedish national breeding program and performed analytical procedures to assess their utility in a range of tasks. SNPs were identified and used for deciphering the genetic structure of the studied population, estimating genomic relationships and implementing an association study for growth-related traits. Missing information and underestimation of heterozygosity in the low coverage set were limiting factors in genetic diversity and genomic relationship analyses, where high coverage performed notably better. On the other hand, the high coverage dataset proved to be valuable when it comes to identifying loci that are associated with phenotypic traits of interest. In general, both genotyping strategies offer sustainable alternatives to hybridization-based genotyping platforms and show potential for applications in aquaculture selective breeding.

Whole-organism eQTL mapping at cellular resolution with single-cell sequencing

Genetic regulation of gene expression underlies variation in disease risk and other complex traits. The effect of expression quantitative trait loci (eQTLs) varies across cell types; however, the complexity of mammalian tissues makes studying cell-type eQTLs highly challenging. We developed a novel approach in the model nematode Caenorhabditis elegans that uses single-cell RNA sequencing to map eQTLs at cellular resolution in a single one-pot experiment. We mapped eQTLs across cell types in an extremely large population of genetically distinct C. elegans individuals. We found cell-type-specific trans eQTL hotspots that affect the expression of core pathways in the relevant cell types. Finally, we found single-cell-specific eQTL effects in the nervous system, including an eQTL with opposite effects in two individual neurons. Our results show that eQTL effects can be specific down to the level of single cells.

Genomic prediction and genomic heritability of grain yield and its related traits in a safflower genebank collection

Safflower, a minor oilseed crop, is gaining increased attention for food and industrial uses. Safflower genebank collections are an important genetic resource for crop enhancement and future breeding programs. In this study, we investigated the population structure of a safflower collection sourced from the Australian Grain Genebank and assessed the potential of genomic prediction (GP) to evaluate grain yield and related traits using single and multi-site models. Prediction accuracies (PA) of genomic best linear unbiased prediction (GBLUP) from single site models ranged from 0.21 to 0.86 for all traits examined and were consistent with estimated genomic heritability (h² ), which varied from low to moderate across traits. We generally observed a low level of genome × environment interactions (g × E). Multi-site g × E GBLUP models only improved PA for accessions with at least some phenotypes in the training set. We observed that relaxing quality filtering parameters for genotype-by-sequencing (GBS), such as missing genotype call rate, did not affect PA but upwardly biased h² estimation. Our results indicate that GP is feasible in safflower evaluation and is potentially a cost-effective tool to facilitate fast introgression of desired safflower trait variation from genebank germplasm into breeding lines.

Management tools for genetic diversity in an isolated population of the honeybee (Apis mellifera) in New Zealand

Context Beekeepers and honeybee queen breeders alike currently have few tools at their disposal for the management of genetic diversity inside their populations. Pedigree information is often absent, beekeepers cannot afford to genotype selection candidates due to costs, and acquisition of material for genotyping without risk to individual queen bees is difficult. However, in New World honeybee populations where import of additional genetic material to refresh the population is restricted (e.g. Australia) or impossible (e.g. New Zealand), management of genetic diversity is important for population sustainability. While the role of individual beekeepers in maintaining genetic resources becomes crucial under these circumstances, a more holistic approach to the management of genetic diversity is needed to allow for maximum impact of their contribution. Aims The establishment of affordable genotyping methodologies for successful strategies in managing honeybee genetic diversity, as well as the necessary delivery systems for the results to support the beekeeping community by providing interpretation in the context of the wider population. Methods Genotyping-by-sequencing of honeybee samples collected as part of a national survey were used as the basis for assessment of New Zealand’s honeybee genetic diversity and development of a tool with largely self-explanatory outputs that can be used directly by beekeepers. Key results It would appear that New Zealand’s honeybee population is sufficiently diverse to maintain population viability. However, both within regions and within companies, genetic diversity is significantly reduced, especially in the case of specialised queen breeders, indicating that active management will be necessary to achieve long-term sustainability. Conclusions Interactive tools are needed to help beekeepers understand their role in maintaining overall genetic diversity in the honeybee population as well as the potential impact of planned bee movement and queen acquisitions. Regular rounds of voluntary bee sampling can be used as the basis for management decisions without concentrating genotyping load on specific operators and restricting diversity assessments to subpopulations. Implications The described strategy is expected to both improve the outlook of New Zealand’s honeybee population as a whole and facilitate stringent genetic improvement programs by enabling queen breeders to make informed selection decisions and giving beekeepers confidence in the viability of their population.

The Future of Livestock Management: A Review of Real-Time Portable Sequencing Applied to Livestock

Oxford Nanopore Technologies' MinION has proven to be a valuable tool within human and microbial genetics. Its capacity to produce long reads in real time has opened up unique applications for portable sequencing. Examples include tracking the recent African swine fever outbreak in China and providing a diagnostic tool for disease in the cassava plant in Eastern Africa. Here we review the current applications of Oxford Nanopore sequencing in livestock, then focus on proposed applications in livestock agriculture for rapid diagnostics, base modification detection, reference genome assembly and genomic prediction. In particular, we propose a future application: 'crush-side genotyping' for real-time on-farm genotyping for extensive industries such as northern Australian beef production. An initial in silico experiment to assess the feasibility of crush-side genotyping demonstrated promising results. SNPs were called from simulated Nanopore data, that included the relatively high base call error rate that is characteristic of the data, and calling parameters were varied to understand the feasibility of SNP calling at low coverages in a heterozygous population. With optimised genotype calling parameters, over 85% of the 10,000 simulated SNPs were able to be correctly called with coverages as low as 6×. These results provide preliminary evidence that Oxford Nanopore sequencing has potential to be used for real-time SNP genotyping in extensive livestock operations.

Optimizing whole-genomic prediction for autotetraploid blueberry breeding

Blueberry (Vaccinium spp.) is an important autopolyploid crop with significant benefits for human health. Apart from its genetic complexity, the feasibility of genomic prediction has been proven for blueberry, enabling a reduction in the breeding cycle time and increasing genetic gain. However, as for other polyploid crops, sequencing costs still hinder the implementation of genome-based breeding methods for blueberry. This motivated us to evaluate the effect of training population sizes and composition, as well as the impact of marker density and sequencing depth on phenotype prediction for the species. For this, data from a large real breeding population of 1804 individuals were used. Genotypic data from 86,930 markers and three traits with different genetic architecture (fruit firmness, fruit weight, and total yield) were evaluated. Herein, we suggested that marker density, sequencing depth, and training population size can be substantially reduced with no significant impact on model accuracy. Our results can help guide decisions toward resource allocation (e.g., genotyping and phenotyping) in order to maximize prediction accuracy. These findings have the potential to allow for a faster and more accurate release of varieties with a substantial reduction of resources for the application of genomic prediction in blueberry. We anticipate that the benefits and pipeline described in our study can be applied to optimize genomic prediction for other diploid and polyploid species.

QTL Mapping for Resistance to Cankers Induced by Pseudomonas syringae pv. actinidiae (Psa) in a Tetraploid Actinidia chinensis Kiwifruit Population

Polyploidy is a key driver of significant evolutionary changes in plant species. The genus Actinidia (kiwifruit) exhibits multiple ploidy levels, which contribute to novel fruit traits, high yields and resistance to the canker-causing dieback disease incited by Pseudomonas syringae pv. actinidiae (Psa) biovar 3. However, the genetic mechanism for resistance to Psa observed in polyploid kiwifruit is not yet known. In this study we performed detailed genetic analysis of a tetraploid Actinidia chinensis var. chinensis population derived from a cross between a female parent that exhibits weak tolerance to Psa and a highly Psa-resistant male parent. We used the capture-sequencing approach across the whole kiwifruit genome and generated the first ultra-dense maps in a tetraploid kiwifruit population. We located quantitative trait loci (QTLs) for Psa resistance on these maps. Our approach to QTL mapping is based on the use of identity-by-descent trait mapping, which allowed us to relate the contribution of specific alleles from their respective homologues in the male and female parent, to the control of Psa resistance in the progeny. We identified genes in the diploid reference genome whose function is suggested to be involved in plant defense, which underly the QTLs, including receptor-like kinases. Our study is the first to cast light on the genetics of a polyploid kiwifruit and suggest a plausible mechanism for Psa resistance in this species.

Evaluation of the efficiency of genomic versus pedigree predictions for growth and wood quality traits in Scots pine

BACKGROUND

Genomic selection (GS) or genomic prediction is a promising approach for tree breeding to obtain higher genetic gains by shortening time of progeny testing in breeding programs. As proof-of-concept for Scots pine (Pinus sylvestris L.), a genomic prediction study was conducted with 694 individuals representing 183 full-sib families that were genotyped with genotyping-by-sequencing (GBS) and phenotyped for growth and wood quality traits. 8719 SNPs were used to compare different genomic with pedigree prediction models. Additionally, four prediction efficiency methods were used to evaluate the impact of genomic breeding value estimations by assigning diverse ratios of training and validation sets, as well as several subsets of SNP markers.

RESULTS

Genomic Best Linear Unbiased Prediction (GBLUP) and Bayesian Ridge Regression (BRR) combined with expectation maximization (EM) imputation algorithm showed slightly higher prediction efficiencies than Pedigree Best Linear Unbiased Prediction (PBLUP) and Bayesian LASSO, with some exceptions. A subset of approximately 6000 SNP markers, was enough to provide similar prediction efficiencies as the full set of 8719 markers. Additionally, prediction efficiencies of genomic models were enough to achieve a higher selection response, that varied between 50-143% higher than the traditional pedigree-based selection.

CONCLUSIONS

Although prediction efficiencies were similar for genomic and pedigree models, the relative selection response was doubled for genomic models by assuming that earlier selections can be done at the seedling stage, reducing the progeny testing time, thus shortening the breeding cycle length roughly by 50%.

Development and validation of a RAD-Seq target-capture based genotyping assay for routine application in advanced black tiger shrimp (Penaeus monodon) breeding programs

Background

The development of genome-wide genotyping resources has provided terrestrial livestock and crop industries with the unique ability to accurately assess genomic relationships between individuals, uncover the genetic architecture of commercial traits, as well as identify superior individuals for selection based on their specific genetic profile. Utilising recent advancements in de-novo genome-wide genotyping technologies, it is now possible to provide aquaculture industries with these same important genotyping resources, even in the absence of existing genome assemblies. Here, we present the development of a genome-wide SNP assay for the Black Tiger shrimp (Penaeus monodon) through utilisation of a reduced-representation whole-genome genotyping approach (DArTseq).

Results

Based on a single reduced-representation library, 31,262 polymorphic SNPs were identified across 650 individuals obtained from Australian wild stocks and commercial aquaculture populations. After filtering to remove SNPs with low read depth, low MAF, low call rate, deviation from HWE, and non-Mendelian inheritance, 7542 high-quality SNPs were retained. From these, 4236 high-quality genome-wide loci were selected for baits-probe development and 4194 SNPs were included within a finalized target-capture genotype-by-sequence assay (DArTcap). This assay was designed for routine and cost effective commercial application in large scale breeding programs, and demonstrates higher confidence in genotype calls through increased call rate (from 80.2 ± 14.7 to 93.0% ± 3.5%), increased read depth (from 20.4 ± 15.6 to 80.0 ± 88.7), as well as a 3-fold reduction in cost over traditional genotype-by-sequencing approaches.

Conclusion

Importantly, this assay equips the P. monodon industry with the ability to simultaneously assign parentage of communally reared animals, undertake genomic relationship analysis, manage mate pairings between cryptic family lines, as well as undertake advance studies of genome and trait architecture. Critically this assay can be cost effectively applied as P. monodon breeding programs transition to undertaking genomic selection.

Robust genome-wide ancestry inference for heterogeneous datasets: illustrated using the 1,000 genome project with 3D facial images

Estimates of individual-level genomic ancestry are routinely used in human genetics, and related fields. The analysis of population structure and genomic ancestry can yield insights in terms of modern and ancient populations, allowing us to address questions regarding admixture, and the numbers and identities of the parental source populations. Unrecognized population structure is also an important confounder to correct for in genome-wide association studies. However, it remains challenging to work with heterogeneous datasets from multiple studies collected by different laboratories with diverse genotyping and imputation protocols. This work presents a new approach and an accompanying open-source toolbox that facilitates a robust integrative analysis for population structure and genomic ancestry estimates for heterogeneous datasets. We show robustness against individual outliers and different protocols for the projection of new samples into a reference ancestry space, and the ability to reveal and adjust for population structure in a simulated case–control admixed population. Given that visually evident and easily recognizable patterns of human facial characteristics co-vary with genomic ancestry, and based on the integration of three different sources of genome data, we generate average 3D faces to illustrate genomic ancestry variations within the 1,000 Genome project and for eight ancient-DNA profiles, respectively.

Genetic Diversity in Invasive Populations of Argentine Stem Weevil Associated with Adaptation to Biocontrol

Modified, agricultural landscapes are susceptible to damage by insect pests. Biological control of pests is typically successful once a control agent has established, but this depends on the agent’s capacity to co-evolve with the host. Theoretical studies have shown that different levels of genetic variation between the host and the control agent will lead to rapid evolution of resistance in the host. Although this has been reported in one instance, the underlying genetics have not been studied. To address this, we measured the genetic variation in New Zealand populations of the pasture pest, Argentine stem weevil (Listronotus bonariensis), which is controlled with declining effectiveness by a parasitoid wasp, Microctonus hyperodae. We constructed a draft reference genome of the weevil, collected samples from a geographical survey of 10 sites around New Zealand, and genotyped them using a modified genotyping-by-sequencing approach. New Zealand populations of Argentine stem weevil have high levels of heterozygosity and low population structure, consistent with a large effective population size and frequent gene flow. This implies that Argentine stem weevils were able to evolve more rapidly than their biocontrol agent, which reproduces asexually. These findings show that monitoring genetic diversity in biocontrol agents and their targets is critical for long-term success of biological control.

Application of Low Coverage Genotyping by Sequencing in Selectively Bred Arctic Charr (Salvelinus alpinus)

Arctic charr (Salvelinus alpinus) is a species of high economic value for the aquaculture industry, and of high ecological value due to its Holarctic distribution in both marine and freshwater environments. Novel genome sequencing approaches enable the study of population and quantitative genetic parameters even on species with limited or no prior genomic resources. Low coverage genotyping by sequencing (GBS) was applied in a selected strain of Arctic charr in Sweden originating from a landlocked freshwater population. For the needs of the current study, animals from year classes 2013 (171 animals, parental population) and 2017 (759 animals; 13 full sib families) were used as a template for identifying genome wide single nucleotide polymorphisms (SNPs). GBS libraries were constructed using the PstI and MspI restriction enzymes. Approximately 14.5K SNPs passed quality control and were used for estimating a genomic relationship matrix. Thereafter a wide range of analyses were conducted in order to gain insights regarding genetic diversity and investigate the efficiency of the genomic information for parentage assignment and breeding value estimation. Heterozygosity estimates for both year classes suggested a slight excess of heterozygotes. Furthermore, FST estimates among the families of year class 2017 ranged between 0.009 - 0.066. Principal components analysis (PCA) and discriminant analysis of principal components (DAPC) were applied aiming to identify the existence of genetic clusters among the studied population. Results obtained were in accordance with pedigree records allowing the identification of individual families. Additionally, DNA parentage verification was performed, with results in accordance with the pedigree records with the exception of a putative dam where full sib genotypes suggested a potential recording error. Breeding value estimation for juvenile growth through the usage of the estimated genomic relationship matrix clearly outperformed the pedigree equivalent in terms of prediction accuracy (0.51 opposed to 0.31). Overall, low coverage GBS has proven to be a cost-effective genotyping platform that is expected to boost the selection efficiency of the Arctic charr breeding program.

SNP discovery and population structure analysis in Lassi and Marecha camel breeds using a genotyping by sequencing method

Pakistani camels have been classified socio-geographically into 20 breeds, but they have not yet been subjected to substantial selective pressures and the genetic basis for these breeds is not understood. However, it should be possible to distinguish them by use of molecular data. This study investigated the genetic diversity and population structure within and between two major Pakistani camel breeds, Marecha and Lassi. As no SNP array is currently available, we first identified 63 619 SNPs using a genotyping by sequencing approach. After quality control, a panel of 36 926 SNPs was used in the analysis. Population structure was investigated with a principal coordinate analysis as well as a cluster analysis using NetView, and multilocus heterozygosity analysis to explore between- and within-breed genetic variation. In addition, between-breed variation was explored using the fixation index, F_ST . We also compared relationship matrices computed using the VanRaden SNP-based method and a method developed specifically for genotyping by sequencing data. Among the two camel breeds, Lassi showed a lower level of genetic diversity whereas Marecha showed a higher level. As a genotyping platform has not yet been developed for the camel, the SNPs discovered in this study will be useful in future genetic studies in camels.

A comparison of pedigree, genetic and genomic estimates of relatedness for informing pairing decisions in two critically endangered birds: Implications for conservation breeding programmes worldwide

Conservation management strategies for many highly threatened species include conservation breeding to prevent extinction and enhance recovery. Pairing decisions for these conservation breeding programmes can be informed by pedigree data to minimize relatedness between individuals in an effort to avoid inbreeding, maximize diversity and maintain evolutionary potential. However, conservation breeding programmes struggle to use this approach when pedigrees are shallow or incomplete. While genetic data (i.e., microsatellites) can be used to estimate relatedness to inform pairing decisions, emerging evidence indicates this approach may lack precision in genetically depauperate species, and more effective estimates will likely be obtained from genomic data (i.e., thousands of genome-wide single nucleotide polymorphisms, or SNPs). Here, we compare relatedness estimates and subsequent pairing decisions using pedigrees, microsatellites and SNPs from whole-genome resequencing approaches in two critically endangered birds endemic to New Zealand: kakī/black stilt (Himantopus novaezelandiae) and kākāriki karaka/orange-fronted parakeet (Cyanoramphus malherbi). Our findings indicate that SNPs provide more precise estimates of relatedness than microsatellites when assessing empirical parent-offspring and full sibling relationships. Further, our results show that relatedness estimates and subsequent pairing recommendations using PMx are most similar between pedigree- and SNP-based approaches. These combined results indicate that in lieu of robust pedigrees, SNPs are an effective tool for informing pairing decisions, which has important implications for many poorly pedigreed conservation breeding programmes worldwide.

Genetic diversity, demographic history and neo-sex chromosomes in the Critically Endangered Raso lark

Small effective population sizes could expose island species to inbreeding and loss of genetic variation. Here, we investigate factors shaping genetic diversity in the Raso lark, which has been restricted to a single islet for approximately 500 years, with a population size of a few hundred. We assembled a reference genome for the related Eurasian skylark and then assessed diversity and demographic history using RAD-seq data (75 samples from Raso larks and two related mainland species). We first identify broad tracts of suppressed recombination in females, indicating enlarged neo-sex chromosomes. We then show that genetic diversity across autosomes in the Raso lark is lower than in its mainland relatives, but inconsistent with long-term persistence at its current population size. Finally, we find that genetic signatures of the recent population contraction are overshadowed by an ancient expansion and persistence of a very large population until the human settlement of Cape Verde. Our findings show how genome-wide approaches to study endangered species can help avoid confounding effects of genome architecture on diversity estimates, and how present-day diversity can be shaped by ancient demographic events.