A minisatellite and a microsatellite polymorphism within 1.5 kb at the human muscle glycogen phosphorylase (PYGM) locus can be amplified by PCR and have combined informativeness of PIC 0.95.

We sequenced a genomic clone (pMCMP1), previously reported to detect a VNTR polymorphism at the PYGM locus, and found a dinucleotide repeat segment (CA)14(GA)25 and a complex (AT)-repeat-rich segment containing 63 repeats spanning 160 bp. Resolution of PCR-amplified genomic DNA from the (CA)(GA) repeat region on DNA sequencing gels revealed a highly informative polymorphism with alleles differing by 2-bp intervals and ranging in size from 156 to 190 bp. Among three racial groups, a total of 18 alleles were observed. Fourteen alleles were observed in Caucasians (PIC 0.89), 12 alleles in American Blacks (PIC 0.89), and 9 alleles in Pima Indians (PIC 0.73). PCR amplification of the (AT) repeat region and resolution of the products on DNA sequencing gels revealed a complex variable length polymorphism with alleles distributed in size from 367 to 970 bp. Twenty-eight alleles were found in American Blacks (PIC 0.94), 6 alleles in Pima Indians (PIC 0.70), and 11 alleles in Caucasians (PIC 0.71). Comparison of the previously described VNTR RFLP alleles visualized by Southern hybridization to the PCR products described in this report demonstrated that the polymorphism described in both assays was identical. However, a larger number of alleles could be detected from the PCR-amplified products. Combined informativeness, PIC 0.95, for the two polymorphisms was determined from haplotype analysis of 100 Caucasian chromosomes. Therefore, for genotyping purposes, informativeness is maximized from using both polymorphisms.     

A polymorphic (CA)n repeat element maps the human glucokinase gene (GCK) to chromosome 7p.

A compound imperfect dinucleotide repeat element, [CA]4TTTGT[CT]7[CA]9AA[CA]4CCACATA[CA]3, was found approximately 10 kb 3' to the human glucokinase gene (GCK) from analysis of contiguous genomic DNA obtained from a bacteriophage lambda chromosome walk. Direct human genomic sequencing revealed the source of polymorphism to be variable numbers of CT and CA repeats. Altogether six alleles that range in length from +10 to -15 nucleotides compared to the most common (Z) allele have been identified. Alleles Z, Z + 2, and Z + 4 were present in American Blacks, Pima Indians, and Caucasians, with somewhat varied frequencies among the groups. Two alleles, Z + 10 and Z - 15, appear to be unique to American Blacks, while a Z + 6 allele was observed only in the Caucasian population studied. Observed heterozygosity of the polymorphism in the CEPH reference pedigree collection is 44% and the PIC 0.44. The polymorphism is assayed by PCR amplification and resolution of 32P-end-labeled products (ranging in length from 180 to 205 bp) on denaturing polyacrylamide sequencing gels. Using the PCR assay, the human glucokinase gene was physically localized to chromosome 7 in a panel of rodent/human somatic cell lines. Genetic analysis in CEPH pedigrees placed the dinucleotide repeat element, and thereby the human glucokinase gene, on chromosome 7p between TCRG and a RFLP locus D7S57. The glucokinase dinucleotide repeat genetic marker can now be used to assess the role of the glucokinase gene in diabetes by population association studies. In addition, this repeat marker and others flanking it on chromosome 7 can be used in linkage studies with families segregating the disorder.     

A 2-cM genetic linkage map of human chromosome 7p that includes 47 loci.

A new high-resolution genetic linkage map for human chromosome 7p has been constructed. The map is composed of 47 loci (54 polymorphic systems), 19 of which are uniquely placed with odds of at least 1000:1. Four genes are represented, including glucokinase (GCK, ATP:D-hexose-6-phosphotransferase, EC 2.7.1.2) which was mapped via a (CA)n dinucleotide repeat polymorphism. The sex-average map measures 94.4 cM and the male and female maps measure 73.2 and 116.1 cM, respectively. We believe that the genetic map extends nearly the full length of the short arm of chromosome 7 since a centromere marker has been incorporated, and the most distal marker, D7S21, has been cytogenetically localized by in situ hybridization to 7p22-pter. The average marker spacing is 2 cM, and the largest interval between uniquely placed markers is 13 cM (sex-average map). Overall, female recombination was observed to be about 1.5 times that of males, and a statistically significant sex-specific recombination frequency was found for a single interval. The map is based on genotypic data gathered from 40 CEPH reference pedigrees and was constructed using the CRI-MAP program package. The map presented here represents a combined and substantially expanded dataset compared to previously published chromosome 7 maps, and it will serve as a "baseline" genetic map that should prove useful for future efforts to develop a 1-cM map and for construction of a contiguous clone-based physical map for this chromosome.     

