Sequence variation in two lignin biosynthesis genes,cinnamoyl CoA reductase (CCR) and cinnamyl alcohol dehydrogenase 2 (CAD2)

Cinnamoyl CoA reductase (CCR) and cinnamyl alcohol dehydrogenase 2 (CAD2) are genes which may influence variation in lignin content and composition within plants. Sequence variation within these genes may be responsible for changes in enzyme activity and/or specificity, which could cause variation in lignin content or composition. This study examines sequence variation within these two genes in Eucalyptus globulus, an important species used in pulp and paper-making. Twenty-one single nucleotide polymorphisms (SNPs) were identified in the exons of CCR, of which nine were neutral mutations and 12 were missense mutations. Six of the missense mutations affected highly conserved amino acids within the protein sequence of CCR. Eight SNPs were identified in the CAD2 exons, six of which were neutral mutations and two which were missense mutations. One of the missense mutations affected a highly conserved amino acid within the protein sequence. In addition, 32 SNPs were identified in the CCR introns along with four insertion/deletions and two polyA length variation regions. Polymorphism affecting highly conserved amino acids may alter enzyme function and this molecular variation may be linked to variation in lignin profiles. Selecting positive alleles which produce favourable lignin profiles would be advantageous in tree breeding programs.     

Detection and validation of EST‐derived SNPs for poplar leaf rust Melampsora medusae f. sp. deltoidae

We report the discovery, characterization and validation of 118 single nucleotide polymorphisms (SNPs) for poplar leaf rust Melampsora medusae f. sp. deltoidae identified using a gene‐targeted approach in an expressed sequence tag (EST) library. We developed a genotyping assay using the iPLEX? primer extension method for two multiplex assays of 28 and 22 SNPs.     

Connecting genomic patterns of local adaptation and niche suitability in teosintes

The central abundance hypothesis predicts that local adaptation is a function of the distance to the centre of a species’ geographic range. To test this hypothesis, we gathered genomic diversity data from 49 populations, 646 individuals and 33,464 SNPs of two wild relatives of maize, the teosintes Zea mays ssp. parviglumis and Zea. mays. ssp. mexicana. We examined the association between the distance to their climatic and geographic centroids and the enrichment of SNPs bearing signals of adaptation. We identified candidate adaptive SNPs in each population by combining neutrality tests and cline analyses. By applying linear regression models, we found that the number of candidate SNPs is positively associated with niche suitability, while genetic diversity is reduced at the limits of the geographic distribution. Our results suggest that overall, populations located at the limit of the species’ niches are adapting locally. We argue that local adaptation to this limit could initiate ecological speciation processes and facilitate adaptation to global change.     

Evaluation of microsatellite markers in association studies: a search for an immune-related susceptibility gene in sarcoidosis

Association studies using linkage disequilibrium (LD) between candidate loci and nearby markers have been proposed to identify susceptibility genes for complex diseases. We analyzed polymorphisms of microsatellites (MSs) and LD patterns of the regions in which candidate genes related to the Th1 immune response have been annotated and attempted to identify a susceptibility gene for sarcoidosis in a marker-based association study. Nineteen MSs were identified in six Th1-related genes (IFNGR1, IFNGR2, IL12RB1, IL12RB2, STAT1 and STAT4) and then eight were further characterized as useful polymorphic markers. Most of these MSs showed LD with single nucleotide polymorphisms (SNPs) on both 5 and 3 ends of these candidate genes, in which r² values between at least one of the MS marker alleles and the SNPs were higher than 0.1. A significant association with one MS allele near STAT4 was shown and a cluster of SNPs in LD with the MS marker was associated with sarcoidosis. These results suggest that association studies using not only SNPs but also multi-allelic MS within or near candidate loci would be useful markers to search for a disease susceptibility gene, especially in populations with unknown LD structure.     

Clade classification of monolignol biosynthesis gene family members reveals target genes to decrease lignin in Lolium perenne

