Evolution and selection of Rhg1, a copy‐number variant nematode‐resistance locus

The soybean cyst nematode (SCN) resistance locus Rhg1 is a tandem repeat of a 31.2 kb unit of the soybean genome. Each 31.2‐kb unit contains four genes. One allele of Rhg1, Rhg1‐b, is responsible for protecting most US soybean production from SCN. Whole‐genome sequencing was performed, and PCR assays were developed to investigate allelic variation in sequence and copy number of the Rhg1 locus across a population of soybean germplasm accessions. Four distinct sequences of the 31.2‐kb repeat unit were identified, and some Rhg1 alleles carry up to three different types of repeat unit. The total number of copies of the repeat varies from 1 to 10 per haploid genome. Both copy number and sequence of the repeat correlate with the resistance phenotype, and the Rhg1 locus shows strong signatures of selection. Significant linkage disequilibrium in the genome outside the boundaries of the repeat allowed the Rhg1 genotype to be inferred using high‐density single nucleotide polymorphism genotyping of 15 996 accessions. Over 860 germplasm accessions were found likely to possess Rhg1 alleles. The regions surrounding the repeat show indications of non‐neutral evolution and high genetic variability in populations from different geographic locations, but without evidence of fixation of the resistant genotype. A compelling explanation of these results is that balancing selection is in operation at Rhg1.     

The soybean (Glycine max) nodulation‐suppressive CLE peptide,GmRIC1, functions interspecifically in common white bean (Phaseolus vulgaris), but not in a supernodulating line mutated in the receptor PvNARK

Legume plants regulate the number of nitrogen‐fixing root nodules they form via a process called the Autoregulation of Nodulation (AON). Despite being one of the most economically important and abundantly consumed legumes, little is known about the AON pathway of common bean (Phaseolus vulgaris). We used comparative‐ and functional‐genomic approaches to identify central components in the AON pathway of common bean. This includes identifying PvNARK, which encodes a LRR receptor kinase that acts to regulate root nodule numbers. A novel, truncated version of the gene was identified directly upstream of PvNARK, similar to Medicago truncatula, but not seen in Lotus japonicus or soybean. Two mutant alleles of PvNARK were identified that cause a classic shoot‐controlled and nitrate‐tolerant supernodulation phenotype. Homeologous over‐expression of the nodulation‐suppressive CLE peptide‐encoding soybean gene, GmRIC1, abolished nodulation in wild‐type bean, but had no discernible effect on PvNARK‐mutant plants. This demonstrates that soybean GmRIC1 can function interspecifically in bean, acting in a PvNARK‐dependent manner. Identification of bean PvRIC1, PvRIC2 and PvNIC1, orthologues of the soybean nodulation‐suppressive CLE peptides, revealed a high degree of conservation, particularly in the CLE domain. Overall, our work identified four new components of bean nodulation control and a truncated copy of PvNARK, discovered the mutation responsible for two supernodulating bean mutants and demonstrated that soybean GmRIC1 can function in the AON pathway of bean.     

Genetic basis of tristyly in tetraploid Oxalis alpina (Oxalidaceae)

The inheritance of style‐morphs was investigated in tetraploid populations of tristylous Oxalis alpina (Oxalidaceae) to determine if alleles controlling style‐morphs are expressed at duplicated loci. In tetraploid populations, a dominant S allele leads to expression of the short‐styled phenotype at the short/non‐short locus and is epistatic to the M allele at the mid/long locus. The M allele results in expression of the mid‐styled phenotype but only if the S allele is absent. Long‐styled morphs are homozygous recessive at the short and mid loci. Test crosses of many tetraploid short‐styled individuals resulted in segregations of short‐, mid‐ and long‐styled individuals which, because of linkage between the short and mid loci, can only occur with polyploidy and expression of alleles at duplicated loci. Segregation patterns from three crosses suggest the possibility of disomic inheritance via preferential pairing of chromosomes in tetraploid populations of O. alpina. Segregation patterns in the progeny of mid‐styled individuals indicated that only a few individuals had more than one copy of the M allele, despite the potential for accumulation of M alleles via self‐fertilization of partially self‐compatible mid‐styled morphs in some populations. © 2015 The Linnean Society of London, Botanical Journal of the Linnean Society, 2015, 179 , 308–318.     

