Sequence variation and haplotype structure at the putative flowering-time locus COL1 of Brassica nigra

Motivated by a previous study indicating that polymorphism at an indel, Ind2, within the Brassica nigra COL1 gene is significantly associated with flowering time, we searched for evidence of selection in a sample of 41 complete sequences of B. nigra COL1. The within-gene population recombination rate is moderate, and all neutrality tests used in the present study failed to detect departure from the standard neutral model or evidence of selection. The haplotype structure of the 5'-half of the gene is primarily associated with the demographic history of the species and more specifically with the split between European and Ethiopian populations, whereas the structure of the 3'-half reflects the polymorphism at Ind2. This could be the result of selection or a combination of recombination and migration during the history of the sample of sequences. Without additional information on polymorphism in flanking areas, these two alternatives are difficult to tell apart. If selection acted on the gene, we suggest that if the indel itself is not the target of selection, among the polymorphic sites cosegregating with the polymorphism at Ind2, replacement polymorphisms around sites 890 and 1260 are the most likely quantitative trait nucleotides within the gene.     

Naturally occurring indel variation in the Brassica nigra COL1 gene is associated with variation in flowering time

Previous QTL mapping identified a Brassica nigra homolog to Arabidopsis thaliana CO as a candidate gene affecting flowering time in B. nigra. Transformation of an A. thaliana co mutant with two different alleles of the B. nigra CO (Bni COa) homolog, one from an early-flowering B. nigra plant and one from a late one, did not show any differential effect of the two alleles on flowering time. The DNA sequence of the coding region of the two alleles was also identical, showing that nucleotide variation influencing flowering time must reside outside the coding region of Bni COa. In contrast, the nucleotide sequence of the B. nigra COL1 (Bni COL1) gene located 3.5 kb upstream of Bni COa was highly diverged between the alleles from early and late plants. One indel polymorphism in the Bni COL1 coding region, present in several natural populations of B. nigra, displayed a significant association with flowering time within a majority of these populations. These data indicate that a quantitative trait nucleotide (QTN) affecting flowering time is located within or close to the Bni COL1 gene. The intergenic sequence between Bni COL1 and Bni COa displayed a prominent peak of divergence 1 kb downstream of the Bni COL1 coding region. This region could contain regulatory elements for the downstream Bni COa gene. Our data suggest that a naturally occurring QTN for flowering time affects the function or expression of either Bni COL1 or Bni COa.     

Polymorphism and divergence at three duplicate genes in Brassica nigra

The CONSTANS-like gene family has been shown to evolve exceptionally fast in Brassicaceae. In the present study we analyzed sequence polymorphism and divergence of three genes from this family: COL1 (CONSTANS-LIKE 1) and two copies of CO (CONSTANS), COa and COb, in B. nigra. There was a significant fourfold difference in overall nucleotide diversity among the three genes, with BniCOb having twice as much variation as BniCOL1, which in turn was twice as variable as BniCOa. The ratio of nonsynonymous-to-synonymous substitutions (dN/dS) was high for all three genes, confirming previous studies. While we did not detect evidence of selection at BniCOa and BniCOb, there was a significant excess of polymorphic synonymous mutations in a McDonald-Kreitman test comparing COL1 in B. nigra and A. thaliana. This is apparently the result of an increase in selective constraint on COL1 in B. nigra combined with a decrease in A. thaliana. In conclusion, a complex scenario involving both demography and selection seems to have shaped the pattern of polymorphism at the three genes.     

The pattern of insertion/deletion polymorphism in Arabidopsis thaliana

Little is known about variation of nucleotide insertion/deletions (indels) within species. In Arabidopsis thaliana, we investigated indel polymorphism patterns between two genome sequences and among 96 accessions at 1215 loci. Our study identified patterns in the variation of indel density, size, GC content and distribution, and a correlation between indels and substitutions. We found that the GC content in indel sequences was lower than that in non-indel sequences and that indels typically occur in regions with lower GC content. Patterns of indel frequency distribution among populations were more consistent with neutral expectation than substitution patterns. We also found that the local level of substitutions is positively correlated with indel density and negatively correlated with their distance to the closed indel, suggesting that indels play an important role in nucleotide variation.     

Repetitive Indel Markers within the ALMT1 Gene Conditioning Aluminium Tolerance in Wheat (Triticum aestivum L.)

