Natural variation in GmGBP1 promoter affects photoperiod control of flowering time and maturity in soybean

The present study screened for polymorphisms in coding and non‐coding regions of the GmGBP1 gene in 278 soybean accessions with variable maturity and growth habit characteristics under natural field conditions in three different latitudes in China. The results showed that the promoter region was highly diversified compared with the coding sequence of GmGBP1. Five polymorphisms and four haplotypes were closely related to soybean flowering time and maturity through association and linkage disequilibrium analyses. Varieties with the polymorphisms SNP_‐796G, SNP_‐770G, SNP_‐307T, InDel_‐242normal, SNP_353A, or haplotypes Hap‐3 and Hap‐4 showed earlier flowering time and maturity in different environments. The shorter growth period might be largely due to higher GmGBP1 expression levels in soybean that were caused by the TCT‐motif with SNP_‐796G in the promoter. In contrast, the lower expression level of GmGBP1 in soybean caused by RNAi interference of GmGBP1 resulted in a longer growth period under different day lengths. Furthermore, the gene interference of GmGBP1 also caused a reduction in photoperiod response sensitivity (PRS) before flowering in soybean. RNA‐seq analysis on GmGBP1 underexpression in soybean showed that 94 and 30 predicted genes were significantly upregulated and downregulated, respectively. Of these, the diurnal photoperiod‐specific expression pattern of three significant flowering time genes GmFT2a, GmFT5a, and GmFULc also showed constantly lower mRNA levels in GmGBP1‐i soybean than in wild type, especially under short day conditions. Together, the results showed that GmGBP1 functioned as a positive regulator upstream of GmFT2a and GmFT5a to activate the expression of GmFULc to promote flowering on short days.     

Nucleotide variability and haplotype heterogeneity at the porcine fat mass and obesity‐associated (FTO) gene

A large number of studies have confirmed that variants within the fat mass and obesity‐associated (FTO) gene are associated with higher obesity risk in humans. We and others have shown that FTO polymorphisms are associated with fat deposition and related traits in several pig populations, thus confirming the role of this gene in fatness across species. However, some differences observed in different pig populations may be derived, at least in part, from genetic heterogeneity at this locus. Here, we characterise the nucleotide variability and haplotype diversity of the porcine FTO gene in breeds having different predispositions to fat deposition traits. We resequenced 4749 bp of coding and non‐coding regions of the porcine FTO gene in 44 pigs of eight different breeds and identified 27 single nucleotide polymorphisms (SNPs) and four insertions/deletions. A positive Tajima's D‐value (P < 0.10) obtained in Italian Duroc pigs may be compatible with putative balancing selection. From the sequenced pig panel, 20 haplotypes were inferred, some of which clustered according to the breed of origin (Meishan and Italian Duroc). Genetic heterogeneity at this locus could complicate the dissection of the effects of this gene on fat deposition and production traits in pigs. This situation resembles, to some extent, what has been reported in humans, thus making the study of the porcine FTO gene variability especially interesting, as it could be used as a model to understand the complex and elusive role of this gene in human obesity.     

HERC1 polymorphisms: population‐specific variations in haplotype composition

Human HERC1 is one of six HERC proteins and may play an important role in intracellular membrane trafficking. The human HERC1 gene is suggested to have been affected by local positive selection. To assess the global frequency distributions of coding and non‐coding single nucleotide polymorphisms (SNPs) in the HERC1 gene, we developed a new simultaneous genotyping method for four SNPs, and applied this method to investigate 1213 individuals from 12 global populations. The results confirmed remarked differences in the allele and haplotype frequencies between East Asian and non‐East Asian populations. One of the three common haplotypes observed was found to be characteristic of East Asians, who showed a relatively uniform distribution of haplotypes. Information on haplotypes would be useful for testing the function of polymorphisms in the HERC1 gene. This is the first study to investigate the distribution of HERC1 polymorphisms in various populations. Copyright © 2009 John Wiley & Sons, Ltd.     

Association of dopamine‐ and serotonin‐related genes with canine aggression

