Bridge to the future: Important lessons from 20 years of ecosystem observations made by the OzFlux network.

In 2020, the Australian and New Zealand flux research and monitoring network, OzFlux, celebrated its 20th anniversary by reflecting on the lessons learned through two decades of ecosystem studies on global change biology. OzFlux is a network not only for ecosystem researchers, but also for those 'next users' of the knowledge, information and data that such networks provide. Here, we focus on eight lessons across topics of climate change and variability, disturbance and resilience, drought and heat stress and synergies with remote sensing and modelling. In distilling the key lessons learned, we also identify where further research is needed to fill knowledge gaps and improve the utility and relevance of the outputs from OzFlux. Extreme climate variability across Australia and New Zealand (droughts and flooding rains) provides a natural laboratory for a global understanding of ecosystems in this time of accelerating climate change. As evidence of worsening global fire risk emerges, the natural ability of these ecosystems to recover from disturbances, such as fire and cyclones, provides lessons on adaptation and resilience to disturbance. Drought and heatwaves are common occurrences across large parts of the region and can tip an ecosystem's carbon budget from a net CO2 sink to a net CO2 source. Despite such responses to stress, ecosystems at OzFlux sites show their resilience to climate variability by rapidly pivoting back to a strong carbon sink upon the return of favourable conditions. Located in under-represented areas, OzFlux data have the potential for reducing uncertainties in global remote sensing products, and these data provide several opportunities to develop new theories and improve our ecosystem models. The accumulated impacts of these lessons over the last 20 years highlights the value of long-term flux observations for natural and managed systems. A future vision for OzFlux includes ongoing and newly developed synergies with ecophysiologists, ecologists, geologists, remote sensors and modellers.

Chromosome-specific KASP markers for detecting Amblyopyrum muticum segments in wheat introgression lines.

Many wild-relative species are being used in prebreeding programs to increase the genetic diversity of wheat (Triticum aestivum L.). Genotyping tools such as single nucleotide polymorphism (SNP)-based arrays and molecular markers have been widely used to characterize wheat-wild relative introgression lines. However, due to the polyploid nature of the recipient wheat genome, it is difficult to develop SNP-based Kompetitive allele-specific polymerase chain reaction (KASP) markers that are codominant to track the introgressions from the wild species. Previous attempts to develop KASP markers have involved both exome- and polymerase chain reaction (PCR)-amplicon-based sequencing of the wild species. But chromosome-specific KASP assays have been hindered by homoeologous SNPs within the wheat genome. This study involved whole genome sequencing of the diploid wheat wild relative Amblyopyrum muticum (Boiss.) Eig and development of a de novo SNP discovery pipeline that generated ∼38,000 SNPs in unique wheat genome sequences. New assays were designed to increase the density of Am. muticum polymorphic KASP markers. With a goal of one marker per 60 Mbp, 335 new KASP assays were validated as diagnostic for Am. muticum in a wheat background. Together with assays validated in previous studies, 498 well distributed chromosome-specific markers were used to recharacterize previously genotyped wheat-Am. muticum doubled haploid (DH) introgression lines. The chromosome-specific nature of the KASP markers allowed clarification of which wheat chromosomes were involved with recombination events or substituted with Am. muticum chromosomes and the higher density of markers allowed detection of new small introgressions in these DH lines.

The FLUXNET2015 dataset and the ONEFlux processing pipeline for eddy covariance data.

The FLUXNET2015 dataset provides ecosystem-scale data on CO2, water, and energy exchange between the biosphere and the atmosphere, and other meteorological and biological measurements, from 212 sites around the globe (over 1500 site-years, up to and including year 2014). These sites, independently managed and operated, voluntarily contributed their data to create global datasets. Data were quality controlled and processed using uniform methods, to improve consistency and intercomparability across sites. The dataset is already being used in a number of applications, including ecophysiology studies, remote sensing studies, and development of ecosystem and Earth system models. FLUXNET2015 includes derived-data products, such as gap-filled time series, ecosystem respiration and photosynthetic uptake estimates, estimation of uncertainties, and metadata about the measurements, presented for the first time in this paper. In addition, 206 of these sites are for the first time distributed under a Creative Commons (CC-BY 4.0) license. This paper details this enhanced dataset and the processing methods, now made available as open-source codes, making the dataset more accessible, transparent, and reproducible.

