Molecular genetics of coat colour variations in White Galloway and White Park cattle

White Galloway cattle exhibit three different white coat colour phenotypes, that is, well marked, strongly marked and mismarked. However, mating of individuals with the preferred well or strongly marked phenotype also results in offspring with the undesired mismarked and/or even fully black coat colour. To elucidate the genetic background of the coat colour variations in White Galloway cattle, we analysed four coat colour relevant genes: mast/stem cell growth factor receptor (KIT), KIT ligand (KITLG), melanocortin 1 receptor (MC1R) and tyrosinase (TYR). Here, we show that the coat colour variations in White Galloway cattle and White Park cattle are caused by a KIT gene (chromosome 6) duplication and aberrant insertion on chromosome 29 (Cs₂₉) as recently described for colour‐sided Belgian Blue. Homozygous (Cs₂₉/Cs₂₉) White Galloway cattle and White Park cattle exhibit the mismarked phenotype, whereas heterozygous (Cs₂₉/wt₂₉) individuals are either well or strongly marked. In contrast, fully black individuals are characterised by the wild‐type chromosome 29. As known for other cattle breeds, mutations in the MC1R gene determine the red colouring. Our data suggest that the white coat colour variations in White Galloway cattle and White Park cattle are caused by a dose‐dependent effect based on the ploidy of aberrant insertions and inheritance of the KIT gene on chromosome 29.     

The Uyghur population and genetic susceptibility to type 2 diabetes: potential role for variants in CAPN10, APM1 and FUT6 genes

Genome‐wide association studies have successfully identified over 70 loci associated with the risk of type 2 diabetes mellitus (T2DM) in multiple populations of European ancestry. However, the risk attributable to an individual variant is modest and does not yet provide convincing evidence for clinical utility. Association between these established genetic variants and T2DM in general populations is hitherto understudied in the isolated populations, such as the Uyghurs, resident in Hetian, far southern Xinjiang Uyghur Autonomous Region, China. In this case–control study, we genotyped 13 single‐nucleotide polymorphisms (SNPs) at 10 genes associated with diabetes in 130 cases with T2DM and 135 healthy controls of Uyghur, a Chinese minority ethnic group. Three of the 13 SNPs demonstrated significant association with T2DM in the Uyghur population. There were significant differences between the T2DM patients and controls in the risk allele distributions of rs3792267 (CAPN10) (P = 0.002), rs1501299 (APM1) (P = 0.017), and rs3760776 (FUT6) (P = 0.031). Allelic carriers of rs3792267‐A, rs1501299‐T, and rs3760776‐T had a 2.24‐fold [OR (95% CI): 1.35–3.71], 0.59‐fold [OR (95% CI): 0.39–0.91], 0.57‐fold [OR (95% CI): 0.34–0.95] increased risk for T2DM respectively. We further confirmed that the cumulative risk allelic scores calculated from the 13 susceptibility loci for T2DM differed significantly between the T2DM patients and controls (P = 0.001), and the effect of obesity/overweight on T2DM was only observed in the subjects with a combined risk allelic score under a value of 17. This study observed that the SNPs rs3792267 in CAPN10, rs1501299 in APM1, and rs3760776 in FUT6 might serve as potential susceptible biomarkers for T2DM in Uyghurs. The cumulative risk allelic scores of multiple loci with modest individual effects are also significant risk factors in Uyghurs for T2DM, particularly among non‐obese individuals. This is the first investigation having observed/found genetic variations on genetic loci functionally linked with glycosylation associated with the risk of T2DM in a Uyghur population.     

FtsHi1/ARC1 is an essential gene in Arabidopsis that links chloroplast biogenesis and division

The Arabidopsis arc1 (accumulation and replication of chloroplasts 1) mutant has pale seedlings and smaller, more numerous chloroplasts than the wild type. Previous work has suggested that arc1 affects the timing of chloroplast division but does not function directly in the division process. We isolated ARC1 by map‐based cloning and discovered it encodes FtsHi1 (At4g23940), one of several FtsHi proteins in Arabidopsis. These poorly studied proteins resemble FtsH metalloproteases important for organelle biogenesis and protein quality control but are presumed to be proteolytically inactive. FtsHi1 bears a predicted chloroplast transit peptide and localizes to the chloroplast envelope membrane. Phenotypic studies showed that arc1 (hereafter ftsHi1‐1), which bears a missense mutation, is a weak allele of FtsHi1 that disrupts thylakoid development and reduces de‐etiolation efficiency in seedlings, suggesting that FtsHi1 is important for chloroplast biogenesis. Consistent with this finding, transgenic plants suppressed for accumulation of an FtsHi1 fusion protein were often variegated. A strong T‐DNA insertion allele, ftsHi1‐2, caused embryo‐lethality, indicating that FtsHi1 is an essential gene product. A wild‐type FtsHi1 transgene rescued both the chloroplast division and pale phenotypes of ftsHi1‐1 and the embryo‐lethal phenotype of ftsHi1‐2. FtsHi1 overexpression produced a subtle increase in chloroplast size and decrease in chloroplast number in wild‐type plants while suppression led to increased numbers of small chloroplasts, providing new evidence that FtsHi1 negatively influences chloroplast division. Taken together, our analyses reveal that FtsHi1 functions in an essential, envelope‐associated process that may couple plastid development with division.     

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.     

