Genome-wide Heterozygosity Research Articles

Genome report: a draft genome of Alliaria petiolata (garlic mustard) as a model system for invasion genetics.

The emerging field of invasion genetics examines the genetic causes and consequences of biological invasions, but few study systems are available that integrate deep ecological knowledge with genomic tools. Here, we report on the de novo assembly and annotation of a genome for the biennial herb Alliaria petiolata (M. Bieb.) Cavara and Grande (Brassicaceae), which is widespread in Eurasia and invasive across much of temperate North America. Our goal was to sequence and annotate a genome to complement resources available from hundreds of published ecological studies, a global field survey, and hundreds of genetic lines maintained in Germany and Canada. We sequenced a genotype (EFCC3-3-20) collected from the native range near Venice, Italy, and sequenced paired-end and mate pair libraries at ∼70 × coverage. A de novo assembly resulted in a highly continuous draft genome (N50 = 121 Mb; L50 = 2) with 99.7% of the 1.1 Gb genome mapping to scaffolds of at least 50 Kb in length. A total of 64,770 predicted genes in the annotated genome include 99% of plant BUSCO genes and 98% of transcriptome reads. Consistent with previous reports of (auto)hexaploidy in western Europe, we found that almost one-third of BUSCO genes (390/1440) mapped to two or more scaffolds despite <2% genome-wide average heterozygosity. The continuity and gene space quality of our draft assembly will enable molecular and functional genomic studies of A. petiolata to address questions relevant to invasion genetics and conservation strategies.

Runs of homozygosity in killer whale genomes provide a global record of demographic histories.

Runs of homozygosity (ROH) occur when offspring inherit haplotypes that are identical by descent from each parent. Length distributions of ROH are informative about population history; specifically, the probability of inbreeding mediated by mating system and/or population demography. Here, we investigated whether variation in killer whale (Orcinus orca) demographic history is reflected in genome-wide heterozygosity and ROH length distributions, using a global data set of 26genomes representative of geographic and ecotypic variation in this species, and two F1 admixed individuals with Pacific-Atlantic parentage. We first reconstructed demographic history for each population as changes in effective population size through time using the pairwise sequential Markovian coalescent (PSMC) method. We found a subset of populations declined in effective population size during the Late Pleistocene, while others had more stable demography. Genomes inferred to have undergone ancestral declines in effective population size, were autozygous at hundreds of short ROH (<1Mb), reflecting high background relatedness due to coalescence of haplotypes deep within the pedigree. In contrast, longer and therefore younger ROH (>1.5Mb) were found in low latitude populations, and populations of known conservation concern. These include a Scottish killer whale, for which 37.8% of the autosomes were comprised of ROH >1.5Mb in length. The fate of this population, in which only two adult males have been sighted in the past five years, and zero fecundity over the last two decades, may be inextricably linked to its demographic history and consequential inbreeding depression.

The endangered White Sands pupfish (Cyprinodon tularosa) genome reveals low diversity and heterogenous patterns of differentiation.

The White Sands pupfish (Cyprinodon tularosa), endemic to New Mexico in Southwestern North America, is of conservation concern due in part to invasive species, chemical pollution, and groundwater withdrawal. Herein, we developed a draft reference genome and use it to provide biological insights into the evolution and conservation of C. tularosa. We used our assembly to localize microsatellite markers previously used to demarcate evolutionary significant units (ESU), quantified genomic divergence and transposable element profiles between species, and compared C. tularosa genomic diversity related species. Our de novo assembly of PacBio Sequel II error-corrected reads resulted in a 1.08Gb draft genome with a contig N50 of 1.4Mb and 25,260 annotated protein coding genes, including 95% of the expected Actinopterygii conserved complete single-copy orthologues. Many of the C. tularosa microsatellite markers used for conservation assessments fell within, or near, genes and exhibited a pattern of increased heterozygosity near genic areas compared to those in intergenic regions. Nuclear alignments between these two species revealed 193genes contained in rapidly diverging tracts; transposable element profiles were largely concordant and suggest a shared, rapid expansion of LINE and Gypsy elements. Genome-wide heterozygosity was markedly lower in C. tularosa compared to estimates from other related species, probably because of smaller long-term effective population sizes constrained by their isolated and limited habitat. Overall, these inferences provide new insights into C. tularosa that should help inform future management efforts.