Development of Wheat-Aegilops caudata Introgression Lines and Their Characterization Using Genome-Specific KASP Markers.

Aegilops caudata L. [syn. Ae. markgrafii (Greuter) Hammer], is a diploid wild relative of wheat (2n = 2x = 14, CC) and a valuable source for new genetic diversity for wheat improvement. It has a variety of disease resistance factors along with tolerance for various abiotic stresses and can be used for wheat improvement through the generation of genome-wide introgressions resulting in different wheat-Ae. caudata recombinant lines. Here, we report the generation of nine such wheat-Ae. caudata recombinant lines which were characterized using wheat genome-specific KASP (Kompetitive Allele Specific PCR) markers and multi-color genomic in situ hybridization (mcGISH). Of these, six lines have stable homozygous introgressions from Ae. caudata and will be used for future trait analysis. Using cytological techniques and molecular marker analysis of the recombinant lines, 182 KASP markers were physically mapped onto the seven Ae. caudata chromosomes, of which 155 were polymorphic specifically with only one wheat subgenome. Comparative analysis of the physical positions of these markers in the Ae. caudata and wheat genomes confirmed that the former had chromosomal rearrangements with respect to wheat, as previously reported. These wheat-Ae. caudata recombinant lines and KASP markers are useful resources that can be used in breeding programs worldwide for wheat improvement. Additionally, the genome-specific KASP markers could prove to be a valuable tool for the rapid detection and marker-assisted selection of other Aegilops species in a wheat background.

Rapid identification of homozygosity and site of wild relative introgressions in wheat through chromosome-specific KASP genotyping assays.

For future food security, it is important that wheat, one of the most widely consumed crops in the world, can survive the threat of abiotic and biotic stresses. New genetic variation is currently being introduced into wheat through introgressions from its wild relatives. For trait discovery, it is necessary that each introgression is homozygous and hence stable. Breeding programmes rely on efficient genotyping platforms for marker-assisted selection (MAS). Recently, single nucleotide polymorphism (SNP)-based markers have been made available on high-throughput Axiom® SNP genotyping arrays. However, these arrays are inflexible in their design and sample numbers, making their use unsuitable for long-term MAS. SNPs can potentially be converted into Kompetitive allele-specific PCR (KASP™) assays that are comparatively cost-effective and efficient for low-density genotyping of introgression lines. However, due to the polyploid nature of wheat, KASP assays for homoeologous SNPs can have difficulty in distinguishing between heterozygous and homozygous hybrid lines in a backcross population. To identify co-dominant SNPs, that can differentiate between heterozygotes and homozygotes, we PCR-amplified and sequenced genomic DNA from potential single-copy regions of the wheat genome and compared them to orthologous copies from different wild relatives. A panel of 620 chromosome-specific KASP assays have been developed that allow rapid detection of wild relative segments and provide information on their homozygosity and site of introgression in the wheat genome. A set of 90 chromosome-nonspecific assays was also produced that can be used for genotyping introgression lines. These multipurpose KASP assays represent a powerful tool for wheat breeders worldwide.

How do leaf and ecosystem measures of water-use efficiency compare?