In monocots, lignin content has a strong impact on the digestibility of the cell wall fraction. Engineering lignin biosynthesis requires a profound knowledge of the role of paralogues in the multigene families that constitute the monolignol biosynthesis pathway. We applied a bioinformatics approach for genome‐wide identification of candidate genes in Lolium perenne that are likely to be involved in the biosynthesis of monolignols. More specifically, we performed functional subtyping of phylogenetic clades in four multigene families: 4CL, COMT, CAD and CCR. Essential residues were considered for functional clade delineation within these families. This classification was complemented with previously published experimental evidence on gene expression, gene function and enzymatic activity in closely related crops and model species. This allowed us to assign functions to novel identified L. perenne genes, and to assess functional redundancy among paralogues. We found that two 4CL paralogues, two COMT paralogues, three CCR paralogues and one CAD gene are prime targets for genetic studies to engineer developmentally regulated lignin in this species. Based on the delineation of sequence conservation between paralogues and a first analysis of allelic diversity, we discuss possibilities to further study the roles of these paralogues in lignin biosynthesis, including expression analysis, reverse genetics and forward genetics, such as association mapping. We propose criteria to prioritise paralogues within multigene families and certain SNPs within these genes for developing genotyping assays or increasing power in association mapping studies. Although L. perenne was the target of the analyses presented here, this functional subtyping of phylogenetic clades represents a valuable tool for studies investigating monolignol biosynthesis genes in other monocot species.     

Patterns of DNA sequence variation at candidate gene loci in black poplar (<Emphasis Type="Italic">Populus nigra</Emphasis> L.) as revealed by single nucleotide polymorphisms

Black poplar (Populus nigra L.) is an economically and ecologically important tree species and an ideal organism for studies of genetic variation. In the present work, we use a candidate gene approach to infer the patterns of DNA variation in natural populations of this species. A total of 312 single nucleotide polymorphisms (SNPs) are found among 8,056 bp sequenced from nine drought-adaptation and photosynthesis-related gene loci. The median SNP frequency is one site per 26 bp. The average nucleotide diversity is calculated to be θ_W = 0.01074 and π_T = 0.00702, higher values than those observed in P. tremula, P. trichocarpa and most conifer species. Tests of neutrality for each gene reveal a general excess of low-frequency mutations, a greater number of haplotypes than expected and an excess of high-frequency derived variants in P. nigra, which is consistent with previous findings that genetic hitchhiking has occurred in this species. Linkage disequilibrium is low, decaying rapidly from 0.45 to 0.20 or less within a distance of 300 bp, although the declines of r ² are variable among different loci. This is similar to the rate of decay reported in most other tree species. Our dataset is expected to enhance understanding of how evolutionary forces shape genetic variation, and it will contribute to molecular breeding in black poplar.     

Development of SNP markers and haplotype analysis of the candidate gene for rhg1, which confers resistance to soybean cyst nematode in soybean

Soybean cyst nematode (SCN; Heterodera glycines Ichinohe) is one of the most destructive pests in the cultivation of soybean (Glycine max (L.) Merr.) worldwide. Markers based on the SCN resistance gene will enable efficient marker-assisted selection (MAS). We sequenced the candidate gene rhg1 in six resistant and two susceptible soybean genotypes and identified 37 SNPs (single nucleotide polymorphisms) among the sequences, of which 11 were in the coding region. Seven of these 11 SNPs led to changes in the amino acid sequence of the gene. The amino acid sequence we obtained differs from the previously published one by a stretch of 26–27 amino acids. Six codominant allele-specific SNP markers based on agarose gel detection were developed and tested in 70 genotypes, among which occurred only nine different haplotypes. Two neutrality tests (Tajima’s D and Fu and Li’s F) were significant for the six SNP loci in the 70 genotypes, which is consistent with intensive directional selection. A strong LD pattern was detected among five SNPs except 2868T > C. Two SNPs (689C > A and 757C > T) formed one haplotype (689C-757C) that was perfectly associated with SCN resistance. The new allele-specific PCR markers located in the alleged sequence of the rhg1 candidate gene, combined with the microsatellite marker BACR-Satt309, will significantly improve the efficiency of MAS during the development of SCN-resistant cultivars.     

Identification of Single Nucleotide Polymorphisms and Analysis of Linkage Disequilibrium in Different Bamboo Species Using the Candidate Gene Approach

Bamboos are one of the most beautiful and useful plants on Earth. The genetic background and population structure of bamboos are well known, which helps accelerate the process of artificial domestication of bamboo. Partial sequences of six genes involved in nitrogen use efficiency in 32 different bamboo species were analyzed for occurrence of single nucleotide polymorphisms (SNPs). The nucleotide diversity θw and total nucleotide polymorphisms π_T of the sequenced DNA regions was 0.05137 and 0.03332, respectively. Both πnonsyn /πsyn and Ka/Ks values were <1. The nucleotide sequences of these six genes were inferred to be relatively conserved, and the haplotype diversity was relatively high. The results of evolutionary neutrality tests showed that the six genes were in line with neutral evolution, and that the NRT2.1 and AMT2.1 gene sequences may have experienced negative selection. An inter-SNP recombination event at the NRT2.1 gene in the all pooled sample, of all 32 bamboo species was the lowest at 0.0645, whereas the AMT gene recombination events were all >0.1. Estimation and analysis of linkage disequilibrium of five genes revealed that with the increase in nucleotide sequence length, the degree of SNP linkage disequilibrium decreased rapidly. We inferred the population genetic structure of 32 bamboo species based on the SNP loci of six genes with frequencies >18%. 32 bamboo species were divided into five categories, which indicated that the combined population of all bamboo species had obvious multivariate characteristics and was heterogeneous; red (Group 1) and green (Group 2) were the main groups.     

