Pre-existing heterogeneity facilitates development of heteroresistance upon gene acquisition.

Antibiotic resistance causes 1.27 million global deaths annually and is predicted to worsen. Heteroresistance is a form of resistance in which only a minor and unstable subpopulation of cells of a bacterial isolate are resistant to a given antibiotic, and are therefore often undetected by clinical diagnostics. These infrequent and undetected resistant cells can be selected during antibiotic therapy, expand in number, and cause unexplained treatment failures. A major question is how heteroresistance evolves. Here, studying the antibiotic fosfomycin, we report that heteroresistance can develop from a pre-existing state of phenotypic heterogeneity in which an isolate harbors a subpopulation with increased minimum inhibitory concentration (MIC), but below the clinical resistance breakpoint. We call this phenomenon heterosusceptibility and demonstrate that acquisition of a resistance gene, fosA, increases the MIC of the subpopulation beyond the breakpoint, making the isolate heteroresistant. Conversely, deletion of fosA from a heteroresistant isolate led to reduction of the MIC of the resistant subpopulation without a loss of heterogeneity, thus generating heterosusceptibility. A survey of 103 carbapenem-resistant Enterobacterales (CRE) revealed that the Escherichia sp. isolates lacked the fosA gene and uniformly exhibited fosfomycin heterosusceptibility, whereas the Klebsiella and Enterobacter encoded the fosA gene and were almost exclusively heteroresistant. Furthermore, some isolates exhibited heterosusceptibility to other antibiotics, demonstrating that this is a widespread phenomenon. These results highlight a mechanism for the evolution of heteroresistance and suggest that surveillance for heterosusceptibility may facilitate the prediction of impending heteroresistance before it evolves.

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.

1,3-Dioxane-Linked Novel Bacterial Topoisomerase Inhibitors: Expanding Structural Diversity and the Antibacterial Spectrum.

Antibacterial resistance continues its devastation of available therapies. Novel bacterial topoisomerase inhibitors (NBTIs) offer one solution to this critical issue. Two series of amine NBTIs bearing tricyclic DNA-binding moieties as well as amide NBTIs with a bicyclic DNA-binding moiety were synthesized and evaluated against methicillin-resistant Staphylococcus aureus (MRSA). Additionally, these compounds and a series of bicyclic amine analogues displayed high activity against susceptible and drug-resistant Neisseria gonorrhoeae, expanding the spectrum of these dioxane-linked NBTIs.

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.

Colistin Heteroresistance Is Largely Undetected among Carbapenem-Resistant Enterobacterales in the United States.

Heteroresistance is a form of antibiotic resistance where a bacterial strain is comprised of a minor resistant subpopulation and a majority susceptible subpopulation. We showed previously that colistin heteroresistance can mediate the failure of colistin therapy in an in vivo infection model, even for isolates designated susceptible by clinical diagnostics. We sought to characterize the extent of colistin heteroresistance among the highly drug-resistant carbapenem-resistant Enterobacterales (CRE). We screened 408 isolates for colistin heteroresistance. These isolates were collected between 2012 and 2015 in eight U.S. states as part of active surveillance for CRE. Colistin heteroresistance was detected in 10.1% (41/408) of isolates, and it was more common than conventional homogenous resistance (7.1%, 29/408). Most (93.2%, 38/41) of these heteroresistant isolates were classified as colistin susceptible by standard clinical diagnostic testing. The frequency of colistin heteroresistance was greatest in 2015, the last year of the study. This was especially true among Enterobacter isolates, of which specific species had the highest rates of heteroresistance. Among Klebsiella pneumoniae isolates, which were the majority of isolates tested, there was a closely related cluster of colistin-heteroresistant ST-258 isolates found mostly in Georgia. However, cladistic analysis revealed that, overall, there was significant diversity in the genetic backgrounds of heteroresistant K. pneumoniae isolates. These findings suggest that due to being largely undetected in the clinic, colistin heteroresistance among CRE is underappreciated in the United States.IMPORTANCE Heteroresistance is an underappreciated phenomenon that may be the cause of some unexplained antibiotic treatment failures. Misclassification of heteroresistant isolates as susceptible may lead to inappropriate therapy. Heteroresistance to colistin was more common than conventional resistance and was overwhelmingly misclassified as susceptibility by clinical diagnostic testing. Higher proportions of colistin heteroresistance observed in certain Enterobacter species and clustering among heteroresistant Klebsiella pneumoniae strains may inform colistin treatment recommendations. Overall, the rate of colistin nonsusceptibility was more than double the level detected by clinical diagnostics, suggesting that the prevalence of colistin nonsusceptibility among CRE may be higher than currently appreciated in the United States.

Antibiotic-Selected Gene Amplification Heightens Metal Resistance.