Dinucleotide Repeat Loci Contribute Highly Informative Genetic Markers to the Human Chromosome 2 Linkage Map

Microsatellite repeat loci can provide informative markers for genetic linkage. Currently, the human chromosome 2 genetic linkage map has very few highly polymorphic markers. Being such a large chromosome, it will require a large number of informative markers for the dense coverage desired to allow disease genes to be mapped quickly and accurately. Dinucleotide repeat loci from two anonymous chromosome 2 genomic DNA clones were sequenced so that oligonucleotide primers could be designed for amplifying each locus using the polymerase chain reaction (PCR). Five sets of PCR primers were also generated from nucleotide sequences in the GenBank Database of chromosome 2 genes containing dinucleotide repeats. In addition, one PCR primer pair was made that amplifies a restriction fragment length polymorphism on the TNP1 gene (Hoth and Engel, 1991). These markers were placed on the CEPH genetic linkage map by screening the CEPH reference DNA panel with each primer set, combining these data with those of other markers previously placed on the map, and analyzing the combined data set using CRI-MAP and LINKAGE. The microsatellite loci are highly informative markers and the TNP1 locus, as expected, is only moderately informative. A map was constructed with 38 ordered loci (odds 1000:1) spanning 296 cM (male) and 476 cM (female) of chromosome 2 compared with 306 cM (male) and 529 cM (female) for a previous map of 20 markers.     

Construction of a high-density linkage map of Italian ryegrass (Lolium multiflorum Lam) using restriction fragment length polymorphism, amplified fragment length polymorphism, and telomeric repeat associated sequence markers.

To construct a high-density molecular linkage map of Italian ryegrass (Lolium multiflorum Lam), we used a two-way pseudo-testcross F1 population consisting of 82 individuals to analyze three types of markers: restriction fragment length polymorphism markers, which we detected by using genomic probes from Italian ryegrass as well as heterologous anchor probes from other species belonging to the Poaceae family, amplified fragment length polymorphism markers, which we detected by using PstI/MseI primer combinations, and telomeric repeat associated sequence markers. Of the restriction fragment length polymorphism probes that we generated from a PstI genomic library, 74% (239 of 323) of randomly selected probes detected hybridization patterns consistent with single-copy or low-copy genetic locus status in the screening. The 385 (mostly restriction fragment length polymorphism) markers that we selected from the 1226 original markers were grouped into seven linkage groups. The maps cover 1244.4 cM, with an average of 3.7 cM between markers. This information will prove useful for gene targeting, quantitative trait loci mapping, and marker-assisted selection in Italian ryegrass.     

Microsatellite marker polymorphism and mapping in pea (Pisum sativum L.)

This paper aims at providing reliable and cost effective genotyping conditions, level of polymorphism in a range of genotypes and map position of newly developed microsatellite markers in order to promote broad application of these markers as a common set for genetic studies in pea. Optimal PCR conditions were determined for 340 microsatellite markers based on amplification in eight genotypes. Levels of polymorphism were determined for 309 of these markers. Compared to data obtained for other species, levels of polymorphism detected in a panel of eight genotypes were high with a mean number of 3.8 alleles per polymorphic locus and an average PIC value of 0.62, indicating that pea represents a rather polymorphic autogamous species. One of our main objectives was to locate a maximum number of microsatellite markers on the pea genetic map. Data obtained from three different crosses were used to build a composite genetic map of 1,430 cM (Haldane) comprising 239 microsatellite markers. These include 216 anonymous SSRs developed from enriched genomic libraries and 13 SSRs located in genes. The markers are quite evenly distributed throughout the seven linkage groups of the map, with 85% of intervals between the adjacent SSR markers being smaller than 10 cM. There was a good conservation of marker order and linkage group assignment across the three populations. In conclusion, we hope this report will promote wide application of these markers and will allow information obtained by different laboratories worldwide in diverse fields of pea genetics, such as QTL mapping studies and genetic resource surveys, to be easily aligned.Electronic Supplementary Material Supplementary material is available for this article at     