Study of polymorphisms in the promoter region of ovine β-lactoglobulin gene and phylogenetic analysis among the Valle del Belice breed and other sheep breeds considered as ancestors

The aim of this work was to sequence the promoter region of β-lactoglobulin (BLG) gene in four sheep breeds, in order to identify polymorphisms, infer and analyze haplotypes, and phylogenetic relationship among the Valle del Belice breed and the other three breeds considered as ancestors. Sequencing analysis and alignment of the obtained sequences showed the presence of 36 single nucleotide polymorphisms (SNPs) and one deletion. A total of 22 haplotypes found in “best” reconstruction were inferred considering the 37 polymorphic sites identified. Haplotypes were used for the reconstruction of a phylogenetic tree using the Neighbor-Joining algorithm. The number of polymorphisms identified showed high variability within breeds. Analysis of genetic diversity indexes showed that the Sarda breed presented the lowest nucleotide diversity, whereas the Comisana breed presented the highest one. Comparing the nucleotide diversity among breeds, the highest value was obtained between Valle del Belice and Pinzirita breeds, whereas the lowest one was between Valle del Belice and Sarda breeds. Considering that polymorphisms in the promoter region of BLG gene could have a functional role associated with milk composition, the lowest value of nucleotide diversity between Valle del Belice and Sarda breeds may be related to a higher similarity of milk composition of these two breeds compared to the others. Further analyses will be conducted in order to evaluate the possible correlation between the genetic diversity indexes and the BLG content in milk of our breeds.     

Disrupting the cinnamyl alcohol dehydrogenase 1 gene (BdCAD1) leads to altered lignification and improved saccharification in Brachypodium distachyon

Brachypodium distachyon (Brachypodium) has been proposed as a model for grasses, but there is limited knowledge regarding its lignins and no data on lignin‐related mutants. The cinnamyl alcohol dehydrogenase (CAD) genes involved in lignification are promising targets to improve the cellulose‐to‐ethanol conversion process. Down‐regulation of CAD often induces a reddish coloration of lignified tissues. Based on this observation, we screened a chemically induced population of Brachypodium mutants (Bd21–3 background) for red culm coloration. We identified two mutants (Bd4179 and Bd7591), with mutations in the BdCAD1 gene. The mature stems of these mutants displayed reduced CAD activity and lower lignin content. Their lignins were enriched in 8–O–4‐ and 4–O–5‐coupled sinapaldehyde units, as well as resistant inter‐unit bonds and free phenolic groups. By contrast, there was no increase in coniferaldehyde end groups. Moreover, the amount of sinapic acid ester‐linked to cell walls was measured for the first time in a lignin‐related CAD grass mutant. Functional complementation of the Bd4179 mutant with the wild‐type BdCAD1 allele restored the wild‐type phenotype and lignification. Saccharification assays revealed that Bd4179 and Bd7591 lines were more susceptible to enzymatic hydrolysis than wild‐type plants. Here, we have demonstrated that BdCAD1 is involved in lignification of Brachypodium. We have shown that a single nucleotide change in BdCAD1 reduces the lignin level and increases the degree of branching of lignins through incorporation of sinapaldehyde. These changes make saccharification of cells walls pre‐treated with alkaline easier without compromising plant growth.     

Evidence for a large-scale population structure of Arabidopsis thaliana from genome-wide single nucleotide polymorphism markers

Population-based methods for the genetic mapping of adaptive traits and the analysis of natural selection require that the population structure and demographic history of a species are taken into account. We characterized geographic patterns of genetic variation in the model plant Arabidopsis thaliana by genotyping 115 genome-wide single nucleotide polymorphism (SNP) markers in 351 accessions from the whole species range using a matrix-assisted laser desorption/ionization time-of-flight assay, and by sequencing of nine unlinked short genomic regions in a subset of 64 accessions. The observed frequency distribution of SNPs is not consistent with a constant-size neutral model of sequence polymorphism due to an excess of rare polymorphisms. There is evidence for a significant population structure as indicated by differences in genetic diversity between geographic regions. Accessions from Central Asia have a low level of polymorphism and an increased level of genome-wide linkage disequilibrium (LD) relative to accessions from the Iberian Peninsula and Central Europe. Cluster analysis with the structure program grouped Eurasian accessions into K=6 clusters. Accessions from the Iberian Peninsula and from Central Asia constitute distinct populations, whereas Central and Eastern European accessions represent admixed populations in which genomes were reshuffled by historical recombination events. These patterns likely result from a rapid postglacial recolonization of Eurasia from glacial refugial populations. Our analyses suggest that mapping populations for association or LD mapping should be chosen from regional rather than a species-wide sample or identified genetically as sets of individuals with similar average genetic distances. Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