ALMT1 gene encoding a membrane protein that facilitates an aluminium stimulated malate efflux has been characterised and mapped in wheat (Triticum aestivum L.). Here, we have identified molecular markers targeting insertion/deletion (indel) and SSR repeats within intron 3 region of the ALMT1 gene. Both the markers: ALMT1-SSR3a and ALMT1-SSR3b based on repetitive indels, exhibited complete cosegregation with Al tolerance, malate efflux, and a CAPS marker discriminating ALMT1-1 and ALMT1-2 alleles, in a doubled haploid population derived from Diamondbird (Al-tolerant)/Janz (Al-sensitive). A parental screen of 20 diverse wheat genotypes with repetitive indel markers indicated that six allele variants exist at the ALMT1SSR3 locus. Sequence analysis confirmed that these variations were due to indels, copy number of SSR repeats, and base substitution within SSR repeats. The higher level of variation in intron 3 suggests that this genomic region has been constrained by indels, SSR and single nucleotide polymorphisms. Results have proven that repetitive indel markers cosegregating with the Al tolerance locus will be useful for marker assisted selection and population and evolution studies.     

Inheritance and expression patterns of BN28, a low temperature induced gene in Brassica napus, throughout the Brassicaceae.

Molecular genetics is becoming an important tool in the breeding and selection of agronomically important traits. BN28 is a low temperature induced gene in Brassicaceae species. PCR and Southern blot analysis indicate that BN28 is polymorphic in the three diploid genomes: Brassica rapa (AA), Brassica nigra (BB), and Brassica oleracea (CC). Of the allotetraploids, Brassica napus (AACC) is the only species to have inherited homologous genes from both parental genomes. Brassica juncea (AABB) and Brassica carinata (BBCC) have inherited homologues from the AA and CC genomes, respectively, while Sinapsis arvensis (SS) contains a single homologue from the BB genome and Sinapsis alba (dd) appears to be different from all the diploid parents. All species show message induction when exposed to low temperature. However, differences in expression were noticed at the protein level, with silencing occurring in the BB genome at the level of translation. Results suggest that silencing is occurring in diploid species where duplication may not have occurred. Molecular characterization and inheritance of BN28 homologues in the Brassicaceae may play an important role in determining their quantitative function during exposure to low temperature. Key words : Brassicaceae, BN28, inheritance, polymorphism.     

The genomic landscape of short insertion and deletion polymorphisms in the chicken (Gallus gallus) Genome: a high frequency of deletions in tandem duplicates

It is increasingly recognized that insertions and deletions (indels) are an important source of genetic as well as phenotypic divergence and diversity. We analyzed length polymorphisms identified through partial (0.25x) shotgun sequencing of three breeds of domestic chicken made by the International Chicken Polymorphism Map Consortium. A data set of 140,484 short indel polymorphisms in unique DNA was identified after filtering for microsatellite structures. There was a significant excess of tandem duplicates at indel sites, with deletions of a duplicate motif outnumbering the generation of duplicates through insertion. Indel density was lower in microchromosomes than in macrochromosomes, in the Z chromosome than in autosomes, and in 100 bp of upstream sequence, 5'-UTR, and first introns than in intergenic DNA and in other introns. Indel density was highly correlated with single nucleotide polymorphism (SNP) density. The mean density of indels in pairwise sequence comparisons was 1.9 x 10(-4) indel events/bp, approximately 5% the density of SNPs segregating in the chicken genome. The great majority of indels involved a limited number of nucleotides (median 1 bp), with A-rich motifs being overrepresented at indel sites. The overrepresentation of deletions at tandem duplicates indicates that replication slippage in duplicate sequences is a common mechanism behind indel mutation. The correlation between indel and SNP density indicates common effects of mutation and/or selection on the occurrence of indels and point mutations.     

Analysis of SSR dynamics in chloroplast genomes of Brassicaceae family

Simple sequence repeats (SSRs) are present abundantly in most eukaryotic genomes. They affect several cellular processes like chromatin organization, regulation of gene activity, DNA repair, DNA recombination, etc. Though considerable data exists on using nuclear SSRs to infer phylogenetic relationships, the potential of chloroplast microsatellites (cpSSR), in this regard, remains largely unexplored. In the present study we probe various nucleotide repeat motifs (NRMs) / types of SSRs present in chloroplast genomes (cpDNA) of 12 species belonging to Brassicaceae family. NRMs show a non-random distribution in coding and non-coding compartments of cpDNA. As expected, trinucleotide repeats are more common in coding regions while other repeat motifs are prominent in non-coding DNA. Total numbers of SSRs in coding region show little variation between species while considerable variation is exhibited by SSRs in non-coding regions. Finally, we have designed universal primers that yield polymorphic amplicons from all 12 species. Our analysis also suggests that amplicon length polymorphism shows no significant relationship with sequence based phylogeny of SSRs in cpDNA of Brassicaceae family.     