Immunocompetent birds choose larger breeding colonies.

Optimal size of social groups may vary between individuals, depending on their phenotypic traits, such as dominance status, age or personality. Larger social groups often enhance transmission rates of pathogens and should be avoided by individuals with poor immune defences. In contrast, more immunocompetent individuals are expected to take advantage of larger group sizes (e.g. better protection, information transfer) with smaller extra costs from pathogen or parasite pressure. Here, we hypothesized that immunocompetence may be a key determinant of group size choice and tested this hypothesis in a colonial waterbird, the common tern Sterna hirundo. We used a unique experimental framework, where formation of breeding colonies of different sizes was induced under uniform environmental conditions. For this purpose, different-size patches of attractive nesting substrate (artificial floating rafts) were provided at a single site with limited availability of natural nesting habitat. Colony size was identified as the only significant predictor of both innate (natural antibody-mediated complement activation) and adaptive (immunoglobulin concentrations) immunological traits in the common terns, as more immunocompetent birds settled in larger experimental colonies. In contrast, we found no significant associations between colony size and genetic diversity of key pathogen-recognition receptors, toll-like receptors (TLRs) and the Major Histocompatibility Complex (MHC) or genome-wide heterozygosity. We conclude that settlement decisions may be flexible within individuals and, thus, are likely to be primarily determined by the current immunological status, rather than fixed immunogenetic traits. Our study sheds new light on the complex interface between immunity and sociality in animals.

Efficacy and safety of rucaparib in previously treated, locally advanced or metastatic urothelial carcinoma from a phase 2, open-label trial (ATLAS)

BackgroundATLAS evaluated the efficacy and safety of the PARP inhibitor rucaparib in patients with previously treated locally advanced/unresectable or metastatic urothelial carcinoma (UC).MethodsPatients with UC were enrolled independent of tumor homologous recombination deficiency (HRD) status and received rucaparib 600 mg BID. The primary endpoint was investigator-assessed objective response rate (RECIST v1.1) in the intent-to-treat and HRD-positive (loss of genome-wide heterozygosity ≥10%) populations. Key secondary endpoints were progression-free survival (PFS) and safety. Disease control rate (DCR) was defined post-hoc as the proportion of patients with a confirmed complete or partial response (PR), or stable disease lasting ≥16 weeks.ResultsOf 97 enrolled patients, 20 (20.6%) were HRD-positive, 30 (30.9%) HRD-negative, and 47 (48.5%) HRD-indeterminate. Among 95 evaluable patients, there were no confirmed responses. However, reductions in the sum of target lesions were observed, including 6 (6.3%) patients with unconfirmed PR. DCR was 11.6%; median PFS was 1.8 months (95% CI, 1.6–1.9). No relationship was observed between HRD status and efficacy endpoints. Median treatment duration was 1.8 months (range, 0.1–10.1). Most frequent any-grade treatment-emergent adverse events were asthenia/fatigue (57.7%), nausea (42.3%), and anemia (36.1%). Of 64 patients with data from tumor tissue samples, 10 (15.6%) had a deleterious alteration in a DNA damage repair pathway gene, including four with a deleterious BRCA1 or BRCA2 alteration.ConclusionsRucaparib did not show significant activity in unselected patients with advanced UC regardless of HRD status. The safety profile was consistent with that observed in patients with ovarian or prostate cancer.Trial registrationThis trial was registered in ClinicalTrials.gov (NCT03397394). Date of registration: 12 January 2018. This trial was registered in EudraCT (2017–004166-10).

Chromosome-Level Assembly of the Atlantic Silverside Genome Reveals Extreme Levels of Sequence Diversity and Structural Genetic Variation.