Minimum sample sizes for population genomics: an empirical study from an Amazonian plant species

High‐throughput DNA sequencing facilitates the analysis of large portions of the genome in nonmodel organisms, ensuring high accuracy of population genetic parameters. However, empirical studies evaluating the appropriate sample size for these kinds of studies are still scarce. In this study, we use double‐digest restriction‐associated DNA sequencing (ddRADseq) to recover thousands of single nucleotide polymorphisms (SNPs) for two physically isolated populations of Amphirrhox longifolia (Violaceae), a nonmodel plant species for which no reference genome is available. We used resampling techniques to construct simulated populations with a random subset of individuals and SNPs to determine how many individuals and biallelic markers should be sampled for accurate estimates of intra‐ and interpopulation genetic diversity. We identified 3646 and 4900 polymorphic SNPs for the two populations of A. longifolia, respectively. Our simulations show that, overall, a sample size greater than eight individuals has little impact on estimates of genetic diversity within A. longifolia populations, when 1000 SNPs or higher are used. Our results also show that even at a very small sample size (i.e. two individuals), accurate estimates of F_ST can be obtained with a large number of SNPs (≥1500). These results highlight the potential of high‐throughput genomic sequencing approaches to address questions related to evolutionary biology in nonmodel organisms. Furthermore, our findings also provide insights into the optimization of sampling strategies in the era of population genomics.     

The adaptive potential of Populus balsamifera L. to phenology requirements in a warmer global climate

The manner in which organisms adapt to climate change informs a broader understanding of the evolution of biodiversity as well as conservation and mitigation plans. We apply common garden and association mapping approaches to quantify genetic variance and identify loci affecting bud flush and bud set, traits that define a tree's season for height growth, in the boreal forest tree Populus balsamifera L. (balsam poplar). Using data from 478 genotypes grown in each of two common gardens, one near the southern edge and another near the northern edge of P. balsamifera's range, we found that broad‐sense heritability for bud flush and bud set was generally high (H² > 0.5 in most cases), suggesting that abundant genetic variation exists for phenological response to changes in the length of the growing season. To identify the molecular genetic basis of this variation, we genotyped trees for 346 candidate single nucleotide polymorphisms (SNPs) from 27 candidate genes for the CO/FT pathway in poplar. Mixed‐model analyses of variance identified SNPs in 10 genes to be associated with variation in either bud flush or bud set. Multiple SNPs within FRIGIDA were associated with bud flush, whereas multiple SNPs in LEAFY and GIGANTEA 5 were associated with bud set. Although there was strong population structure in stem phenology, the geographic distribution of multilocus association SNP genotypes was widespread except at the most northern populations, indicating that geographic regions may harbour sufficient diversity in functional genes to facilitate adaption to future climatic conditions in many sites.     

Brown‐midrib maize (bm1) – a mutation affecting the cinnamyl alcohol dehydrogenase gene

Brown-midrib (bm) mutants of maize have modified lignin of reddish-brown colour. Although four independent bm loci are known, only one of the mutant genes has been previously identified. We report here that maize bm1, one of the less characterised mutants, shows severely reduced CAD activity in lignified tissues, resulting in the production of a modified lignin. Both the total lignin content and the structure of the polymer are altered by the mutation. We further describe the isolation and characterisation of the maize CAD cDNA and mapping of the CAD gene. CAD maps very closely to the known location of bm1 and co-segregates with the bm1 locus in two independent recombinant inbred populations. These data strongly support the premise that maize bm1 directly affects expression of the CAD gene.     

Association mapping of morphological traits in wild and captive zebra finches: reliable within,but not between populations