Comparative mapping between Arabidopsis thaliana and Brassica nigra indicates that Brassica genomes have evolved through extensive genome replication accompanied by chromosome fusions and frequent rearrangements.

Chromosome organization and evolution in the Brassicaceae family was studied using comparative linkage mapping. A total of 160 mapped Arabidopsis thaliana DNA fragments identified 284 homologous loci covering 751 cM in Brassica nigra. The data support that modern diploid Brassica species are descended from a hexaploid ancestor, and that the A. thaliana genome is similar in structure and complexity to those of each of the hypothetical diploid progenitors of the proposed hexaploid. Thus, the Brassica lineage probably went through a triplication after the divergence of the lineages leading to A. thaliana and B. nigra. These duplications were also accompanied by an exceptionally high rate of chromosomal rearrangements. The average length of conserved segments between A. thaliana and B. nigra was estimated at 8 cM. This estimate corresponds to approximately 90 rearrangements since the divergence of the two species. The estimated rate of chromosomal rearrangements is higher than any previously reported data based on comparative mapping. Despite the large number of rearrangements, fine-scale comparative mapping between model plant A. thaliana and Brassica crops is likely to result in the identification of a large number of genes that affect important traits in Brassica crops.     

Development of genome-wide insertion and deletion markers for maize,based on next-generation sequencing data

Background

Insertions and deletions (indels) are the most abundant form of structural variation in all genomes. Indels have been increasingly recognized as an important source of molecular markers due to high-density occurrence, cost-effectiveness, and ease of genotyping. Coupled with developments in bioinformatics, next-generation sequencing (NGS) platforms enable the discovery of millions of indel polymorphisms by comparing the whole genome sequences of individuals within a species.

Results

A total of 1,973,746 unique indels were identified in 345 maize genomes, with an overall density of 958.79 indels/Mbp, and an average allele number of 2.76, ranging from 2 to 107. There were 264,214 indels with polymorphism information content (PIC) values greater than or equal to 0.5, accounting for 13.39 % of overall indels. Of these highly polymorphic indels, we designed primer pairs for 83,481 and 29,403 indels with major allele differences (i.e. the size difference between the most and second most frequent alleles) greater than or equal to 3 and 8 bp, respectively, based on the differing resolution capabilities of gel electrophoresis. The accuracy of our indel markers was experimentally validated, and among 100 indel markers, average accuracy was approximately 90 %. In addition, we also validated the polymorphism of the indel markers. Of 100 highly polymorphic indel markers, all had polymorphisms with average PIC values of 0.54.

Conclusions

The maize genome is rich in indel polymorphisms. Intriguingly, the level of polymorphism in genic regions of the maize genome was higher than that in intergenic regions. The polymorphic indel markers developed from this study may enhance the efficiency of genetic research and marker-assisted breeding in maize.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1797-5) contains supplementary material, which is available to authorized users.     

Incorporating indels as phylogenetic characters: Impact for interfamilial relationships within Arctoidea (Mammalia: Carnivora)

Insertion and deletion events (indels) provide a suite of markers with enormous potential for molecular phylogenetics. Using many more indel characters than those in previous studies, we here for the first time address the impact of indel inclusion on the phylogenetic inferences of Arctoidea (Mammalia: Carnivora). Based on 6843 indel characters from 22 nuclear intron loci of 16 species of Arctoidea, our analyses demonstrate that when the indels were not taken into consideration, the monophyly of Ursidae and Pinnipedia tree and the monophyly of Pinnipedia and Musteloidea tree were both recovered, whereas inclusion of indels by using three different indel coding schemes give identical phylogenetic tree topologies supporting the monophyly of Ursidae and Pinnipedia. Our work brings new perspectives on the previously controversial placements among Arctoidea families, and provides another example demonstrating the importance of identifying and incorporating indels in the phylogenetic analyses of introns. In addition, comparison of indel incorporation methods revealed that the three indel coding methods are all advantageous over treating indels as missing data, given that incorporating indels produces consistent results across methods. This is the first report of the impact of different indel coding schemes on phylogenetic reconstruction at the family level in Carnivora, which indicates that indels should be taken into account in the future phylogenetic analyses.     

Recent degeneration of an old duplicated flowering time gene in Brassica nigra

