Characterizing a century of genetic diversity and contemporary antigenic diversity of N1 neuraminidase in influenza A virus from North American swine.

Influenza A viruses (IAVs) of the H1N1 classical swine lineage became endemic in North American swine following the 1918 pandemic. Additional human-to-swine transmission events after 1918, and a spillover of H1 viruses from wild birds in Europe, potentiated a rapid increase in genomic diversity via reassortment between introductions and the endemic classical swine lineage. To determine mechanisms affecting reassortment and evolution, we conducted a phylogenetic analysis of N1 and paired HA swine IAV genes in North America between 1930 and 2020. We described fourteen N1 clades within the N1 Eurasian avian lineage (including the N1 pandemic clade), the N1 classical swine lineage, and the N1 human seasonal lineage. Seven N1 genetic clades had evidence for contemporary circulation. To assess antigenic drift associated with N1 genetic diversity, we generated a panel of representative swine N1 antisera and quantified the antigenic distance between wild-type viruses using enzyme-linked lectin assays and antigenic cartography. Within the N1 genes, the antigenic similarity was variable and reflected shared evolutionary history. Sustained circulation and evolution of N1 genes in swine had resulted in a significant antigenic distance between the N1 pandemic clade and the classical swine lineage. Between 2010 and 2020, N1 clades and N1-HA pairings fluctuated in detection frequency across North America, with hotspots of diversity generally appearing and disappearing within 2 years. We also identified frequent N1-HA reassortment events (n = 36), which were rarely sustained (n = 6) and sometimes also concomitant with the emergence of new N1 genetic clades (n = 3). These data form a baseline from which we can identify N1 clades that expand in range or genetic diversity that may impact viral phenotypes or vaccine immunity and subsequently the health of North American swine.

Swine-to-Ferret Transmission of Antigenically Drifted Contemporary Swine H3N2 Influenza A Virus Is an Indicator of Zoonotic Risk to Humans.

Human-to-swine transmission of influenza A (H3N2) virus occurs repeatedly and plays a critical role in swine influenza A virus (IAV) evolution and diversity. Human seasonal H3 IAVs were introduced from human-to-swine in the 1990s in the United States and classified as 1990.1 and 1990.4 lineages; the 1990.4 lineage diversified into 1990.4.A-F clades. Additional introductions occurred in the 2010s, establishing the 2010.1 and 2010.2 lineages. Human zoonotic cases with swine IAV, known as variant viruses, have occurred from the 1990.4 and 2010.1 lineages, highlighting a public health concern. If a variant virus is antigenically drifted from current human seasonal vaccine (HuVac) strains, it may be chosen as a candidate virus vaccine (CVV) for pandemic preparedness purposes. We assessed the zoonotic risk of US swine H3N2 strains by performing phylogenetic analyses of recent swine H3 strains to identify the major contemporary circulating genetic clades. Representatives were tested in hemagglutination inhibition assays with ferret post-infection antisera raised against existing CVVs or HuVac viruses. The 1990.1, 1990.4.A, and 1990.4.B.2 clade viruses displayed significant loss in cross-reactivity to CVV and HuVac antisera, and interspecies transmission potential was subsequently investigated in a pig-to-ferret transmission study. Strains from the three lineages were transmitted from pigs to ferrets via respiratory droplets, but there were differential shedding profiles. These data suggest that existing CVVs may offer limited protection against swine H3N2 infection, and that contemporary 1990.4.A viruses represent a specific concern given their widespread circulation among swine in the United States and association with multiple zoonotic cases.

Interspecies Transmission from Pigs to Ferrets of Antigenically Distinct Swine H1 Influenza A Viruses with Reduced Reactivity to Candidate Vaccine Virus Antisera as Measures of Relative Zoonotic Risk.