Human‐directed canine aggression was studied using 50 aggressive and 81 non‐aggressive dogs. We examined 62 single nucleotide polymorphisms (SNPs) occurring in or in the close vicinity of 16 neurotransmitter‐related genes. Allelic associations with aggression were identified for DRD1, HTR1D, HTR2C and SLC6A1. Risk or protective haplotypes for aggressive behaviour based on 2–5 SNPs were identified. The frequency of aggressive dogs varied significantly between the haplotypes within loci and the odds ratios of aggression in dogs with risk haplotypes compared with protective haplotypes varied from 4.4 (HTR2C) to 9.0 (SLC6A1). A risk haplotype across the neurotransmitter receptor gene HTR1D harboured a non‐synonymous SNP with a potential effect on protein function. We identified no haplotypes in complete association with the recorded phenotypes, supporting a complex inheritance of aggression.     

KIT and melanoma predisposition in pigs: sequence variants and association analysis

KIT mutations have been detected in different cancer subtypes, including melanoma. The gene also has been extensively studied in farm animals for its prominent role in coat color. The present work aimed at detecting KIT variants in a porcine model of cutaneous melanoma, the melanoblastoma‐bearing Libechov Minipig (MeLiM). By sequencing exons and intron borders, 36 SNPs and one indel were identified. Of 10 coding SNPs, three were non‐synonymous mutations, likely to affect the protein conformation. A promising variant, located in exon 19 (p.Val870Ala), was genotyped in a MeLiM × Duroc cross, and an association analysis was conducted on several melanoma‐related traits. This variant showed a significant association with melanoma development, tumor ulceration and cutaneous invasion. In conclusion, although the KIT gene would not be a major causal gene for melanoma development in pig, its genetic variation could be influencing this trait.     

Chromosome-level wild Hevea brasiliensis genome provides new tools for genomic-assisted breeding and valuable loci to elevate rubber yield

The rubber tree (Hevea brasiliensis) is grown in tropical regions and is the major source of natural rubber. Using traditional breeding approaches, the latex yield has increased by sixfold in the last century. However, the underlying genetic basis of rubber yield improvement is largely unknown. Here, we present a high-quality, chromosome-level genome sequence of the wild rubber tree, the first report on selection signatures and a genome-wide association study (GWAS) of its yield traits. Population genomic analysis revealed a moderate population divergence between the Wickham clones and wild accessions. Interestingly, it is suggestive that H. brasiliensis and six relatives of the Hevea genus might belong to the same species. The selective sweep analysis found 361 obvious signatures in the domesticated clones associated with 245 genes. In a 15-year field trial, GWAS identified 155 marker–trait associations with latex yield, in which 326 candidate genes were found. Notably, six genes related to sugar transport and metabolism, and four genes related to ethylene biosynthesis and signalling are associated with latex yield. The homozygote frequencies of the causal nonsynonymous SNPs have been greatly increased under selection, which may have contributed to the fast latex yield improvement during the short domestication history. Our study provides insights into the genetic basis of the latex yield trait and has implications for genomic-assisted breeding by offering valuable resources in this new domesticated crop.     

Haplotype hitchhiking promotes trait coselection in Brassica napus

Local haplotype patterns surrounding densely spaced DNA markers with significant trait associations can reveal information on selective sweeps and genome diversity associated with important crop traits. Relationships between haplotype and phenotype diversity, coupled with analysis of gene content in conserved haplotype blocks, can provide insight into coselection for nonrelated traits. We performed genome‐wide analysis of haplotypes associated with the important physiological and agronomic traits leaf chlorophyll and seed glucosinolate content, respectively, in the major oilseed crop species Brassica napus. A locus on chromosome A01 showed opposite effects on leaf chlorophyll content and seed glucosinolate content, attributed to strong linkage disequilibrium (LD) between orthologues of the chlorophyll biosynthesis genes EARLY LIGHT‐INDUCED PROTEIN and CHLOROPHYLL SYNTHASE, and the glucosinolate synthesis gene ATP SULFURYLASE 1. Another conserved haplotype block, on chromosome A02, contained a number of chlorophyll‐related genes in LD with orthologues of the key glucosinolate biosynthesis genes METHYLTHIOALKYMALATE SYNTHASE‐LIKE 1 and 3. Multigene haplogroups were found to have a significantly greater contribution to variation for chlorophyll content than haplotypes for any single gene, suggesting positive effects of additive locus accumulation. Detailed reanalysis of population substructure revealed a clade of ten related accessions exhibiting high leaf chlorophyll and low seed glucosinolate content. These accessions each carried one of the above‐mentioned haplotypes from A01 or A02, generally in combination with further chlorophyll‐associated haplotypes from chromosomes A05 and/or C05. The phenotypic rather than pleiotropic correlations between leaf chlorophyll content index and seed GSL suggest that LD may have led to inadvertent coselection for these two traits.     