The levels and distribution of standing genetic variation in a genome can provide a wealth of insights about the adaptive potential, demographic history, and genome structure of a population or species. As structural variants are increasingly associated with traits important for adaptation and speciation, investigating both sequence and structural variation is essential for wholly tapping this potential. Using a combination of shotgun sequencing, 10x Genomics linked reads and proximity-ligation data (Chicago and Hi-C), we produced and annotated a chromosome-level genome assembly for the Atlantic silverside (Menidia menidia)—an established ecological model for studying the phenotypic effects of natural and artificial selection—and examined patterns of genomic variation across two individuals sampled from different populations with divergent local adaptations. Levels of diversity varied substantially across each chromosome, consistently being highly elevated near the ends (presumably near telomeric regions) and dipping to near zero around putative centromeres. Overall, our estimate of the genome-wide average heterozygosity in the Atlantic silverside is among the highest reported for a fish, or any vertebrate (1.32–1.76% depending on inference method and sample). Furthermore, we also found extreme levels of structural variation, affecting ∼23% of the total genome sequence, including multiple large inversions (> 1 Mb and up to 12.6 Mb) associated with previously identified haploblocks showing strong differentiation between locally adapted populations. These extreme levels of standing genetic variation are likely associated with large effective population sizes and may help explain the remarkable adaptive divergence among populations of the Atlantic silverside.

Moving beyond assumptions: Polyploidy and environmental effects explain a geographical parthenogenesis scenario in European plants.

Geographical parthenogenesis (GP) describes the phenomenon that apomicts tend to have larger distribution areas and/or occur at higher altitudes or latitudes compared to sexual relatives. However, the complex effects of genome-wide heterozygosity, ploidy, reproduction mode (sexual vs. apomictic), and environment shaping GP of plants are still not well understood. We ascertained ploidy and reproduction mode by flow cytometry of 221 populations, and added genomic RADseq data (maximum 33,165 loci) of 80 taxa of the Ranunculusauricomus polyploid plant complex in temperate Europe. We observed 7% mainly diploid sexual, 28% facultative apomictic (mean sexuality 7.1%), and 65% obligate apomictic populations. Sexuals occupied a more southern, smaller distribution area, whereas apomicts expanded their range to higher latitudes. Within the complex, we detected three main genetic clusters and highly reticulate relationships. A genetically-informed path analysis using GLMMs revealed several significant relationships. Sexuality of populations (percent of sexual seeds) was higher in diploids compared to polyploids, associated with more petals, and similar between forests and open habitats. In contrast to other apomictic plant complexes, sexuality was mainly positively correlated to solar radiation and isothermality, which fits the southern distribution. We found up to three times higher heterozygosity in polyploids compared to diploids, and generally more heterozygous individuals in forests compared with open habitats. Interestingly, we revealed a previously unknown positive association between heterozygosity and temperature seasonality, suggesting a higher resistance of polyploids to more extreme climatic conditions. We provide empirical evidence for intrinsic and extrinsic factors shaping the GP pattern in a polyploid plant complex.

Testing a Local Inbreeding Hypothesis as a Cause of Observed Antler Characteristics in Managed Populations of White-Tailed Deer

The increased use of antler restrictions by state game agencies has led to a focus on antlers by the hunting public, particularly the potential for an association between genetics and antler characteristics. We analyzed microsatellite data from 1231 male white-tailed deer (Odocoileus virginianus) from three states (Oklahoma, Mississippi, and Texas) within USA to determine if genetic relatedness, internal relatedness (IR), homozygosity weighted by locus (HL), or correlations among uniting gametes (Fis) influenced total antler points, antler score, non-typical points or antler malformations. Within each location, deer in the lower and upper quartile intervals for number of antler points and score were unrelated (95% CI included 0 or was &lt;0) and relatively heterozygous for four measures of inbreeding. Antler score and points were positively influenced by age but negatively influenced by IR and HL, except for antler score in Mississippi. Relatedness, HL, IR and Fis did not differ between groups of deer with and without antler malformations. Perceived differences in antler quality do not appear to be affected by heterozygosity or a result of close inbreeding because we found deer were unrelated and measures of inbreeding and genome-wide heterozygosity were not correlated with antler characteristics.