During the last decade, endemic swine H1 influenza A viruses (IAV) from six different genetic clades of the hemagglutinin gene caused zoonotic infections in humans. The majority of zoonotic events with swine IAV were restricted to a single case with no subsequent transmission. However, repeated introduction of human-seasonal H1N1, continual reassortment between endemic swine IAV, and subsequent drift in the swine host resulted in highly diverse swine IAV with human-origin genes that may become a risk to the human population. To prepare for the potential of a future swine-origin IAV pandemic in humans, public health laboratories selected candidate vaccine viruses (CVV) for use as vaccine seed strains. To assess the pandemic risk of contemporary US swine H1N1 or H1N2 strains, we quantified the genetic diversity of swine H1 HA genes, and identified representative strains from each circulating clade. We then characterized the representative swine IAV against human seasonal vaccine and CVV strains using ferret antisera in hemagglutination inhibition assays (HI). HI assays revealed that 1A.3.3.2 (pdm09) and 1B.2.1 (delta-2) demonstrated strong cross reactivity to human seasonal vaccines or CVVs. However, swine IAV from three clades that represent more than 50% of the detected swine IAVs in the USA showed significant reduction in cross-reactivity compared to the closest CVV virus: 1A.1.1.3 (alpha-deletion), 1A.3.3.3-clade 3 (gamma), and 1B.2.2.1 (delta-1a). Representative viruses from these three clades were further characterized in a pig-to-ferret transmission model and shown to exhibit variable transmission efficiency. Our data prioritize specific genotypes of swine H1N1 and H1N2 to further investigate in the risk they pose to the human population.

A Chromosome-Level Genome of the Camphor Tree and the Underlying Genetic and Climatic Factors for Its Top-Geoherbalism.

Camphor tree [Cinnamomum camphora (L.) J. Presl], a species in the magnoliid family Lauraceae, is known for its rich volatile oils and is used as a medical cardiotonic and as a scent in many perfumed hygiene products. Here, we present a high-quality chromosome-scale genome of C. camphora with a scaffold N50 of 64.34 Mb and an assembled genome size of 755.41 Mb. Phylogenetic inference revealed that the magnoliids are a sister group to the clade of eudicots and monocots. Comparative genomic analyses identified two rounds of ancient whole-genome duplication (WGD). Tandem duplicated genes exhibited a higher evolutionary rate, a more recent evolutionary history and a more clustered distribution on chromosomes, contributing to the production of secondary metabolites, especially monoterpenes and sesquiterpenes, which are the principal essential oil components. Three-dimensional analyses of the volatile metabolites, gene expression and climate data of samples with the same genotype grown in different locations showed that low temperature and low precipitation during the cold season modulate the expression of genes in the terpenoid biosynthesis pathways, especially TPS genes, which facilitates the accumulation of volatile compounds. Our study lays a theoretical foundation for policy-making regarding the agroforestry applications of camphor tree.

octoFLUshow: an Interactive Tool Describing Spatial and Temporal Trends in the Genetic Diversity of Influenza A Virus in U.S. Swine.

Influenza A virus (IAV) is passively surveilled in swine in the United States through a U.S. Department of Agriculture administered surveillance system. We present an interactive Web tool to visualize and explore trends in the genetic and geographic diversity of IAV derived from the surveillance system.

De novo assembly, annotation, and comparative analysis of 26 diverse maize genomes.

We report de novo genome assemblies, transcriptomes, annotations, and methylomes for the 26 inbreds that serve as the founders for the maize nested association mapping population. The number of pan-genes in these diverse genomes exceeds 103,000, with approximately a third found across all genotypes. The results demonstrate that the ancient tetraploid character of maize continues to degrade by fractionation to the present day. Excellent contiguity over repeat arrays and complete annotation of centromeres revealed additional variation in major cytological landmarks. We show that combining structural variation with single-nucleotide polymorphisms can improve the power of quantitative mapping studies. We also document variation at the level of DNA methylation and demonstrate that unmethylated regions are enriched for cis-regulatory elements that contribute to phenotypic variation.

SequelTools: a suite of tools for working with PacBio Sequel raw sequence data.

BACKGROUND: PacBio sequencing is an incredibly valuable third-generation DNA sequencing method due to very long read lengths, ability to detect methylated bases, and its real-time sequencing methodology. Yet, hitherto no tool was available for analyzing the quality of, subsampling, and filtering PacBio data. RESULTS: Here we present SequelTools, a command-line program containing three tools: Quality Control, Read Subsampling, and Read Filtering. The Quality Control tool quickly processes PacBio Sequel raw sequence data from multiple SMRTcells producing multiple statistics and publication-quality plots describing the quality of the data including N50, read length and count statistics, PSR, and ZOR. The Read Subsampling tool allows the user to subsample reads by one or more of the following criteria: longest subreads per CLR or random CLR selection. The Read Filtering tool provides options for normalizing data by filtering out certain low-quality scraps reads and/or by minimum CLR length. SequelTools is implemented in bash, R, and Python using only standard libraries and packages and is platform independent. CONCLUSIONS: SequelTools is a program that provides the only free, fast, and easy-to-use quality control tool, and the only program providing this kind of read subsampling and read filtering for PacBio Sequel raw sequence data, and is available at https://github.com/ISUgenomics/SequelTools .

