Major Impacts of Widespread Structural Variation on Gene Expression and Crop Improvement in Tomato.

Structural variants (SVs) underlie important crop improvement and domestication traits. However, resolving the extent, diversity, and quantitative impact of SVs has been challenging. We used long-read nanopore sequencing to capture 238,490 SVs in 100 diverse tomato lines. This panSV genome, along with 14 new reference assemblies, revealed large-scale intermixing of diverse genotypes, as well as thousands of SVs intersecting genes and cis-regulatory regions. Hundreds of SV-gene pairs exhibit subtle and significant expression changes, which could broadly influence quantitative trait variation. By combining quantitative genetics with genome editing, we show how multiple SVs that changed gene dosage and expression levels modified fruit flavor, size, and production. In the last example, higher order epistasis among four SVs affecting three related transcription factors allowed introduction of an important harvesting trait in modern tomato. Our findings highlight the underexplored role of SVs in genotype-to-phenotype relationships and their widespread importance and utility in crop improvement.

Defining Genetic Variation in Widely Used Congenic and Backcrossed Mouse Models Reveals Varied Regulation of Genes Important for Immune Responses.

Genetic variation influences how the genome is interpreted in individuals and in mouse strains used to model immune responses. We developed approaches to utilize next-generation sequencing datasets to identify sequence variation in genes and enhancer elements in congenic and backcross mouse models. We defined genetic variation in the widely used B6-CD45.2 and B6.SJL-CD45.1 congenic model, identifying substantial differences in SJL genetic content retained in B6.SJL-CD45.1 strains on the basis of the vendor source of the mice. Genes encoding PD-1, CD62L, Bcl-2, cathepsin E, and Cxcr4 were within SJL genetic content in at least one vendor source of B6.SJL-CD45.1 mice. SJL genetic content affected enhancer elements, gene regulation, protein expression, and amino acid content in CD4+ T helper 1 cells, and mice infected with influenza showed reduced expression of Cxcr4 on B6.SJL-CD45.1 T follicular helper cells. These findings provide information on experimental variables and aid in creating approaches that account for genetic variables.

Genetic load: genomic estimates and applications in non-model animals.

Genetic variation, which is generated by mutation, recombination and gene flow, can reduce the mean fitness of a population, both now and in the future. This 'genetic load' has been estimated in a wide range of animal taxa using various approaches. Advances in genome sequencing and computational techniques now enable us to estimate the genetic load in populations and individuals without direct fitness estimates. Here, we review the classic and contemporary literature of genetic load. We describe approaches to quantify the genetic load in whole-genome sequence data based on evolutionary conservation and annotations. We show that splitting the load into its two components - the realized load (or expressed load) and the masked load (or inbreeding load) - can improve our understanding of the population genetics of deleterious mutations.

Island-specific evolution of a sex-primed autosome in a sexual planarian.

The sexual strain of the planarian Schmidtea mediterranea, indigenous to Tunisia and several Mediterranean islands, is a hermaphrodite1,2. Here we isolate individual chromosomes and use sequencing, Hi-C3,4 and linkage mapping to assemble a chromosome-scale genome reference. The linkage map reveals an extremely low rate of recombination on chromosome 1. We confirm suppression of recombination on chromosome 1 by genotyping individual sperm cells and oocytes. We show that previously identified genomic regions that maintain heterozygosity even after prolonged inbreeding make up essentially all of chromosome 1. Genome sequencing of individuals isolated in the wild indicates that this phenomenon has evolved specifically in populations from Sardinia and Corsica. We find that most known master regulators5-13 of the reproductive system are located on chromosome 1. We used RNA interference14,15 to knock down a gene with haplotype-biased expression, which led to the formation of a more pronounced female mating organ. On the basis of these observations, we propose that chromosome 1 is a sex-primed autosome primed for evolution into a sex chromosome.

An Adversarial DNA N6-Methyladenine-Sensor Network Preserves Polycomb Silencing.