Next‐generation sequencing (NGS)‐based identification of induced mutations in a doubly mutagenized tomato (Solanum lycopersicum) population

The identification of mutations in targeted genes has been significantly simplified by the advent of TILLING (Targeting Induced Local Lesions In Genomes), speeding up the functional genomic analysis of animals and plants. Next‐generation sequencing (NGS) is gradually replacing classical TILLING for mutation detection, as it allows the analysis of a large number of amplicons in short durations. The NGS approach was used to identify mutations in a population of Solanum lycopersicum (tomato) that was doubly mutagenized by ethylmethane sulphonate (EMS). Twenty‐five genes belonging to carotenoids and folate metabolism were PCR‐amplified and screened to identify potentially beneficial alleles. To augment efficiency, the 600‐bp amplicons were directly sequenced in a non‐overlapping manner in Illumina MiSeq, obviating the need for a fragmentation step before library preparation. A comparison of the different pooling depths revealed that heterozygous mutations could be identified up to 128‐fold pooling. An evaluation of six different software programs (camba , crisp , gatk unified genotyper , lofreq , snver and vipr ) revealed that no software program was robust enough to predict mutations with high fidelity. Among these, crisp and camba predicted mutations with lower false discovery rates. The false positives were largely eliminated by considering only mutations commonly predicted by two different software programs. The screening of 23.47 Mb of tomato genome yielded 75 predicted mutations, 64 of which were confirmed by Sanger sequencing with an average mutation density of 1/367 Kb. Our results indicate that NGS combined with multiple variant detection tools can reduce false positives and significantly speed up the mutation discovery rate.     

Generation of chromosomal deletions in dicotyledonous plants employing a user‐friendly genome editing toolkit

Genome editing facilitated by Cas9‐based RNA‐guided nucleases (RGNs) is becoming an increasingly important and popular technique for reverse genetics in both model and non‐model species. So far, RGNs were mainly applied for the induction of point mutations, and one major challenge consists in the detection of genome‐edited individuals from a mutagenized population. Also, point mutations are not appropriate for functional dissection of non‐coding DNA. Here, the multiplexing capacity of a newly developed genome editing toolkit was exploited for the induction of inheritable chromosomal deletions at six different loci in Nicotiana benthamiana and Arabidopsis. In both species, the preferential formation of small deletions was observed, suggesting reduced efficiency with increasing deletion size. Importantly, small deletions (<100 bp) were detected at high frequencies in N. benthamiana T₀ and Arabidopsis T₂ populations. Thus, targeting of small deletions by paired nucleases represents a simple approach for the generation of mutant alleles segregating as size polymorphisms in subsequent generations. Phenotypically selected deletions of up to 120 kb occurred at low frequencies in Arabidopsis, suggesting larger population sizes for the discovery of valuable alleles from addressing gene clusters or non‐coding DNA for deletion by programmable nucleases.     

Genetic control of α‐farnesene production in apple fruit and its role in fungal pathogenesis