High diversity, inbreeding and a dynamic Pleistocene demographic history revealed by African buffalo genomes

Genomes retain records of demographic changes and evolutionary forces that shape species and populations. Remnant populations of African buffalo (Syncerus caffer) in South Africa, with varied histories, provide an opportunity to investigate signatures left in their genomes by past events, both recent and ancient. Here, we produce 40 low coverage (7.14×) genome sequences of Cape buffalo (S. c. caffer) from four protected areas in South Africa. Genome-wide heterozygosity was the highest for any mammal for which these data are available, while differences in individual inbreeding coefficients reflected the severity of historical bottlenecks and current census sizes in each population. PSMC analysis revealed multiple changes in Ne between approximately one million and 20 thousand years ago, corresponding to paleoclimatic changes and Cape buffalo colonisation of southern Africa. The results of this study have implications for buffalo management and conservation, particularly in the context of the predicted increase in aridity and temperature in southern Africa over the next century as a result of climate change.

What evolutionary processes maintain MHC IIꞵ diversity within and among populations of stickleback?

Major Histocompatibility Complex (MHC) genes code for proteins that recognize foreign protein antigens to initiate T‐cell‐mediated adaptive immune responses. They are often the most polymorphic genes in vertebrate genomes. How evolution maintains this diversity remains of debate. Three main hypotheses seek to explain the maintenance of MHC diversity by invoking pathogen‐mediated selection: heterozygote advantage, frequency‐dependent selection, and fluctuating selection across landscapes or through time. Here, we use a large‐scale field parasite survey in a stickleback metapopulation to test predictions derived from each of these hypotheses. We identify over 1000 MHC IIβ variants (alleles spanning paralogous genes) and find that many of them covary positively or negatively with parasite load, suggesting that these genes contribute to resistance or susceptibility. However, despite our large sample‐size, we find no evidence for the widely cited stabilizing selection on MHC heterozygosity, in which individuals with an intermediate number of MHC variants have the lowest parasite burden. Nor do we observe a rare‐variant advantage, or widespread fluctuating selection across populations. In contrast, we find that MHC diversity is best predicted by neutral genome‐wide heterozygosity and between‐population genomic divergence, suggesting neutral processes are important in shaping the pattern of metapopulation MHC diversity. Thus, although MHC IIβ is highly diverse and relevant to the type and intensity of macroparasite infection in these populations of stickleback, the main models of MHC evolution still provide little explanatory power in this system.

Evaluating the correlation between genome-wide diversity and the release of plastic phenotypic variation in experimental translocations to novel natural environments.

Phenotypic reaction norms are often shaped and constrained by selection and are important for allowing organisms to respond to environmental change. However, selection cannot constrain reaction norms for environmental conditions that populations have not experienced. Consequently, cryptic neutral genetic variation for the reaction norm can accumulate such that a release of phenotypic variation occurs upon exposure to novel14 conditions. Most genomic diversity behaves as if functionally neutral. Therefore, genome-wide diversity metrics may correlate with levels of cryptic genetic variation and, as a result, exhibit a positive relationship with a release of phenotypic variation in novel environments. To test this hypothesis, we conducted translocations of juvenile brook trout (Salvelinus fontinalis) from 12 populations to novel uninhabited ponds that represented a gradient of environmental conditions. We assessed reaction norms for morphological traits (body size and four morphometric relative warps) across pond environmental gradients and evaluated the effect of genome-wide heterozygosity on phenotypic variability. All traits displayed plastic reaction norms. Overall, we found some evidence that a release of phenotypic variation consistent with cryptic genetic variation can occur in novel environmental conditions. However, the extent to which this release correlated with average genome-wide diversity was limited to only one of five traits examined. Our results suggest a limited link between genomic diversity26 and the accumulation of cryptic genetic variation in reaction norms. Similarly, reaction norms were constrained for many of the morphological traits examined. Past conditions may have constrained reaction norms in the putatively novel environments despite significant deviations from contemporary source population habitat. Additionally, as a generalist colonizing species brook trout may exhibit plastic phenotypes across a wide range of environmental conditions.