Soybean cyst nematode resistance in soybean is independent of the <Emphasis Type="Italic">Rhg4</Emphasis> locus <Emphasis Type="Italic">LRR-RLK</Emphasis> gene

To test the function of candidate genes in soybean for resistance to the soybean cyst nematode (SCN), a large collection of EMS-mutants from the SCN-resistant soybean cultivar “Forrest” was developed for Targeting Induced Local Lesions IN Genomes (TILLING). Additionally, due to the complexity of the soybean genome, an integrated set of genomic and genetic analysis tools was employed to complement the TILLING approach. The efficiency of this integrated set of tools was tested using a candidate soybean gene for resistance to SCN, encoding a leucine-rich repeat receptor-like kinase (LRR-RLK) that was identified by map-based cloning at the Rhg4 locus. The Rhg4 locus is one of the major quantitative trait loci controlling soybean resistance against SCN race 3 (HG type 0) in cv. Forrest, but the gene(s) sequence for resistance remains to be determined. Using TILLING, a Forrest mutant containing a nonsense mutation in the LRR domain of the candidate resistance protein was identified and confirmed; however, the SCN-resistant phenotype of the mutant was not altered. Haplotyping and EcoTILLING of recombinant inbred lines along with complementation analysis corroborated the TILLING result and ruled out the possibility of functional redundancy by a second copy of the LRR-RLK gene identified in the soybean genome. This study validates the use of TILLING, in combination with an integrated set of genomic tools, as an efficient means of testing candidate genes for SCN resistance in soybean.     

An efficient method for measuring copy number variation applied to improvement of nematode resistance in soybean

Copy number variation (CNV) is implicated in important traits in multiple crop plants, but can be challenging to genotype using conventional methods. The Rhg1 locus of soybean, which confers resistance to soybean cyst nematode (SCN), is a CNV of multiple 31.2‐kb genomic units each containing four genes. Reliable, high‐throughput methods to quantify Rhg1 and other CNVs for selective breeding were developed. The CNV genotyping assay described here uses a homeologous gene copy within the paleopolyploid soybean genome to provide the internal control for a single‐tube TaqMan copy number assay. Using this assay, CNV in breeding populations can be tracked with high precision. We also show that extensive CNV exists within Fayette, a released, inbred SCN‐resistant soybean cultivar with a high copy number at Rhg1 derived from a single donor parent. Copy number at Rhg1 is therefore unstable within a released variety over a relatively small number of generations. Using this assay to select for individuals with altered copy number, plants were obtained with both increased copy number and increased SCN resistance relative to control plants. Thus, CNV genotyping technologies can be used as a new type of marker‐assisted selection to select for desirable traits in breeding populations, and to control for undesirable variation within cultivars.     

Allopolyploidy in the Wintergreen Group of tribe Gaultherieae (Ericaceae) inferred from low‐copy nuclear genes

DNA sequence data from the low‐copy nuclear genes waxy (GBSSI) and leafy were compared with plastid and ITS sequence data from prior studies to reconstruct phylogenetic relationships in the Wintergreen Group of tribe Gaultherieae (Ericaceae). We conducted phylogenetic analysis with 102 species that includes representatives of all 15 major clades previously identified within the Wintergreen Group and that together span its circum‐Pacific distribution. Results yielded two distinct homeologous copies of waxy for two of the clades, each in widely separated parts of the tree. It also yielded two copies of leafy for one of the clades; only one copy of leafy was found for the other clade, but it was placed in the same major clade as its waxy counterpart and well away from its placement in a prior plastid analysis. A combined four‐locus (waxy, leafy, ITS and plastid data) phylogenetic analysis of all available relevant data placed the copies of each of the clades in two distinct positions in the phylogeny with strong overall statistical support. In combination with evidence from morphology, reproductive biology and cytology, the results suggest that these clades arose through allopolyploid hybridization between lineages deep in the phylogeny but relatively close geographically. This finding confirms previous assumptions that hybridization has played an important role in the evolution of the Gaultherieae.     