The dinucleotide (CA) repeat polymorphism of estrogen receptor beta but not the dinucleotide (TA) repeat polymorphism of estrogen receptor alpha is associated with venous ulceration

Venous ulcers are the predominant form of chronic wound in the elderly, accounting for around 70% of all cases. The steroid sex hormone estrogen plays a crucial role in normal human skin maintenance and during cutaneous wound repair following injury. Estrogen can reverse age-related impaired wound healing by dampening the inflammatory response and increasing matrix deposition at the wound site. The molecular actions of estrogen are mediated through two nuclear sex steroid hormone receptors, estrogen receptor alpha (ERalpha) and beta (ERbeta). We have conducted a case-control study to investigate whether dinucleotide repeat polymorphisms in the estrogen receptor genes are associated with venous ulceration in the UK Caucasian population. Genomic fragments containing the ERalpha dinucleotide (TA)(n) repeat polymorphism or the ERbeta dinucleotide (CA)(n) repeat polymorphism were amplified by polymerase chain reaction in subject DNA samples and genotyped according to fragment length by capillary electrophoresis. There was no evidence to suggest that the TA repeat polymorphism of ERalpha was associated with venous ulceration. However, the CA*18 allele of the ERbeta CA repeat polymorphism was significantly associated with venous ulceration (n = 120, OR = 1.8, 95% CI = 1.1-2.8, P = 0.02). When the CA repeats alleles were grouped together into either low (L < or = 18) or high (H > 18) numbers of CA repeats, the low (L) repeat allele was significantly associated with venous ulceration (OR = 1.5, 95% CI = 1.0-2.2, P = 0.03). Our results show that a specific ERbeta variant is associated with impaired healing in the elderly, predisposing individuals to venous ulceration.     

The first genetic linkage map of Luohanguo (Siraitia grosvenorii ) based on ISSR and SRAP markers

In this study, the first genetic map of Luohanguo (Siraitia grosvenorii (Swingle) C. Jeffrey) was constructed with 150 F? population individuals using inter-simple sequence repeat (ISSR) and sequence-related amplified polymorphism (SRAP) markers. A total of 100 ISSRs and 196 SRAP primer combinations generated 51 and 222 polymorphic markers, respectively. Among the 273 markers obtained, 199 markers (29 ISSRs and 170 SRAPs) were mapped to 25 linkage groups. The map covered 1463.3 cM with a mean map distance of 7.35 cM between adjacent markers and a maximum map distance of 52.6 cM between two markers. The markers were distributed randomly in 25 groups except for minor clusters in the distal region of linkage groups. All 25 linkage groups consisted of 2-36 loci ranging in length from 19.5 to 152.6 cM and accounted for 59.8% of the total map distance. This map provides reference information for future molecular breeding work on Luohanguo.     

Molecular mapping of the Oregon Wolfe Barleys: a phenotypically polymorphic doubled-haploid population

A phenotypically polymorphic barley (Hordeum vulgare L.) mapping population was developed using morphological marker stocks as parents. Ninety-four doubled-haploid lines were derived for genetic mapping from an F₁ using the Hordeum bulbosum system. A linkage map was constructed using 12 morphological markers, 87 restriction fragment length polymorphism (RFLP), five random amplified polymorphic DNA (RAPD), one sequence-tagged site (STS), one intron fragment length polymorphism (IFLP), 33 simple sequence repeat (SSR), and 586 amplified fragment length polymorphism (AFLP) markers. The genetic map spanned 1,387 cM with an average density of one marker every 1.9 cM. AFLP markers tended to cluster on centromeric regions and were more abundant on chromosome 1 (7H). RAPD markers showed a high level of segregation distortion, 54% compared with the 26% observed for AFLP markers, 27% for SSR markers, and 18% for RFLP markers. Three major regions of segregation distortion, based on RFLP and morphological markers, were located on chromosomes 2 (2H), 3 (3H), and 7 (5H). Segregation distortion may indicate that preferential gametic selection occurred during the development of the doubled-haploid lines. This may be due to the extreme phenotypes determined by alleles at morphological trait loci of the dominant and recessive parental stocks. Several molecular markers were found to be closely linked to morphological loci. The linkage map reported herein will be useful in integrating data on quantitative traits with morphological variants and should aid in map-based cloning of genes controlling morphological traits. Received: 23 August 2000 / Accepted: 15 December 2000     