Insight into the population structure of hardhead silverside,Atherinomorus stipes (Teleostei: Atherinidae), in Belize and the Florida Keys using nd2

Little is known about the natural history, biology, and population genetic structure of the Hardhead Silverside, Atherinomorus stipes, a small schooling fish found around islands throughout the Caribbean. Our field observations of A. stipes in the cays of Belize and the Florida Keys found that populations tend to be in close association with the shoreline in mangrove habitats. Due to this potential island‐based population structuring, A. stipes represents an ideal system to examine questions about gene flow and isolation by distance at different geographic scales. For this study, the mitochondrial gene nd2 was amplified from 394 individuals collected from seven different Belizean Cays (N = 175) and eight different Floridian Keys (N = 219). Results show surprisingly high haplotype diversity both within and between island‐groups, as well as a high prevalence of unique haplotypes within each island population. The results are consistent with models that require gene flow among populations as well as in situ evolution of rare haplotypes. There was no evidence for an isolation by distance model. The nd2 gene tree consists of two well‐supported monophyletic groups: a Belizean‐type clade and a Floridian‐type clade, indicating potential species‐level differentiation.     

Expression of the sucrose transporter 3 (<Emphasis Type="Italic">HbSUT3</Emphasis>) in rubber tree and its relation to latex yield

Sucrose is not only a precursor of rubber biosynthesis, but it also participates in other latex regeneration between repeated tapping. Sucrose transporter (HbSUT) genes play key roles in sinks and sources of the rubber biosynthesis pathway. HbSUT3 was dominant in expression in latex among the SUT member genes. Therefore, the expression level of HbSUT3 in latex is a potential indicator distinguishing between high- and low-yielding clones. The aim of this research was to assess the potential of this gene in selective breeding to improve latex yield, from the correlation of HbSUT3 expression with latex yield. Four high-yielding clones were sampled in this study and compared with the common RRIM 600 clone as the baseline, paired by the field. Among the putative full-length cDNAs of Hevea sucrose transporters available in the NCBI database, only HbSUT3 was detected. The HbSUT3 gene was overexpressed in the four rubber clones relative to the control, and the NK1 clone had the highest expression level. The expression level of HbSUT3 correlated positively with latex yield but negatively with the sucrose content of latex. Gene expression analysis of rubber seedlings indicated that the bark had higher expression of this gene than the leaves, and the levels correlated with latex yields of these clones. These data provided new candidate selection criteria for use in the early selection of high-yielding rubber clones, necessary for rapid cycle breeding programs.     

Effect of swine leukocyte antigen‐DQA gene variation on diarrhea in Large White,Landrace, and Duroc piglets

Piglet diarrhea is one of the most common factors that affects the benefits of the swine industry. Although recent studies have shown that exon 2 of SLA‐DQA is associated with piglet resistance to diarrhea, contributions of genetic variation in the additional exon coding regions of this gene remain unclear. Here, we investigated variation in exons 1, 3 and 4 of the SLA‐DQA gene and evaluated their effects on diarrheal infection in 425 suckling piglets. No variation was identified in exon 1. In exon 3, there were eight alleles detected, generated by 14 single nucleotide polymorphisms (SNPs) and three nucleotide deletions, eight SNPs being newly identified. Four allele sequences and three SNPs were identified in exon 4, only one SNP being newly identified. Statistical analysis showed that the genotypes of exon 3 are significantly associated with piglet diarrhea; indeed, genotypes DQA*wb01/wb02 and wb04/wb05 are clearly associated with resistance to piglet diarrhea, as they have the lowest probabilities of infection (P < 0.05). However, no significant association was found between the genotypes of exon 4 and diarrhea (P > 0.05). These results provide important new information concerning the level of genetic diversity at the SLA‐DQA locus and suggest that further genetic association studies of piglet diarrhea resistance should include analyses of both exons 2 and 3 of this locus.     

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.     