Effect of sequence depth and length in long-read assembly of the maize inbred NC358.

Improvements in long-read data and scaffolding technologies have enabled rapid generation of reference-quality assemblies for complex genomes. Still, an assessment of critical sequence depth and read length is important for allocating limited resources. To this end, we have generated eight assemblies for the complex genome of the maize inbred line NC358 using PacBio datasets ranging from 20 to 75 × genomic depth and with N50 subread lengths of 11-21 kb. Assemblies with ≤30 × depth and N50 subread length of 11 kb are highly fragmented, with even low-copy genic regions showing degradation at 20 × depth. Distinct sequence-quality thresholds are observed for complete assembly of genes, transposable elements, and highly repetitive genomic features such as telomeres, heterochromatic knobs, and centromeres. In addition, we show high-quality optical maps can dramatically improve contiguity in even our most fragmented base assembly. This study provides a useful resource allocation reference to the community as long-read technologies continue to mature.

Molecular Evidence for Functional Divergence and Decay of a Transcription Factor Derived from Whole-Genome Duplication in Arabidopsis thaliana.

Functional divergence between duplicate transcription factors (TFs) has been linked to critical events in the evolution of land plants and can result from changes in patterns of expression, binding site divergence, and/or interactions with other proteins. Although plant TFs tend to be retained post polyploidization, many are lost within tens to hundreds of million years. Thus, it can be hypothesized that some TFs in plant genomes are in the process of becoming pseudogenes. Here, we use a pair of salt tolerance-conferring transcription factors, DWARF AND DELAYED FLOWERING1 (DDF1) and DDF2, that duplicated through paleopolyploidy 50 to 65 million years ago, as examples to illustrate potential mechanisms leading to duplicate retention and loss. We found that the expression patterns of Arabidopsis thaliana (At)DDF1 and AtDDF2 have diverged in a highly asymmetric manner, and AtDDF2 has lost most inferred ancestral stress responses. Consistent with promoter disablement, the AtDDF2 promoter has fewer predicted cis-elements and a methylated repetitive element. Through comparisons of AtDDF1, AtDDF2, and their Arabidopsis lyrata orthologs, we identified significant differences in binding affinities and binding site preference. In particular, an AtDDF2-specific substitution within the DNA-binding domain significantly reduces binding affinity. Cross-species analyses indicate that both AtDDF1 and AtDDF2 are under selective constraint, but among A. thaliana accessions, AtDDF2 has a higher level of nonsynonymous nucleotide diversity compared with AtDDF1. This may be the result of selection in different environments or may point toward the possibility of ongoing functional decay despite retention for millions of years after gene duplication.

Retained duplicate genes in green alga Chlamydomonas reinhardtii tend to be stress responsive and experience frequent response gains.

BACKGROUND: Green algae belong to a group of photosynthetic organisms that occupy diverse habitats, are closely related to land plants, and have been studied as sources of food and biofuel. Although multiple green algal genomes are available, a global comparative study of algal gene families has not been carried out. To investigate how gene families and gene expression have evolved, particularly in the context of stress response that have been shown to correlate with gene family expansion in multiple eukaryotes, we characterized the expansion patterns of gene families in nine green algal species, and examined evolution of stress response among gene duplicates in Chlamydomonas reinhardtii. RESULTS: Substantial variation in domain family sizes exists among green algal species. Lineage-specific expansion of families occurred throughout the green algal lineage but inferred gene losses occurred more often than gene gains, suggesting a continuous reduction of algal gene repertoire. Retained duplicates tend to be involved in stress response, similar to land plant species. However, stress responsive genes tend to be pseudogenized as well. When comparing ancestral and extant gene stress response state, we found that response gains occur in 13% of duplicate gene branches, much higher than 6% in Arabidopsis thaliana. CONCLUSION: The frequent gains of stress response among green algal duplicates potentially reflect a high rate of innovation, resulting in a species-specific gene repertoire that contributed to adaptive response to stress. This could be further explored towards deciphering the mechanism of stress response, and identifying suitable green algal species for oil production.

The origin and evolution of maize in the Southwestern United States.