A PCR-Based Genetic Map for Human Chromosome 3

Oligonucleotide primers for 125 simple sequence repeat microsatellite-based genetic markers have been assayed by polymerase chain reaction (PCR) in the CEPH reference family panel. These microsatellites include 101 dinucleotide repeats as well as 24 new tetranucleotide repeats. The average heterozygosity of this marker set was 72.4%. Genetic data were analyzed with the genetic mapping package LINKAGE. A subset of these microsatellite markers define a set of 56 uniquely ordered loci (>1000:1 against local inversion) that span 271 cM. Sixty-seven additional loci were tightly linked to markers on the uniquely ordered map, but could not be ordered with such high precision. These markers were positioned by CMAP into confidence intervals. One hundred thirteen of the microsatellite markers were also tested on a chromosome 3 framework somatic cell hybrid panel that divides this chromosome into 23 cytogenetically defined regions, integrating the genetic and physical maps of this chromosome. The high density, high heterozygosity, and PCR format of this genetically and physically mapped set of markers will accelerate the mapping and positional cloning of new chromosome 3 genes.     

The construction of a genetic linkage map of non-heading Chinese cabbage （Brassica campestris ssp. chinensis Makino）

Non-heading Chinese cabbage （Brassica carnpestris ssp. chinensis Makino） is one of the most important vegetables in eastern China. A genetic linkage map was constructed using 127 doubled haploid （DH） lines, and the DH population was derived from a commercial hybrid ＂Hanxiao＂ （lines SW-13 x L-118）. Out of the 614 polyrnorphic markers, 43.49% were not assigned to any of the linkage groups （LGs）. Chi-square tests showed that 42.67% markers were distorted from expected Mendelian segregation ratios, and the direction of distorted segregation was mainly toward the paternal parent L-118. After sequentially removing the markers that had an interval distance smaller than 1 cM from the upper marker, the overall quality of the linkage map was increased. Two hundred and sixty-eight molecular markers were mapped into 10 LGs, which were anchored to the corresponding chromosome of the B. rapa reference map based on com- mon simple sequence repeat （SSR） markers. The map covers 973.38 cM of the genome and the average interval distance between markers was 3.63 cM. The number of markers on each LG ranged from 18 （R08） to 64 （R07）, with an average interval distance within a single LG from 1.70 cM （R07） to 6.71 cM （R06）. Among these mapped markers, 169 were sequence-related amplified polymorphism （SRAP） molecular markers, 50 were SSR markers and 49 were random amplification polymorphic DNA （RAPD） markers. With further saturation to the LG9 the current map offers a genetic tool for loci analysis for important agronomic traits.     

Genetic mapping of 14 short tandem repeat polymorphisms on human chromosome 22

We have constructed a linkage map of 14 short tandem repeat polymorphisms (11 with heterozygosity > 70%) on the long arm of human chromosome 22 using 23 non-CEPH pedigrees. Twelve of the markers could be positioned uniquely with a likelihood of at least 1,000:1, and distributed at an average distance of 6.62 cM (range 1.5–16.1 cM). The sex-combined map covers a total of 79.6 cM, the female map 93.2 cM and the male map 64.6 cM. Based on comparisons between physical maps and other genetic maps, we estimate that our map covers 70%–80% of the chromosome. The map integrates markers from previous genetic maps and uniquely positions one marker (D22S307). Data from physical mapping on the location of four genetic markers correlates well with our linkage map, and provides information on an additional marker (D22S315). This map will facilitate high resolution mapping of additional polymorphic loci and disease genes on chromosome 22, and act as a reference for building and verifying physical maps.     

A (CA)n repeat polymorphism for the human skeletal muscle alpha-actinin gene ACTN2 and its localization on the linkage map of chromosome 1.

A CA dinucleotide repeat polymorphism has been identified for the skeletal muscle alpha-actinin gene ACTN2. The observed heterozygosity is 44% (predicted heterozygosity 50%, PIC 0.47). This polymorphic marker has been localized between D1S74 and D1S103 on the multipoint linkage map of chromosome 1 at a position 44.4 cM from the most distal marker D1S68 at 1 qter.     