Unique haplotypes in ant‐attended aphids and widespread haplotypes in non‐attended aphids

Aphid species within the genus Tuberculatus Mordvilko (Hemiptera: Aphididae) exhibit a variety of interactions with ants, ranging from close associations to non‐attendance. A previous study indicated that despite wing possession, ant‐attended Tuberculatus species exhibited low dispersal rates compared with non‐attended species. This study examined if presence or absence of mutualistic interactions and habitat continuity of host plants affected intraspecific genetic diversity and genetic differentiation in mitochondrial DNA cytochrome oxidase I (COI) sequences. Sympatric ant‐attended Tuberculatus quercicola (Matsumura) (Hemiptera: Aphididae) and non‐attended Tuberculatus paiki Hille Ris Lambers (Hemiptera: Aphididae) were collected from the daimyo oak Quercus dentata Thunberg (Fagales: Fagaceae) in Japan and examined for haplotype variability. Seventeen haplotypes were identified in 568 T. quercicola individuals representing 23 populations and seven haplotypes in 425 T. paiki representing 19 populations. Haplotype diversity, which indicates the mean number of differences between all pairs of haplotypes in the sample, and nucleotide diversity were higher in T. quercicola than T. paiki. Analysis of molecular variance (AMOVA) showed higher genetic differentiation among populations within groups of T. quercicola (39.8%) than T. paiki (22.6%). The effects of attendant ant species on genetic differentiation in T. quercicola were not distinguishable from geographic factors. Despite low dispersal rates, host plant habitat continuity might facilitate widespread dispersal of a T. quercicola haplotype in Hokkaido. These results suggested that following T. quercicola colonization, gene flow among populations was limited, resulting in genetic drift within populations. However, frequent T. paiki dispersal is clearly evident by low genetic differentiation among populations within groups, resulting in lower haplotype diversity.     

Novel Y‐chromosome polymorphisms in Chinese domestic yak

Y‐chromosome‐specific haplotypes (Y‐haplotypes) constructed using single nucleotide polymorphisms (Y‐SNPs) in the MSY (male‐specific region of the Y‐chromosome) are valuable in population genetic studies. But sequence variants in the yak MSY region have been poorly characterized so far. In this study, we screened a total of 16 Y‐chromosome‐specific gene segments from the ZFY, SRY, UTY, USP9Y, AMELY and OFD1Y genes to identify Y‐SNPs in domestic yaks. Six novel Y‐SNPs distributed in the USP9Y (g.223C>T), UTY19 (g.158A>C and g.169C>T), AMELY2 (g.261C>T), OFD1Y9 (g.165A>G) and SRY4 (g.104G>A) loci, which can define three Y‐haplotypes (YH1, YH2 and YH3) in yaks, were discovered. YH1 was the dominant and presumably most ancient haplotype based on the comparison of UTY19 locus with other bovid species. Interestingly, we found informative UTY19 markers (g.158A>C and g.169C>T) that can effectively distinguish the three yak Y‐haplotypes. The nucleotide diversity was 1.7 × 10^?4 ± 0.3 × 10^?4, indicating rich Y‐chromosome diversity in yaks. We identified two highly divergent lineages (YH1 and YH2 vs. YH3) that share similar frequencies (YH1 +  YH2: 0.82–0.89, YH3: 0.11–0.18) among all three populations. In agreement with previous mtDNA studies, we supported the hypothesis that the two highly divergent lineages (YH1 and YH2 vs. YH3) derived from a single gene pool, which can be explained by the reunion of at least two paternal populations with the divergent lineages already accumulated before domestication. We estimated a divergence time of 408 110 years between the two divergent lineages, which is consistent with the data from mitochondrial DNA in yaks.     

