凡纳滨对虾分子标记筛选、连锁图谱构建和QTL定位

利用AFLP和微卫星标记，以凡纳滨对虾F1全同胞家系为作图群体，构建了凡纳滨对虾的雌性和雄性遗传连锁图谱并对生长相关性状体长和体重进行了QTL分析。利用经过筛选的108对AFLP引物组合，对亲本和94个子代个体进行了分离分析。共得到2041个多态AFLP标记（1:1分离）。平均每个引物组合产生20个多态片段。有826个AFLP标记偏离孟德尔遗传（40.5%, P<0.05）。所筛选的100个微卫星标记中有30个在家系中有作图信息，分别有24个和20个位点可以用于母本和父本的作图。对父母本β-1，3-葡聚糖结合蛋白(BGBP)，脂多糖和葡聚糖结合蛋白(LGBP)和蜕皮抑制激素（MIH）等基因片段进行序列分析得到一些可用于构建连锁图谱的SNP标记。 对所有的分离标记进行了连锁分析，分别绘制了凡纳滨对虾的雌性和雄性连锁图谱。雌性连锁图谱的框架图有319个遗传标记组成，分布于45个连锁群，其中AFLP标记300个，微卫星位点18个，性别标记一个，连锁群长度从29.5cM到260.0cM, 每个连锁群含4－16个标记。图谱长度为4134.4 cM, 各连锁群平均图距在7.6到25.9 cM之间，总平均图距为15.1 cM。有267个标记（含14个微卫星）整合到雄性框架图上。 雄性框架图也含45个连锁群，长度从14.2cM到161.1cM, 图谱长度为3220.9 cM, 各连锁群平均图距在4.1到25.5 cM之间，总平均图距为14.5 cM。 作图群体所有个体的性别均作为标记整合到雌雄分离信息中，在94个F1个体中，54个为雌虾，40个是雄虾。在雌性图谱第29连连锁群上，性别与三个微卫星标记紧密连锁（v1f148, v145f120, v95f166）, 图距分别为6.6, 8.6 和8.6cM, 相应LOD值分别为17.8 14.3 和 16.4。 与雌性图相对应的是，雄性图谱未发现任何与性别连锁的标记。推测凡纳滨对虾的性别决定机制可能是ZW型，其中雌性为异配性别。 体长和体重均显示出连续变异的特点，显示这些与生长相关的性状都是典型的数量性状或多基因遗传。体重和体长符合正态分布且两性状之间存在着显著的相关性（P<0.001），Pearson 相关系数为0.95。在凡纳滨对虾雌性图谱和雄性图谱上，共定位了6个与体长和体重等生长性状相关的QTL，可解释的表型变异从15.1％到21.3％。共检测到3个与体长相关的QTL，包括两个正向效应的QTL和一个负向的QTL，分别定位到雌性的第6连锁群和雄性的第4和11 连锁群上。体重也检测到3个QTL，其中一个QTL加性效应是正向的，分别位于雌性的第6和7连锁群及雄性的第11连锁群上。QTL主要集中在雌性图谱的Fc6和Fc7及雄性图的Mc4和Mc11上。体长和体重的QTL定位结果中各有两个是十分相似的，其相邻的标记位点完全一样，只是相应的LOD值略有差异，还各存在两性状特异的一个QTL。大多QTL与相邻标记之间的距离只有0.1cM(仅Wfc7的距离较大)，为进一步的精细定位奠定了基础。 分子标记筛选、遗传图谱的构建及生长相关性状的QTL定位为我们下一步从事分子标记辅助选择育种,QTL精细定位和比较基因作图打下基础，并最终推动凡纳滨对虾的遗传改良。

Genetic mapping of size-related quantitative trait loci (QTL) in the bay scallop (Argopecten irradians) using AFLP and microsatellite markers

