Genome-Wide Association of Carbon and Nitrogen Metabolism in the Maize Nested Association Mapping Population

Carbon (C) and nitrogen (N) metabolism are critical to plant growth and development and are at the basis of crop yield and adaptation. We performed high-throughput metabolite analyses on over 12,000 samples from the nested association mapping population to identify genetic variation in C and N metabolism in maize (Zea mays ssp. mays). All samples were grown in the same field and used to identify natural variation controlling the levels of 12 key C and N metabolites, namely chlorophyll a, chlorophyll b, fructose, fumarate, glucose, glutamate, malate, nitrate, starch, sucrose, total amino acids, and total protein, along with the first two principal components derived from them. Our genome-wide association results frequently identified hits with single-gene resolution. In addition to expected genes such as invertases, natural variation was identified in key C₄ metabolism genes, including carbonic anhydrases and a malate transporter. Unlike several prior maize studies, extensive pleiotropy was found for C and N metabolites. This integration of field-derived metabolite data with powerful mapping and genomics resources allows for the dissection of key metabolic pathways, providing avenues for future genetic improvement.Carbon (C) and nitrogen (N) metabolism are the basis for life on Earth. The production, balance, and tradeoffs of C and N metabolism are critical to all plant growth, yield, and local adaptation (Coruzzi and Bush, 2001; Coruzzi et al., 2007). In plants, there is a critical balance between the tissues that are producing energy (sources) and those using it (sinks), as the identities and locations of these vary through time and developmental stage (Smith et al., 2004). While a great deal of research has focused on the key genes and proteins involved in these processes (Wang et al., 1993; Kim et al., 2000; Takahashi et al., 2009), relatively little is known about the natural variation within a species that fine-tunes these processes in individual plants.In addition, a key aspect of core C metabolism involves the nature of plant photosynthesis. While the majority of plants use standard C₃ photosynthetic pathways, some, including maize (Zea mays) and many other grasses, use C₄ photosynthesis to concentrate CO₂ in bundle sheath cells to avoid wasteful photorespiration (Sage, 2004). Under some conditions (such as drought or high temperatures), C₄ photosynthesis is much more efficient than C₃ photosynthesis. Since these conditions are expected to become more prevalent in the near future due to climate change, various research groups are working to convert C₃ crop species to C₄ metabolism in order to boost crop production and food security (Sage and Zhu, 2011). Beyond this, better understanding of both C₃ and C₄ metabolic pathways will aid efforts to breed crops for superior yield, N-use efficiency, and other traits important for global food production.In the last two decades, quantitative trait locus (QTL) mapping, first with linkage analysis and later with association mapping, has been used to dissect C and N metabolism in several species, including Arabidopsis (Arabidopsis thaliana; Mitchell-Olds and Pedersen, 1998; Keurentjes et al., 2008; Lisec et al., 2008; Sulpice et al., 2009), tomato (Solanum lycopersicum; Schauer et al., 2006), and maize (Hirel et al., 2001; Limami et al., 2002; Zhang et al., 2006, 2010a, 2010b). These studies identified key genetic regions underlying variation in core C and N metabolism, many of which include candidate genes known to be involved in these processes.Previous studies of genetic variation for C and N metabolism are limited by the fact that they identified trait loci only through linkage mapping in artificial families or through association mapping across populations of unrelated individuals. Linkage mapping benefits from high statistical power due to many individuals sharing the same genotype at any given location, but it suffers from low resolution due to the limited number of generations (and hence recombination events) since the initial founders. Association mapping, in turn, enjoys high resolution due to the long recombination histories of natural populations but suffers from low power, since most genotypes occur in only a few individuals. In addition, many of these studies focused on C and N in artificial settings (e.g. greenhouses or growth chambers) instead of field conditions, running the risk that important genetic loci could be missed if the conditions do not include important (and potentially unknown) natural environmental variables.To address these issues and improve our understanding of C and N metabolism in maize, we used a massive and diverse germplasm resource, the maize nested association mapping (NAM) population (Buckler et al., 2009; McMullen et al., 2009), to evaluate genetic variation underlying the accumulation of 12 targeted metabolites in maize leaf tissue under field conditions. This population was formed by mating 25 diverse maize lines to the reference line, B73, and creating a 200-member biparental family from each of these crosses. The entire 5,000-member NAM population thus combines the strengths of both linkage and association mapping (McMullen et al., 2009), and it has been used to identify QTLs for important traits such as flowering time (Buckler et al., 2009), disease resistance (Kump et al., 2011; Poland et al., 2011), and plant architecture (Tian et al., 2011; Peiffer et al., 2013). Most importantly, this combination of power and resolution frequently resolves associations down to the single-gene level, even when using field-based data.The metabolites we profiled are key indicators of photosynthesis, respiration, glycolysis, and protein and sugar metabolism in the plant (Sulpice et al., 2009). By taking advantage of a robotized metabolic phenotyping platform (Gibon et al., 2004), we performed more than 100,000 assays across 12,000 samples, with two independent samples per experimental plot. Raw data and the best linear unbiased predictors (BLUPs) of these data were included as part of a study of general functional variation in maize (Wallace et al., 2014), but, to our knowledge, this is the first in-depth analysis of these metabolic data. We find strong correlations among several of the metabolites, and we also find extensive pleiotropy among the different traits. Many of the top QTLs are also near or within candidate genes relating to C and N metabolism, thus identifying targets for future breeding and selection. These results provide a powerful resource for those working with core C and N metabolism in plants and for improving maize performance in particular.     