Exclusion of linkage between familial Mediterranean fever and the human serum amyloid A (SAA) gene cluster

Summary We studied the relationship between the autosomal recessive trait familial Mediterranean fever (FMF) and the serum amyloid A (SAA) genes by comparing alleles of a highly polymorphic dinucleotide repeat and a conventional restriction fragment length polymorphism (RFLP) in the SAA gene cluster in Israeli FMF kindreds. By haplotype analysis, our data indicate a minimum crossover frequency of 22% between the SAA gene marker and FMF. By conventional linkage analysis this eliminates a minimum of 10.4cM including and surrounding the SAA gene cluster as the site of the FMF mutation although SAA proteins are prominent physiologic markers of the acute attacks.     

Molecular mapping of the male-sterile, female-sterile mutant gene (st8) in soybean

Soybean male-sterile, female-sterile mutant genes have been identified by genetic and cytological studies. The St8 gene has been identified as an asynaptic mutation resulting in male and female sterility. This mutant gene was derived from a gene-tagging study using the soybean w4-mutable line. In this report we identified the genetic map position of st8 via restriction fragment length polymorphism (RFLP) and simple sequence repeat (SSR) markers. The St8 gene mutation was located between RFLP marker E107 and SSR markers Satt132, Sct_065, and Satt414 on molecular linkage group J and linked to each by 7.8 cM and 3.4 cM, respectively.     

Angiotensin-converting enzyme gene (ACE) insertion/deletion polymorphism in Mexican populations

The angiotensin-converting enzyme gene (ACE) insertion/deletion polymorphism was determined in 211 Mexican healthy individuals belonging to different Mexican ethnic groups (98 Mestizos, 64 Teenek, and 49 Nahuas). ACE polymorphism differed among Mexicans with a high frequency of the D allele and the D/D genotype in Mexican Mestizos. The D/D genotype was absent in Teenek and present in only one Nahua individual (2.0%). When comparisons were made, we observed that Caucasian, African, and Asian populations presented the highest frequencies of the D allele, whereas Amerindian (Teenek and Pima) and Australian Aboriginals showed the highest frequencies of the I allele. The distribution of I/D genotype was heterogeneous in all populations: Australian Aboriginals presented the lowest frequency (4.9%), whereas Nahuas presented the highest (73.4%). The present study shows the frequencies of a polymorphism not analyzed previously in Mexican populations and establishes that this polymorphism distinguishes the Amerindian populations of other groups. On the other hand, since ACE alleles have been associated with genetic susceptibility to developing cardiovascular diseases and hypertension, knowledge of the distribution of these alleles could help to define the true significance of ACE polymorphism as a genetic susceptibility marker in the Amerindian populations.     

An integrative genetic linkage map of winter wheat (<Emphasis Type="Italic">Triticum aestivum</Emphasis> L.)

We constructed a genetic linkage map based on a cross between two Swiss winter wheat (Triticum aestivum L.) varieties, Arina and Forno. Two-hundred and forty F₅ single-seed descent (SSD)-derived lines were analysed with 112 restriction fragment length polymorphism (RFLP) anonymous probes, 18 wheat cDNA clones coding for putative stress or defence-related proteins and 179 simple-sequence repeat (SSR) primer-pairs. The 309 markers revealed 396 segregating loci. Linkage analysis defined 27 linkage groups that could all be assigned to chromosomes or chromosome arms. The resulting genetic map comprises 380 loci and spans 3,086 cM with 1,131 cM for the A genome, 920 cM for the B genome and 1,036 cM for the D genome. Seventeen percent of the loci showed a significant (P < 0.05) deviation from a 1:1 ratio, most of them in favour of the Arina alleles. This map enabled the mapping of QTLs for resistance against several fungal diseases such as Stagonospora glume blotch, leaf rust and Fusarium head blight. It will also be very useful for wheat genetic mapping, as it combines RFLP and SSR markers that were previously located on separate maps. S. Paillard and T. Schnurbusch contributed equally to the work     