Terpenes are important compounds in plant trophic interactions. A meta‐analysis of GC‐MS data from a diverse range of apple (Malus × domestica) genotypes revealed that apple fruit produces a range of terpene volatiles, with the predominant terpene being the acyclic branched sesquiterpene (E,E)‐α‐farnesene. Four quantitative trait loci (QTLs) for α‐farnesene production in ripe fruit were identified in a segregating ‘Royal Gala’ (RG) × ‘Granny Smith’ (GS) population with one major QTL on linkage group 10 co‐locating with the MdAFS1 (α‐farnesene synthase‐1) gene. Three of the four QTLs were derived from the GS parent, which was consistent with GC‐MS analysis of headspace and solvent‐extracted terpenes showing that cold‐treated GS apples produced higher levels of (E,E)‐α‐farnesene than RG. Transgenic RG fruit downregulated for MdAFS1 expression produced significantly lower levels of (E,E)‐α‐farnesene. To evaluate the role of (E,E)‐α‐farnesene in fungal pathogenesis, MdAFS1 RNA interference transgenic fruit and RG controls were inoculated with three important apple post‐harvest pathogens [Colletotrichum acutatum, Penicillium expansum and Neofabraea alba (synonym Phlyctema vagabunda)]. From results obtained over four seasons, we demonstrate that reduced (E,E)‐α‐farnesene is associated with decreased disease initiation rates of all three pathogens. In each case, the infection rate was significantly reduced 7 days post‐inoculation, although the size of successful lesions was comparable with infections on control fruit. These results indicate that (E,E)‐α‐farnesene production is likely to be an important factor involved in fungal pathogenesis in apple fruit.     

Polymorphisms in the AKT1 and AKT2 genes and oesophageal squamous cell carcinoma risk in an Eastern Chinese population

Ethnic Han Chinese are at high risk of developing oesophageal squamous cell carcinoma (ESCC). Aberrant activation of the AKT signalling pathway is involved in many cancers, including ESCC. Some single nucleotide polymorphisms (SNPs) in genes involved in this pathway may contribute to ESCC susceptibility. We selected five potentially functional SNPs in AKT1 (rs2494750, rs2494752 and rs10138277) and AKT2 (rs7254617 and rs2304186) genes and investigated their associations with ESCC risk in 1117 ESCC cases and 1096 controls in an Eastern Chinese population. None of individual SNPs exhibited an association with ESCC risk. However, the combined analysis of three AKT1 SNPs suggested that individuals carrying one of AKT1 variant genotypes had a decreased ESCC risk [adjusted odds ratio (OR) = 0.60, 95% CI = 0.42–0.87]. Further stratified analysis found that AKT1 rs2294750 SNP was associated with significantly decreased ESCC risk among women (adjusted OR = 0.63, 95% CI = 0.43–0.94) and non‐drinkers (OR = 0.79, 95% CI = 0.64–0.99). Similar protective effects on women (adjusted OR = 0.56, 95% CI = 0.37–0.83) and non‐drinker (adjusted OR = 0.75, 95% CI = 0.60–0.94) were also observed for the combined genotypes of AKT1 SNPs. Consistently, logistic regression analysis indicated significant gene–gene interactions among three AKT1 SNPs (P < 0.015). A three‐AKT1 SNP haplotype (C‐A‐C) showed a significant association with a decreased ESCC risk (adjusted OR = 0.70, 95% CI = 0.52–0.94). Multifactor dimensionality reduction analysis confirmed a high‐order gene–environment interaction in ESCC risk. Overall, we found that three AKT1 SNPs might confer protection against ESCC risk; nevertheless, these effects may be dependent on other risk factors. Our results provided evidence of important gene–environment interplay in ESCC carcinogenesis.     

The genetics of brown coat color and white spotting in domestic yaks (Bos grunniens)

Domestic yaks (Bos grunniens) exhibit two major coat color variations: a brown vs. wild‐type black pigmentation and a white spotting vs. wild‐type solid color pattern. The genetic basis for these variations in color and distribution remains largely unknown and may be complicated by a breeding history involving hybridization between yaks and cattle. Here, we investigated 92 domestic yaks from China using a candidate gene approach. Sequence variations in MC1R, PMEL and TYRP1 were surveyed in brown yaks; TYRP1 was unassociated with the coloration and excluded. Recessive mutations from MC1R, or p.Gln34*, p.Met73Leu and possibly p.Arg142Pro, are reported in bovids for the first time and accounted for approximately 40% of the brown yaks in this study. The remaining 60% of brown individuals correlated with a cattle‐derived deletion mutation from PMEL (p.Leu18del) in a dominant manner. Degrees of white spotting found in yaks vary from color sidedness and white face, to completely white. After examining the candidate gene KIT, we suggest that color‐sided and all‐white yaks are caused by the serial translations of KIT (Cs₆ or Cs₂₉) as reported for cattle. The white‐faced phenotype in yaks is associated with the KIT haplotype S^wf. All KIT mutations underlying the serial phenotypes of white spotting in yaks are identical to those in cattle, indicating that cattle are the likely source of white spotting in yaks. Our results reveal the complex genetic origins of domestic yak coat color as either native in yaks through evolution and domestication or as introduced from cattle through interspecific hybridization.     