Chromosome-Scale Assembly and Annotation of the Macadamia Genome (Macadamia integrifolia HAES 741).

Macadamia integrifolia is a representative of the large basal eudicot family Proteaceae and the main progenitor species of the Australian native nut crop macadamia. Since its commercialisation in Hawaii fewer than 100 years ago, global production has expanded rapidly. However, genomic resources are limited in comparison to other horticultural crops. The first draft assembly of M. integrifolia had good coverage of the functional gene space but its high fragmentation has restricted its use in comparative genomics and association studies. Here we have generated an improved assembly of cultivar HAES 741 (4,094 scaffolds, 745 Mb, N50 413 kb) using a combination of Illumina paired and PacBio long read sequences. Scaffolds were anchored to 14 pseudo-chromosomes using seven genetic linkage maps. This assembly has improved contiguity and coverage, with >120 Gb of additional sequence. Following annotation, 34,274 protein-coding genes were predicted, representing 90% of the expected gene content. Our results indicate that the macadamia genome is repetitive and heterozygous. The total repeat content was 55% and genome-wide heterozygosity, estimated by read mapping, was 0.98% or an average of one SNP per 102 bp. This is the first chromosome-scale genome assembly for macadamia and the Proteaceae. It is expected to be a valuable resource for breeding, gene discovery, conservation and evolutionary genomics.

A demonstration of conservation genomics for threatened species management.

As species extinction rates increase, genomics provides a powerful tool to support intensive management of threatened species. We use the Tasmanian devil (Sarcophilus harrisii) to demonstrate how conservation genomics can be implemented in threatened species management. We conducted whole genome sequencing (WGS) of 25 individuals from the captive breeding programme and reduced-representation sequencing (RRS) of 98 founders of the same programme. A subset of the WGS samples was also sequenced by RRS, allowing us to directly compare genome-wide heterozygosity with estimates from RRS data. We found good congruence in interindividual variation and gene-ontology classifications between the two data sets, indicating that our RRS data reflect the genome well. We also attempted genome-wide association studies with both data sets (regarding breeding success), but the genomic data suffered from small sample size, while the RRS data suffered from lack of precision, highlighting a key trade-off in the design of conservation genomic research. Nevertheless, we identified a number of candidate genes that may be associated with variation in breeding success. Individual heterozygosity, as measured by WGS or RRS, was not associated with breeding success in captivity but was negatively associated with litter sizes of breeding females in the RRS data set. Our findings enable conservation managers to have confidence in RRS data while understanding its limitations, and provide avenues for further investigation into which processes underlie variation in breeding success in captive Tasmanian devils. We caution, however, that deep functional insights using RRS may be impaired by a lack of precision, especially when marker density is low.

Genome-wide insights into introgression and its consequences for genome-wide heterozygosity in the Mytilus species complex across Europe.

The three mussel species comprising the Mytilus complex are widespread across Europe and readily hybridize when they occur in sympatry, resulting in a mosaic of populations with varying genomic backgrounds. Two of these species, M. edulis and M. galloprovincialis, are extensively cultivated across Europe, with annual production exceeding 230,000 tonnes. The third species, M. trossulus, is considered commercially damaging as hybridization with this species results in weaker shells and poor meat quality. We therefore used restriction site associated DNA sequencing to generate high‐resolution insights into the structure of the Mytilus complex across Europe and to resolve patterns of introgression. Inferred species distributions were concordant with the results of previous studies based on smaller numbers of genetic markers, with M. edulis and M. galloprovincialis predominating in northern and southern Europe respectively, while introgression between these species was most pronounced in northern France and the Shetland Islands. We also detected traces of M. trossulus ancestry in several northern European populations, especially around the Baltic and in northern Scotland. Finally, genome‐wide heterozygosity, whether quantified at the population or individual level, was lowest in M. edulis, intermediate in M. galloprovincialis, and highest in M. trossulus, while introgression was positively associated with heterozygosity in M. edulis but negatively associated with heterozygosity in M. galloprovincialis. Our study will help to inform mussel aquaculture by providing baseline information on the genomic backgrounds of different Mytilus populations across Europe and by elucidating the effects of introgression on genome‐wide heterozygosity, which is known to influence commercially important traits such as growth, viability, and fecundity in mussels.