The increasing frequency of antibiotic resistance poses myriad challenges to modern medicine. Environmental survival of multidrug-resistant bacteria in health care facilities, including hospitals, creates reservoirs for transmission of these difficult to treat pathogens. To prevent bacterial colonization, these facilities deploy an array of infection control measures, including bactericidal metals on surfaces, as well as implanted devices. Although antibiotics are routinely used in these health care environments, it is unknown whether and how antibiotic exposure affects metal resistance. We identified a multidrug-resistant Enterobacter clinical isolate that displayed heteroresistance to the antibiotic colistin, where only a minor fraction of cells within the population resist the drug. When this isolate was grown in the presence of colistin, a 9-kb DNA region was duplicated in the surviving resistant subpopulation, but surprisingly, was not required for colistin heteroresistance. Instead, the amplified region included a three-gene locus (ncrABC) that conferred resistance to the bactericidal metal, nickel. ncrABC expression alone was sufficient to confer nickel resistance to E. coli K-12. Due to its selection for the colistin-resistant subpopulation harboring the duplicated 9-kb region that includes ncrABC, colistin treatment led to enhanced nickel resistance. Taken together, these data suggest that the use of antibiotics may inadvertently promote enhanced resistance to antimicrobial metals, with potentially profound implications for bacterial colonization and transmission in the health care environment.IMPORTANCE To inhibit bacterial transmission and infection, health care facilities use bactericidal metal coatings to prevent colonization of surfaces and implanted devices. In these environments, antibiotics are commonly used, but their effect on metal resistance is unclear. The data described here reveal that exposure of a human isolate of Enterobacter cloacae to a last-line antibiotic, colistin, resulted in a DNA amplification that does not confer antibiotic resistance but instead facilitates resistance to the toxic metal nickel. This highlights a novel aspect of antibiotic and metal interplay. Concerningly, these data suggest the use of antibiotics could in some cases promote bacterial survival and colonization in the health care environment and ultimately increase transmission and infection of patients.

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.

Antibiotic combinations that exploit heteroresistance to multiple drugs effectively control infection.

Antibiotic-resistant bacteria are a significant threat to human health, with one estimate suggesting they will cause 10 million worldwide deaths per year by 2050, surpassing deaths due to cancer1. Because new antibiotic development can take a decade or longer, it is imperative to effectively use currently available drugs. Antibiotic combination therapy offers promise for treating highly resistant bacterial infections, but the factors governing the sporadic efficacy of such regimens have remained unclear. Dogma suggests that antibiotics ineffective as monotherapy can be effective in combination2. Here, using carbapenem-resistant Enterobacteriaceae (CRE) clinical isolates, we reveal the underlying basis for the majority of effective combinations to be heteroresistance. Heteroresistance is a poorly understood mechanism of resistance reported for different classes of antibiotics3-6 in which only a subset of cells are phenotypically resistant7. Within an isolate, the subpopulations resistant to different antibiotics were distinct, and over 88% of CRE isolates exhibited heteroresistance to multiple antibiotics ('multiple heteroresistance'). Combinations targeting multiple heteroresistance were efficacious, whereas those targeting homogenous resistance were ineffective. Two pan-resistant Klebsiella isolates were eradicated by combinations targeting multiple heteroresistance, highlighting a rational strategy to identify effective combinations that employs existing antibiotics and could be clinically implemented immediately.

Thiol Starvation Induces Redox-Mediated Dysregulation of Escherichia coli Biofilm Components.

A hallmark of bacterial biofilms is the production of an extracellular matrix (ECM) that encases and protects the community from environmental stressors. Biofilm formation is an integral portion of the uropathogenic Escherichia coli (UPEC) life cycle. Approximately 2% of UPEC isolates are cysteine auxotrophs. Here, we investigated how cysteine homeostasis impacted UPEC UTI89 strain biofilm formation and, specifically, the production of the ECM components curli and cellulose. Cysteine auxotrophs produced less cellulose and slightly more curli compared to wild-type (WT) strains, and cysteine auxotrophs formed smooth, nonrugose colonies. Cellulose production was restored in cysteine auxotrophs when YfiR was inactivated. YfiR is a redox-sensitive regulator of the diguanylate cyclase, YfiN. The production of curli, a temperature-regulated appendage, was independent of temperature in UTI89 cysteine auxotrophs. In a screen of UPEC isolates, we found that ∼60% of UPEC cysteine auxotrophs produced curli at 37°C, but only ∼2% of cysteine prototrophic UPEC isolates produced curli at 37°C. Interestingly, sublethal concentrations of amdinocillin and trimethoprim-sulfamethoxazole inhibited curli production, whereas strains auxotrophic for cysteine continued to produce curli even in the presence of amdinocillin and trimethoprim-sulfamethoxazole. The dysregulation of ECM components and resistance to amdinocillin in cysteine auxotrophs may be linked to hyperoxidation, since the addition of exogenous cysteine or glutathione restored WT biofilm phenotypes to mutants unable to produce cysteine and glutathione.IMPORTANCE Uropathogenic Escherichia coli (UPEC) bacteria are the predominant causative agent of urinary tract infections (UTIs). UTIs account for billions of dollars of financial burden annually to the health care industry in the United States. Biofilms are an important aspect of the UPEC pathogenesis cascade and for the establishment of chronic infections. Approximately 2% of UPEC isolates from UTIs are cysteine auxotrophs, yet there is relatively little known about the biofilm formation of UPEC cysteine auxotrophs. Here we show that cysteine auxotrophs have dysregulated biofilm components due to a change in the redox state of the periplasm. Additionally, we show the relationship between cysteine auxotrophs, biofilms, and antibiotics frequently used to treat UTIs.