The origin of maize (Zea mays mays) in the US Southwest remains contentious, with conflicting archaeological data supporting either coastal(1-4) or highland(5,6) routes of diffusion of maize into the United States. Furthermore, the genetics of adaptation to the new environmental and cultural context of the Southwest is largely uncharacterized(7). To address these issues, we compared nuclear DNA from 32 archaeological maize samples spanning 6,000 years of evolution to modern landraces. We found that the initial diffusion of maize into the Southwest about 4,000 years ago is likely to have occurred along a highland route, followed by gene flow from a lowland coastal maize beginning at least 2,000 years ago. Our population genetic analysis also enabled us to differentiate selection during domestication for adaptation to the climatic and cultural environment of the Southwest, identifying adaptation loci relevant to drought tolerance and sugar content.

Consequences of Whole-Genome Triplication as Revealed by Comparative Genomic Analyses of the Wild Radish Raphanus raphanistrum and Three Other Brassicaceae Species.

Polyploidization events are frequent among flowering plants, and the duplicate genes produced via such events contribute significantly to plant evolution. We sequenced the genome of wild radish (Raphanus raphanistrum), a Brassicaceae species that experienced a whole-genome triplication event prior to diverging from Brassica rapa. Despite substantial gene gains in these two species compared with Arabidopsis thaliana and Arabidopsis lyrata, ∼70% of the orthologous groups experienced gene losses in R. raphanistrum and B. rapa, with most of the losses occurring prior to their divergence. The retained duplicates show substantial divergence in sequence and expression. Based on comparison of A. thaliana and R. raphanistrum ortholog floral expression levels, retained radish duplicates diverged primarily via maintenance of ancestral expression level in one copy and reduction of expression level in others. In addition, retained duplicates differed significantly from genes that reverted to singleton state in function, sequence composition, expression patterns, network connectivity, and rates of evolution. Using these properties, we established a statistical learning model for predicting whether a duplicate would be retained postpolyploidization. Overall, our study provides new insights into the processes of plant duplicate loss, retention, and functional divergence and highlights the need for further understanding factors controlling duplicate gene fate.

MAKER-P: a tool kit for the rapid creation, management, and quality control of plant genome annotations.

We have optimized and extended the widely used annotation engine MAKER in order to better support plant genome annotation efforts. New features include better parallelization for large repeat-rich plant genomes, noncoding RNA annotation capabilities, and support for pseudogene identification. We have benchmarked the resulting software tool kit, MAKER-P, using the Arabidopsis (Arabidopsis thaliana) and maize (Zea mays) genomes. Here, we demonstrate the ability of the MAKER-P tool kit to automatically update, extend, and revise the Arabidopsis annotations in light of newly available data and to annotate pseudogenes and noncoding RNAs absent from The Arabidopsis Informatics Resource 10 build. Our results demonstrate that MAKER-P can be used to manage and improve the annotations of even Arabidopsis, perhaps the best-annotated plant genome. We have also installed and benchmarked MAKER-P on the Texas Advanced Computing Center. We show that this public resource can de novo annotate the entire Arabidopsis and maize genomes in less than 3 h and produce annotations of comparable quality to those of the current The Arabidopsis Information Resource 10 and maize V2 annotation builds.

Evolution of TPS20-related terpene synthases influences chemical diversity in the glandular trichomes of the wild tomato relative Solanum habrochaites.

A systematic screen of volatile terpene production in the glandular trichomes of 79 accessions of Solanum habrochaites was conducted and revealed the presence of 21 mono- and sesquiterpenes that exhibit a range of qualitative and quantitative variation. Hierarchical clustering identified distinct terpene phenotypic modules with shared patterns of terpene accumulation across accessions. Several terpene modules could be assigned to previously identified terpene synthase (TPS) activities that included members of the TPS-e/f subfamily that utilize the unusual cis-prenyl diphosphate substrates neryl diphosphate and 2z,6z-farnesyl diphosphate. DNA sequencing and in vitro enzyme activity analysis of TPS-e/f members from S. habrochaites identified three previously unassigned enzyme activities that utilize these cisoid substrates. These produce either the monoterpenes α-pinene and limonene, or the sesquiterpene 7-epizingiberene, with the in vitro analyses that recapitulated the trichome chemistry found in planta. Comparison of the distribution of S. habrochaites accessions with terpene content revealed a strong preference for the presence of particular TPS20 alleles at distinct geographic locations. This study reveals that the unusually high intra-specific variation of volatile terpene synthesis in glandular trichomes of S. habrochaites is due at least in part to evolution at the TPS20 locus.