Detection of induced mutations in CaFAD2 genes by next‐generation sequencing leading to the production of improved oil composition in Crambe abyssinica

Crambe abyssinica is a hexaploid oil crop for industrial applications. An increase of erucic acid (C22:1) and reduction of polyunsaturated fatty acid (PUFA) contents in crambe oil is a valuable improvement. An increase in oleic acid (C18:1), a reduction in PUFA and possibly an increase in C22:1 can be obtained by down‐regulating the expression of fatty acid desaturase2 genes (CaFAD2), which code for the enzyme that converts C18:1 into C18:2. We conducted EMS‐mutagenesis in crambe, followed by Illumina sequencing, to screen mutations in three expressed CaFAD2 genes. Two novel analysis strategies were used to detect mutation sites. In the first strategy, mutation detection targeted specific sequence motifs. In the second strategy, every nucleotide position in a CaFAD2 fragment was tested for the presence of mutations. Seventeen novel mutations were detected in 1100 one‐dimensional pools (11 000 individuals) in three expressed CaFAD2 genes, including non‐sense mutations and mis‐sense mutations in CaFAD2‐C1, ‐C2 and ‐C3. The homozygous non‐sense mutants for CaFAD2‐C3 resulted in a 25% higher content of C18:1 and 25% lower content of PUFA compared to the wild type. The mis‐sense mutations only led to small changes in oil composition. Concluding, targeted mutation detection using NGS in a polyploid was successfully applied and it was found that a non‐sense mutation in even a single CaFAD2 gene can lead to changes in crambe oil composition. Stacking the mutations in different CaFAD2 may gain additional changes in C18:1 and PUFA contents.     

Development of high‐density SNP genotyping arrays for white spruce (Picea glauca) and transferability to subtropical and nordic congeners

High‐density SNP genotyping arrays can be designed for any species given sufficient sequence information of high quality. Two high‐density SNP arrays relying on the Infinium iSelect technology (Illumina) were designed for use in the conifer white spruce (Picea glauca). One array contained 7338 segregating SNPs representative of 2814 genes of various molecular functional classes for main uses in genetic association and population genetics studies. The other one contained 9559 segregating SNPs representative of 9543 genes for main uses in population genetics, linkage mapping of the genome and genomic prediction. The SNPs assayed were discovered from various sources of gene resequencing data. SNPs predicted from high‐quality sequences derived from genomic DNA reached a genotyping success rate of 64.7%. Nonsingleton in silico SNPs (i.e. a sequence polymorphism present in at least two reads) predicted from expressed sequenced tags obtained with the Roche 454 technology and Illumina GAII analyser resulted in a similar genotyping success rate of 71.6% when the deepest alignment was used and the most favourable SNP probe per gene was selected. A variable proportion of these SNPs was shared by other nordic and subtropical spruce species from North America and Europe. The number of shared SNPs was inversely proportional to phylogenetic divergence and standing genetic variation in the recipient species, but positively related to allele frequency in P. glauca natural populations. These validated SNP resources should open up new avenues for population genetics and comparative genetic mapping at a genomic scale in spruce species.     

Development and preliminary evaluation of a genomewide single nucleotide polymorphisms resource generated by RAD‐seq for the small yellow croaker (Larimichthys polyactis)

Recent advances in high‐throughput sequencing technologies have offered the possibility to generate genomewide sequence data to delineate previously unidentified genetic structure, obtain more accurate estimates of demographic parameters and to evaluate potential adaptive divergence. Here, we identified 27 556 single nucleotide polymorphisms for the small yellow croaker (Larimichthys polyactis) using restriction‐site‐associated DNA (RAD) sequencing of 24 individuals from two populations. Significant sources of genetic variation were identified, with an average nucleotide diversity (π) of 0.00105 ± 0.000425 across individuals, and long‐term effective population size was thus estimated to range between 26 172 and 261 716. According to the results, no differentiation between the two populations was detected based on the SNP data set of top quality score per contig or neutral loci. However, the two analysed populations were highly differentiated based on SNP data set of both top F_ST value per contig and the outlier SNPs. Moreover, local adaptation was highlighted by an F_ST‐based outlier tests implemented in LOSITAN and a total of 538 potentially locally selected SNPs were identified. blast2go annotation of contigs containing the outlier SNPs yielded hits for 37 (66%) of 56 significant blastx matches. Candidate genes for local adaptation constituted a wide array of biological functions, including cellular response to oxidative stress, actin filament binding, ion transmembrane transport and synapse assembly. The generated SNP resources in this study provided a valuable tool for future population genetics and genomics studies of L. polyactis.     