The Catabolite Repressor Protein-Cyclic AMP Complex Regulates csgD and Biofilm Formation in Uropathogenic Escherichia coli.

The extracellular matrix protects Escherichia coli from immune cells, oxidative stress, predation, and other environmental stresses. Production of the E. coli extracellular matrix is regulated by transcription factors that are tuned to environmental conditions. The biofilm master regulator protein CsgD upregulates curli and cellulose, the two major polymers in the extracellular matrix of uropathogenic E. coli (UPEC) biofilms. We found that cyclic AMP (cAMP) regulates curli, cellulose, and UPEC biofilms through csgD The alarmone cAMP is produced by adenylate cyclase (CyaA), and deletion of cyaA resulted in reduced extracellular matrix production and biofilm formation. The catabolite repressor protein (CRP) positively regulated csgD transcription, leading to curli and cellulose production in the UPEC isolate, UTI89. Glucose, a known inhibitor of CyaA activity, blocked extracellular matrix formation when added to the growth medium. The mutant strains ΔcyaA and Δcrp did not produce rugose biofilms, pellicles, curli, cellulose, or CsgD. Three putative CRP binding sites were identified within the csgD-csgB intergenic region, and purified CRP could gel shift the csgD-csgB intergenic region. Additionally, we found that CRP binded upstream of kpsMT, which encodes machinery for K1 capsule production. Together our work shows that cAMP and CRP influence E. coli biofilms through transcriptional regulation of csgD IMPORTANCE The catabolite repressor protein (CRP)-cyclic AMP (cAMP) complex influences the transcription of ∼7% of genes on the Escherichia coli chromosome (D. Zheng, C. Constantinidou, J. L. Hobman, and S. D. Minchin, Nucleic Acids Res 32:5874-5893, 2004, https://dx.doi.org/10.1093/nar/gkh908). Glucose inhibits E. coli biofilm formation, and ΔcyaA and Δcrp mutants show impaired biofilm formation (D. W. Jackson, J.W. Simecka, and T. Romeo, J Bacteriol 184:3406-3410, 2002, https://dx.doi.org/10.1128/JB.184.12.3406-3410.2002). We determined that the cAMP-CRP complex regulates curli and cellulose production and the formation of rugose and pellicle biofilms through csgD Additionally, we propose that cAMP may work as a signaling compound for uropathogenic E. coli (UPEC) to transition from the bladder lumen to inside epithelial cells for intracellular bacterial community formation through K1 capsule regulation.

The Biology of the Escherichia coli Extracellular Matrix.

Escherichia coli is one of the world's best-characterized organisms, because it has been extensively studied for over a century. However, most of this work has focused on E. coli grown under laboratory conditions that do not faithfully simulate its natural environments. Therefore, the historical perspectives on E. coli physiology and life cycle are somewhat skewed toward experimental systems that feature E. coli growing logarithmically in a test tube. Typically a commensal bacterium, E. coli resides in the lower intestines of a slew of animals. Outside of the lower intestine, E. coli can adapt and survive in a very different set of environmental conditions. Biofilm formation allows E. coli to survive, and even thrive, in environments that do not support the growth of planktonic populations. E. coli can form biofilms virtually everywhere: in the bladder during a urinary tract infection, on in-dwelling medical devices, and outside of the host on plants and in the soil. The E. coli extracellular matrix (ECM), primarily composed of the protein polymer named curli and the polysaccharide cellulose, promotes adherence to organic and inorganic surfaces and resistance to desiccation, the host immune system, and other antimicrobials. The pathways that govern E. coli biofilm formation, cellulose production, and curli biogenesis will be discussed in this article, which concludes with insights into the future of E. coli biofilm research and potential therapies.

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.