The terrestrial carbon and water cycles are intimately linked: the carbon cycle is driven by photosynthesis, while the water balance is dominated by transpiration, and both fluxes are controlled by plant stomatal conductance. The ratio between these fluxes, the plant water-use efficiency (WUE), is a useful indicator of vegetation function. WUE can be estimated using several techniques, including leaf gas exchange, stable isotope discrimination, and eddy covariance. Here we compare global compilations of data for each of these three techniques. We show that patterns of variation in WUE across plant functional types (PFTs) are not consistent among the three datasets. Key discrepancies include the following: leaf-scale data indicate differences between needleleaf and broadleaf forests, but ecosystem-scale data do not; leaf-scale data indicate differences between C3 and C4 species, whereas at ecosystem scale there is a difference between C3 and C4 crops but not grasslands; and isotope-based estimates of WUE are higher than estimates based on gas exchange for most PFTs. Our study quantifies the uncertainty associated with different methods of measuring WUE, indicates potential for bias when using WUE measures to parameterize or validate models, and indicates key research directions needed to reconcile alternative measures of WUE.

Genetic Variation Controlling Wrinkled Seed Phenotypes in Pisum: How Lucky Was Mendel?

One of the traits studied by Mendel in pea (Pisum sativum L.) was the wrinkled-seeded phenotype, and the molecular basis for a mutation underlying this phenotype was discovered in the 1990s. Although the starch-branching enzyme gene mutation identified at the genetic locus r is most likely to be that in seeds available to Mendel in the mid-1800s, it has remained an open question as to whether or not additional natural mutations in this gene exist within Pisum germplasm collections. Here, we explore this question and show that all but two wrinkled-seeded variants in one such collection correspond to either the mutant allele described previously for the r locus or a mutation at a second genetic locus, rb, affecting the gene encoding the large subunit of Adenosine diphosphoglucose (ADP-glucose) pyrophosphorylase; the molecular basis for the rb mutation is described here. The genetic basis for the phenotype of one (JI 2110) of the two lines which are neither r nor rb has been studied in crosses with a round-seeded variant (JI 281); for which extensive genetic marker data were expected. In marked contrast to the trait studied by Mendel and the rb phenotype; the data suggest that the wrinkled-seeded phenotype in JI 2110 is maternally determined, controlled by two genetic loci, and the extent to which it is manifested is very sensitive to the environment. Metabolite analysis of the cotyledons of JI 2110 revealed a profile for sucrose and sucrose-derived compounds that was more similar to that of wild-type round-seeded, than that of wrinkled-seeded r, pea lines. However, the metabolite profile of the seed coat (testa) of JI 2110 was distinct from that of other round-seeded genotypes tested which, together with analysis of recombinant inbred progeny lines, suggests an explanation for the seed phenotype.

Eliminating anti-nutritional plant food proteins: the case of seed protease inhibitors in pea.

Several classes of seed proteins limit the utilisation of plant proteins in human and farm animal diets, while plant foods have much to offer to the sustainable intensification of food/feed production and to human health. Reduction or removal of these proteins could greatly enhance seed protein quality and various strategies have been used to try to achieve this with limited success. We investigated whether seed protease inhibitor mutations could be exploited to enhance seed quality, availing of induced mutant and natural Pisum germplasm collections to identify mutants, whilst acquiring an understanding of the impact of mutations on activity. A mutant (TILLING) resource developed in Pisum sativum L. (pea) and a large germplasm collection representing Pisum diversity were investigated as sources of mutations that reduce or abolish the activity of the major protease inhibitor (Bowman-Birk) class of seed protein. Of three missense mutations, predicted to affect activity of the mature trypsin / chymotrypsin inhibitor TI1 protein, a C77Y substitution in the mature mutant inhibitor abolished inhibitor activity, consistent with an absolute requirement for the disulphide bond C77-C92 for function in the native inhibitor. Two further classes of mutation (S85F, E109K) resulted in less dramatic changes to isoform or overall inhibitory activity. The alternative strategy to reduce anti-nutrients, by targeted screening of Pisum germplasm, successfully identified a single accession (Pisum elatius) as a double null mutant for the two closely linked genes encoding the TI1 and TI2 seed protease inhibitors. The P. elatius mutant has extremely low seed protease inhibitory activity and introgression of the mutation into cultivated germplasm has been achieved. The study provides new insights into structure-function relationships for protease inhibitors which impact on pea seed quality. The induced and natural germplasm variants identified provide immediate potential for either halving or abolishing the corresponding inhibitory activity, along with associated molecular markers for breeding programmes. The potential for making large changes to plant protein profiles for improved and sustainable food production through diversity is illustrated. The strategy employed here to reduce anti-nutritional proteins in seeds may be extended to allergens and other seed proteins with negative nutritional effects. Additionally, the novel variants described for pea will assist future studies of the biological role and health-related properties of so-called anti-nutrients.