Neutralising antibodies for Mayaro virus in Pantanal, Brazil

The Pantanal hosts diverse wildlife species and therefore is a hotspot for arbovirus studies in South America. A serosurvey for Mayaro virus (MAYV), eastern (EEEV), western (WEEV) and Venezuelan (VEEV) equine encephalitis viruses was conducted with 237 sheep, 87 free-ranging caimans and 748 equids, including 37 collected from a ranch where a neurologic disorder outbreak had been recently reported. Sera were tested for specific viral antibodies using plaque-reduction neutralisation test. From a total of 748 equids, of which 264 were immunised with vaccine composed of EEEV and WEEV and 484 had no history of immunisation, 10 (1.3%) were seropositive for MAYV and two (0.3%) for VEEV using criteria of a ≥ 4-fold antibody titre difference. Among the 484 equids without history of immunisation, 48 (9.9%) were seropositive for EEEV and four (0.8%) for WEEV using the same criteria. Among the sheep, five were sero- positive for equine encephalitis alphaviruses, with one (0.4%) for EEEV, one (0.4%) for WEEV and three (1.3%) for VEEV. Regarding free-ranging caimans, one (1.1%) and three (3.4%), respectively, had low titres for neutralising antibodies to VEEV and undetermined alphaviruses. The neurological disorder outbreak could not be linked to the alphaviruses tested. Our findings represent strong evidence that MAYV and all equine encephalitis alphaviruses circulated in the Pantanal.     

Gut Symbionts from Distinct Hosts Exhibit Genotoxic Activity via Divergent Colibactin Biosynthesis Pathways

Secondary metabolites produced by nonribosomal peptide synthetase (NRPS) or polyketide synthase (PKS) pathways are chemical mediators of microbial interactions in diverse environments. However, little is known about their distribution, evolution, and functional roles in bacterial symbionts associated with animals. A prominent example is colibactin, a largely unknown family of secondary metabolites produced by Escherichia coli via a hybrid NRPS-PKS biosynthetic pathway that inflicts DNA damage upon eukaryotic cells and contributes to colorectal cancer and tumor formation in the mammalian gut. Thus far, homologs of this pathway have only been found in closely related Enterobacteriaceae, while a divergent variant of this gene cluster was recently discovered in a marine alphaproteobacterial Pseudovibrio strain. Herein, we sequenced the genome of Frischella perrara PEB0191, a bacterial gut symbiont of honey bees and identified a homologous colibactin biosynthetic pathway related to those found in Enterobacteriaceae. We show that the colibactin genomic island (GI) has conserved gene synteny and biosynthetic module architecture across F. perrara, Enterobacteriaceae, and the Pseudovibrio strain. Comparative metabolomics analyses of F. perrara and E. coli further reveal that these two bacteria produce related colibactin pathway-dependent metabolites. Finally, we demonstrate that F. perrara, like E. coli, causes DNA damage in eukaryotic cells in vitro in a colibactin pathway-dependent manner. Together, these results support that divergent variants of the colibactin biosynthetic pathway are widely distributed among bacterial symbionts, producing related secondary metabolites and likely endowing its producer with functional capabilities important for diverse symbiotic associations.     

The HIV-1 Rev response element (RRE) adopts alternative conformations that promote different rates of virus replication