Adenine N6 methylation in DNA (6mA) is widespread among bacteria and phage and is detected in mammalian genomes, where its function is largely unexplored. Here we show that 6mA deposition and removal are catalyzed by the Mettl4 methyltransferase and Alkbh4 dioxygenase, respectively, and that 6mA accumulation in genic elements corresponds with transcriptional silencing. Inactivation of murine Mettl4 depletes 6mA and causes sublethality and craniofacial dysmorphism in incross progeny. We identify distinct 6mA sensor domains of prokaryotic origin within the MPND deubiquitinase and ASXL1, a component of the Polycomb repressive deubiquitinase (PR-DUB) complex, both of which act to remove monoubiquitin from histone H2A (H2A-K119Ub), a repressive mark. Deposition of 6mA by Mettl4 triggers the proteolytic destruction of both sensor proteins, preserving genome-wide H2A-K119Ub levels. Expression of the bacterial 6mA methyltransferase Dam, in contrast, fails to destroy either sensor. These findings uncover a native, adversarial 6mA network architecture that preserves Polycomb silencing.

Detecting inbreeding depression in structured populations.

Measuring inbreeding and its consequences on fitness is central for many areas in biology including human genetics and the conservation of endangered species. However, there is no consensus on the best method, neither for quantification of inbreeding itself nor for the model to estimate its effect on specific traits. We simulated traits based on simulated genomes from a large pedigree and empirical whole-genome sequences of human data from populations with various sizes and structures (from the 1,000 Genomes project). We compare the ability of various inbreeding coefficients ([Formula: see text]) to quantify the strength of inbreeding depression: allele-sharing, two versions of the correlation of uniting gametes which differ in the weight they attribute to each locus and two identical-by-descent segments-based estimators. We also compare two models: the standard linear model and a linear mixed model (LMM) including a genetic relatedness matrix (GRM) as random effect to account for the nonindependence of observations. We find LMMs give better results in scenarios with population or family structure. Within the LMM, we compare three different GRMs and show that in homogeneous populations, there is little difference among the different [Formula: see text] and GRM for inbreeding depression quantification. However, as soon as a strong population or family structure is present, the strength of inbreeding depression can be most efficiently estimated only if i) the phenotypes are regressed on [Formula: see text] based on a weighted version of the correlation of uniting gametes, giving more weight to common alleles and ii) with the GRM obtained from an allele-sharing relatedness estimator.

Purging due to self-fertilization does not prevent accumulation of expansion load.

As species expand their geographic ranges, colonizing populations face novel ecological conditions, such as new environments and limited mates, and suffer from evolutionary consequences of demographic change through bottlenecks and mutation load accumulation. Self-fertilization is often observed at species range edges and, in addition to countering the lack of mates, is hypothesized as an evolutionary advantage against load accumulation through increased homozygosity and purging. We study how selfing impacts the accumulation of genetic load during range expansion via purging and/or speed of colonization. Using simulations, we disentangle inbreeding effects due to demography versus due to selfing and find that selfers expand faster, but still accumulate load, regardless of mating system. The severity of variants contributing to this load, however, differs across mating system: higher selfing rates purge large-effect recessive variants leaving a burden of smaller-effect alleles. We compare these predictions to the mixed-mating plant Arabis alpina, using whole-genome sequences from refugial outcrossing populations versus expanded selfing populations. Empirical results indicate accumulation of expansion load along with evidence of purging in selfing populations, concordant with our simulations, suggesting that while purging is a benefit of selfing evolving during range expansions, it is not sufficient to prevent load accumulation due to range expansion.

Conservation management strategy impacts inbreeding and mutation load in scimitar-horned oryx.