Allele frequencies in the VRN-A1, VRN-B1 and VRN-D1 vernalization response and PPD-B1 and PPD-D1 photoperiod sensitivity genes, and their effects on heading in a diverse set of wheat cultivars (Triticum aestivum L.).

Heading of cereals is determined by complex genetic and environmental factors in which genes responsible for vernalization and photoperiod sensitivity play a decisive role. Our aim was to use diagnostic molecular markers to determine the main allele types in VRN-A1, VRN-B1, VRN-D1, PPD-B1 and PPD-D1 in a worldwide wheat collection of 683 genotypes and to investigate the effect of these alleles on heading in the field. The dominant VRN-A1, VRN-B1 and VRN-D1 alleles were present at a low frequency. The PPD-D1a photoperiod-insensitive allele was carried by 57 % of the cultivars and was most frequent in Asian and European cultivars. The PPD-B1 photoperiod-insensitive allele was carried by 22 % of the genotypes from Asia, America and Europe. Nine versions of the PPD-B1-insensitive allele were identified based on gene copy number and intercopy structure. The allele compositions in PPD-D1, PPD-B1 and VRN-D1 significantly influenced heading and together explained 37.5 % of the phenotypic variance. The role of gene model increased to 39.1 % when PPD-B1 intercopy structure was taken into account instead of overall PPD-B1 type (sensitive vs. insensitive). As a single component, PPD-D1 had the most important role (28.0 % of the phenotypic variance), followed by PPD-B1 (12.3 % for PPD-B1_overall, and 15.1 % for PPD-B1_intercopy) and VRN-D1 (2.2 %). Significant gene interactions were identified between the marker alleles within PPD-B1 and between VRN-D1 and the two PPD1 genes. The earliest heading genotypes were those with the photoperiod-insensitive allele in PPD-D1 and PPD-B1, and with the spring allele for VRN-D1 and the winter alleles for VRN-A1 and VRN-B1. This combination could only be detected in genotypes from Southern Europe and Asia. Late-heading genotypes had the sensitivity alleles for both PPD1 genes, regardless of the allelic composition of the VRN1 genes. There was a 10-day difference in heading between the earliest and latest groups under field conditions.

Genetic distinctiveness of the Herdwick sheep breed and two other locally adapted hill breeds of the UK.

There is considerable interest in locally adapted breeds of livestock as reservoirs of genetic diversity that may provide important fitness traits for future use in agriculture. In marginal areas, these animals contribute to food security and extract value from land unsuitable for other systems of farming. In England, close to 50% of the national sheep flock is farmed on grassland designated as disadvantaged areas for agricultural production. Many of these areas are in the uplands, where some native breeds of sheep continue to be commercially farmed only in highly localised geographical regions to which they are adapted. This study focuses on three of these breeds, selected for their adaptation to near identical environments and their geographical concentration in regions close to one another. Our objective has been to use retrotyping, microsatellites and single nucleotide polymorphisms to explore the origins of the breeds and whether, despite their similar adaptations and proximity, they are genetically distinctive. We find the three breeds each have a surprisingly different pattern of retrovirus insertions into their genomes compared with one another and with other UK breeds. Uniquely, there is a high incidence of the R0 retrotype in the Herdwick population, characteristic of a primitive genome found previously in very few breeds worldwide and none in the UK mainland. The Herdwick and Rough Fells carry two rare retroviral insertion events, common only in Texels, suggesting sheep populations in the northern uplands have a historical association with the original pin-tail sheep of Texel Island. Microsatellite data and analyses of SNPs associated with RXFP2 (horn traits) and PRLR (reproductive performance traits) also distinguished the three breeds. Significantly, an SNP linked to TMEM154, a locus controlling susceptibility to infection by Maedi-Visna, indicated that all three native hill breeds have a lower than average risk of infection to the lentivirus.