We constructed genetic linkage maps for the bay scallop Argopecten irradians using AFLP and microsatellite markers and conducted composite interval mapping (CIM) of body-size-related traits. Three hundred seventeen AFLP and 10 microsatellite markers were used for map construction. The female parent map contained 120 markers in 15 linkage groups, spanning 479.6 cM with an average interval of 7.0 cM. The male parent map had 190 markers in 17 linkage groups, covering 883.8 cM at 7.2 cM per marker. The observed coverage was 70.4% for the female and 81.1% for the male map. Markers that were distorted toward the same direction were closely linked to each other on the genetic maps, suggesting the presence of genes important for survival. Six size-related traits, shell length, shell height, shell width, total weight, soft tissue weight, and shell weight, were measured for QTL mapping. The size data were significantly correlated with each other. We subjected the data, log transformed firstly, to a principle component analysis and use the first principle component for QTL mapping. CIM analysis revealed one significant QTL (LOD=2.69, 1000 permutation, P<0.05) in linkage group 3 on the female parent map. The mapping of size-related QTL in this study raises the possibility of improving the growth of bay scallops through marker-assisted selection. (c) 2007 Published by Elsevier B.V.

Two genes on A/J chromosome 18 are associated with susceptibility to Staphylococcus aureus infection by combined microarray and QTL analyses.

Although it has recently been shown that A/J mice are highly susceptible to Staphylococcus aureus sepsis as compared to C57BL/6J, the specific genes responsible for this differential phenotype are unknown. Using chromosome substitution strains (CSS), we found that loci on chromosomes 8, 11, and 18 influence susceptibility to S. aureus sepsis in A/J mice. We then used two candidate gene selection strategies to identify genes on these three chromosomes associated with S. aureus susceptibility, and targeted genes identified by both gene selection strategies. First, we used whole genome transcription profiling to identify 191 (56 on chr. 8, 100 on chr. 11, and 35 on chr. 18) genes on our three chromosomes of interest that are differentially expressed between S. aureus-infected A/J and C57BL/6J. Second, we identified two significant quantitative trait loci (QTL) for survival post-infection on chr. 18 using N(2) backcross mice (F(1) [C18A]xC57BL/6J). Ten genes on chr. 18 (March3, Cep120, Chmp1b, Dcp2, Dtwd2, Isoc1, Lman1, Spire1, Tnfaip8, and Seh1l) mapped to the two significant QTL regions and were also identified by the expression array selection strategy. Using real-time PCR, 6 of these 10 genes (Chmp1b, Dtwd2, Isoc1, Lman1, Tnfaip8, and Seh1l) showed significantly different expression levels between S. aureus-infected A/J and C57BL/6J. For two (Tnfaip8 and Seh1l) of these 6 genes, siRNA-mediated knockdown of gene expression in S. aureus-challenged RAW264.7 macrophages induced significant changes in the cytokine response (IL-1 beta and GM-CSF) compared to negative controls. These cytokine response changes were consistent with those seen in S. aureus-challenged peritoneal macrophages from CSS 18 mice (which contain A/J chromosome 18 but are otherwise C57BL/6J), but not C57BL/6J mice. These findings suggest that two genes, Tnfaip8 and Seh1l, may contribute to susceptibility to S. aureus in A/J mice, and represent promising candidates for human genetic susceptibility studies.

Identificación de QTL para dormición y análisis molecular de la interacción ABA-GAs en cariopses de sorgo granífero