De novo mutations across 1,465 diverse genomes reveal mutational insights and reductions in the Amish founder population

De novo mutations (DNMs), or mutations that appear in an individual despite not being seen in their parents, are an important source of genetic variation whose impact is relevant to studies of human evolution, genetics, and disease. Utilizing high-coverage whole-genome sequencing data as part of the Trans-Omics for Precision Medicine (TOPMed) Program, we called 93,325 single-nucleotide DNMs across 1,465 trios from an array of diverse human populations, and used them to directly estimate and analyze DNM counts, rates, and spectra. We find a significant positive correlation between local recombination rate and local DNM rate, and that DNM rate explains a substantial portion (8.98 to 34.92%, depending on the model) of the genome-wide variation in population-level genetic variation from 41K unrelated TOPMed samples. Genome-wide heterozygosity does correlate with DNM rate, but only explains <1% of variation. While we are underpowered to see small differences, we do not find significant differences in DNM rate between individuals of European, African, and Latino ancestry, nor across ancestrally distinct segments within admixed individuals. However, we did find significantly fewer DNMs in Amish individuals, even when compared with other Europeans, and even after accounting for parental age and sequencing center. Specifically, we found significant reductions in the number of C→A and T→C mutations in the Amish, which seem to underpin their overall reduction in DNMs. Finally, we calculated near-zero estimates of narrow sense heritability (h2), which suggest that variation in DNM rate is significantly shaped by nonadditive genetic effects and the environment.

Joint Estimates of Heterozygosity and Runs of Homozygosity for Modern and Ancient Samples.

Both the total amount and the distribution of heterozygous sites within individual genomes are informative about the genetic diversity of the population they belong to. Detecting true heterozygous sites in ancient genomes is complicated by the generally limited coverage achieved and the presence of post-mortem damage inflating sequencing errors. Additionally, large runs of homozygosity found in the genomes of particularly inbred individuals and of domestic animals can skew estimates of genome-wide heterozygosity rates. Current computational tools aimed at estimating runs of homozygosity and genome-wide heterozygosity levels are generally sensitive to such limitations. Here, we introduce ROHan, a probabilistic method which substantially improves the estimate of heterozygosity rates both genome-wide and for genomic local windows. It combines a local Bayesian model and a Hidden Markov Model at the genome-wide level and can work both on modern and ancient samples. We show that our algorithm outperforms currently available methods for predicting heterozygosity rates for ancient samples. Specifically, ROHan can delineate large runs of homozygosity (at megabase scales) and produce a reliable confidence interval for the genome-wide rate of heterozygosity outside of such regions from modern genomes with a depth of coverage as low as 5–6× and down to 7–8× for ancient samples showing moderate DNA damage. We apply ROHan to a series of modern and ancient genomes previously published and revise available estimates of heterozygosity for humans, chimpanzees and horses.

Too much of a good thing? Finding the most informative genetic data set to answer conservation questions.

Molecular markers are a useful tool allowing conservation and population managers to shed light on genetic processes affecting threatened populations. However, as technological advancements in molecular techniques continue to evolve, conservationists are frequently faced with new genetic markers, each with nuanced variation in their characteristics as well as advantages and disadvantages for informing various questions. We used a well-studied population of Tasmanian devils (Sarcophilus harrisii) from Maria Island, Tasmania, to illustrate the issues associated with combining multiple genetic data sets and to help answer a question posed by many population managers: which data set will provide the most precise and accurate estimates of the population processes we are trying to measure? We analysed individual heterozygosity (as internal relatedness, IR) of 96 individuals, calculated using four genetic marker types (putatively neutral microsatellites, major histocompatibility complex-linked microsatellites, reduced representation sequencing, and candidate region resequencing). We found no correlation in IR values across marker types, suggesting that various genetic markers reflect different aspects of genomic diversity. In addition, some marker types were more informative than others for conservation decision-making. Reduced representation sequencing provided the highest precision (lowest error) for estimating population-level genetic diversity, and most closely reflected genome-wide heterozygosity both theoretically and empirically. Within the conservation context, our results highlight important considerations when choosing a molecular technique for wildlife genetics.