Exploiting a fast neutron mutant genetic resource in Pisum sativum (pea) for functional genomics.

A fast neutron (FN)-mutagenised population was generated in Pisum sativum L. (pea) to enable the identification and isolation of genes underlying traits and processes. Studies of several phenotypic traits have clearly demonstrated the utility of the resource by associating gene deletions with phenotype followed by functional tests exploiting additional mutant sources, from both induced and natural variant germplasm. For forward genetic screens, next generation sequencing methodologies provide an opportunity for identifying genes associated with deletions rapidly and systematically. The application of rapid reverse genetic screens of the fast neutron mutant pea population supports conclusions on the frequency of deletions based on phenotype alone. These studies also suggest that large deletions affecting one or more loci can be non-deleterious to the pea genome, yielding mutants that could not be obtained by other means. Deletion mutants affecting genes associated with seed metabolism and storage are providing unique opportunities to identify the products of complex and related gene families, and to study the downstream consequences of such deletions.

Copy number variation affecting the Photoperiod-B1 and Vernalization-A1 genes is associated with altered flowering time in wheat (Triticum aestivum).

The timing of flowering during the year is an important adaptive character affecting reproductive success in plants and is critical to crop yield. Flowering time has been extensively manipulated in crops such as wheat (Triticum aestivum L.) during domestication, and this enables them to grow productively in a wide range of environments. Several major genes controlling flowering time have been identified in wheat with mutant alleles having sequence changes such as insertions, deletions or point mutations. We investigated genetic variants in commercial varieties of wheat that regulate flowering by altering photoperiod response (Ppd-B1 alleles) or vernalization requirement (Vrn-A1 alleles) and for which no candidate mutation was found within the gene sequence. Genetic and genomic approaches showed that in both cases alleles conferring altered flowering time had an increased copy number of the gene and altered gene expression. Alleles with an increased copy number of Ppd-B1 confer an early flowering day neutral phenotype and have arisen independently at least twice. Plants with an increased copy number of Vrn-A1 have an increased requirement for vernalization so that longer periods of cold are required to potentiate flowering. The results suggest that copy number variation (CNV) plays a significant role in wheat adaptation.

The genetic map of Artemisia annua L. identifies loci affecting yield of the antimalarial drug artemisinin.

Artemisinin is a plant natural product produced by Artemisia annua and the active ingredient in the most effective treatment for malaria. Efforts to eradicate malaria are increasing demand for an affordable, high-quality, robust supply of artemisinin. We performed deep sequencing on the transcriptome of A. annua to identify genes and markers for fast-track breeding. Extensive genetic variation enabled us to build a detailed genetic map with nine linkage groups. Replicated field trials resulted in a quantitative trait loci (QTL) map that accounts for a significant amount of the variation in key traits controlling artemisinin yield. Enrichment for positive QTLs in parents of new high-yielding hybrids confirms that the knowledge and tools to convert A. annua into a robust crop are now available.

Tendril-less regulates tendril formation in pea leaves.