El objetivo de esta tesis fue contribuir a la dilucidación de los mecanismos moleculares y genéticos que participan en la expresión de la dormición de semillas de cereales, utilizando Sorghum bicolor (L.) Moench como sistema modelo. Para ello se utilizaron dos aproximaciones complementarias: la identificación de QTL para el carácter dormición y la evaluación de la ocurrencia de interacciones in vitro entre componentes de la señalización del ácido abscísico (ABA) y el catabolismo de las giberelinas (GAs), candidatos a tener un rol importante durante la expresión de la dormición en granos de sorgo inmaduros (i.e. antes de madurez fisiológica). Los resultados obtenidos permitieron identificar tres QTL (qDOR-5; qDOR-9 y qDOR-10) que explican una proporción de la variabilidad que se observa en el patrón de expresión de dormición de granos de sorgo maduros (i.e. después de madurez fisiológica). Un análisis in silico de las secuencias abarcadas por estos QTL mostró que ninguno ellos incluye genes considerados como candidatos para dormición de sorgo. En ese sentido, esta tesis aportó nuevas regiones genómicas que contienen genes hasta ahora desconocidos, que serían importantes en la expresión del carácter dormición en granos maduros. Por otra parte, los análisis de unión in vitro realizados mostraron que las proteínas SbABI4 y SbABI5 (componentes de la señalización del ABA) pueden interactuar de manera específica con el ABRC (complejo de respuesta al ABA) del promotor del gen SbGA2ox3, responsable de la degradación de giberelinas activas. Este mecanismo de cross-talk ABA-GAs podría ser uno de los responsables del mantenimiento de la dormición en cariopses inmaduros resistentes al brotado pre-cosecha. Más aún, el ABRC del promotor de SbGA2ox3, involucrado en las interacciones, se encontró además en los promotores de genes GA2ox de otras especies monocotiledóneas como Brachypodium y arroz (Oryza sativa), pero no así en las dicotiledóneas analizadas, sugiriendo que el cross-talk ABA-GAs podría tener lugar en otras especies además de sorgo. Los resultados de esta tesis en forma conjunta aportaron nuevas evidencias acerca del rol preponderante que tienen ciertas regiones del genoma o genes puntuales en la expresión de la dormición tanto en granos maduros como inmaduros de sorgo granífero.

Identificación de QTL para dormición y análisis molecular de la interacción ABA-GAs en cariopses de sorgo granífero

El objetivo de esta tesis fue contribuir a la dilucidación de los mecanismos moleculares y genéticos que participan en la expresión de la dormición de semillas de cereales, utilizando Sorghum bicolor (L.)Moench como sistema modelo. Para ello se utilizaron dos aproximaciones complementarias: la identificación de QTL para el carácter dormición y la evaluación de la ocurrencia de interacciones in vitro entre componentes de la señalización del ácido abscísico (ABA)y el catabolismo de las giberelinas (GAs), candidatos a tener un rol importante durante la expresión de la dormición en granos de sorgo inmaduros (i.e. antes de madurez fisiológica). Los resultados obtenidos permitieron identificar tres QTL (qDOR-5; qDOR-9 y qDOR-10)que explican una proporción de la variabilidad que se observa en el patrón de expresión de dormición de granos de sorgo maduros (i.e. después de madurez fisiológica). Un análisis in silico de las secuencias abarcadas por estos QTL mostró que ninguno ellos incluye genes considerados como candidatos para dormición de sorgo. En ese sentido, esta tesis aportó nuevas regiones genómicas que contienen genes hasta ahora desconocidos, que serían importantes en la expresión del carácter dormición en granos maduros. Por otra parte, los análisis de unión in vitro realizados mostraron que las proteínas SbABI4 y SbABI5 (componentes de la señalización del ABA)pueden interactuar de manera específica con el ABRC (complejo de respuesta al ABA)del promotor del gen SbGA2ox3, responsable de la degradación de giberelinas activas. Este mecanismo de cross-talk ABA-GAs podría ser uno de los responsables del mantenimiento de la dormición en cariopses inmaduros resistentes al brotado pre-cosecha. Más aún, el ABRC del promotor de SbGA2ox3, involucrado en las interacciones, se encontró además en los promotores de genes GA2ox de otras especies monocotiledóneas como Brachypodium y arroz (Oryza sativa), pero no así en las dicotiledóneas analizadas, sugiriendo que el cross-talk ABA-GAs podría tener lugar en otras especies además de sorgo. Los resultados de esta tesis en forma conjunta aportaron nuevas evidencias acerca del rol preponderante que tienen ciertas regiones del genoma o genes puntuales en la expresión de la dormición tanto en granos maduros como inmaduros de sorgo granífero.