The HIV Rev protein forms a complex with a 351 nucleotide sequence present in unspliced and incompletely spliced human immunodeficiency virus (HIV) mRNAs, the Rev response element (RRE), to recruit the cellular nuclear export receptor Crm1 and Ran-GTP. This complex facilitates nucleo-cytoplasmic export of these mRNAs. The precise secondary structure of the HIV-1 RRE has been controversial, since studies have reported alternative structures comprising either four or five stem-loops. The published structures differ only in regions that lie outside of the primary Rev binding site. Using in-gel SHAPE, we have now determined that the wt NL4-3 RRE exists as a mixture of both structures. To assess functional differences between these RRE ‘conformers’, we created conformationally locked mutants by site-directed mutagenesis. Using subgenomic reporters, as well as HIV replication assays, we demonstrate that the five stem-loop form of the RRE promotes greater functional Rev/RRE activity compared to the four stem-loop counterpart.     

Effect of older age on treatment decisions and outcomes among patients with traumatic spinal cord injury

Background:

Older people are at increased risk of traumatic spinal cord injury from falls. We evaluated the impact of older age (≥ 70 yr) on treatment decisions and outcomes.

Methods:

We identified patients with traumatic spinal cord injury for whom consent and detailed data were available from among patients recruited (2004–2013) at any of the 31 acute care and rehabilitation hospitals participating in the Rick Hansen Spinal Cord Injury Registry. Patients were assessed by age group (< 70 v. ≥ 70 yr). The primary outcome was the rate of acute surgical treatment. We used bivariate and multivariate regression models to assess patient and injury-related factors associated with receiving surgical treatment and with the timing of surgery after arrival to a participating centre.

Results:

Of the 1440 patients included in our study cohort, 167 (11.6%) were 70 years or older at the time of injury. Older patients were more likely than younger patients to be injured by falling (83.1% v. 37.4%; p < 0.001), to have a cervical injury (78.0% v. 61.6%; p = 0.001), to have less severe injuries on admission (American Spinal Injury Association Impairment Scale grade C or D: 70.5% v. 46.9%; p < 0.001), to have a longer stay in an acute care hospital (median 35 v. 28 d; p < 0.005) and to have a higher in-hospital mortality (4.2% v. 0.6%; p < 0.001). Multivariate analysis did not show that age of 70 years or more at injury was associated with a decreased likelihood of surgical treatment (adjusted odds ratio [OR] 0.48, 95% confidence interval [CI] 0.22–1.07). An unplanned sensitivity analysis with different age thresholds showed that a threshold of 65 years was associated with a decreased chance of surgical treatment (OR 0.39, 95% CI 0.19–0.80). Older patients who underwent surgical treatment had a significantly longer wait time from admission to surgery than younger patients (37 v. 19 h; p < 0.001).

Interpretation:

We found chronological age to be a factor influencing treatment decisions but not at the 70-year age threshold that we had hypothesized. Older patients waited longer for surgery and had a substantially higher in-hospital mortality despite having less severe injuries than younger patients. Further research into the link between treatment delays and outcomes among older patients could inform surgical guideline development.Globally there has been an epidemiologic shift in the age of patients who sustain a traumatic spinal cord injury.^1–3 Although most people who have traumatic spinal cord injuries are 16–30 years old, there has been a progressive increase in the number who are over 70. The average age at injury has increased from 29 to 40 years.⁴ By 2032, patients over 70 are predicted to account for most patients with new traumatic spinal cord injuries.⁵ This change is attributed in part to aging baby boomers. It is unknown whether the management and outcomes of these older patients differ compared with younger patients.Older patients typically have more comorbid conditions, including cardiovascular disease, respiratory disorders, cerebrovascular disease and dementia, which are thought to increase their risk of perioperative adverse events.⁶ The use of anticoagulants for cardiac and cerebrovascular indications can delay timely surgical interventions. Older patients are also at increased risk of postoperative and medication-related adverse events, such as delirium.⁷ As a direct consequence of this perceived risk of perioperative adverse events and ambiguity about the optimal treatment for spinal cord injury in older patients, surgeons may deliberate for some time before making a clear therapeutic decision, they may choose nonoperative treatment,⁸ or they may delay the surgical treatment in an effort to optimize the patient’s condition medically.Given the increasing incidence of traumatic spinal cord injury in older adults, and the potential for differences in treatment among older and younger patients, we evaluated the impact of age on treatment decisions and outcomes among patients with traumatic spinal cord injury. We hypothesized that surgical management would differ at an age threshold of 70 years.