In an age of habitat loss and overexploitation, small populations, both captive and wild, are increasingly facing the effects of isolation and inbreeding. Genetic management has therefore become a vital tool for ensuring population viability. However, little is known about how the type and intensity of intervention shape the genomic landscape of inbreeding and mutation load. We address this using whole-genome sequence data of the scimitar-horned oryx (Oryx dammah), an iconic antelope that has been subject to contrasting management strategies since it was declared extinct in the wild. We show that unmanaged populations are enriched for long runs of homozygosity (ROH) and have significantly higher inbreeding coefficients than managed populations. Additionally, despite the total number of deleterious alleles being similar across management strategies, the burden of homozygous deleterious genotypes was consistently higher in unmanaged groups. These findings emphasize the risks associated with deleterious mutations through multiple generations of inbreeding. As wildlife management strategies continue to diversify, our study reinforces the importance of maintaining genome-wide variation in vulnerable populations and has direct implications for one of the largest reintroduction attempts in the world.

Increased homozygosity due to endogamy results in fitness consequences in a human population.

Recessive alleles have been shown to directly affect both human Mendelian disease phenotypes and complex traits. Pedigree studies also suggest that consanguinity results in increased childhood mortality and adverse health phenotypes, presumably through penetrance of recessive mutations. Here, we test whether the accumulation of homozygous, recessive alleles decreases reproductive success in a human population. We address this question among the Namibian Himba, an endogamous agro-pastoralist population, who until very recently practiced natural fertility. Using a sample of 681 individuals, we show that Himba exhibit elevated levels of "inbreeding," calculated as the fraction of the genome in runs of homozygosity (FROH). Many individuals contain multiple long segments of ROH in their genomes, indicating that their parents had high kinship coefficients. However, we do not find evidence that this is explained by first-cousin consanguinity, despite a reported social preference for cross-cousin marriages. Rather, we show that elevated haplotype sharing in the Himba is due to a bottleneck, likely in the past 60 generations. We test whether increased recessive mutation load results in observed fitness consequences by assessing the effect of FROH on completed fertility in a cohort of postreproductive women (n = 69). We find that higher FROH is significantly associated with lower fertility. Our data suggest a multilocus genetic effect on fitness driven by the expression of deleterious recessive alleles, especially those in long ROH. However, these effects are not the result of consanguinity but rather elevated background identity by descent.

Functional genomic diversity is correlated with neutral genomic diversity in populations of an endangered rattlesnake.

Theory predicts that genetic erosion in small, isolated populations of endangered species can be assessed using estimates of neutral genetic variation, yet this widely used approach has recently been questioned in the genomics era. Here, we leverage a chromosome-level genome assembly of an endangered rattlesnake (Sistrurus catenatus) combined with whole genome resequencing data (N = 110 individuals) to evaluate the relationship between levels of genome-wide neutral and functional diversity over historical and future timescales. As predicted, we found positive correlations between genome-wide estimates of neutral genetic diversity (π) and inferred levels of adaptive variation and an estimate of inbreeding mutation load, and a negative relationship between neutral diversity and an estimate of drift mutation load. However, these correlations were half as strong for projected future levels of neutral diversity based on contemporary effective population sizes. Broadly, our results confirm that estimates of neutral genetic diversity provide an accurate measure of genetic erosion in populations of a threatened vertebrate. They also provide nuance to the neutral-functional diversity controversy by suggesting that while these correlations exist, anthropogenetic impacts may have weakened these associations in the recent past and into the future.

Genome sequences and population genomics provide insights into the demographic history, inbreeding, and mutation load of two 'living fossil' tree species of Dipteronia.

'Living fossils', that is, ancient lineages of low taxonomic diversity, represent an exceptional evolutionary heritage, yet we know little about how demographic history and deleterious mutation load have affected their long-term survival and extinction risk. We performed whole-genome sequencing and population genomic analyses on Dipteronia sinensis and D. dyeriana, two East Asian Tertiary relict trees. We found large-scale genome reorganizations and identified species-specific genes under positive selection that are likely involved in adaptation. Our demographic analyses suggest that the wider-ranged D. sinensis repeatedly recovered from population bottlenecks over late Tertiary/Quaternary periods of adverse climate conditions, while the population size of the narrow-ranged D. dyeriana steadily decreased since the late Miocene, especially after the Last Glacial Maximum (LGM). We conclude that the efficient purging of deleterious mutations in D. sinensis facilitated its survival and repeated demographic recovery. By contrast, in D. dyeriana, increased genetic drift and reduced selection efficacy, due to recent severe population bottlenecks and a likely preponderance of vegetative propagation, resulted in fixation of strongly deleterious mutations, reduced fitness, and continuous population decline, with likely detrimental consequences for the species' future viability and adaptive potential. Overall, our findings highlight the significant impact of demographic history on levels of accumulation and purging of putatively deleterious mutations that likely determine the long-term survival and extinction risk of Tertiary relict trees.