Resolving candidate genes of mouse skeletal muscle QTL via RNA-Seq and expression network analyses

BACKGROUND:

We have recently identified a number of Quantitative Trait Loci (QTL) contributing to the 2-fold muscle weight difference between the LG/J and SM/J mouse strains and refined their confidence intervals. To facilitate nomination of the candidate genes responsible for these differences we examined the transcriptome of the tibialis anterior (TA) muscle of each strain by RNA-Seq.

13,726 genes were expressed in mouse skeletal muscle. Intersection of a set of 1061 differentially expressed transcripts with a mouse muscle Bayesian Network identified a coherent set of differentially expressed genes that we term the LG/J and SM/J Regulatory Network (LSRN). The integration of the QTL, transcriptome and the network analyses identified eight key drivers of the LSRN (Kdr, Plbd1, Mgp, Fah, Prss23, 2310014F06Rik, Grtp1, Stk10) residing within five QTL regions, which were either polymorphic or differentially expressed between the two strains and are strong candidates for quantitative trait genes (QTGs) underlying muscle mass. The insight gained from network analysis including the ability to make testable predictions is illustrated by annotating the LSRN with knowledge-based signatures and showing that the SM/J state of the network corresponds to a more oxidative state. We validated this prediction by NADH tetrazolium reductase staining in the TA muscle revealing higher oxidative potential of the SM/J compared to the LG/J strain (p<0.03).

Thus, integration of fine resolution QTL mapping, RNA-Seq transcriptome information and mouse muscle Bayesian Network analysis provides a novel and unbiased strategy for nomination of muscle QTGs.

Identification de QTL à pression artérielle dans le chromosome 18 du rat et analyse des gènes candidats Adrb2 et Nedd4l associés à l’hypertension essentielle

L’hypertension constitue un facteur majeur de risque de maladies cardiovasculaires et touche à un pourcentage important de la population humaine. Il s’agit d’une maladie complexe étant donné son caractère multifactoriel. La régulation de la pression artérielle (PA) est sous contrôle de plusieurs gènes, appelés loci pour traits quantitatifs ou QTL, et elle est le résultat de leur interaction. Étant donné que la PA est un trait quantitatif sous contrôle de plusieurs composantes comme les facteurs génétiques et environnementaux, l’étude de l’hypertension est limitée chez les populations humaines. Ainsi la stratégie principale pour l’étude de cette maladie est l’identification de QTL à PA chez des souches congéniques de rat construites à partir des lignées hyper- et normotendues; à savoir les souches Dahl salt-sensitive [1] et Lewis, respectivement. Des études précédentes dans notre laboratoire ont localisé trois QTL à PA au niveau du chromosome 18 chez le rat. Au cours de ce projet, de nouvelles sous-souches ont été construites afin de raffiner la cartographie de ces QTL. Ainsi, les C18QTL1, C18QTL3 et C18QTL4 ont été définis. Des analyses moléculaires ont été effectuées sur deux gènes candidats pour le C18QTL3; à savoir, Adrb2 et Nedd4l associés précédemment à l’hypertension. La comparaison des résultats de séquençage des régions régulatrices et codantes de ces deux gènes, ainsi que leur analyse d’expression par qRT-PCR chez les souches contrastantes DSS et Lewis, n’ont pas montré de différence significative pouvant expliquer la variation du phénotype observé. Des études plus poussées devront être effectuées sur ces deux gènes et, le cas échéant, l’analyse d’autres gènes contenus dans le C18QTL3 devra être entamée afin d’identifier le gène responsable de ce QTL.

Identification et regroupement de QTL influençant la pression artérielle en modules épistatiques et analyse de deux gènes candidats chez la souche Dahl Salt-Sensitive