Targeting novel chemical and constitutive primed metabolites against Plectosphaerella cucumerina

Priming is a physiological state for protection of plants against a broad range of pathogens, and is achieved through stimulation of the plant immune system. Various stimuli, such as beneficial microbes and chemical induction, activate defense priming. In the present study, we demonstrate that impairment of the high‐affinity nitrate transporter 2.1 (encoded by NRT2.1) enables Arabidopsis to respond more quickly and strongly to Plectosphaerella cucumerina attack, leading to enhanced resistance. The Arabidopsis thaliana mutant lin1 (affected in NRT2.1) is a priming mutant that displays constitutive resistance to this necrotroph, with no associated developmental or growth costs. Chemically induced priming by β–aminobutyric acid treatment, the constitutive priming mutant ocp3 and the constitutive priming present in the lin1 mutant result in a common metabolic profile within the same plant–pathogen interactions. The defense priming significantly affects sugar metabolism, cell‐wall remodeling and shikimic acid derivatives levels, and results in specific changes in the amino acid profile and three specific branches of Trp metabolism, particularly accumulation of indole acetic acid, indole‐3–carboxaldehyde and camalexin, but not the indolic glucosinolates. Metabolomic analysis facilitated identification of three metabolites in the priming fingerprint: galacturonic acid, indole‐3–carboxylic acid and hypoxanthine. Treatment of plants with the latter two metabolites by soil drenching induced resistance against P. cucumerina, demonstrating that these compounds are key components of defense priming against this necrotrophic fungus. Here we demonstrate that indole‐3–carboxylic acid induces resistance by promoting papillae deposition and H₂O₂ production, and that this is independent of PR1, VSP2 and PDF1.2 priming.     

RNA‐Seq bulked segregant analysis enables the identification of high‐resolution genetic markers for breeding in hexaploid wheat

The identification of genetic markers linked to genes of agronomic importance is a major aim of crop research and breeding programmes. Here, we identify markers for Yr15, a major disease resistance gene for wheat yellow rust, using a segregating F₂ population. After phenotyping, we implemented RNA sequencing (RNA‐Seq) of bulked pools to identify single‐nucleotide polymorphisms (SNP) associated with Yr15. Over 27 000 genes with SNPs were identified between the parents, and then classified based on the results from the sequenced bulks. We calculated the bulk frequency ratio (BFR) of SNPs between resistant and susceptible bulks, selecting those showing sixfold enrichment/depletion in the corresponding bulks (BFR > 6). Using additional filtering criteria, we reduced the number of genes with a putative SNP to 175. The 35 SNPs with the highest BFR values were converted into genome‐specific KASP assays using an automated bioinformatics pipeline (PolyMarker) which circumvents the limitations associated with the polyploid wheat genome. Twenty‐eight assays were polymorphic of which 22 (63%) mapped in the same linkage group as Yr15. Using these markers, we mapped Yr15 to a 0.77‐cM interval. The three most closely linked SNPs were tested across varieties and breeding lines representing UK elite germplasm. Two flanking markers were diagnostic in over 99% of lines tested, thus providing a reliable haplotype for marker‐assisted selection in these breeding programmes. Our results demonstrate that the proposed methodology can be applied in polyploid F₂ populations to generate high‐resolution genetic maps across target intervals.     

Development of SNP‐genotyping arrays in two shellfish species

Use of SNPs has been favoured due to their abundance in plant and animal genomes, accompanied by the falling cost and rising throughput capacity for detection and genotyping. Here, we present in vitro (obtained from targeted sequencing) and in silico discovery of SNPs, and the design of medium‐throughput genotyping arrays for two oyster species, the Pacific oyster, Crassostrea gigas, and European flat oyster, Ostrea edulis. Two sets of 384 SNP markers were designed for two Illumina GoldenGate arrays and genotyped on more than 1000 samples for each species. In each case, oyster samples were obtained from wild and selected populations and from three‐generation families segregating for traits of interest in aquaculture. The rate of successfully genotyped polymorphic SNPs was about 60% for each species. Effects of SNP origin and quality on genotyping success (Illumina functionality Score) were analysed and compared with other model and nonmodel species. Furthermore, a simulation was made based on a subset of the C. gigas SNP array with a minor allele frequency of 0.3 and typical crosses used in shellfish hatcheries. This simulation indicated that at least 150 markers were needed to perform an accurate parental assignment. Such panels might provide valuable tools to improve our understanding of the connectivity between wild (and selected) populations and could contribute to future selective breeding programmes.