Solanum pennellii (LA5240) backcross inbred lines (BILs) for high resolution mapping in tomato.

Wild species are an invaluable source of new traits for crop improvement. Over the years, the tomato community bred cultivated lines that carry introgressions from different species of the tomato tribe to facilitate trait discovery and mapping. The next phase in such projects is to find the genes that drive the identified phenotypes. This can be achieved by genotyping a few thousand individuals resulting in fine mapping that can potentially identify the causative gene. To couple trait discovery and fine mapping, we are presenting large, recombination-rich, Backcross Inbred Line (BIL) populations involving an unexplored accession of the wild, green-fruited species Solanum pennellii (LA5240; the 'Lost' Accession) with two modern tomato inbreds: LEA, determinate, and TOP, indeterminate. The LEA and TOP BILs are in BC2F6-8 generation and include 1400 and 500 lines, respectively. The BILs were genotyped with 5000 SPET markers, showing that in the euchromatic regions there was one recombinant every 17-18 Kb while in the heterochromatin a recombinant every 600-700 Kb (TOP and LEA respectively). To gain perspective on the topography of recombination we compared five independent members of the Self-pruning gene family with their respective neighboring genes; based on PCR markers, in all cases we found recombinants. Further mapping analysis of two known morphological mutations that segregated in the BILs (self-pruning and hairless) showed that the maximal delimited intervals were 73 Kb and 210 Kb, respectively, and included the known causative genes. The 'Lost'_BILs provide a solid framework to study traits derived from a drought-tolerant wild tomato.

Wright's Hierarchical F-Statistics.

This perspective article offers a meditation on FST and other quantities developed by Sewall Wright to describe the population structure, defined as any departure from reproduction through random union of gametes. Concepts related to the F-statistics draw from studies of the partitioning of variation, identity coefficients, and diversity measures. Relationships between the first two approaches have recently been clarified and unified. This essay addresses the third pillar of the discussion: Nei's GST and related measures. A hierarchy of probabilities of identity-by-state provides a description of the relationships among levels of a structured population with respect to genetic diversity. Explicit expressions for the identity-by-state probabilities are determined for models of structured populations undergoing regular inbreeding and recurrent mutation. Levels of genetic diversity within and between subpopulations reflect mutation as well as migration. Accordingly, indices of the population structure are inherently locus-specific, contrary to the intentions of Wright. Some implications of this locus-specificity are explored.

Selfing Promotes Spread and Introgression of Segregation Distorters in Hermaphroditic Plants.

Segregation distorters (SDs) are genetic elements that distort the Mendelian segregation ratio to favor their own transmission and are able to spread even when they incur fitness costs on organisms carrying them. Depending on the biology of the host organisms and the genetic architecture of the SDs, the population dynamics of SDs can be highly variable. Inbreeding is considered an effective mechanism for inhibiting the spread of SDs in populations, and can evolve as a defense mechanism against SDs in some systems. However, we show that inbreeding in the form of selfing in fact promotes the spread of SDs acting as pollen killers in a toxin-antidote system in hermaphroditic plants by two mechanisms: (i) By reducing the effective recombination rate between killer and antidote loci in the two-locus system and (ii) by increasing the proportion of SD alleles in individual flowers, rather than in the general gene-pool. We also show that in rice (Oryza sativa L.), a typical hermaphroditic plant, all molecularly characterized SDs associated with pollen killing were involved in population hybridization and have introgressed across different species. Paradoxically, these loci, which are associated with hybrid incompatibility and can be thought of as Bateson-Dobzhansky-Muller incompatibility loci are expected to reduce gene-flow between species, in fact cross species boundaries more frequently than random loci, and may act as important drivers of introgression.