Structural and functional analysis of nitrogenase genes from the broad-host-range Rhizobium strain ANU240

The genes encoding the structural components of nitrogenase, nifH, nifD and nifK, from the fast-growing, broad-host-range Rhizobium strain ANU240 have been identified and characterized. They are duplicated and linked in an operon nifHDK in both copies. Sequence analysis of the nifH gene from each copy, together with partial sequence analysis of the nifD and nifK genes, and restriction endonuclease analysis suggested that the duplication is precise. Comparison of the Fe-protein sequence from strain ANU240 with that from other nitrogen-fixing organisms revealed that, despite its broad host range and certain physiological properties characteristic of Bradyrhizobium strains, ANU240 is more closely related to the narrow-host-range Rhizobium strains than to the broad-host-range Bradyrhizobium strains. The promoter regions of both copies of the nif genes contain the consensus sequence characteristic of nif promoters, and functional analysis of the two promoters suggested that both nif operons are transcribed in nodules.     

The evolutionary fate of heterogeneous gene duplications: A precarious overdominant equilibrium between environment,sublethality and complementation

Gene duplications occur at a high rate. Although most appear detrimental, some homogeneous duplications (identical gene copies) can be selected for beneficial increase in produced proteins. Heterogeneous duplications, which combine divergent alleles of a single locus, are seldom studied due to the paucity of empirical data. We investigated their role in an ongoing adaptive process at the ace‐1 locus in Culex pipiens mosquitoes. We assessed the worldwide diversity of the ace‐1 alleles (single‐copy, susceptible S and insecticide‐resistant R, and duplicated D that pair one S and one R copy), analysed their phylogeography and measured their fitness to understand their early dynamics using population genetics models. It provides a coherent and comprehensive evolutionary scenario. We show that D alleles are present in most resistant populations and display a higher diversity than R alleles (27 vs. 4). Most appear to result from independent unequal crossing‐overs between local single‐copy alleles, suggesting a recurrent process. Most duplicated alleles have a limited geographic distribution, probably resulting from their homozygous sublethality (HS phenotype). In addition, heterozygotes carrying different HS D alleles showed complementation, indicating different recessive lethal mutations. Due to mosaic insecticide control practices, balancing selection (overdominance) plays a key role in the early dynamics of heterogeneous duplicated alleles; it also favours a high local polymorphism of HS D alleles in natural populations (overdominance reinforced by complementation). Overall, our study shows that the evolutionary fate of heterogeneous duplications (and their long‐term role) depends on finely balanced selective pressures due to the environment and to their genomic structure.     

Molecular Evidence for Asymmetric Evolution of Sister Duplicated Blocks after Cereal Polyploidy

Polyploidy (genome duplication) is thought to have contributed to the evolution of the eukaryotic genome, but complex genome structures and massive gene loss during evolution has complicated detection of these ancestral duplication events. The major factors determining the fate of duplicated genes are currently unclear, as are the processes by which duplicated genes evolve after polyploidy. Fine-scale analysis between homologous regions may allow us to better understand post-polyploidy evolution. Here, using gene-by-gene and gene-by-genome strategies, we identified the S5 region and four homologous regions within the japonica genome. Additional phylogenomic analyses of the comparable duplicated blocks indicate that four successive duplication events gave rise to these five regions, allowing us to propose a model for this local chromosomal evolution. According to this model, gene loss may play a major role in post-duplication genetic evolution at the segmental level. Moreover, we found molecular evidence that one of the sister duplicated blocks experienced more gene loss and a more rapid evolution subsequent to two recent duplication events. Given that these two recent duplication events were likely involved in polyploidy, this asymmetric evolution (gene loss and gene divergence) may be one possible mechanism accounting for the diploidization at the segmental level. Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1007/s11103-005-4414-1     

