Pine cones that remain closed and retain seeds until fire causes the cones to open (cone serotiny) represent a key adaptive trait in a variety of pine species. In lodgepole pine, there is substantial ...geographical variation in serotiny across the Rocky Mountain region. This variation in serotiny has evolved as a result of geographically divergent selection, with consequences that extend to forest communities and ecosystems. An understanding of the genetic architecture of this trait is of interest owing to the wide‐reaching ecological consequences of serotiny and also because of the repeated evolution of the trait across the genus. Here, we present and utilize an inexpensive and time‐effective method for generating population genomic data. The method uses restriction enzymes and PCR amplification to generate a library of fragments that can be sequenced with a high level of multiplexing. We obtained data for more than 95 000 single nucleotide polymorphisms across 98 serotinous and nonserotinous lodgepole pines from three populations. We used a Bayesian generalized linear model (GLM) to test for an association between genotypic variation at these loci and serotiny. The probability of serotiny varied by genotype at 11 loci, and the association between genotype and serotiny at these loci was consistent in each of the three populations of pines. Genetic variation across these 11 loci explained 50% of the phenotypic variation in serotiny. Our results provide a first genome‐wide association map of serotiny in pines and demonstrate an inexpensive and efficient method for generating population genomic data.
See also the Perspective by Eckert and Dyer
• Arbuscular mycorrhizal (AM) fungi form endosymbioses with most plants, and they themselves are hosts for Mollicutes/Mycoplasma-related endobacteria (MRE). Despite their significance, genomic ...information for AM fungi and their MRE are relatively sparse, which hinders our understanding of their biology and evolution.
• We assembled the genomes of the AM fungus Diversispora epigaea (formerly Glomus versiforme) and its MRE and performed comparative genomics and evolutionary analyses.
• The D. epigaea genome showed a pattern of substantial gene duplication and differential evolution of gene families, including glycosyltransferase family 25, whose activities are exclusively lipopolysaccharide biosynthesis. Genes acquired by horizontal transfer from bacteria possibly function in defense against foreign DNA or viruses. The MRE population was diverse, with multiple genomes displaying characteristics of differential evolution and encoding many MRE-specific genes as well as genes of AM fungal origin. Gene family expansion in D. epigaea may enhance adaptation to both external and internal environments, such as expansion of kinases for signal transduction upon external stimuli and expansion of nucleoside salvage pathway genes potentially for competition with MRE, whose genomes lack purine and pyrimidine biosynthetic pathways.
• Collectively, this metagenome provides high-quality references and begins to reveal the diversity within AM fungi and their MRE.
Meiotic crossovers (COs) are not uniformly distributed across the genome. Factors affecting this phenomenon are not well understood. Although many species exhibit large differences in CO numbers ...between sexes, sex-specific aspects of CO landscape are particularly poorly elucidated. Here, we conduct high-resolution CO mapping in maize. Our results show that CO numbers as well as their overall distribution are similar in male and female meioses. There are, nevertheless, dissimilarities at local scale. Male and female COs differ in their locations relative to transcription start sites in gene promoters and chromatin marks, including nucleosome occupancy and tri-methylation of lysine 4 of histone H3 (H3K4me3). Our data suggest that sex-specific factors not only affect male-female CO number disparities but also cause fine differences in CO positions. Differences between male and female CO landscapes indicate that recombination has distinct implications for population structure and gene evolution in male and in female meioses.
Oomycete plant pathogens are difficult to control and routine genetic research is challenging. A major problem is instability of isolates. Here we characterize >600 field and single zoospore isolates ...of Phytophthora capsici for inheritance of mating type, sensitivity to mefenoxam, chromosome copy number and heterozygous allele frequencies. The A2 mating type was highly unstable with 26% of 241 A2 isolates remaining A2. The A1 mating type was stable. Isolates intermediately resistant to mefenoxam produced fully resistant single-spore progeny. Sensitive isolates remained fully sensitive. Genome re-sequencing of single zoospore isolates revealed extreme aneuploidy; a phenomenon dubbed Dynamic Extreme Aneuploidy (DEA). DEA is characterized by the asexual inheritance of diverse intra-genomic combinations of chromosomal ploidy ranging from 2N to 3N and heterozygous allele frequencies that do not strictly correspond to ploidy. Isolates sectoring on agar media showed dramatically altered heterozygous allele frequencies. DEA can explain the rapid increase of advantageous alleles (e.g. drug resistance), mating type switches and copy neutral loss of heterozygosity (LOH). Although the mechanisms driving DEA are unknown, it can play an important role in adaptation and evolution and seriously hinders all aspects of P. capsici research.
Of 7028 disorders with suspected Mendelian inheritance, 1139 are recessive and have an established molecular basis. Although individually uncommon, Mendelian diseases collectively account for ~20% of ...infant mortality and ~10% of pediatric hospitalizations. Preconception screening, together with genetic counseling of carriers, has resulted in remarkable declines in the incidence of several severe recessive diseases including Tay-Sachs disease and cystic fibrosis. However, extension of preconception screening to most severe disease genes has hitherto been impractical. Here, we report a preconception carrier screen for 448 severe recessive childhood diseases. Rather than costly, complete sequencing of the human genome, 7717 regions from 437 target genes were enriched by hybrid capture or microdroplet polymerase chain reaction, sequenced by next-generation sequencing (NGS) to a depth of up to 2.7 gigabases, and assessed with stringent bioinformatic filters. At a resultant 160x average target coverage, 93% of nucleotides had at least 20x coverage, and mutation detection/genotyping had ~95% sensitivity and ~100% specificity for substitution, insertion/deletion, splicing, and gross deletion mutations and single-nucleotide polymorphisms. In 104 unrelated DNA samples, the average genomic carrier burden for severe pediatric recessive mutations was 2.8 and ranged from 0 to 7. The distribution of mutations among sequenced samples appeared random. Twenty-seven percent of mutations cited in the literature were found to be common polymorphisms or misannotated, underscoring the need for better mutation databases as part of a comprehensive carrier testing strategy. Given the magnitude of carrier burden and the lower cost of testing compared to treating these conditions, carrier screening by NGS made available to the general population may be an economical way to reduce the incidence of and ameliorate suffering associated with severe recessive childhood disorders.
Medicago truncatula is a model for investigating legume genetics, including the genetics and evolution of legume–rhizobia symbiosis. We used whole-genome sequence data to identify and characterize ...sequence polymorphisms and linkage disequilibrium (LD) in a diverse collection of 26 M. truncatula accessions. Our analyses reveal that M. truncatula harbors both higher diversity and less LD than soybean (Glycine max) and exhibits patterns of LD and recombination similar to Arabidopsis thaliana. The population-scaled recombination rate is approximately one-third of the mutation rate, consistent with expectations for a species with a high selfing rate. Linkage disequilibrium, however, is not extensive, and therefore, the low recombination rate is likely not a major constraint to adaptation. Nucleotide diversity in 100-kb windows was negatively correlated with gene density, which is expected if diversity is shaped by selection acting against slightly deleterious mutations. Among putative coding regions, members of four gene families harbor significantly higher diversity than the genome-wide average. Three of these families are involved in resistance against pathogens; one of these families, the nodule-specific, cysteine-rich gene family, is specific to the galegoid legumes and is involved in control of rhizobial differentiation. The more than 3 million SNPs that we detected, approximately one-half of which are present in more than one accession, are a valuable resource for genome-wide association mapping of genes responsible for phenotypic diversity in legumes, especially traits associated with symbiosis and nodulation.
Previous studies exploring sequence variation in the model legume, Medicago truncatula, relied on mapping short reads to a single reference. However, read-mapping approaches are inadequate to examine ...large, diverse gene families or to probe variation in repeat-rich or highly divergent genome regions. De novo sequencing and assembly of M. truncatula genomes enables near-comprehensive discovery of structural variants (SVs), analysis of rapidly evolving gene families, and ultimately, construction of a pan-genome.
Genome-wide synteny based on 15 de novo M. truncatula assemblies effectively detected different types of SVs indicating that as much as 22% of the genome is involved in large structural changes, altogether affecting 28% of gene models. A total of 63 million base pairs (Mbp) of novel sequence was discovered, expanding the reference genome space for Medicago by 16%. Pan-genome analysis revealed that 42% (180 Mbp) of genomic sequences is missing in one or more accession, while examination of de novo annotated genes identified 67% (50,700) of all ortholog groups as dispensable - estimates comparable to recent studies in rice, maize and soybean. Rapidly evolving gene families typically associated with biotic interactions and stress response were found to be enriched in the accession-specific gene pool. The nucleotide-binding site leucine-rich repeat (NBS-LRR) family, in particular, harbors the highest level of nucleotide diversity, large effect single nucleotide change, protein diversity, and presence/absence variation. However, the leucine-rich repeat (LRR) and heat shock gene families are disproportionately affected by large effect single nucleotide changes and even higher levels of copy number variation.
Analysis of multiple M. truncatula genomes illustrates the value of de novo assemblies to discover and describe structural variation, something that is often under-estimated when using read-mapping approaches. Comparisons among the de novo assemblies also indicate that different large gene families differ in the architecture of their structural variation.
Parasite host switches may trigger disease emergence, but prehistoric host ranges are often unknowable. Lymphatic filariasis and loiasis are major human diseases caused by the insect-borne filarial ...nematodes Brugia, Wuchereria and Loa. Here we show that the genomes of these nematodes and seven tropical bird lineages exclusively share a novel retrotransposon, AviRTE, resulting from horizontal transfer (HT). AviRTE subfamilies exhibit 83-99% nucleotide identity between genomes, and their phylogenetic distribution, paleobiogeography and invasion times suggest that HTs involved filarial nematodes. The HTs between bird and nematode genomes took place in two pantropical waves, >25-22 million years ago (Myr ago) involving the Brugia/Wuchereria lineage and >20-17 Myr ago involving the Loa lineage. Contrary to the expectation from the mammal-dominated host range of filarial nematodes, we hypothesize that these major human pathogens may have independently evolved from bird endoparasites that formerly infected the global breadth of avian biodiversity.
Meiosis is a critical process in the reproduction and life cycle of flowering plants in which homologous chromosomes pair, synapse, recombine and segregate. Understanding meiosis will not only ...advance our knowledge of the mechanisms of genetic recombination, but also has substantial applications in crop improvement. Despite the tremendous progress in the past decade in other model organisms (e.g., Saccharomyces cerevisiae and Drosophila melanogaster), the global identification of meiotic genes in flowering plants has remained a challenge due to the lack of efficient methods to collect pure meiocytes for analyzing the temporal and spatial gene expression patterns during meiosis, and for the sensitive identification and quantitation of novel genes.
A high-throughput approach to identify meiosis-specific genes by combining isolated meiocytes, RNA-Seq, bioinformatic and statistical analysis pipelines was developed. By analyzing the studied genes that have a meiosis function, a pipeline for identifying meiosis-specific genes has been defined. More than 1,000 genes that are specifically or preferentially expressed in meiocytes have been identified as candidate meiosis-specific genes. A group of 55 genes that have mitochondrial genome origins and a significant number of transposable element (TE) genes (1,036) were also found to have up-regulated expression levels in meiocytes.
These findings advance our understanding of meiotic genes, gene expression and regulation, especially the transcript profiles of MGI genes and TE genes, and provide a framework for functional analysis of genes in meiosis.
Mosquitoes rely on RNA interference (RNAi) as their primary defense against viral infections. To this end, the combination of RNAi and invertebrate cell culture systems has become an invaluable tool ...in studying virus-vector interactions. Nevertheless, a recent study failed to detect an active RNAi response to West Nile virus (WNV) infection in C6/36 (Aedes albopictus) cells, a mosquito cell line frequently used to study arthropod-borne viruses (arboviruses). Therefore, we sought to determine if WNV actively evades the host's RNAi response or if C6/36 cells have a dysfunctional RNAi pathway. C6/36 and Drosophila melanogaster S2 cells were infected with WNV (Flaviviridae), Sindbis virus (SINV, Togaviridae) and La Crosse virus (LACV, Bunyaviridae) and total RNA recovered from cell lysates. Small RNA (sRNA) libraries were constructed and subjected to high-throughput sequencing. In S2 cells, virus-derived small interfering RNAs (viRNAs) from all three viruses were predominantly 21 nt in length, a hallmark of the RNAi pathway. However, in C6/36 cells, viRNAs were primarily 17 nt in length from WNV infected cells and 26-27 nt in length in SINV and LACV infected cells. Furthermore, the origin (positive or negative viral strand) and distribution (position along viral genome) of S2 cell generated viRNA populations was consistent with previously published studies, but the profile of sRNAs isolated from C6/36 cells was altered. In total, these results suggest that C6/36 cells lack a functional antiviral RNAi response. These findings are analogous to the type-I interferon deficiency described in Vero (African green monkey kidney) cells and suggest that C6/36 cells may fail to accurately model mosquito-arbovirus interactions at the molecular level.