Measuring the Efficiency of Purging by non-random Mating in Human Populations.

Human populations harbor a high concentration of deleterious genetic variants. Here, we tested the hypothesis that non-random mating practices affect the distribution of these variants, through exposure in the homozygous state, leading to their purging from the population gene pool. To do so, we produced whole-genome sequencing data for two pairs of Asian populations exhibiting different alliance rules and rates of inbreeding, but with similar effective population sizes. The results show that populations with higher rates of inbred matings do not purge deleterious variants more efficiently. Purging therefore has a low efficiency in human populations, and different mating practices lead to a similar mutational load.

Whole-genome resequencing of temporally stratified samples reveals substantial loss of haplotype diversity in the highly inbred Scandinavian wolf population.

Genetic drift can dramatically change allele frequencies in small populations and lead to reduced levels of genetic diversity, including loss of segregating variants. However, there is a shortage of quantitative studies of how genetic diversity changes over time in natural populations, especially on genome-wide scales. Here, we analyzed whole-genome sequences from 76 wolves of a highly inbred Scandinavian population, founded by only one female and two males, sampled over a period of 30 yr. We obtained chromosome-level haplotypes of all three founders and found that 10%-24% of their diploid genomes had become lost after about 20 yr of inbreeding (which approximately corresponds to five generations). Lost haplotypes spanned large genomic regions, as expected from the amount of recombination during this limited time period. Altogether, 160,000 SNP alleles became lost from the population, which may include adaptive variants as well as wild-type alleles masking recessively deleterious alleles. Although not sampled, we could indirectly infer that the two male founders had megabase-sized runs of homozygosity and that all three founders showed significant haplotype sharing, meaning that there were on average only 4.2 unique haplotypes in the six copies of each autosome that the founders brought into the population. This violates the assumption of unrelated founder haplotypes often made in conservation and management of endangered species. Our study provides a novel view of how whole-genome resequencing of temporally stratified samples can be used to visualize and directly quantify the consequences of genetic drift in a small inbred population.

Estimation of inbreeding and kinship coefficients via latent identity-by-descent states.

MOTIVATION: Estimating the individual inbreeding coefficient and pairwise kinship is an important problem in human genetics (e.g. in disease mapping) and in animal and plant genetics (e.g. inbreeding design). Existing methods, such as sample correlation-based genetic relationship matrix, KING, and UKin, are either biased, or not able to estimate inbreeding coefficients, or produce a large proportion of negative estimates that are difficult to interpret. This limitation of existing methods is partly due to failure to explicitly model inbreeding. Since all humans are inbred to various degrees by virtue of shared ancestries, it is prudent to account for inbreeding when inferring kinship between individuals. RESULTS: We present "Kindred," an approach that estimates inbreeding and kinship by modeling latent identity-by-descent states that accounts for all possible allele sharing-including inbreeding-between two individuals. Kindred used non-negative least squares method to fit the model, which not only increases computation efficiency compared to the maximum likelihood method, but also guarantees non-negativity of the kinship estimates. Through simulation, we demonstrate the high accuracy and non-negativity of kinship estimates by Kindred. By selecting a subset of SNPs that are similar in allele frequencies across different continental populations, Kindred can accurately estimate kinship between admixed samples. In addition, we demonstrate that the realized kinship matrix estimated by Kindred is effective in reducing genomic control values via linear mixed model in genome-wide association studies. Finally, we demonstrate that Kindred produces sensible heritability estimates on an Australian height dataset. AVAILABILITY AND IMPLEMENTATION: Kindred is implemented in C with multi-threading. It takes vcf file or stream as input and works seamlessly with bcftools. Kindred is freely available at https://github.com/haplotype/kindred.

Purging of deleterious burden in the endangered Iberian lynx.