Allopolyploidization in Cucumis contributes to delayed leaf maturation with repression of redundant homoeologous genes

The important role of polyploidy in plant evolution is widely recognized. However, many questions remain to be explored to address how polyploidy affects the phenotype of the plant. To shed light on the phenotypic and molecular impacts of allopolyploidy, we investigated the leaf development of a synthesized allotetraploid (Cucumis × hytivus), with an emphasis on chlorophyll development. Delayed leaf maturation was identified in C. × hytivus, based on delayed leaf expansion, initial chlorophyll deficiency in the leaves and disordered sink‐source transition. Anatomical observations also revealed disturbed chloroplast development in C. ×hytivus. The determination of chlorophyll biosynthesis intermediates suggested that the chlorophyll biosynthesis pathway of C. × hytivus is blocked at the site at which uroporphyrinogen III is catalysed to coproporphyrinogen III. Three chlorophyll biosynthesis‐related genes, HEMA1, HEME2 and POR, were significantly repressed in C. × hytivus. Sequence alignment showed both synonymous and non‐synonymous substitutions in the HEMA1, HEME2 and POR genes of the parents. Cloning of the chlorophyll biosynthetic genes suggested the retention of homoeologs. In addition, a chimeric clone of the HEMA1 gene that consisted of homologous genes from the parents was identified in C. × hytivus. Overall, our results showed that allopolyploidization in Cucumis has resulted in disturbed chloroplast development and reduced chlorophyll biosynthesis caused by the repressed expression of duplicated homologous genes, which further led to delayed leaf maturation in the allotetraploid, C. × hytivus. The preferential retention/loss of certain types of genes and non‐reciprocal homoeologous recombination were also supported in the present study, which provides new insights into the impact of allopolyploidy.     

The Na(+) cycle in Acetobacterium woodii: identification and characterization of a Na(+) translocating F(1)F(0)-ATPase with a mixed oligomer of 8 and 16 kDa proteolipids

The homoacetogenic bacterium Acetobacterium woodii relies on a sodium ion current across its cytoplasmic membrane for energy-dependent reactions. The sodium ion potential is established by a yet to be identified primary, electrogenic pump connected to the Wood-Ljungdahl pathway. Reactions possibly involved in Na(+) export are discussed. The electrochemical sodium ion potential generated is used to drive endergonic reactions such as flagellar rotation and ATP synthesis. Biochemical and molecular data identified the Na(+)-ATPase of A. woodii as a typical member of the F(1)F(0) class of ATPases. Its catalytic properties and the hypothetical sodium ion binding site in subunit c are discussed. The encoding genes were cloned and, surprisingly, the atp operon was shown to contain multiple copies of genes encoding subunit c. Two copies encode identical 8 kDa proteolipids, and a third copy arose by duplication and subsequent fusion of two genes. Furthermore, the duplicated subunit c does not contain the ion binding site in hair pin two. Biochemical and molecular data revealed that all three copies of subunit c constitute a mixed oligomer. The evolution of the structure and function of subunit c in ATPases from eucarya, bacteria, and archaea is discussed.     

Parallel origins of duplications and the formation of pseudogenes in mitochondrial DNA from parthenogenetic lizards (Heteronotia binoei; Gekkonidae)