Is genomic diversity a useful proxy for census population size? Evidence from a species-rich community of desert lizards.

Species abundance data are critical for testing ecological theory, but obtaining accurate empirical estimates for many taxa is challenging. Proxies for species abundance can help researchers circumvent time and cost constraints that are prohibitive for long-term sampling. Under simple demographic models, genetic diversity is expected to correlate with census size, such that genome-wide heterozygosity may provide a surrogate measure of species abundance. We tested whether nucleotide diversity is correlated with long-term estimates of abundance, occupancy and degree of ecological specialization in a diverse lizard community from arid Australia. Using targeted sequence capture, we obtained estimates of genomic diversity from 30 species of lizards, recovering an average of 5,066 loci covering 3.6Mb of DNA sequence per individual. We compared measures of individual heterozygosity to a metric of habitat specialization to investigate whether ecological preference exerts a measurable effect on genetic diversity. We find that heterozygosity is significantly correlated with species abundance and occupancy, but not habitat specialization. Demonstrating the power of genomic sampling, the correlation between heterozygosity and abundance/occupancy emerged from considering just one or two individuals per species. However, genetic diversity does no better at predicting abundance than a single day of traditional sampling in this community. We conclude that genetic diversity is a useful proxy for regional-scale species abundance and occupancy, but a large amount of unexplained variation in heterozygosity suggests additional constraints or a failure of ecological sampling to adequately capture variation in true population size.

Heterozygosity and asymmetry: Ectodysplasin as a form of genetic stress in marine threespine stickleback.

Genome-wide heterozygosity has long been hypothesized to play a role in buffering organisms against developmental perturbations, potentially resulting in increased symmetry. If true, this could in part explain the maintenance of standing genetic variation in wild populations. Marine threespine sticklebacks (Gasterosteus aculeatus) were sampled across their eastern Pacific coastal distribution from Alaska to California and variations in asymmetry for both structural and nonstructural armor traits (lateral plates) were assessed. Structural plates consistently showed less asymmetry than nonstructural plates, but standardized measures of heterozygosity were not correlated with the extent of asymmetry expressed by a fish. Fish that were heterozygous for the major-effect gene controlling lateral plate variation (Ectodysplasin) had higher occurrences of asymmetry, even when the individuals were phenotypically fully plated. Collectively, this suggests that heterozygosity at a major-effect locus can have a greater impact on asymmetry than heterozygosity sampled across the genome.

A male-expressed rice embryogenic trigger redirected for asexual propagation through seeds.

The molecular pathways that trigger the initiation of embryogenesis after fertilization in flowering plants, and prevent its occurrence without fertilization, are not well understood1. Here we show in rice (Oryza sativa) that BABY BOOM1 (BBM1), a member of the AP2 family2 of transcription factors that is expressed in sperm cells, has a key role in this process. Ectopic expression of BBM1 in the egg cell is sufficient for parthenogenesis, which indicates that a single wild-type gene can bypass the fertilization checkpoint in the female gamete. Zygotic expression of BBM1 is initially specific to the male allele but is subsequently biparental, and this is consistent with its observed auto-activation. Triple knockout of the genes BBM1, BBM2 and BBM3 causes embryo arrest and abortion, which are fully rescued by male-transmitted BBM1. These findings suggest that the requirement for fertilization in embryogenesis is mediated by male-genome transmission of pluripotency factors. When genome editing to substitute mitosis for meiosis (MiMe)3,4 is combined with the expression of BBM1 in the egg cell, clonal progeny can be obtained that retain genome-wide parental heterozygosity. The synthetic asexual-propagation trait is heritable through multiple generations of clones. Hybrid crops provide increased yields that cannot be maintained by their progeny owing to genetic segregation. This work establishes the feasibility of asexual reproduction in crops, and could enable the maintenance of hybrids clonally through seed propagation5,6.