Identifying causal genetic variants underlying heritable phenotypic variation is a long‐standing goal in evolutionary genetics. We previously identified several quantitative trait loci (QTL) for five morphological traits in a captive population of zebra finches (Taeniopygia guttata) by whole‐genome linkage mapping. We here follow up on these studies with the aim to narrow down on the quantitative trait variants (QTN) in one wild and three captive populations. First, we performed an association study using 672 single nucleotide polymorphisms (SNPs) within candidate genes located in the previously identified QTL regions in a sample of 939 wild‐caught zebra finches. Then, we validated the most promising SNP–phenotype associations (n = 25 SNPs) in 5228 birds from four populations. Genotype–phenotype associations were generally weak in the wild population, where linkage disequilibrium (LD) spans only short genomic distances. In contrast, in captive populations, where LD blocks are large, apparent SNP effects on morphological traits (i.e. associations) were highly repeatable with independent data from the same population. Most of those SNPs also showed significant associations with the same trait in other captive populations, but the direction and magnitude of these effects varied among populations. This suggests that the tested SNPs are not the causal QTN but rather physically linked to them, and that LD between SNPs and causal variants differs between populations due to founder effects. While the identification of QTN remains challenging in nonmodel organisms, we illustrate that it is indeed possible to confirm the location and magnitude of QTL in a population with stable linkage between markers and causal variants.     

Implication of genetic variants near TMEM18, BCDIN3D/FAIM2, and MC4R with coronary artery disease and obesity in Chinese: a angiography-based study

Coronary artery disease (CAD) is multifactorial disease which occurs as a result of the interaction of genetic and environmental factors. Obesity is an independent risk factor for cardiovascular disease. Recent genome-wide association studies have identified several genes associated with obesity in Europeans. We wondered whether these genetic variants were associated with CAD. Three single nucleotide polymorphisms (SNPs) rs7561317 near TMEM18, rs7138803 near BCDIN3D/FAIM2 and rs12970134 near MC4R were examined in 930 Han Chinese subjects based on coronary angiography, using polymerase chain reaction (PCR) followed by restriction fragment length polymorphism (RFLP) analysis. There were no significant differences in genotypes and allele distributions of three SNPs between CAD and CAD-free groups. The AA genotype of SNP rs12970134 near MC4R was associated to obesity both in CAD group and CAD-free group in Han Chinese population (P < 0.001, OR = 2.96, 95% CI 2.01–3.73; and P = 0.003, OR = 2.59, 95% CI 1.86–3.19, respectively). Our observations suggest that the polymorphism rs12970134 near MC4R may be associated to the risk of obesity in Han Chinese population.     

Analysis of dyslexia candidate genes in the Raine cohort representing the general Australian population

Several genes have been suggested as dyslexia candidates. Some of these candidate genes have been recently shown to be associated with literacy measures in sample cohorts derived from the general population. Here, we have conducted an association study in a novel sample derived from the Australian population (the Raine cohort) to further investigate the role of dyslexia candidate genes. We analysed markers, previously reported to be associated with dyslexia, located within the MRPL19/C2ORF3, KIAA0319, DCDC2 and DYX1C1 genes in a sample of 520 individuals and tested them for association with reading and spelling measures. Association signals were detected for several single nucleotide polymorphisms (SNPs) within DYX1C1 with both the reading and spelling tests. The high linkage disequilibrium (LD) we observed across the DYX1C1 gene suggests that the association signal might not be refined by further genetic mapping.     

Multilocus OCA2 genotypes specify human iris colors

Human iris color is a quantitative, multifactorial phenotype that exhibits quasi-Mendelian inheritance. Recent studies have shown that OCA2 polymorphism underlies most of the natural variability in human iris pigmentation but to date, only a few associated polymorphisms in this gene have been described. Herein, we describe an iris color score (C) for quantifying iris melanin content in-silico and undertake a more detailed survey of the OCA2 locus (n = 271 SNPs). In 1,317 subjects, we confirmed six previously described associations and identified another 27 strongly associated with C that were not explained by continental population stratification (OR 1.5–17.9, P = 0.03 to <0.001). Haplotype analysis with respect to these 33 SNPs revealed six haplotype blocks and 11 hap-tags within these blocks. To identify genetic features for best-predicting iris color, we selected sets of SNPs by parsing P values among possible combinations and identified four discontinuous and non-overlapping sets across the LD blocks (p-Selected SNP sets). In a second, partially overlapping sample of 1,072, samples with matching diplotypes comprised of these p-Selected OCA2 SNPs exhibited a rate of C concordance of 96.3% (n = 82), which was significantly greater than that obtained from randomly selected samples (62.6%, n = 246, P<0.0001). In contrast, the rate of C concordance using diplotypes comprised of the 11 identified hap-tags was only 83.7%, and that obtained using diplotypes comprised of all 33 SNPs organized as contiguous sets along the locus (defined by the LD block structure) was only 93.3%. These results confirm that OCA2 is the major human iris color gene and suggest that using an empirical database-driven system, genotypes from a modest number of SNPs within this gene can be used to accurately predict iris melanin content from DNA. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.