Gene and genome duplications play a major role in the evolution of plant species. The Brassica nigra genome is highly replicated as a result of ancient polyploidization events. Two copies of the flowering time gene CONSTANS (COa and COb) have been identified in B. nigra, and previous studies showed that COa is functional. In the present study, the polymorphism of 92 COb alleles sampled in seven populations was analyzed. Both polymorphism and recombination levels were elevated and varied strongly among populations and 8% of COb alleles exhibit apparently disabling mutations. Sequence data, however, do not provide unambiguous support for the presence of relaxed selective constraint on COb as compared to known functional CO genes. On the one hand, some of the disabling mutations reached high-frequency arguing for a loss of function but, on the other hand, the ratio of nonsynonymous to synonymous nucleotide polymorphism and diversity is low and similar to that observed in other B. nigra CO and CO-like genes, supporting the conservation of some function. We also showed that COb is still transcribed. Finally, the flowering time of Arabidopsis thaliana co mutant plants transformed with COb alleles with and without apparent disabling mutations was similar. We propose that COb was retained for a long period after duplication, but a recent fixation of a detrimental mutation, possibly as an effect of a bottleneck, resulted in its nonfunctionalization. We also speculate as to the presence of subsequent selection for rapid degeneration of the gene.     

Phylogenetic analysis of <Emphasis Type="Italic">Brassiceae</Emphasis> based on the nucleotide sequences of the <Emphasis Type="Italic">S</Emphasis>-locus related gene, <Emphasis Type="Italic">SLR1</Emphasis>

Nucleotide sequences of orthologs of the S-locus related gene, SLR1, in 20 species of Brassicaceae were determined and compared with the previously reported SLR1 sequences of six species. Identities of deduced amino-acid sequences with Brassica oleracea SLR1 ranged from 66.0% to 97.6%, and those with B. oleracea SRK and SLR2 were less than 62% and 55%, respectively. In multiple alignment of deduced amino-acid sequences, the 180-190th amino-acid residues from the initial methionine were highly variable, this variable region corresponding to hypervariable region I of SLG and SRK. A phylogenetic tree based on the deduced amino-acid sequences showed a close relationship of SLR1 orthologs of species in the Brassicinae and Raphaninae. Brassica nigra SLR1 was found to belong to the same clade as Sinapis arvensis and Diplotaxis siifolia, while the sequences of the other Brassica species belonged to another clade together with B. oleracea and Brassica rapa. The phylogenetic tree was similar to previously reported trees constructed using the data of electrophoretic band patterns of chloroplast DNA, though minor differences were found. Based on synonymous substitution rates in SLR1, the diversification time of SLR1 orthologs between species in the Brassicinae was estimated. The evolution and function of SLR1 and the phylogenetic relationship of Brassiceae plants are discussed.     

Pattern and rate of indel evolution inferred from whole chloroplast intergenic regions in sugarcane, maize and rice.

Microstructural changes such as insertions and deletions (=indels) are a major driving force in the evolution of non-coding DNA sequences. To better understand the mechanisms by which indel mutations arise, as well as the molecular evolution of non-coding regions, the number and pattern of indels and nucleotide substitutions were compared in the whole chloroplast genomes. Comparisons were made for a total of over 38 kb non-coding DNA sequences from 126 intergenic regions in two data sets representing species with different divergence times: sugarcane and maize and Oryza sativa var. indica and japonica. The main findings of this study are: (i) Approximately half of all indels are single nucleotide indels. This observation agrees with previous studies in various organisms. (ii) The distribution and number of indels was different between two data sets, and different patterns were observed for tandem repeat and non-repeat indels. (iii) Distribution pattern of tandem repeat indels showed statistically significant bias towards A/T-rich. (iv) The rate of indel mutation was estimated to be approximately 0.8 +/- 0.04 x 10(-9) per site per year, which was similar to previous estimates in other organisms. (v) The frequencies of nucleotide substitutions and indels were significantly lower in inverted repeat (IR).     

Patterns and relative rates of nucleotide and insertion/deletion evolution at six chloroplast intergenic regions in new world species of the Lecythidaceae