Summary Analysis of mitochondrial DNAs (mtDNAs) from parthenogenetic lizards of theHeteronotia binoei complex with restriction enzymes revealed an 5-kb addition present in all 77 individuals. Cleavage site mapping suggested the presence of a direct tandem duplication spanning the 16S and 12S rRNA genes, the control region and most, if not all, of the gene for the subunit 1 of NADH dehydrogenase (ND1). The location of the duplication was confirmed by Southern hybridization. A restriction enzyme survey provided evidence for modifications to each copy of the duplicated sequence, including four large deletions. Each gene affected by a deletion was complemented by an intact version in the other copy of the sequence, although for one gene the functional copy was heteroplasmic for another deletion. Sequencing of a fragment from one copy of the duplication which encompassed the tRNA^leu(UUR) and parts of the 16S rRNA and ND1 genes, revealed mutations expected to disrupt function. Thus, evolution subsequent to the duplication event has resulted in mitochondrial pseudogenes. The presence of duplications in all of these parthenogens, but not among representatives of their maternal sexual ancestors, suggests that the duplications arose in the parthenogenetic form. This provides the second instance inH. binoei of mtDNA duplication associated with the transition from sexual to parthenogenetic reproduction. The increased incidence of duplications in parthenogenetic lizards may be caused by errors in mtDNA replication due to either polyploidy or hybridity of their nuclear genomes.     

SPO11.2 is essential for programmed double‐strand break formation during meiosis in bread wheat (Triticum aestivum L.)

Meiotic recombination is initiated by formation of DNA double‐strand breaks (DSBs). This involves a protein complex that includes in plants the two similar proteins, SPO11‐1 and SPO11‐2. We analysed the sequences of SPO11‐2 in hexaploid bread wheat (Triticum aestivum), as well as in its diploid and tetraploid progenitors. We investigated its role during meiosis using single, double and triple mutants. The three homoeologous SPO11‐2 copies of hexaploid wheat exhibit high nucleotide and amino acid similarities with those of the diploids, tetraploids and Arabidopsis. Interestingly, however, two nucleotides deleted in exon‐2 of the A copy lead to a premature stop codon and suggest that it encodes a non‐functional protein. Remarkably, the mutation was absent from the diploid A‐relative Triticum urartu, but present in the tetraploid Triticum dicoccoides and in different wheat cultivars indicating that the mutation occurred after the first polyploidy event and has since been conserved. We further show that triple mutants with all three copies (A, B, D) inactivated are sterile. Cytological analyses of these mutants show synapsis defects, accompanied by severe reductions in bivalent formation and numbers of DMC1 foci, thus confirming the essential role of TaSPO11‐2 in meiotic recombination in wheat. In accordance with its 2‐nucleotide deletion in exon‐2, double mutants for which only the A copy remained are also sterile. Notwithstanding, some DMC1 foci remain visible in this mutant, suggesting a residual activity of the A copy, albeit not sufficient to restore fertility.     

Regulation of cyclic peptide biosynthesis and pathogenicity in Cochliobolus carbonum by TOXEp, a novel protein with a bZIP basic DNA-binding motif and four ankyrin repeats

HC-toxin is an epoxide-containing cyclic tetrapeptide that is a critical virulence determinant in the pathogenic interaction between the filamentous fungus Cochliobolus carbonum and maize. HC-toxin exerts a potent cytostatic effect on plant and animal cells by inhibiting histone deacetylase. The biosynthesis of HC-toxin by C. carbonum is controlled by a complex genetic locus, TOX2, that contains multiple, duplicated copies of genes encoding export and biosynthetic enzymes. A new gene in the TOX2 complex, TOXE, has now been isolated. Mutation of TOXE by targeted gene disruption has no effect on growth and sporulation but abolishes HC-toxin production and pathogenicity. TOXE is required for the expression of three genes with a known or putative role in HC-toxin production, but is not required for expression of HTS1, which encodes the large, multifunctional peptide synthetase that is the central enzyme in HC-toxin biosynthesis. At its N-terminus, TOXEp has a bZIP basic DNA binding domain, but it does not contain any discernible leucine zipper or helix-loop-helix. At its carboxy terminus, TOXEp contains four ankyrin repeats. In having these two common regulatory motifs in a single polypeptide, TOXEp appears to represent a novel class of regulatory protein. TOXE is present only in HC-toxin-producing (Tox2⁺) isolates of C. carbonum. Most Tox2⁺ isolates have two copies; in strain SB111, one copy of TOXE is on the same 3.5-Mb chromosome that contains all of the other genes known to be involved in HC-toxin biosynthesis, and the second copy of TOXE is on a 0.7-Mb chromosome. Received: 20 April 1998 / Accepted: 21 September 1998     