Tendrils are contact-sensitive, filamentous organs that permit climbing plants to tether to their taller neighbors. Tendrilled legume species are grown as field crops, where the tendrils contribute to the physical support of the crop prior to harvest. The homeotic tendril-less (tl) mutation in garden pea (Pisum sativum), identified almost a century ago, transforms tendrils into leaflets. In this study, we used a systematic marker screen of fast neutron-generated tl deletion mutants to identify Tl as a Class I homeodomain leucine zipper (HDZIP) transcription factor. We confirmed the tendril-less phenotype as loss of function by targeting induced local lesions in genomes (TILLING) in garden pea and by analysis of the tendril-less phenotype of the t mutant in sweet pea (Lathyrus odoratus). The conversion of tendrils into leaflets in both mutants demonstrates that the pea tendril is a modified leaflet, inhibited from completing laminar development by Tl. We provide evidence to show that lamina inhibition requires Unifoliata/LEAFY-mediated Tl expression in organs emerging in the distal region of the leaf primordium. Phylogenetic analyses show that Tl is an unusual Class I HDZIP protein and that tendrils evolved either once or twice in Papilionoid legumes. We suggest that tendrils arose in the Fabeae clade of Papilionoid legumes through acquisition of the Tl gene.

A high-density microsatellite consensus map for bread wheat (Triticum aestivum L.).

A microsatellite consensus map was constructed by joining four independent genetic maps of bread wheat. Three of the maps were F(1)-derived, doubled-haploid line populations and the fourth population was 'Synthetic' x 'Opata', an F(6)-derived, recombinant-inbred line population. Microsatellite markers from different research groups including the Wheat Microsatellite Consortium, GWM, GDM, CFA, CFD, and BARC were used in the mapping. A sufficient number of common loci between genetic maps, ranging from 52 to 232 loci, were mapped on different populations to facilitate joining the maps. Four genetic maps were developed using MapMaker V3.0 and JoinMap V3.0. The software CMap, a comparative map viewer, was used to align the four maps and identify potential errors based on consensus. JoinMap V3.0 was used to calculate marker order and recombination distances based on the consensus of the four maps. A total of 1,235 microsatellite loci were mapped, covering 2,569 cM, giving an average interval distance of 2.2 cM. This consensus map represents the highest-density public microsatellite map of wheat and is accompanied by an allele database showing the parent allele sizes for every marker mapped. This enables users to predict allele sizes in new breeding populations and develop molecular breeding and genomics strategies.

A microarray-based high throughput molecular marker genotyping method: the tagged microarray marker (TAM) approach.

A microarray-based method has been developed for scoring thousands of DNAs for a co-dominant molecular marker on a glass slide. The approach was developed to detect insertional polymorphism of transposons and works well with single nucleotide polymorphism (SNP) markers. Biotin- terminated allele-specific PCR products are spotted unpurified onto streptavidin-coated glass slides and visualised by hybridisation of fluorescent detector oligonucleotides to tags attached to the allele- specific PCR primers. Two tagged primer oligonucleotides are used per locus and each tag is detected by hybridisation to a concatameric DNA probe labelled with multiple fluorochromes.

Toward positional cloning of Vgt1, a QTL controlling the transition from the vegetative to the reproductive phase in maize.

Vgt1 (Vegetative to generative transition 1) is a quantitative trait locus (QTL) for flowering time in maize (Zea mays L.). Vgt1 was initially mapped in a ca. 5-cM interval on chromosome bin 8.05, using a set of near-isogenic lines (NILs) in the genetic background of the late dent line N28, with the earliness allele introgressed from the early variety Gaspé Flint. A new large mapping population was produced by crossing N28 and one early NIL with a ca. 6-cM long Gaspé Flint introgression at the Vgt1 region. Using PCR-based assays at markers flanking Vgt1, 69 segmental NILs homozygous for independent crossovers near the QTL were developed. When the NILs were tested in replicated field trials for days to pollen shed (DPS) and plant node number (ND), the QTL followed a Mendelian segregation. Using bulk segregant analysis and AFLP profiling, 17 AFLP markers linked to the QTL region were identified. Statistical analysis indicated a substantial coincidence of the effects of Vgt1 on both DPS and ND. Vgt1 was mapped at ca. 0.3 cM from an AFLP marker. As compared to DPS, the higher heritability of ND allowed for a more accurate assessment of the effects of Vgt1. The feasibility of the positional cloning of Vgt1 is discussed.