Deep sequencing of a candidate region harboring the SOX9 gene for the canine XX disorder of sex development

A disorder of sex development (DSD) in dogs with female sex chromosomes (78, XX), a lack of the SRY gene and the presence of testes or ovotestes is commonly diagnosed in numerous breeds. The molecular background of DSD is not fully recognized but has been linked to the copy number variation in the region harboring the SOX9 gene. We applied a genome‐wide association study and targeted next‐generation sequencing techniques to compare DSD and normal female dogs. The genome‐wide association study did not indicate a significant chromosome region. Targeted next‐generation sequencing of a 1.5‐Mb region on canine chromosome 9 harboring the SOX9 gene revealed two putatively DSD‐associated copy number variations 355 kb upstream and 691 kb downstream of SOX9, four blocks of low polymorphism and two blocks of an elevated heterozygosity. An initial next‐generation sequencing analysis showed an association with two SNPs, but validation in larger cohorts did not confirm this result. We identified a large homologous fragment (over 243.8 kb), named hfMAGI2, located upstream of SOX9, that overlaps a known copy number variation region. It shows a high sequence similarity with the 5′ flanking region of the MAGI2 gene located on canine chromosome 18 that encodes a protein involved in ovary formation during early embryonic development. Our study showed that the identified copy number variation region located upstream of the SOX9 gene contains potential regulatory sequences (long non‐coding RNA and hfMAGI2) and led to the assumption that a multiplication of this element may alter expression of the SOX9 gene, triggering the DSD phenotype.     

A genome‐wide association study in Caucasian women suggests the involvement of HLA genes in the severity of facial solar lentigines

Solar lentigines are a common feature of sun‐induced skin ageing. Little is known, however, about the genetic factors contributing to their development. In this genome‐wide association study, we aimed to identify genetic loci associated with solar lentigines on the face in 502 middle‐aged French women. Nine SNPs, gathered in two independent blocks on chromosome 6, exhibited a false discovery rate below 25% when looking for associations with the facial lentigine score. The first block, in the 6p22 region, corresponded to intergenic SNPs and also exhibited a significant association with forehead lentigines (P = 1.37 × 10^?8). The second block, within the 6p21 HLA region, was associated with decreased HLA‐C expression according to several eQTL databases. Interestingly, these SNPs were also in high linkage disequilibrium with the HLA‐C*0701 allele (r² = 0.95). We replicated an association recently found by GWAS in the IRF4 gene. Finally, a complementary study on 44 selected candidate SNPs revealed novel associations in the MITF gene. Overall, our results point to several mechanisms involved in the severity of facial lentigines, including HLA/immunity and the melanogenesis pathway.     

Novel Chloroplast Microsatellite (cpSSR) Markers for Genetic Diversity Assessment of Cultivated and Wild Hevea Rubber

Understanding the genetic diversity and population structure of Hevea is of great importance for managing its conservation and for utilization of rubber genetic resources. In this study, we investigated the genetic diversity and population structure of eight populations of Hevea rubber genotypes from Malaysia (MY), India (IN), Sri Lanka (LK), Indonesia (ID), France (FR), Thailand (TH), Brazil (BR), and China (CN), in addition to individual primary clones, using 10 nuclear and 11 polymorphic novel chloroplast microsatellite markers (nSSRs and cpSSRs, respectively). The BR population exhibited the greatest genetic diversity (H _e ) for both nSSR (0.841) and haploid genetic diversity (cpSSR; 0.207) markers. Bayesian analysis was applied to infer genetic structure based on nSSR data using STRUCTURE software, and a neighbor-joining tree was used to construct an unrooted phylogram based on Nei’s genetic distance that clustered these Hevea genotypes into three groups. Six polymorphic cpSSR markers produced 13 alleles and eight chlorotypes. Seven chlorotypes, A, B, C, D, E, F, and H were detected among Brazilian populations from Acre (AC), Rondônia (RO), and Mato Grosso (MT) locations. The G and H chlorotypes were found in Asiatic genotypes or Wickham clones and only one genotype, 15 AC-F-7 38-62 from AC location. These results provide valuable data for in situ or ex situ conservation and utilization of germplasm collections for breeding programs.