Diploidization and genome size change in allopolyploids is associated with differential dynamics of low‐ and high‐copy sequences

Recent advances have highlighted the ubiquity of whole‐genome duplication (polyploidy) in angiosperms, although subsequent genome size change and diploidization (returning to a diploid‐like condition) are poorly understood. An excellent system to assess these processes is provided by Nicotiana section Repandae, which arose via allopolyploidy (approximately 5 million years ago) involving relatives of Nicotiana sylvestris and Nicotiana obtusifolia. Subsequent speciation in Repandae has resulted in allotetraploids with divergent genome sizes, including Nicotiana repanda and Nicotiana nudicaulis studied here, which have an estimated 23.6% genome expansion and 19.2% genome contraction from the early polyploid, respectively. Graph‐based clustering of next‐generation sequence data enabled assessment of the global genome composition of these allotetraploids and their diploid progenitors. Unexpectedly, in both allotetraploids, over 85% of sequence clusters (repetitive DNA families) had a lower abundance than predicted from their diploid relatives; a trend seen particularly in low‐copy repeats. The loss of high‐copy sequences predominantly accounts for the genome downsizing in N. nudicaulis. In contrast, N. repanda shows expansion of clusters already inherited in high copy number (mostly chromovirus‐like Ty3/Gypsy retroelements and some low‐complexity sequences), leading to much of the genome upsizing predicted. We suggest that the differential dynamics of low‐ and high‐copy sequences reveal two genomic processes that occur subsequent to allopolyploidy. The loss of low‐copy sequences, common to both allopolyploids, may reflect genome diploidization, a process that also involves loss of duplicate copies of genes and upstream regulators. In contrast, genome size divergence between allopolyploids is manifested through differential accumulation and/or deletion of high‐copy‐number sequences.     

Use of Illumina sequencing to identify transposon insertions underlying mutant phenotypes in high‐copy Mutator lines of maize

High‐copy transposons have been effectively exploited as mutagens in a variety of organisms. However, their utility for phenotype‐driven forward genetics has been hampered by the difficulty of identifying the specific insertions responsible for phenotypes of interest. We describe a new method that can substantially increase the throughput of linking a disrupted gene to a known phenotype in high‐copy Mutator (Mu) transposon lines in maize. The approach uses the Illumina platform to obtain sequences flanking Mu elements in pooled, bar‐coded DNA samples. Insertion sites are compared among individuals of suitable genotype to identify those that are linked to the mutation of interest. DNA is prepared for sequencing by mechanical shearing, adapter ligation, and selection of DNA fragments harboring Mu flanking sequences by hybridization to a biotinylated oligonucleotide corresponding to the Mu terminal inverted repeat. This method yields dense clusters of sequence reads that tile approximately 400 bp flanking each side of each heritable insertion. The utility of the approach is demonstrated by identifying the causal insertions in four genes whose disruption blocks chloroplast biogenesis at various steps: thylakoid protein targeting (cpSecE), chloroplast gene expression (polynucleotide phosphorylase and PTAC12), and prosthetic group attachment (HCF208/CCB2). This method adds to the tools available for phenotype‐driven Mu tagging in maize, and could be adapted for use with other high‐copy transposons. A by‐product of the approach is the identification of numerous heritable insertions that are unrelated to the targeted phenotype, which can contribute to community insertion resources.