Insertions and deletions (indels) in chloroplast noncoding regions are common genetic markers to estimate population structure and gene flow, although relatively little is known about indel evolution among recently diverged lineages such as within plant families. Because indel events tend to occur nonrandomly along DNA sequences, recurrent mutations may generate homoplasy for indel haplotypes. This is a potential problem for population studies, because indel haplotypes may be shared among populations after recurrent mutation as well as gene flow. Furthermore, indel haplotypes may differ in fitness and therefore be subject to natural selection detectable as rate heterogeneity among lineages. Such selection could contribute to the spatial patterning of cpDNA haplotypes, greatly complicating the interpretation of cpDNA population structure. This study examined both nucleotide and indel cpDNA variation and divergence at six noncoding regions (psbB-psbH, atpB-rbcL, trnL-trnH, rpl20-5'rps12, trnS-trnG, and trnH-psbA) in 16 individuals from eight species in the Lecythidaceae and a Sapotaceae outgroup. We described patterns of cpDNA changes, assessed the level of indel homoplasy, and tested for rate heterogeneity among lineages and regions. Although regression analysis of branch lengths suggested some degree of indel homoplasy among the most divergent lineages, there was little evidence for indel homoplasy within the Lecythidaceae. Likelihood ratio tests applied to the entire phylogenetic tree revealed a consistent pattern rejecting a molecular clock. Tajima's 1D and 2D tests revealed two taxa with consistent rate heterogeneity, one showing relatively more and one relatively fewer changes than other taxa. In general, nucleotide changes showed more evidence of rate heterogeneity than did indel changes. The rate of evolution was highly variable among the six cpDNA regions examined, with the trnS-trnG and trnH-psbA regions showing as much as 10% and 15% divergence within the Lecythidaceae. Deviations from rate homogeneity in the two taxa were constant across cpDNA regions, consistent with lineage-specific rates of evolution rather than cpDNA region-specific natural selection. There is no evidence that indels are more likely than nucleotide changes to experience homoplasy within the Lecythidaceae. These results support a neutral interpretation of cpDNA indel and nucleotide variation in population studies within species such as Corythophora alta.     

Pervasive indels and their evolutionary dynamics after the fish-specific genome duplication

Insertions and deletions (indels) in protein-coding genes are important sources of genetic variation. Their role in creating new proteins may be especially important after gene duplication. However, little is known about how indels affect the divergence of duplicate genes. We here study thousands of duplicate genes in five fish (teleost) species with completely sequenced genomes. The ancestor of these species has been subject to a fish-specific genome duplication (FSGD) event that occurred approximately 350 Ma. We find that duplicate genes contain at least 25% more indels than single-copy genes. These indels accumulated preferentially in the first 40 my after the FSGD. A lack of widespread asymmetric indel accumulation indicates that both members of a duplicate gene pair typically experience relaxed selection. Strikingly, we observe a 30-80% excess of deletions over insertions that is consistent for indels of various lengths and across the five genomes. We also find that indels preferentially accumulate inside loop regions of protein secondary structure and in regions where amino acids are exposed to solvent. We show that duplicate genes with high indel density also show high DNA sequence divergence. Indel density, but not amino acid divergence, can explain a large proportion of the tertiary structure divergence between proteins encoded by duplicate genes. Our observations are consistent across all five fish species. Taken together, they suggest a general pattern of duplicate gene evolution in which indels are important driving forces of evolutionary change.     

Cross-species amplification of mitochondrial DNA sequence-tagged-site markers in conifers: the nature of polymorphism and variation within and among species in Picea

Primers previously developed to amplify specific non-coding regions of the mitochondrial genome in Angiosperms, and new primers for additional non-coding mtDNA regions, were tested for their ability to direct DNA amplification in 12 conifer taxa and to detect sequence-tagged-site (STS) polymorphisms within and among eight species in Picea. Out of 12 primer pairs, nine were successful at amplifying mtDNA in most of the taxa surveyed. In conifers, indels and substitutions were observed for several loci, allowing them to distinguish between families, genera and, in some cases, between species within genera. In Picea, interspecific polymorphism was detected for four loci, while intraspecific variation was observed for three of the mtDNA regions studied. One of these (SSU rRNA V1 region) exhibited indel polymorphisms, and the two others ( nad1 intron b/c and nad5 intron1) revealed restriction differences after digestion with Sau3AI (PCR-RFLP). A fourth locus, the nad4L- orf25 intergenic region, showed a multibanding pattern for most of the spruce species, suggesting a possible gene duplication. Maternal inheritance, expected for mtDNA in conifers, was observed for all polymorphic markers except the intergenic region nad4L- orf25. Pooling of the variation observed with the remaining three markers resulted in two to six different mtDNA haplotypes within the different species of Picea. Evidence for intra-genomic recombination was observed in at least two taxa. Thus, these mitotypes are likely to be more informative than single-locus haplotypes. They should be particularly useful for the study of biogeography and the dynamics of hybrid zones.