A genetic map of melon (Cucumis melo L.) based on amplified fragment length polymorphism (AFLP) markers

 Genetic maps facilitate the study of genome structure and evolution, and the identification of monogenic traits or Mendelian components of quantitative traits. We evaluated 228 RAPD, microsatellite and AFLP markers for linkage analysis in melon (Cucumis melo L.) varieties MR-1 (resistant to Fusarium wilt, powdery and downy mildews) and Ananas Yokneum (AY; susceptible to these diseases) and constructed a detailed genetic map. The mapping population consisted of 66 backcross progenies derived from AY×(MR-1×AY). Despite a relatively low level of polymorphism in the species, AFLP markers were found to be more efficient in mapping the melon genome than RAPD or microsatellite markers. The map contains 197 AFLPs, six RAPDs and one microsatellite marker assigned to 14 major and six minor linkage groups, and covers 1942 cM with the average distance between adjacent markers of approximately 10 cM. The maximum distance allowed between markers is 27.5 cM. About 11% of the intervals (20 out of 173) are over 20 cM (but less than 27.5 cM). The map has immediate utility for identifying markers linked to disease resistance genes that are suitable for marker-assisted breeding. The use of microsatellite markers for integration with other maps is also discussed. Received: 12 March 1997 / Accepted: 20 May 1997     

Optimization of linkage mapping strategy and construction of a high-density American lotus linkage map

Background

Lotus is a diploid plant with agricultural, medicinal, and ecological significance. Genetic linkage maps are fundamental resources for genome and genetic study, and also provide molecular markers for breeding in agriculturally important species. Genotyping by sequencing revolutionized genetic mapping, the restriction-site associated DNA sequencing (RADseq) allowed rapid discovery of thousands of SNPs markers, and a crucial aspect of the sequence based mapping strategy is the reference sequences used for marker identification.

Results

We assessed the effectiveness of linkage mapping using three types of references for scoring markers: the unmasked genome, repeat masked genome, and gene models. Overall, the repeat masked genome produced the optimal genetic maps. A high-density genetic map of American lotus was constructed using an F₁ population derived from a cross between Nelumbo nucifera ‘China Antique’ and N. lutea ‘AL1’. A total of 4,098 RADseq markers were used to construct the American lotus ‘AL1’ genetic map, and 147 markers were used to construct the Chinese lotus ‘China Antique’ genetic map. The American lotus map has 9 linkage groups, and spans 494.3 cM, with an average distance of 0.7 cM between adjacent markers. The American lotus map was used to anchor scaffold sequences in the N. nucifera ‘China Antique’ draft genome. 3,603 RADseq markers anchored 234 individual scaffold sequences into 9 megascaffolds spanning 67% of the 804 Mb draft genome.

Conclusions

Among the unmasked genome, repeat masked genome and gene models, the optimal reference sequences to call RADseq markers for map construction is repeat masked genome. This high density genetic map is a valuable resource for genomic research and crop improvement in lotus.     

Genetic maps of SSR and SRAP markers in diploid orchardgrass (Dactylis glomerata L.) using the pseudo-testcross strategy

Orchardgrass (Dactylis glomerata L.) is one of the most important cool-season forage grasses commonly grown throughout the temperate regions of the world. The objective of this work was to construct a diploid (2n = 2x = 14) orchardgrass genetic linkage map useful as a framework for basic genetic studies and plant breeding. A combination of simple sequence repeat (SSR) and sequence-related amplified polymorphism (SRAP) molecular markers were used for map construction. The linkage relationships among 164 SSRs and 108 SRAPs, assayed in a pseudo-testcross F1 segregating population generated from a cross between two diploid parents, were used to construct male (01996) and female (YA02-103) parental genetic maps. The paternal genetic map contains 90 markers (57 SSRs and 33 SRAPs) over 9 linkage groups (LGs), and the maternal genetic map is composed of 87 markers (54 SSRs and 33 SRAPs) assembled over 10 LGs. The total map distance of the male map is 866.7 centimorgans (cM), representing 81% genome coverage, whereas the female map spans 772.0 cM, representing 75% coverage. The mean map distance between markers is 9.6 cM in the male map and 8.9 cM in the female map. About 14% of the markers remained unassigned. The level of segregation distortion observed in this cross was 15%. Homology between the two maps was established between five LGs of the male map and five LGs of the female map using 10 bridging markers. The information presented in this study establishes a foundation for extending genetic mapping in this species, serves as a framework for mapping quantitative trait loci (QTLs), and provides basic information for future molecular breeding studies.