L’hypertension essentielle étant un facteur majeur de morbidité, la compréhension de son l’étiologie est prépondérante. Ainsi, la découverte de nouvelles composantes ou mécanismes de régulation de la PA par l’identification de QTL et l’étude de leurs interactions s’avère une approche prometteuse. L’utilisation de souches congéniques de rats pour l’étude de l’hypertension est une stratégie payante puisqu’elle permet de masquer les effets de l’environnement, tout en gardant le caractère polygénique de la PA. Longtemps conçu comme un trait issu de l’accumulation des effets minimes des QTL, la PA est régulée par une architecture basée sur l’existence d’interactions épistatiques. L’analyse par paires de QTL individuels a permis d’établir une modularité dans l’organisation des QTL chez le rat Dahl Salt-sensitive en fonction de la présence ou de l’absence d’une interaction épistatique entre eux. Ainsi, deux modules épistatiques ont été établis; EM1 et EM2 où tous les QTL appartenant à EM1 sont épistatiques entre eux et agissent de façon additive avec les membres de EM2. Des hiérarchies dans la régulation peuvent alors être révélées si les QTL d’un même EM ont des effets opposés. L’identification de la nature moléculaire des candidats C18QTL4/Hdhd2 et C18QTL3/Tcof1, membres du EM1, et de l’interaction épistatique entre ces deux QTL, a permis, en plus, d’élucider une régulation séquentielle au sein du module. Hdhd2 pourrait agir en amont de Tcof1 et réguler ce dernier par une modification post-traductionnelle. Cette interaction est la première évidence expérimentale de la prédiction des relations entre QTL, phénomène établi par leur modularisation. Le dévoilement du fonctionnement de l’architecture génétique à la base du contrôle de la PA et la découverte des gènes responsables des QTL permettrait d’élargir les cibles thérapeutiques et donc de développer des traitements antihypertenseurs plus efficaces.

Use of synteny to identify candidate genes underlying QTL controlling stomatal traits in faba bean (Vicia faba L.)

Key message We have identified QTLs for stomatal characteristics on chromosome II of faba bean by applying SNPs derived from M. truncatula , and have identified candidate genes within these QTLs using synteny between the two species. Abstract Faba bean (Vicia faba L.) is a valuable food and feed crop worldwide, but drought often limits its production, and its genome is large and poorly mapped. No information is available on the effects of genomic regions and genes on drought adaptation characters such as stomatal characteristics in this species, but the synteny between the sequenced model legume, Medicago truncatula, and faba bean can be used to identify candidate genes. A mapping population of 211 F5 recombinant inbred lines (Mélodie/2 × ILB 938/2) were phenotyped to identify quantitative trait loci (QTL) affecting stomatal morphology and function, along with seed weight, under well-watered conditions in a climate-controlled glasshouse in 2013 and 2014. Canopy temperature (CT) was evaluated in 2013 under water-deficit (CTd). In total, 188 polymorphic single nucleotide polymorphisms (SNPs), developed from M. truncatula genome data, were assigned to nine linkage groups that covered ~928 cM of the faba bean genome with an average inter-marker distance of 5.8 cM. 15 putative QTLs were detected, of which eight (affecting stomatal density, length and conductance and CT) co-located on chromosome II, in the vicinity of a possible candidate gene—a receptor-like protein kinase found in the syntenic interval of M. truncatula chromosome IV. A ribose-phosphate pyrophosphokinase from M. truncatula chromosome V, postulated as a possible candidate gene for the QTL for CTd, was found some distance away in the same chromosome. These results demonstrate that genomic information from M. truncatula can successfully be translated to the faba bean genome.

The identification of QTL controlling ergot sclerotia size in hexaploid wheat implicates a role for the Rht dwarfing alleles