SignificanceThe dynamics of deleterious variation under contrasting demographic scenarios remain poorly understood in spite of their relevance in evolutionary and conservation terms. Here we apply a genomic approach to study differences in the burden of deleterious alleles between the endangered Iberian lynx (Lynx pardinus) and the widespread Eurasian lynx (Lynx lynx). Our analysis unveils a significantly lower deleterious burden in the former species that should be ascribed to genetic purging, that is, to the increased opportunities of selection against recessive homozygotes due to the inbreeding caused by its smaller population size, as illustrated by our analytical predictions. This research provides theoretical and empirical evidence on the evolutionary relevance of genetic purging under certain demographic conditions.

Genetic diversity analysis of Inner Mongolia cashmere goats (Erlangshan subtype) based on whole genome re-sequencing.

BACKGROUND: Inner Mongolia cashmere goat (IMCG), renowned for its superior cashmere quality, is a Chinese indigenous goat breed that has been developed through natural and artificial selection over a long period. However, recently, the genetic resources of IMCGs have been significantly threatened by the introduction of cosmopolitan goat breeds and the absence of adequate breed protection systems. RESULTS: In order to assess the conservation effectiveness of IMCGs and efficiently preserve and utilize the purebred germplasm resources, this study analyzed the genetic diversity, kinship, family structure, and inbreeding of IMCGs utilizing resequencing data from 225 randomly selected individuals analyzed using the Plink (v.1.90), GCTA (v.1.94.1), and R (v.4.2.1) software. A total of 12,700,178 high-quality SNPs were selected through quality control from 34,248,064 SNP sites obtained from 225 individuals. The average minor allele frequency (MAF), polymorphic information content (PIC), and Shannon information index (SHI) were 0.253, 0.284, and 0.530, respectively. The average observed heterozygosity (Ho) and the average expected heterozygosity (He) were 0.355 and 0.351, respectively. The analysis of the identity by state distance matrix and genomic relationship matrix has shown that most individuals' genetic distance and genetic relationship are far away, and the inbreeding coefficient is low. The family structure analysis identified 10 families among the 23 rams. A total of 14,109 runs of homozygosity (ROH) were identified in the 225 individuals, with an average ROH length of 1014.547 kb. The average inbreeding coefficient, calculated from ROH, was 0.026 for the overall population and 0.027 specifically among the 23 rams, indicating a low level of inbreeding within the conserved population. CONCLUSIONS: The IMCGs exhibited moderate polymorphism and a low level of kinship with inbreeding occurring among a limited number of individuals. Simultaneously, it is necessary to prevent the loss of bloodline to guarantee the perpetuation of the IMCGs' germplasm resources.

Whole-genome resequencing provides insights into the diversity and adaptation to desert environment in Xinjiang Mongolian cattle.

BACKGROUND: Xinjiang Mongolian cattle is an indigenous breed that inhabits the Taklimakan Desert and is characterized by its small body size. However, the genomic diversity, origin, and genetic basis underlying the adaptation to the desert environment have been poorly studied. RESULTS: We analyzed patterns of Xinjiang Mongolian cattle genetic variation by sequencing 20 genomes together with seven previously sequenced genomes and comparing them to the 134 genomes of nine representative breeds worldwide. Among the breeds of Bos taurus, we found the highest nucleotide diversity (0.0024) associated with the lower inbreeding coefficient (2.0110-6), the lowest linkage disequilibrium (r2 = 0.3889 at distance of 10 kb), and the highest effective population size (181 at 20 generations ago) in Xinjiang Mongolian cattle. The genomic diversity pattern could be explained by a limited introgression of Bos indicus genes. More importantly, similarly to desert-adapted camel and same-habitat sheep, we also identified signatures of selection including genes, GO terms, and/or KEGG pathways controlling water reabsorption and osmoregulation, metabolic regulation and energy balance, as well as small body size in Xinjiang Mongolian cattle. CONCLUSIONS: Our results imply that Xinjiang Mongolian cattle might have acquired distinct genomic diversity by virtue of the introgression of Bos indicus, which helps understand the demographic history. The identification of selection signatures can provide novel insights into the genomic basis underlying the adaptation of Xinjiang Mongolian cattle to the desert environment.