The fungal pathogen Claviceps purpurea infects ovaries of a broad range of temperate grasses and cereals, including hexaploid wheat, causing a disease commonly known as ergot. Sclerotia produced in place of seed carry a cocktail of harmful alkaloid compounds that result in a range of symptoms in humans and animals, causing ergotism. Following a field assessment of C. purpurea infection in winter wheat, two varieties ‘Robigus’ and ‘Solstice’ were selected which consistently produced the largest differential effect on ergot sclerotia weights. They were crossed to produce a doubled haploid mapping population, and a marker map, consisting of 714 genetic loci and a total length of 2895 cM was produced. Four ergot reducing QTL were identified using both sclerotia weight and size as phenotypic parameters; QCp.niab.2A and QCp.niab.4B being detected in the wheat variety ‘Robigus’, and QCp.niab.6A and QCp.niab.4D in the variety ‘Solstice’. The ergot resistance QTL QCp.niab.4B and QCp.niab.4D peaks mapped to the same markers as the known reduced height (Rht) loci on chromosomes 4B and 4D, Rht-B1 and Rht-D1, respectively. In both cases, the reduction in sclerotia weight and size was associated with the semi-dwarfing alleles, Rht-B1b from ‘Robigus’ and Rht-D1b from ‘Solstice’. Two-dimensional, two-QTL scans identified significant additive interactions between QTL QCp.niab.4B and QCp.niab.4D, and between QCp.niab.2A and QCp.niab.4B when looking at sclerotia size, but not between QCp.niab.2A and QCp.niab.4D. The two plant height QTL, QPh.niab.4B and QPh.niab.4D, which mapped to the same locations as QCp.niab.4B and QCp.niab.4D, also displayed significant genetic interactions.

Fine mapping and replication of QTL in outbred chicken advanced intercross lines

Background: Linkage mapping is used to identify genomic regions affecting the expression of complex traits. However, when experimental crosses such as F2 populations or backcrosses are used to map regions containing a Quantitative Trait Locus (QTL), the size of the regions identified remains quite large, i.e. 10 or more Mb. Thus, other experimental strategies are needed to refine the QTL locations. Advanced Intercross Lines (AIL) are produced by repeated intercrossing of F2 animals and successive generations, which decrease linkage disequilibrium in a controlled manner. Although this approach is seen as promising, both to replicate QTL analyses and fine-map QTL, only a few AIL datasets, all originating from inbred founders, have been reported in the literature. Methods: We have produced a nine-generation AIL pedigree (n = 1529) from two outbred chicken lines divergently selected for body weight at eight weeks of age. All animals were weighed at eight weeks of age and genotyped for SNP located in nine genomic regions where significant or suggestive QTL had previously been detected in the F2 population. In parallel, we have developed a novel strategy to analyse the data that uses both genotype and pedigree information of all AIL individuals to replicate the detection of and fine-map QTL affecting juvenile body weight. Results: Five of the nine QTL detected with the original F2 population were confirmed and fine-mapped with the AIL, while for the remaining four, only suggestive evidence of their existence was obtained. All original QTL were confirmed as a single locus, except for one, which split into two linked QTL. Conclusions: Our results indicate that many of the QTL, which are genome-wide significant or suggestive in the analyses of large intercross populations, are true effects that can be replicated and fine-mapped using AIL. Key factors for success are the use of large populations and powerful statistical tools. Moreover, we believe that the statistical methods we have developed to efficiently study outbred AIL populations will increase the number of organisms for which in-depth complex traits can be analyzed.

qtl.outbred: Interfacing outbred line cross data with the R/qtl mapping software

Background qtl.outbred is an extendible interface in the statistical environment, R, for combining quantitative trait loci (QTL) mapping tools. It is built as an umbrella package that enables outbred genotype probabilities to be calculated and/or imported into the software package R/qtl. Findings Using qtl.outbred, the genotype probabilities from outbred line cross data can be calculated by interfacing with a new and efficient algorithm developed for analyzing arbitrarily large datasets (included in the package) or imported from other sources such as the web-based tool, GridQTL. Conclusion qtl.outbred will improve the speed for calculating probabilities and the ability to analyse large future datasets. This package enables the user to analyse outbred line cross data accurately, but with similar effort than inbred line cross data.

Genetic parameters for body weight, carcass chemical composition and yield in a broiler-layer cross developed for QTL mapping

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)